From 486a188b074c0c4f3ad771f536159e1bd2045fe2 Mon Sep 17 00:00:00 2001 From: Tavian Barnes Date: Wed, 6 Oct 2010 23:00:58 -0400 Subject: Add Bogacki-Shampine method. --- tests/BS23.cpp | 39 +++++++++++++++++++++++++++++++++++++++ 1 file changed, 39 insertions(+) create mode 100644 tests/BS23.cpp (limited to 'tests/BS23.cpp') diff --git a/tests/BS23.cpp b/tests/BS23.cpp new file mode 100644 index 0000000..15abd0c --- /dev/null +++ b/tests/BS23.cpp @@ -0,0 +1,39 @@ +#include "vZ.hpp" +#include +#include +#include +#include + +// y' = y (y == C*exp(t)) +double +f(double t, double y) +{ + return y; +} + +int +main() +{ + vZ::BS23Integrator integrator(f); + integrator.tol(1e-6).y(1.0).x(0.0).h(0.03); + + integrator.integrate(2.0); + + double actual = integrator.y(); + double expected = std::exp(2.0); + + std::cout << std::setprecision(10) + << "Numerical: " << actual << std::endl + << "Expected: " << expected << std::endl + << "h: " << integrator.h() << std::endl + << "rejections: " << integrator.rejections() << std::endl; + + double error = std::fabs(expected - actual)/expected; + if (error > 4.8e-6) { + std::cerr << "Error: " << 100.0*error << "%" << std::endl; + return EXIT_FAILURE; + } else { + std::cout << "Error: " << 100.0*error << "%" << std::endl; + return EXIT_SUCCESS; + } +} -- cgit v1.2.3