ode.hpp

/* ode.hpp */
/* $Id$ */
#ifndef __ode_hpp
#define __ode_hpp

#include <vector>
#include <utility>

namespace ode
{
  template <typename T, typename X, typename F> class Solver1D
  {
  public:
    
    virtual ~Solver1D(){;}
        
    bool operator()(const T &t, 
            const T &dt, 
            F &f, 
            X &x,
            X &dx
            ) const 
    {return step(t, dt, f, x, dx);}
    
  protected:
    
    virtual bool step(const T &t,  
              const T &dt, 
              F &f,        
              X &x,        
              X &dx        
              ) const = 0;
  };
  
  template <typename T, typename X, typename F> class Solver
  {
  public:
    
    virtual ~Solver(){;}
    
    bool operator()(const T &t, 
            const T &dt, 
            const std::vector<F> &f, 
            std::vector<X> &x,
            std::vector<X> &dx
            ) const 
    {return step(t, dt, f, x, dx);}
        
  protected:
    
    virtual bool step(const T &t,              
              const T &dt,             
              const std::vector<F> &f, 
              std::vector<X> &x,       
              std::vector<X> &dx       
              ) const = 0;
  };
};

#endif /* ode.hpp */