diff options
| author |  Antonio Sanchez <cantonios@google.com> | 2021-04-07 15:08:39 -0700 |
|---|---|---|
| committer | Antonio Sanchez <cantonios@google.com> | 2021-04-07 15:08:39 -0700 |
| commit | fcb5106c6e16599eeabadf7d82a465d52229698f (patch) | |
| tree | 5ecc5884c857db7f26bc7394025b5d3d2ca5d962 /Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h | |
| parent | 6197ce1a3544f52e881a112ede9679d84fafebf3 (diff) | |
Scaled epsilon the wrong way.
Should have been 0.5 to widen the bounds, since this is inverse
precision. Setting to 0.5, however, leads to many more failing
tests at Google, so reverting to 1 for now.
Diffstat (limited to 'Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h')
| -rw-r--r-- | Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h index a6f006c5f..59e59644e 100644 --- a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +++ b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h @@ -508,7 +508,7 @@ ComputationInfo computeFromTridiagonal_impl(DiagType& diag, SubDiagType& subdiag typedef typename DiagType::RealScalar RealScalar; const RealScalar considerAsZero = (std::numeric_limits<RealScalar>::min)(); - const RealScalar precision_inv = RealScalar(2)/NumTraits<RealScalar>::epsilon(); + const RealScalar precision_inv = RealScalar(1)/NumTraits<RealScalar>::epsilon(); while (end>0) { for (Index i = start; i<end; ++i) { |