9#ifndef Tempus_StepperTrapezoidal_decl_hpp
10#define Tempus_StepperTrapezoidal_decl_hpp
12#include "Tempus_config.hpp"
13#include "Tempus_StepperImplicit.hpp"
97 const Teuchos::RCP<Thyra::NonlinearSolverBase<Scalar> >& solver,
99 std::string ICConsistency,
100 bool ICConsistencyCheck,
101 bool zeroInitialGuess,
109 virtual Teuchos::RCP<StepperTrapezoidalAppAction<Scalar> >
getAppAction()
const
137 virtual Scalar
getAlpha(
const Scalar dt)
const {
return Scalar(2.0)/dt; }
139 virtual Scalar
getBeta (
const Scalar )
const {
return Scalar(1.0); }
143 virtual void describe(Teuchos::FancyOStream & out,
144 const Teuchos::EVerbosityLevel verbLevel)
const;
147 virtual bool isValidSetup(Teuchos::FancyOStream & out)
const;
164template <
typename Scalar>
185 xDotDot = Teuchos::null;
187 Thyra::V_StVpStV(xDot.ptr(),
s_,*x,-
s_,*
xOld_);
188 Thyra::V_VpStV (xDot.ptr(),*xDot,Scalar(-1.0),*
xDotOld_);
199 Teuchos::RCP<const Thyra::VectorBase<Scalar> >
xOld_;
200 Teuchos::RCP<const Thyra::VectorBase<Scalar> >
xDotOld_;
206template<
class Scalar>
207Teuchos::RCP<StepperTrapezoidal<Scalar> >
210 Teuchos::RCP<Teuchos::ParameterList> pl);
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
Thyra Base interface for implicit time steppers.
Application Action for StepperTrapezoidal.
Time-derivative interface for Trapezoidal method.
Teuchos::RCP< const Thyra::VectorBase< Scalar > > xOld_
Teuchos::RCP< const Thyra::VectorBase< Scalar > > xDotOld_
virtual ~StepperTrapezoidalTimeDerivative()
Destructor.
virtual void initialize(Scalar s, Teuchos::RCP< const Thyra::VectorBase< Scalar > > xOld, Teuchos::RCP< const Thyra::VectorBase< Scalar > > xDotOld)
StepperTrapezoidalTimeDerivative(Scalar s, Teuchos::RCP< const Thyra::VectorBase< Scalar > > xOld, Teuchos::RCP< const Thyra::VectorBase< Scalar > > xDotOld)
Constructor.
virtual void compute(Teuchos::RCP< const Thyra::VectorBase< Scalar > > x, Teuchos::RCP< Thyra::VectorBase< Scalar > > xDot, Teuchos::RCP< Thyra::VectorBase< Scalar > > xDotDot=Teuchos::null)
Compute the time derivative.
Trapezoidal method time stepper.
virtual bool isExplicit() const
virtual void takeStep(const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory)
Take the specified timestep, dt, and return true if successful.
virtual Scalar getAlpha(const Scalar dt) const
Return alpha = d(xDot)/dx.
virtual void setInitialConditions(const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory)
Set the initial conditions and make them consistent.
virtual bool isMultiStepMethod() const
virtual OrderODE getOrderODE() const
virtual bool isValidSetup(Teuchos::FancyOStream &out) const
virtual bool isExplicitImplicit() const
virtual Teuchos::RCP< StepperTrapezoidalAppAction< Scalar > > getAppAction() const
virtual bool isOneStepMethod() const
virtual void setAppAction(Teuchos::RCP< StepperTrapezoidalAppAction< Scalar > > appAction)
virtual void setUseFSAL(bool a)
StepperTrapezoidal()
Default constructor.
virtual Scalar getOrder() const
virtual void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
virtual bool isImplicit() const
virtual Scalar getOrderMax() const
virtual Scalar getOrderMin() const
virtual Teuchos::RCP< Tempus::StepperState< Scalar > > getDefaultStepperState()
Get a default (initial) StepperState.
virtual Scalar getBeta(const Scalar) const
Return beta = d(x)/dx.
Teuchos::RCP< StepperTrapezoidalAppAction< Scalar > > stepperTrapAppAction_
void setUseFSALTrueOnly(bool a)
This interface defines the time derivative connection between an implicit Stepper and WrapperModelEva...
Teuchos::RCP< StepperTrapezoidal< Scalar > > createStepperTrapezoidal(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model, Teuchos::RCP< Teuchos::ParameterList > pl)
Nonmember constructor - ModelEvaluator and ParameterList.
@ FIRST_ORDER_ODE
Stepper integrates first-order ODEs.