diff options
| author |  Chen-Pang He <jdh8@ms63.hinet.net> | 2013-06-24 23:57:57 +0800 |
|---|---|---|
| committer | Chen-Pang He <jdh8@ms63.hinet.net> | 2013-06-24 23:57:57 +0800 |
| commit | b9543ce237ec4ca346f9088fadfb2c75bc40c9f1 (patch) | |
| tree | 838564542f52c23c22f4df3566ef2ce15d2c88e9 /unsupported | |
| parent | b9fc9d8f32749b86bcd7a9b65bd0859e570976a3 (diff) | |
Matrix square root can process 0 eigenvalue.
Diffstat (limited to 'unsupported')
| -rw-r--r-- | unsupported/Eigen/src/MatrixFunctions/MatrixSquareRoot.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixSquareRoot.h b/unsupported/Eigen/src/MatrixFunctions/MatrixSquareRoot.h index c1a04036b..b48ea9d46 100644 --- a/unsupported/Eigen/src/MatrixFunctions/MatrixSquareRoot.h +++ b/unsupported/Eigen/src/MatrixFunctions/MatrixSquareRoot.h @@ -94,7 +94,7 @@ void MatrixSquareRootQuasiTriangular<MatrixType>::computeDiagonalPartOfSqrt(Matr const Index size = m_A.rows(); for (Index i = 0; i < size; i++) { if (i == size - 1 || T.coeff(i+1, i) == 0) { - eigen_assert(T(i,i) > 0); + eigen_assert(T(i,i) >= 0); sqrtT.coeffRef(i,i) = sqrt(T.coeff(i,i)); } else { |