test_embeddedrungekutta.hpp
f71ccdc1dc33fixpoint

test_embeddedrungekutta.hpp
View Options

			/*-------------------------------------------------------

			- Module : tests/odeint
			- File : test_embeddedrungekutta.hpp
			- Author : Fabien Georget

			Copyright (c) 2014, Fabien Georget, Princeton University

			---------------------------------------------------------*/

			#include "cxxtest/TestSuite.h"
			#include "odeint/runge_kutta_step.hpp"


			using namespace specmicp::odeint;

			void rhs1(double x, const double& yin, double& yout)
			{
			yout = -2*yin;
			}

			class TestEmbeddedRungeKutta : public CxxTest::TestSuite
			{
			public:

			void test_simpletest()
			{
			EmbeddedRungeKuttaStep45<1> algo(rhs1, CashKarp45);
			double yout, yerr;
			double tf = 0.01;
			algo.run_step(1, -2, 0, tf, yout, yerr);
			TS_ASSERT_DELTA(yout, std::exp(-2*tf), 1e-3);


			tf = 0.05;
			algo.run_step(1, -2, 0, tf, yout, yerr);
			TS_ASSERT_DELTA(yout, std::exp(-2*tf), 1e-2);

			}
			};