TAPsODESolverRungeKutta4.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsODESolverRungeKutta4.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (07/10/2005)
UPDATE          (08/18/2010)
******************************************************************************/
#ifndef TAPs_ODE_SOLVER_RUNGE_KUTTA_4_HPP
#define TAPs_ODE_SOLVER_RUNGE_KUTTA_4_HPP

#include "TAPsODESolver.hpp"

BEGIN_NAMESPACE_TAPs__Simulation
//=============================================================================
//-----------------------------------------------------------------------------
template <typename T, typename T_SET = T>
class ODESolverRungeKutta4 : public ODESolver<T, T_SET>
{
public:

    ODESolverRungeKutta4 ();

    virtual ~ODESolverRungeKutta4 ();

    virtual void SetSize ( int vectorSize );    // set data size

    virtual void Solve ( 
        VectorSet<T_SET> &x0,   // i/p: x(t_0) an initial state vector
        VectorSet<T_SET> &xEnd, // o/p: x(t_1) an end state vector
        T t0, T t1,             // i/p: starting and ending times
        DerivFunc dxdt,         // i/p: a derivative function
        void *userData          // o/p: array of user data
    );

    //virtual void Solve_ForPtrData ( 
    //  VectorSetPtr<T_SET,T> &x0,      // i/p: x(t_0) an initial state vector
    //  VectorSetPtr<T_SET,T> &xEnd,    // o/p: x(t_1) an end state vector
    //  T t0, T t1,                     // i/p: starting and ending times
    //  DerivFunc_ForPtrData dxdt,      // i/p: a derivative function
    //  void *userData                  // o/p: array of user data
    //);

private:
    // temp VectorSets used by the RungeKutta4 ode solver
    VectorSet<T_SET> vK1, vK2, vK3, vK4;
};
//-----------------------------------------------------------------------------
//=============================================================================
END_NAMESPACE_TAPs__Simulation
//-----------------------------------------------------------------------------
// Include definition if TAPs_USE_EXPORT is not defined
//#if !defined( TAPs_USE_EXPORT )
    #include "TAPsODESolverRungeKutta4.cpp"
//#endif
//-----------------------------------------------------------------------------
#endif
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----