From 1d2bcd60d3e0cb8649ba5d2ae8672457849656ce Mon Sep 17 00:00:00 2001 From: Tavian Barnes Date: Fri, 10 Jun 2011 01:21:32 -0600 Subject: Make the test DE care about x. --- tests/Heun.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'tests/Heun.cpp') diff --git a/tests/Heun.cpp b/tests/Heun.cpp index 6eca969..f1b4364 100644 --- a/tests/Heun.cpp +++ b/tests/Heun.cpp @@ -23,11 +23,11 @@ #include #include -// y' = y (y == C*exp(x)) +// y' = x*y (y == C*exp(x^2/2)) double f(double x, double y) { - return y; + return x*y; } int @@ -49,7 +49,7 @@ main() << "Iterations: " << integrator.iterations() << std::endl; double error = std::abs(expected - actual)/std::abs(expected); - if (error > 1.4e-4 || !std::isfinite(error)) { + if (error > 2.7e-4 || !std::isfinite(error)) { std::cerr << "Error: " << 100.0*error << "%" << std::endl; return EXIT_FAILURE; } else { -- cgit v1.2.3