numerics/row_method.h

// -*- c++ -*-

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

#ifndef _MATHTL_ROW_METHOD_H
#define _MATHTL_ROW_METHOD_H

#include <algebra/triangular_matrix.h>
#include <algebra/vector.h>
#include <numerics/w_method.h>

namespace MathTL
{
  template <class VECTOR>
  class ROWMethod
    : public WMethod<VECTOR>, public WMethodStageEquationHelper<VECTOR>
  {
  public:
    ROWMethod(const typename WMethod<VECTOR>::Method method);

    virtual ~ROWMethod() {}

    void solve_W_stage_equation(const AbstractIVP<VECTOR>* ivp,
                                const double t,
                                const VECTOR& v,
                                const double alpha,
                                const VECTOR& y,
                                const double tolerance,
                                VECTOR& x) const
    {
      // use the exact jacobian
      ivp->solve_ROW_stage_equation(t, v, alpha, y, tolerance, x);
    }

    void approximate_ft(const AbstractIVP<VECTOR>* ivp,
                        const double t,
                        const VECTOR& v,
                        const double tolerance,
                        VECTOR& result) const
    {
      // use the exact derivative f_t
      ivp->evaluate_ft(t, v, tolerance, result);
    }
  };
}

#include <numerics/row_method.cpp>

#endif