Leverage triangular square root in matrix log.

2026-04-10 11:34:33 +08:00 · 2011-08-25 07:42:32 +01:00
parent c01ed935dd
commit 7ee084f82f
1 changed files with 3 additions and 1 deletions
--- a/unsupported/Eigen/src/MatrixFunctions/MatrixLogarithm.h
+++ b/unsupported/Eigen/src/MatrixFunctions/MatrixLogarithm.h
@@ -158,7 +158,9 @@ void MatrixLogarithmAtomic<MatrixType>::computeBig(const MatrixType& A, MatrixTy
 	break;
      ++numberOfExtraSquareRoots;
    }
-    T = T.sqrt();
+    MatrixType sqrtT;
+    MatrixSquareRootTriangular<MatrixType>(T).compute(sqrtT);
+    T = sqrtT;
    ++numberOfSquareRoots;
  }