From 39617c503faf08de1324516adc8ebb7150c218ca Mon Sep 17 00:00:00 2001 From: Tavian Barnes Date: Wed, 6 Oct 2010 17:57:48 -0400 Subject: Add RK4. --- tests/Makefile.am | 4 +++- tests/RK4.cpp | 25 +++++++++++++++++++++++++ 2 files changed, 28 insertions(+), 1 deletion(-) create mode 100644 tests/RK4.cpp (limited to 'tests') diff --git a/tests/Makefile.am b/tests/Makefile.am index 5f94aab..edef1ae 100644 --- a/tests/Makefile.am +++ b/tests/Makefile.am @@ -20,8 +20,10 @@ INCLUDES = -I$(top_srcdir)/src check_PROGRAMS = Euler-test \ - Midpoint-test + Midpoint-test \ + RK4-test TESTS = $(check_PROGRAMS) Euler_test_SOURCES = Euler.cpp Midpoint_test_SOURCES = Midpoint.cpp +RK4_test_SOURCES = RK4.cpp diff --git a/tests/RK4.cpp b/tests/RK4.cpp new file mode 100644 index 0000000..8e0ca0a --- /dev/null +++ b/tests/RK4.cpp @@ -0,0 +1,25 @@ +#include "vZ.hpp" +#include +#include +#include + +// y' = y (y == C*exp(t)) +double +f(double t, double y) +{ + return y; +} + +int +main() +{ + vZ::RK4Integrator integrator(f); + integrator.y(1.0).x(0.0).h(0.04); + + integrator.integrate(2.0); + + std::cout << integrator.y() << std::endl + << std::exp(2.0) << std::endl; + + return EXIT_SUCCESS; +} -- cgit v1.2.3