Eigenvalues 3x3 matrix
Submitted by tho..@..gmx.de
Assigned to Nobody
Link to original bugzilla bug (#1014)
Created attachment 575
There is floating point comparison issue with the routine in SelfAdjointEigenSolver.h (Version 3.2.4) for fixed-size 3x3 matrices.
A tolerance is created as follows:
Scalar safeNorm2 = Eigen::NumTraits<Scalar>::epsilon();
In an older version (3.2.2), if followed the line
safeNorm2 *= safeNorm2;
which has been omitted. For matrices which are close to a multiple of the identity matrix, this change becomes crucial in some comparisons.
Attachment 575, "Test code":