diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2015-05-12 18:38:43 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2015-05-12 18:38:43 +0200 |
| commit | ef817306252f33e844eb789ef850778de84f47cb (patch) | |
| tree | de57e87243e46cda16d9780a6ec79f6cacb53b4b /Eigen | |
| parent | a605a1d7df3bd94e78d23b779c85660046f29829 (diff) | |
Ignore denormal numbers in selfadjoint eigensolver.
Diffstat (limited to 'Eigen')
| -rw-r--r-- | Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h | 7 |
1 files changed, 5 insertions, 2 deletions
diff --git a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h index 8240b9465..fff331d33 100644 --- a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +++ b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h @@ -474,11 +474,14 @@ ComputationInfo computeFromTridiagonal_impl(DiagType& diag, SubDiagType& subdiag Index end = n-1; Index start = 0; Index iter = 0; // total number of iterations - + + typedef typename DiagType::RealScalar RealScalar; + const RealScalar considerAsZero = (std::numeric_limits<RealScalar>::min)(); + while (end>0) { for (Index i = start; i<end; ++i) - if (internal::isMuchSmallerThan(abs(subdiag[i]),(abs(diag[i])+abs(diag[i+1])))) + if (internal::isMuchSmallerThan(abs(subdiag[i]),(abs(diag[i])+abs(diag[i+1]))) || abs(subdiag[i]) <= considerAsZero) subdiag[i] = 0; // find the largest unreduced block |