TAPsODESolverMidpoint.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsODESolverMidpoint.cpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (07/10/2005)
UPDATE          (08/18/2010)
******************************************************************************/
#include "TAPsODESolverMidpoint.hpp"
// Using Inclusion Model (i.e. definitions are included in declarations)
//                       (this name.cpp is included in name.hpp)
// Each friend is defined directly inside its declaration.


BEGIN_NAMESPACE_TAPs__Simulation
//=============================================================================
//-----------------------------------------------------------------------------
template <typename T, typename T_SET>
ODESolverMidpoint<T, T_SET>::ODESolverMidpoint () : ODESolver<T, T_SET>()   {}
//-------------------------------------------------------------------
template <typename T, typename T_SET>
ODESolverMidpoint<T, T_SET>::~ODESolverMidpoint ()  {}
//-------------------------------------------------------------------
template <typename T, typename T_SET>
void ODESolverMidpoint<T, T_SET>::SetSize ( int vectorSize )
{ 
    vStart.resize( vectorSize );
    vMid.resize( vectorSize );
}
//-------------------------------------------------------------------
template <typename T, typename T_SET>
void ODESolverMidpoint<T, T_SET>::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
)
{
    bool restartNotRequired = true;
    T h = t1 - t0;
    // Midpoint Method
    restartNotRequired = dxdt( t0, x0, vStart, userData );
    assert( restartNotRequired );
    restartNotRequired = dxdt( t0 + h/2, x0 + h/2*vStart, vMid, userData );
    assert( restartNotRequired );
    // Update the state
    xEnd = x0 + h * vMid;
}
//-------------------------------------------------------------------
//template <typename T, typename T_SET>
//void ODESolverMidpoint<T, T_SET>::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
//)
//{
//}
//-----------------------------------------------------------------------------
//=============================================================================
END_NAMESPACE_TAPs__Simulation
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----