runge_kutta.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: runge_kutta.h 1348 2013-08-16 00:45:07Z secanell $
//
//    Copyright (C) 2006, 2007 by Guido Kanschat
//
//    This file is subject to QPL and may not be  distributed
//    without copyright and license information. Please refer
//    to the file deal.II/doc/license.html for the  text  and
//    further information on this license.
//
//---------------------------------------------------------------------------

#ifndef __deal2__appframe__runge_kutta_h
#define __deal2__appframe__runge_kutta_h

#include <appframe/base.h>
#include <appframe/timestep_control.h>

namespace dealii
{
  class ParameterHandler;
  template<typename T> class BlockVector;
}

namespace AppFrame
{
  class RungeKuttaExplicit :
      public ApplicationCopy
  {
    public:
      RungeKuttaExplicit (ApplicationBase& app);
      
      virtual void solve(FEVector& start, const FEVectors& rhs);
      virtual void declare_parameters (ParameterHandler& param);
      virtual void initialize (ParameterHandler& param);
      
      void _initialize (ParameterHandler& param);

      void set_scheme(const std::string name);

      DeclException1(ExcUnknownScheme, std::string,
                     << "A scheme with name " << arg1 << " is not known");
      
    protected:
      TimestepControl control;

    private:
      void resize_vectors(unsigned int n);
      std::vector<std::vector<double> > a;
      std::vector<double> b;
      std::vector<double> c;

  };
  
}

#endif