diff options
author | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2016-09-16 14:24:47 +0200 |
---|---|---|
committer | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2016-09-16 14:24:47 +0200 |
commit | ce3557ca69742af477546d031d644a6dab1ff614 (patch) | |
tree | 7f099d3c4dc0dfd7b012af7a51c6a2d6ab05e72b /Eigen/src | |
parent | 6edd2e2851ca9080a43afa5ee64031a92750efdd (diff) |
Make makeHouseholder more stable for cases where real(c0) is not very small (but the rest is).
Diffstat (limited to 'Eigen/src')
-rw-r--r-- | Eigen/src/Householder/Householder.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/Householder/Householder.h b/Eigen/src/Householder/Householder.h index 4c1f499a1..31f804a99 100644 --- a/Eigen/src/Householder/Householder.h +++ b/Eigen/src/Householder/Householder.h @@ -77,7 +77,7 @@ void MatrixBase<Derived>::makeHouseholder( Scalar c0 = coeff(0); const RealScalar tol = (std::numeric_limits<RealScalar>::min)(); - if(tailSqNorm <= tol && numext::abs2(numext::imag(c0))<=tol) + if(tailSqNorm <= tol && numext::abs2(c0)<=tol) { tau = RealScalar(0); beta = numext::real(c0); |