odesolverinterface.hh

#ifndef DUNE_FEM_SOLVER_ODESOLVERINTERFACE_HH
#define DUNE_FEM_SOLVER_ODESOLVERINTERFACE_HH
 
#include <iostream>
 
//- dune-common includes
#include <dune/common/exceptions.hh>
 
#include <dune/fem/io/parameter.hh>
#include <dune/fem/misc/mpimanager.hh>
#include <dune/fem/operator/common/spaceoperatorif.hh>
 
namespace DuneODE
{
  using namespace Dune;
  using namespace Fem;
 
  template <class DestinationImp>
  class OdeSolverInterface
  {
  protected:
    OdeSolverInterface () {}
 
    struct Monitor
    {
      double odeSolveTime_;
      double operatorTime_;
      double error_;
 
      std::size_t numberOfElements_;
 
      int newtonIterations_;
      int linearSolverIterations_;
      int maxNewtonIterations_;
      int maxLinearSolverIterations_;
      int spaceOperatorCalls_;
 
      Monitor() { reset(); }
 
      // reset all counters
      void reset()
      {
        odeSolveTime_ = 0;
        operatorTime_ = 0;
        error_ = 0;
        numberOfElements_ = 0;
        newtonIterations_ = 0;
        linearSolverIterations_ = 0;
        maxNewtonIterations_ = 0;
        maxLinearSolverIterations_ = 0;
        spaceOperatorCalls_ = 0;
      }
    };
 
  public:
    typedef Monitor MonitorType;
 
    typedef DestinationImp DestinationType;
 
    virtual ~OdeSolverInterface () {}
 
    virtual void initialize(const DestinationType& arg) = 0;
 
    virtual void solve(DestinationType& u)
    {
      MonitorType monitor;
      solve( u, monitor );
    }
 
    virtual void solve ( DestinationType &u, MonitorType &monitor ) = 0;
 
    virtual void description(std::ostream&) const = 0;
  };
 
 
} // namespace DuneODE
 
#endif // #ifndef DUNE_FEM_SOLVER_ODESOLVERINTERFACE_HH