ODEIntegrators.h

Go to the documentation of this file.

#ifndef ODEINTEGRATORS_H
#define ODEINTEGRATORS_H
 
#include <vector>
 
namespace ODEIntegrators {
 
class AbstractODEIntegrator
{
   public:
    virtual std::vector<double> get_initial_array() const = 0;
 
    virtual void pre_step_callback() = 0;
 
    virtual void post_deriv_callback() = 0;
 
    virtual void post_step_callback(double t, double h, std::vector<double>& x) = 0;
 
    virtual bool premature_termination() = 0;
 
    virtual void derivs(double t, std::vector<double>& x, std::vector<double>& f) = 0;
};
 
bool AdaptiveRK54(AbstractODEIntegrator& ode, double tmin, double tmax, double hmin, double hmax, double eps_allowed, double step_relax);
 
}  // namespace ODEIntegrators
 
#endif