GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-14 20:56:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

Ignore denormal numbers in selfadjoint eigensolver.

Browse Source
This commit is contained in:
Gael Guennebaud 2015-05-12 18:38:43 +02:00
parent a605a1d7df
commit ef81730625
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h
View File

@ -474,11 +474,14 @@ ComputationInfo computeFromTridiagonal_impl(DiagType& diag, SubDiagType& subdiag
Index end = n-1; Index end = n-1;
Index start = 0; Index start = 0;
Index iter = 0; // total number of iterations Index iter = 0; // total number of iterations
typedef typename DiagType::RealScalar RealScalar;
const RealScalar considerAsZero = (std::numeric_limits<RealScalar>::min)();
while (end>0) while (end>0)
{ {
for (Index i = start; i<end; ++i) 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; subdiag[i] = 0;
// find the largest unreduced block // find the largest unreduced block