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

Matrix square root can process 0 eigenvalue.

Browse Source
This commit is contained in:
Chen-Pang He 2013-06-24 23:57:57 +08:00
parent b9fc9d8f32
commit b9543ce237
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
unsupported/Eigen/src/MatrixFunctions/MatrixSquareRoot.h
View File

@ -94,7 +94,7 @@ void MatrixSquareRootQuasiTriangular<MatrixType>::computeDiagonalPartOfSqrt(Matr
const Index size = m_A.rows(); const Index size = m_A.rows();
for (Index i = 0; i < size; i++) { for (Index i = 0; i < size; i++) {
if (i == size - 1 || T.coeff(i+1, i) == 0) { 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)); sqrtT.coeffRef(i,i) = sqrt(T.coeff(i,i));
} }
else { else {