DynamicalSystem.hxx

Go to the documentation of this file.

#ifndef _DynamicalSystem_
#define _DynamicalSystem_ 

#include <math.h>
#include <vector>
#include "MureHeaders.hxx" 


using namespace std; 

//-----------------------------------------------------------------------------//


      

//________________________________________________________________________

class DynamicalSystem
{
 public :

    DynamicalSystem();                                
    DynamicalSystem(const DynamicalSystem & DS);    
    virtual ~DynamicalSystem();                    
    
    int GetNumberOfEquationSize(){return fNVar;}    
    void SetNumberOfEquationSize(int n){fNVar=n;}   



    void SetPrecision(double eps=1e-5){fPrecision=eps;} 
    

    void SetForbidNegativeValue(){fIsNegativeValueAllowed=false;}
    

    void RungeKutta(double *YStart, double t1, double t2);
    

    virtual void BuildEqns(double t, double *Y, double *dYdt){} 
        
    

    
 protected :

    void RK4(double *y, double *dydx, double x, double h, double *yout);    

    void AdaptStepSize(double *y, double *dydx, double *x, double htry, double eps, double *yscal, double *hdid, double *hnext); 

    int      fNVar;      
    double  fPrecision; 
    double  fHestimate; 
    double  fHmin;      
    double  fMaxHdid;   
    double  fMinHdid;   
    bool    fIsNegativeValueAllowed; 
}; 

#endif