diff options
| author |  Jitse Niesen <jitse@maths.leeds.ac.uk> | 2010-03-01 12:05:57 +0000 |
|---|---|---|
| committer | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2010-03-01 12:05:57 +0000 |
| commit | 2d7bd1ec9124ec4e1145321626426ca7ea2e6a3b (patch) | |
| tree | 6f0ed6614dbb5aa1394c63f752ae6fcdacb5de5f /unsupported/test/matrix_exponential.cpp | |
| parent | f1f3c30ddc0e957a0165ae197d6c61b0ee9f5cf2 (diff) | |
Make MatrixFunctions tests more robust.
* Use absolute error instead of relative error.
* Test on well-conditioned matrices.
* Do not repeat the same test g_repeat times (bug fix).
* Correct diagnostic output in matrix_exponential.cpp .
Diffstat (limited to 'unsupported/test/matrix_exponential.cpp')
| -rw-r--r-- | unsupported/test/matrix_exponential.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/unsupported/test/matrix_exponential.cpp b/unsupported/test/matrix_exponential.cpp index 86e942edb..61f30334d 100644 --- a/unsupported/test/matrix_exponential.cpp +++ b/unsupported/test/matrix_exponential.cpp @@ -133,7 +133,7 @@ void randomTest(const MatrixType& m, double tol) m1 = MatrixType::Random(rows, cols); m2 = ei_matrix_function(m1, expfn) * ei_matrix_function(-m1, expfn); - std::cout << "randomTest: error funm = " << relerr(identity, m2 * m3); + std::cout << "randomTest: error funm = " << relerr(identity, m2); VERIFY(identity.isApprox(m2, static_cast<RealScalar>(tol))); m2 = ei_matrix_exponential(m1) * ei_matrix_exponential(-m1); |