diff options
author | Gael Guennebaud <g.gael@free.fr> | 2016-07-23 17:52:31 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2016-07-23 17:52:31 +0200 |
commit | 1b0353c659891c4cad61dddbe60c5e1245a85eac (patch) | |
tree | aeb1ba234dfbc286cbcf67740546a266dd1da8d6 /Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h | |
parent | 72744d93efda9d2394d2b0ababf2bd654c05d047 (diff) |
Fix misuse of dummy_precesion in eigenvalues solvers
Diffstat (limited to 'Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h')
-rw-r--r-- | Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h index f96fe698e..9f0583fd2 100644 --- a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +++ b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h @@ -492,11 +492,12 @@ ComputationInfo computeFromTridiagonal_impl(DiagType& diag, SubDiagType& subdiag typedef typename DiagType::RealScalar RealScalar; const RealScalar considerAsZero = (std::numeric_limits<RealScalar>::min)(); + const RealScalar precision = RealScalar(2)*NumTraits<RealScalar>::epsilon(); while (end>0) { for (Index i = start; i<end; ++i) - if (internal::isMuchSmallerThan(abs(subdiag[i]),(abs(diag[i])+abs(diag[i+1]))) || abs(subdiag[i]) <= considerAsZero) + if (internal::isMuchSmallerThan(abs(subdiag[i]),(abs(diag[i])+abs(diag[i+1])),precision) || abs(subdiag[i]) <= considerAsZero) subdiag[i] = 0; // find the largest unreduced block |