[eigen] Permutation times Sparse matrix |
[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]
When I execute the program below I expect to get L= 1 0 0 1 But I get L = 0 0 0 0 # include <Eigen/SparseCore> # include <iostream> int main() { Eigen::SparseMatrix<double> L(2,2); L.insert(0,0) = 1.0; L.insert(1,1) = 1.0; Eigen::PermutationMatrix<2,2> P; P.setIdentity(); L = P * L; std::cout << "L =\n" << Eigen::Matrix<double,2,2>(L) << "\n"; return 0; }
Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |