summaryrefslogtreecommitdiffstats
path: root/tests/RK4.cpp
diff options
context:
space:
mode:
authorTavian Barnes <tavianator@gmail.com>2011-06-10 01:21:32 -0600
committerTavian Barnes <tavianator@gmail.com>2011-06-10 01:21:32 -0600
commit1d2bcd60d3e0cb8649ba5d2ae8672457849656ce (patch)
tree0fdec95e8d38eab2bd922dfdada6731823e5a37b /tests/RK4.cpp
parentd03ab99aa999544403652e0739c9bea3b8b7835e (diff)
downloadvz-1d2bcd60d3e0cb8649ba5d2ae8672457849656ce.tar.xz
Make the test DE care about x.
Diffstat (limited to 'tests/RK4.cpp')
-rw-r--r--tests/RK4.cpp6
1 files changed, 3 insertions, 3 deletions
diff --git a/tests/RK4.cpp b/tests/RK4.cpp
index 9edeb1c..f442a71 100644
--- a/tests/RK4.cpp
+++ b/tests/RK4.cpp
@@ -23,11 +23,11 @@
#include <iostream>
#include <iomanip>
-// 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 > 4.2e-8 || !std::isfinite(error)) {
+ if (error > 2.2e-7 || !std::isfinite(error)) {
std::cerr << "Error: " << 100.0*error << "%" << std::endl;
return EXIT_FAILURE;
} else {