---

stepctrl.h

Go to the documentation of this file.

// FILE stepctrl.h 

#ifndef  INCLUDED_STEPCTRL_H
#define  INCLUDED_STEPCTRL_H


#include "odenum.h"
#include "datarepr.h"


class STEP_CTRL : private NAME
{
  
public:
  
  NAME::GetName;
  
  STEP_CTRL( const string & name );
  
  void SetStepSize( double h );
  
  double GetStepSize () const;

  virtual double PredictFirstStep( const PtrODENumeric ODE, 
                                 PtrDataRepr Data, 
                                 int order ) = 0;
  
  virtual bool PredictStep( double & hPred, double h, 
                            const PtrODENumeric    ODE, 
                            const INTERVAL_VECTOR  & LocErr, 
                            const INTERVAL_VECTOR  & Ytight, 
                            int order ) = 0;
  
  double CompMinStepSize( double t, double T );  
    
  double CompNextPoint( double h, double t, double T );
  
  double IfLastStep( double t, double h, double T );
  
  virtual void Init( const PtrODENumeric ODE ) = 0;  

  virtual ~STEP_CTRL() {};
  
protected:

  STEP_CTRL() { hCur = 0; }

private:

  double hCur; 
  
};


typedef STEP_CTRL* PtrStepCtrl;


inline STEP_CTRL :: STEP_CTRL( const string & name )
  : NAME( name )
{ 
  ;
}


inline void STEP_CTRL :: SetStepSize( double h ) 
{ 
  assert( h != 0 );  
  hCur = fabs(h); 
}


inline double STEP_CTRL :: GetStepSize() const 
{ 
  return hCur; 
}


#endif 

---