timestepping/newmark.hh

Go to the documentation of this file.

// newmark strategy

#ifndef newmark_hh
#define newmark_hh

#include "basics/typedefs.hh"
#include "function/vector.hh"
#include "strategy.hh"

namespace concepts {
  template<typename F>
  class Operator;
}

namespace timestepping {

  using concepts::Real;

  class TimeVector;

  // *************************************************************** Newmark **

  class Newmark : public TimeStepStrategy {
  public:
    Newmark(concepts::SolverFabric<Real>& fabric,
            concepts::Operator<Real>& D2,
            concepts::Operator<Real>& D0,
            timestepping::TimeVector& trhs,
            const concepts::Vector<Real>& Y0,
            const concepts::Vector<Real>& Z0, 
            Real dt,
            Real beta = 0.25, Real gamma = 0.5);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual void next();
  private:
    concepts::Operator<Real> &D2_, &D0_;
    timestepping::TimeVector &trhs_;
    concepts::Vector<Real> Yn1_, Yn2_;   // Yn1 = Y(t_{n-1}), Yn2 =  Y(t_{n-2})
    Real beta_, gamma_;
    std::auto_ptr<concepts::Operator<Real> > lhs_;
  };

}

#endif // newmark_hh