TAPsODESolverEuler.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsODESolverEuler.cpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (07/10/2005)
UPDATE          (08/18/2010)
******************************************************************************/
#include "TAPsODESolverEuler.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>
ODESolverEuler<T, T_SET>::ODESolverEuler () : ODESolver<T, T_SET>() {}
//-------------------------------------------------------------------
template <typename T, typename T_SET>
ODESolverEuler<T, T_SET>::~ODESolverEuler ()    {}
//-------------------------------------------------------------------
template <typename T, typename T_SET>
void ODESolverEuler<T, T_SET>::SetSize ( int vectorSize )
{ vStart.resize( vectorSize ); }
//-------------------------------------------------------------------
template <typename T, typename T_SET>
void ODESolverEuler<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;
    // Euler's method
    restartNotRequired = dxdt( t0, x0, vStart, userData );
    assert( restartNotRequired );
    // Update the state
    xEnd = x0 + h * vStart;
}
//-------------------------------------------------------------------
//template <typename T, typename T_SET>
//void ODESolverEuler<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----+----