RK4OneStepTempl< T, N > Class Template Reference

#include <RK4OneStepTempl.h>

Public Types
typedef T	Scalar
typedef RKSmallVector< T, N >	Vector

Public Member Functions
Vector	operator() (Scalar startPar, const Vector &startState, const RKDerivative< T, N > &deriv, Scalar step) const

Detailed Description

template<typename T, int N>
class RK4OneStepTempl< T, N >

Definition at line 9 of file RK4OneStepTempl.h.

Member Typedef Documentation

template<typename T, int N>

typedef T RK4OneStepTempl< T, N >::Scalar

Definition at line 11 of file RK4OneStepTempl.h.

template<typename T, int N>

typedef RKSmallVector<T, N> RK4OneStepTempl< T, N >::Vector

Definition at line 12 of file RK4OneStepTempl.h.

Member Function Documentation

template<typename T, int N>

Vector RK4OneStepTempl< T, N >::operator() ( Scalar  startPar, const Vector &  startState, const RKDerivative< T, N > &  deriv, Scalar  step  ) const

inline

Definition at line 14 of file RK4OneStepTempl.h.

References relval_steps::k2, and mps_fire::result.

                                                                                                                   {
    // cout << "RK4OneStepTempl: starting from " << startPar << startState << endl;

    Vector k1 = step * deriv(startPar, startState);
    Vector k2 = step * deriv(startPar + step / 2, startState + k1 / 2);
    Vector k3 = step * deriv(startPar + step / 2, startState + k2 / 2);
    Vector k4 = step * deriv(startPar + step, startState + k3);

    Vector result = startState + k1 / 6 + k2 / 3 + k3 / 3 + k4 / 6;

    // cout << "RK4OneStepTempl: result for step " << step << " is " << result << endl;

    return result;
  }