Created attachment 31 [details] test case 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. Best regards, JL
-- GitLab Migration Automatic Message -- This bug has been migrated to gitlab.com's GitLab instance and has been closed from further activity. You can subscribe and participate further through the new bug through this link to our GitLab instance: https://gitlab.com/libeigen/eigen/issues/107.