NOTE: I tried to file the issue below as a (potential) Eigen bug in Bugzilla,
but I never received the promised confirmation mail for my registration.
More Bugzilla problems still?
Dear friends,
With today's trunk version of Eigen, it appears that left-multiplication of a
ColMajor sparse matrix and right-multiplication of a RowMajor sparse
matrix with a DiagonalWrapper fail. On gcc-4.8.3 I get a segmentation
fault or a bus error when passing -DNDEBUG *AND* -O[1,2} is used.
Otherwise, everything works fine
I reduced this to the code below. (I apologize if i mis-understand Eigen
and that code is in error. I am an Eigen-novice still...)
Regards, Jan.
// sparse_test.cpp:
#include "Eigen/Core"
#include "Eigen/Sparse"
#include <iostream>
typedef double value_type;
typedef Eigen::Matrix<value_type, Eigen::Dynamic, 1> vector_type;
typedef Eigen::SparseMatrix<value_type, Eigen::ColMajor> cm_matrix_type;
typedef Eigen::SparseMatrix<value_type, Eigen::RowMajor> rm_matrix_type;
int main()
{
cm_matrix_type cmA(1,1);
cmA.insert(0,0) = 1;
rm_matrix_type rmA(1,1);
rmA.insert(0,0) = 1;
vector_type d(1);
d[0] = 2;
// OK:
std::cout << ( cmA*d.asDiagonal() ).eval();
std::cout << ( d.asDiagonal()*rmA ).eval();
// with gcc-4.8.3:
// - Bus error or segmentation fault with -DNDEBUG and -O{1,2}
std::cout << ( rmA*d.asDiagonal() ).eval();
std::cout << ( d.asDiagonal()*cmA ).eval();
return 0;
}
jan@gum15: gcc --version
gcc (SUSE Linux) 4.8.3 20140627 [gcc-4_8-branch revision 212064]
jan@gum15: g++ -DNDEBUG -O2 -I /home/jan/local/eigen/include/eigen3
sparse_test.cpp && ./a.out
2
2
Segmentation fault (core dumped)
(Omit -DNDEBUG *or* -O2 and everything works like a charm.)