numerics/ivp.h

// -*- c++ -*-

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

#ifndef _MATHTL_IVP_H
#define _MATHTL_IVP_H

#include <geometry/point.h>

namespace MathTL
{
  template <unsigned int DIM>
  class IVP
  {
  public:
    Point<DIM> u0;

    virtual ~IVP ();
    
    virtual void apply_f(const double t, const Point<DIM>& v,
                         Point<DIM>& result) const = 0;

    virtual void apply_ft(const double t, const Point<DIM>& v,
                          Point<DIM>& result) const = 0;

    virtual void solve_jacobian(const double t, const Point<DIM>& v, const double alpha,
                                Point<DIM>& result) const = 0;
  };

  template <class VECTOR>
  class AbstractIVP
  {
  public:
    VECTOR u0;

    virtual ~AbstractIVP() = 0;

    virtual void evaluate_f(const double t,
                            const VECTOR& v,
                            const double tolerance,
                            VECTOR& result) const = 0;

    virtual void evaluate_ft(const double t,
                             const VECTOR& v,
                             const double tolerance,
                             VECTOR& result) const = 0;

    virtual void solve_ROW_stage_equation(const double t,
                                          const VECTOR& v,
                                          const double alpha,
                                          const VECTOR& y,
                                          const double tolerancs,
                                          VECTOR& result) const = 0;
  };
}

#include <numerics/ivp.cpp>

#endif