mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2024-12-21 07:19:46 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix underflow in 3x3 tridiagonalization

Browse Source
This commit is contained in:
Gael Guennebaud 2015-06-23 14:54:31 +02:00
parent 18c9d155f3
commit b0d08869a9
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
Eigen/src/Eigenvalues/Tridiagonalization.h
View File

@ -464,9 +464,10 @@ struct tridiagonalization_inplace_selector<MatrixType,3,false>
static void run(MatrixType& mat, DiagonalType& diag, SubDiagonalType& subdiag, bool extractQ) static void run(MatrixType& mat, DiagonalType& diag, SubDiagonalType& subdiag, bool extractQ)
{ {
using std::sqrt; using std::sqrt;
const RealScalar tol = (std::numeric_limits<RealScalar>::min)();
diag[0] = mat(0,0); diag[0] = mat(0,0);
RealScalar v1norm2 = numext::abs2(mat(2,0)); RealScalar v1norm2 = numext::abs2(mat(2,0));
if(v1norm2 == RealScalar(0)) if(v1norm2 <= tol)
{ {
diag[1] = mat(1,1); diag[1] = mat(1,1);
diag[2] = mat(2,2); diag[2] = mat(2,2);