numerics/one_step_scheme.h

// -*- c++ -*-

// +--------------------------------------------------------------------+
// | This file is part of MathTL - the Mathematical Template Library    |
// |                                                                    |
// | Copyright (c) 2002-2009                                            |
// | Thorsten Raasch, Manuel Werner                                     |
// +--------------------------------------------------------------------+

#ifndef _MATHTL_ONESTEPSCHEME_H
#define _MATHTL_ONESTEPSCHEME_H

#include <list>
#include <map>
#include <numerics/ivp.h>

namespace MathTL
{
  template <class VECTOR>
  class OneStepScheme
  {
  public:
    virtual ~OneStepScheme() = 0;

    virtual void increment(const AbstractIVP<VECTOR>* ivp,
                           const double t_m, const VECTOR& u_m,
                           const double tau,
                           VECTOR& u_mplus1,
                           VECTOR& error_estimate,
                           const double tolerance = 1e-2) const = 0;

    virtual int order() const = 0;
  };

  template <class VECTOR>
  class IVPSolution
  {
  public:
    std::list<double> t;
    std::list<VECTOR> u;
  };

  template <class VECTOR>
  void solve_IVP(const AbstractIVP<VECTOR>* ivp,
                 const OneStepScheme<VECTOR>* scheme,
                 const double T,
                 const double atol,
                 const double rtol,
                 const double q,
                 const double tau_max,
                 IVPSolution<VECTOR>& result);
}

#include <numerics/one_step_scheme.cpp>

#endif