aboutsummaryrefslogtreecommitdiffhomepage
path: root/unsupported/test/matrix_exponential.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'unsupported/test/matrix_exponential.cpp')
-rw-r--r--unsupported/test/matrix_exponential.cpp13
1 files changed, 6 insertions, 7 deletions
diff --git a/unsupported/test/matrix_exponential.cpp b/unsupported/test/matrix_exponential.cpp
index a5b40adde..6150439c5 100644
--- a/unsupported/test/matrix_exponential.cpp
+++ b/unsupported/test/matrix_exponential.cpp
@@ -61,7 +61,7 @@ void test2dRotation(double tol)
std::cout << "test2dRotation: i = " << i << " error funm = " << relerr(C, B);
VERIFY(C.isApprox(B, static_cast<T>(tol)));
- ei_matrix_exponential(angle*A, &C);
+ C = ei_matrix_exponential(angle*A);
std::cout << " error expm = " << relerr(C, B) << "\n";
VERIFY(C.isApprox(B, static_cast<T>(tol)));
}
@@ -86,7 +86,7 @@ void test2dHyperbolicRotation(double tol)
std::cout << "test2dHyperbolicRotation: i = " << i << " error funm = " << relerr(C, B);
VERIFY(C.isApprox(B, static_cast<T>(tol)));
- ei_matrix_exponential(A, &C);
+ C = ei_matrix_exponential(A);
std::cout << " error expm = " << relerr(C, B) << "\n";
VERIFY(C.isApprox(B, static_cast<T>(tol)));
}
@@ -110,7 +110,7 @@ void testPascal(double tol)
std::cout << "testPascal: size = " << size << " error funm = " << relerr(C, B);
VERIFY(C.isApprox(B, static_cast<T>(tol)));
- ei_matrix_exponential(A, &C);
+ C = ei_matrix_exponential(A);
std::cout << " error expm = " << relerr(C, B) << "\n";
VERIFY(C.isApprox(B, static_cast<T>(tol)));
}
@@ -137,10 +137,9 @@ void randomTest(const MatrixType& m, double tol)
std::cout << "randomTest: error funm = " << relerr(identity, m2 * m3);
VERIFY(identity.isApprox(m2 * m3, static_cast<RealScalar>(tol)));
- ei_matrix_exponential(m1, &m2);
- ei_matrix_exponential(-m1, &m3);
- std::cout << " error expm = " << relerr(identity, m2 * m3) << "\n";
- VERIFY(identity.isApprox(m2 * m3, static_cast<RealScalar>(tol)));
+ m2 = ei_matrix_exponential(m1) * ei_matrix_exponential(-m1);
+ std::cout << " error expm = " << relerr(identity, m2) << "\n";
+ VERIFY(identity.isApprox(m2, static_cast<RealScalar>(tol)));
}
}