From 084dc63b4ccfcc9a83a12973505af74a8bc32839 Mon Sep 17 00:00:00 2001
From: Jitse Niesen <jitse@maths.leeds.ac.uk>
Date: Mon, 22 Jul 2013 13:56:15 +0100
Subject: Clean-up of MatrixSquareRoot.

---
 unsupported/test/matrix_power.cpp | 8 +++-----
 1 file changed, 3 insertions(+), 5 deletions(-)

(limited to 'unsupported/test/matrix_power.cpp')

diff --git a/unsupported/test/matrix_power.cpp b/unsupported/test/matrix_power.cpp
index 849e4287b..4c4cac509 100644
--- a/unsupported/test/matrix_power.cpp
+++ b/unsupported/test/matrix_power.cpp
@@ -100,8 +100,6 @@ template<typename MatrixType>
 void testSingular(MatrixType m, double tol)
 {
   const int IsComplex = NumTraits<typename internal::traits<MatrixType>::Scalar>::IsComplex;
-  typedef typename internal::conditional< IsComplex, MatrixSquareRootTriangular<MatrixType>,
-          MatrixSquareRootQuasiTriangular<MatrixType> >::type SquareRootType;
   typedef typename internal::conditional<IsComplex, TriangularView<MatrixType,Upper>, const MatrixType&>::type TriangularType;
   typename internal::conditional< IsComplex, ComplexSchur<MatrixType>, RealSchur<MatrixType> >::type schur;
   MatrixType T;
@@ -116,13 +114,13 @@ void testSingular(MatrixType m, double tol)
     processTriangularMatrix<MatrixType>::run(m, T, U);
     MatrixPower<MatrixType> mpow(m);
 
-    SquareRootType(T).compute(T);
+    T = T.sqrt();
     VERIFY(mpow(0.5).isApprox(U * (TriangularType(T) * U.adjoint()), tol));
 
-    SquareRootType(T).compute(T);
+    T = T.sqrt();
     VERIFY(mpow(0.25).isApprox(U * (TriangularType(T) * U.adjoint()), tol));
 
-    SquareRootType(T).compute(T);
+    T = T.sqrt();
     VERIFY(mpow(0.125).isApprox(U * (TriangularType(T) * U.adjoint()), tol));
   }
 }
-- 
cgit v1.2.3