From 7ee084f82febace669fae5334a13d3465aeac7d4 Mon Sep 17 00:00:00 2001
From: Jitse Niesen <jitse@maths.leeds.ac.uk>
Date: Thu, 25 Aug 2011 07:42:32 +0100
Subject: Leverage triangular square root in matrix log.

---
 unsupported/Eigen/src/MatrixFunctions/MatrixLogarithm.h | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

(limited to 'unsupported/Eigen/src/MatrixFunctions/MatrixLogarithm.h')

diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixLogarithm.h b/unsupported/Eigen/src/MatrixFunctions/MatrixLogarithm.h
index 4500483fb..90afb59ff 100644
--- 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;
   }
 
-- 
cgit v1.2.3