Eigenvalues 3x3 matrix

Submitted by tho..@..gmx.de

Assigned to Nobody

Link to original bugzilla bug (#1014)
Version: 3.2

Description

Created attachment 575
Test code

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":
eigenValues3x3.cpp

Blocking

#814 (closed)