strategy.hh

Go to the documentation of this file.

 
#ifndef timestpStrg_hh
#define timestpStrg_hh
 
#include "basics/typedefs.hh"
#include "operator/compositions.hh"
#include "function/vector.hh"
#include "operator/solverFabric.hh"
 
namespace timestepping {
 
  using concepts::Real;
  
  class TimeStepping;
 
  // ****************************************************** TimeStepStrategy **
 
  class TimeStepStrategy : public concepts::OutputOperator {
  public:
    TimeStepStrategy(concepts::SolverFabric<Real>& fabric,
                     concepts::Operator<Real>* lhs, 
                     const concepts::Vector<Real>& sol, 
                     const concepts::Vector<Real>& rhs,
                     Real dt) : 
      solver_(fabric(*lhs)), 
      liCo_(lhs),
      sol_(sol),
      rhs_(rhs),
      dt_(dt),
      t_(0) {}
  protected:
    friend class TimeStepping;
    virtual void next() = 0;
    std::unique_ptr<concepts::Operator<Real> > solver_;
    concepts::Operator<Real>* liCo_;
    concepts::Vector<Real> sol_;       // sol_ = Y(t_n)
    concepts::Vector<Real> rhs_;       // rhs_ = rhs(t_n)
    Real dt_;                          // dt_ = t_n - t_{n-1}, 
    Real t_;                           // t_ = t_n
  };
 
}
 
#endif // timestpStrg_hh