From 1f25636c64f0a7c07fa268c0d66c084948b2dd15 Mon Sep 17 00:00:00 2001 From: Tavian Barnes Date: Thu, 7 Oct 2010 01:47:25 -0400 Subject: Fix order parameter for RKF45. --- src/vZ/RKF45.hpp | 2 +- tests/RKF45.cpp | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/vZ/RKF45.hpp b/src/vZ/RKF45.hpp index 5d40b88..5c558b7 100644 --- a/src/vZ/RKF45.hpp +++ b/src/vZ/RKF45.hpp @@ -45,7 +45,7 @@ namespace vZ typedef typename GenericAdaptiveIntegrator::Function Function; GenericRKF45Integrator(Function f) - : GenericAdaptiveIntegrator(f, 3, s_a, s_b, s_bStar) { } + : GenericAdaptiveIntegrator(f, 5, s_a, s_b, s_bStar) { } ~GenericRKF45Integrator() { } private: diff --git a/tests/RKF45.cpp b/tests/RKF45.cpp index e2764cf..a3cd259 100644 --- a/tests/RKF45.cpp +++ b/tests/RKF45.cpp @@ -30,7 +30,7 @@ main() << "Rejections: " << integrator.rejections() << std::endl; double error = std::fabs(expected - actual)/expected; - if (error > 1.5e-6) { + if (error > 1.7e-6) { std::cerr << "Error: " << 100.0*error << "%" << std::endl; return EXIT_FAILURE; } else { -- cgit v1.2.3