[eigen] Bug in SelfAdjointEigenSolver<>, or am I missing something important?

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

Hello all,

I attach a little test program. It basically takes a square symmetric
matrix and computes its eigen{vectors,values} with
SelfAdjointEigenSolver<>, then reconstruct the original matrix with
EVEC * diag(EVALS) * EVEC^T.

The point is, that works for a ColMajor ordering, but not with a
RowMajor, so, for these two tests:

	myTest< Eigen::Matrix<double,4,4,ColMajor> >();
	myTest< Eigen::Matrix<double,4,4,RowMajor> >();

it fails in the latter.

Mail converted by MHonArc 2.6.19+ http://listengine.tuxfamily.org/