diff options
author | Gael Guennebaud <g.gael@free.fr> | 2016-04-13 22:49:51 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2016-04-13 22:49:51 +0200 |
commit | feef39e2d1bfd2a703ff9125b60e899802f0c3d9 (patch) | |
tree | 4b2bb42fbf657966250a9d557fb422fa43f4d9ec /Eigen | |
parent | f4e12272f15ecde2751d1be615527925b19dee0e (diff) |
Fix underflow in JacoviSVD's complex to real preconditioner
Diffstat (limited to 'Eigen')
-rw-r--r-- | Eigen/src/SVD/JacobiSVD.h | 12 |
1 files changed, 9 insertions, 3 deletions
diff --git a/Eigen/src/SVD/JacobiSVD.h b/Eigen/src/SVD/JacobiSVD.h index bf5ff48c3..88bc0688e 100644 --- a/Eigen/src/SVD/JacobiSVD.h +++ b/Eigen/src/SVD/JacobiSVD.h @@ -368,9 +368,15 @@ struct svd_precondition_2x2_block_to_be_real<MatrixType, QRPreconditioner, true> if(n==0) { - z = abs(work_matrix.coeff(p,q)) / work_matrix.coeff(p,q); - work_matrix.row(p) *= z; - if(svd.computeU()) svd.m_matrixU.col(p) *= conj(z); + // make sure firt column is zero (deflation) + work_matrix.coeffRef(p,p) = work_matrix.coeffRef(q,p) = Scalar(0); + if(work_matrix.coeff(p,q)!=Scalar(0)) + { + // work_matrix.coeff(p,q) can be zero if work_matrix.coeff(q,p) is not zero but small enough to underflow when computing n + z = abs(work_matrix.coeff(p,q)) / work_matrix.coeff(p,q); + work_matrix.row(p) *= z; + if(svd.computeU()) svd.m_matrixU.col(p) *= conj(z); + } if(work_matrix.coeff(q,q)!=Scalar(0)) { z = abs(work_matrix.coeff(q,q)) / work_matrix.coeff(q,q); |