numeric.h

Go to the documentation of this file.

#ifndef __NUMERIC_H__
#define __NUMERIC_H__

#include <vector>
#include <exception>

class StateLengthError : public std::exception
{
public:
    StateLengthError();
    virtual ~StateLengthError() throw();
    const char* what() const throw();
};

typedef std::vector<double> StateVector;

StateVector operator+(const StateVector&, const StateVector&);  
StateVector operator-(const StateVector&, const StateVector&);  
StateVector operator*(const StateVector&, double);  
StateVector operator*(double, const StateVector&);  
StateVector operator/(const StateVector&, double);  
double abs(StateVector);    

class DiffEq
{
public:
    DiffEq();
    virtual ~DiffEq();
    

    virtual StateVector derivative(StateVector v) = 0;
};

class Integrator
{
protected:
    double stepSize;    
    double initStepSize;
public:

    Integrator(double stepSize = 0.01);
    virtual ~Integrator();
    

    virtual StateVector next(StateVector state, DiffEq* equation, double step = 0.0) = 0;
    
    void setStepSize(double);
    void resetStepSize();
    
    double getStepSize();
};

#endif