			#ifndef __SCHEME_H__
			#define __SCHEME_H__


			// runge kutta order 4
			// x = position
			// v = speed
			// t = time
			// a = acceleration function, owned by a Base class
			// returns the delta of position dx
			//
			template<class T, class Base>
			T runge_kutta4(const T &args, double dt, double t, T (Base::f)(const T&, double) const, Base base)
			{
			#define ACC(x, t) (base->*f)(x, t)

			T k1 = dt * ACC(args, t);
			T k2 = dt * ACC(args + k1 / 2, t + dt / 2);
			T k3 = dt * ACC(args + k2 / 2, t + dt / 2);
			T k4 = dt * ACC(args + k3, t + dt);

			return args + (k1 + 2 * k2 + 2 * k3 + k4) / 6;
			}

			#endif

rSTICAZZI/physnum/rap4/include/scheme.h
596b593deee7master

scheme.h
View Options