---

vodesolv.h

Go to the documentation of this file.

// FILE vodesolv.h 


#ifndef  INCLUDED_VODESOLV_H
#define  INCLUDED_VODESOLV_H


#include "orderctrl.h"
#include "initencl.h"
#include "tgtencl.h"


class   VODE_SOLVER;

typedef VODE_SOLVER* PtrVODESolver;


class VODE_SOLVER : public ODE_SOLVER
{
public:
  
  VODE_SOLVER( const PtrODENumeric  ODE , 
               const PtrDataRepr    Data, 
               const PtrOrderCtrl   Order,
               const PtrStepCtrl    Step,
               const PtrInitEncl    Init,
               const PtrTightEncl   Tight );
  
  
  void IntegrateTo( double tEnd, bool Continue = false );
  
  void SetOrderCtrl ( const PtrOrderCtrl Order );  
  
  void SetInitEncl  ( const PtrInitEncl  Init  );  
  
  void SetTightEncl ( const PtrTightEncl Tight );  

  PtrTightEncl GetTightEncl () const;
      
  virtual ~VODE_SOLVER();   

  
protected:
  
  void InitSolver ();


protected:
  
  PtrOrderCtrl    OrderControl; 
  
  PtrInitEncl     InitEncl;   
  
  PtrTightEncl    TightEncl;    
  
  INTERVAL_VECTOR LocalError;
  
  INTERVAL_VECTOR Ytight;
  
  INTERVAL_VECTOR Yinit;
  
private:
  
  void ResizeSpace( int n );
  
};   



inline VODE_SOLVER :: VODE_SOLVER( const PtrODENumeric  ODE, 
                                   const PtrDataRepr    Data, 
                                   const PtrOrderCtrl   Order, 
                                   const PtrStepCtrl    Step, 
                                   const PtrInitEncl    Init, 
                                   const PtrTightEncl   Tight )
  
  : ODE_SOLVER(ODE,Data,Step),  OrderControl(Order),  
    InitEncl(Init), TightEncl(Tight)
{
  assert( NotNull(Order) );
  assert( NotNull(Init)  );
  assert( NotNull(Tight) );
}


inline void VODE_SOLVER :: SetOrderCtrl ( const PtrOrderCtrl  Order )
{ 
  assert(Order);
  DELETE(OrderControl);
  OrderControl = Order; 
}


inline void VODE_SOLVER :: SetInitEncl  ( const PtrInitEncl Init )
{ 
  assert(Init);
  DELETE(InitEncl);
  InitEncl = Init; 
}


inline void VODE_SOLVER :: SetTightEncl ( const PtrTightEncl Tight ) 
{ 
  assert(Tight);
  DELETE(TightEncl);
  TightEncl = Tight; 
}

inline PtrTightEncl VODE_SOLVER :: GetTightEncl () const
{
  assert(NotNull(TightEncl));
  return TightEncl;
}


#endif

---