aboutsummaryrefslogtreecommitdiffhomepage
path: root/unsupported/test/matrix_power.cpp
diff options
context:
space:
mode:
authorGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2013-07-22 13:56:15 +0100
committerGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2013-07-22 13:56:15 +0100
commit084dc63b4ccfcc9a83a12973505af74a8bc32839 (patch)
tree75b3adbabc05897ff38f446e1c6177e273d8b767 /unsupported/test/matrix_power.cpp
parent463343fb37cbd76e2527c560557a4d27a3d29c80 (diff)
Clean-up of MatrixSquareRoot.
Diffstat (limited to 'unsupported/test/matrix_power.cpp')
-rw-r--r--unsupported/test/matrix_power.cpp8
1 files changed, 3 insertions, 5 deletions
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));
}
}