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

[Jose Luis Blanco <joseluisblancoc@xxxxxxxxx>]

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.