Tempus  Version of the Day
Time Integration
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
Tempus::StepperDIRK_1StageTheta< Scalar > Class Template Reference

SDIRK 1 Stage Theta. More...

#include <Tempus_StepperRKButcherTableau.hpp>

Inheritance diagram for Tempus::StepperDIRK_1StageTheta< Scalar >:
Tempus::StepperDIRK< Scalar > Tempus::StepperImplicit< Scalar > Tempus::StepperRKBase< Scalar > Tempus::Stepper< Scalar > Tempus::Stepper< Scalar >

Public Member Functions

 StepperDIRK_1StageTheta ()
 Default constructor. More...
 
 StepperDIRK_1StageTheta (const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel, const Teuchos::RCP< StepperRKObserverComposite< Scalar > > &obs, const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &solver, bool useFSAL, std::string ICConsistency, bool ICConsistencyCheck, bool useEmbedded, bool zeroInitialGuess, Scalar theta=Scalar(0.5))
 
 StepperDIRK_1StageTheta (const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel, const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &solver, bool useFSAL, std::string ICConsistency, bool ICConsistencyCheck, bool useEmbedded, bool zeroInitialGuess, const Teuchos::RCP< StepperRKAppAction< Scalar > > &stepperRKAppAction, Scalar theta=Scalar(0.5))
 
void setTheta (Scalar theta)
 
Scalar getTheta ()
 
std::string getDescription () const
 
virtual bool getICConsistencyCheckDefault () const
 
Teuchos::RCP< const
Teuchos::ParameterList > 
getValidParameters () const
 
- Public Member Functions inherited from Tempus::StepperDIRK< Scalar >
std::vector< Teuchos::RCP
< Thyra::VectorBase< Scalar > > > & 
getStageXDot ()
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > & 
getXTilde ()
 
virtual Scalar getAlpha (const Scalar dt) const
 Return alpha = d(xDot)/dx. More...
 
virtual Scalar getBeta (const Scalar) const
 Return beta = d(x)/dx. More...
 
virtual bool isValidSetup (Teuchos::FancyOStream &out) const
 
virtual void setObserver (Teuchos::RCP< StepperObserver< Scalar > > obs=Teuchos::null)
 Set Observer. More...
 
virtual Teuchos::RCP
< StepperObserver< Scalar > > 
getObserver () const
 Get Observer. More...
 
virtual void initialize ()
 Initialize after construction and changing input parameters. More...
 
virtual void setInitialConditions (const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory)
 Set the initial conditions and make them consistent. More...
 
virtual void setResetInitialGuess (bool reset_guess)
 Set parameter so that the initial guess is reset at the beginning of each timestep. More...
 
virtual bool getResetInitialGuess () const
 
virtual void takeStep (const Teuchos::RCP< SolutionHistory< Scalar > > &solutionHistory)
 Take the specified timestep, dt, and return true if successful. More...
 
virtual Teuchos::RCP
< Tempus::StepperState< Scalar > > 
getDefaultStepperState ()
 Get a default (initial) StepperState. More...
 
virtual bool isExplicit () const
 
virtual bool isImplicit () const
 
virtual bool isExplicitImplicit () const
 
virtual bool isOneStepMethod () const
 
virtual bool isMultiStepMethod () const
 
virtual OrderODE getOrderODE () const
 
void getValidParametersBasicDIRK (Teuchos::RCP< Teuchos::ParameterList > pl) const
 
virtual void describe (Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
 
virtual void setUseEmbedded (bool a)
 Use embedded if avialable. More...
 
virtual bool getUseEmbedded () const
 
virtual bool getUseEmbeddedDefault () const
 
- Public Member Functions inherited from Tempus::StepperImplicit< Scalar >
virtual void setModel (const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel)
 
virtual Teuchos::RCP< const
Thyra::ModelEvaluator< Scalar > > 
getModel ()
 
virtual Teuchos::RCP< const
WrapperModelEvaluator< Scalar > > 
getWrapperModel ()
 
virtual void setDefaultSolver ()
 
virtual void setSolver (Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > solver)
 Set solver. More...
 
virtual Teuchos::RCP
< Thyra::NonlinearSolverBase
< Scalar > > 
getSolver () const
 Get solver. More...
 
const Thyra::SolveStatus< Scalar > solveImplicitODE (const Teuchos::RCP< Thyra::VectorBase< Scalar > > &x)
 Solve problem using x in-place. (Needs to be deprecated!) More...
 
const Thyra::SolveStatus< Scalar > solveImplicitODE (const Teuchos::RCP< Thyra::VectorBase< Scalar > > &x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &xDot, const Scalar time, const Teuchos::RCP< ImplicitODEParameters< Scalar > > &p)
 Solve implicit ODE, f(x, xDot, t, p) = 0. More...
 
void evaluateImplicitODE (Teuchos::RCP< Thyra::VectorBase< Scalar > > &f, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &xDot, const Scalar time, const Teuchos::RCP< ImplicitODEParameters< Scalar > > &p)
 Evaluate implicit ODE residual, f(x, xDot, t, p). More...
 
virtual void setInitialGuess (Teuchos::RCP< const Thyra::VectorBase< Scalar > > initialGuess)
 Pass initial guess to Newton solver (only relevant for implicit solvers) More...
 
virtual void setZeroInitialGuess (bool zIG)
 Set parameter so that the initial guess is set to zero (=True) or use last timestep (=False). More...
 
virtual bool getZeroInitialGuess () const
 
virtual Scalar getInitTimeStep (const Teuchos::RCP< SolutionHistory< Scalar > > &) const
 
- Public Member Functions inherited from Tempus::Stepper< Scalar >
virtual void setNonConstModel (const Teuchos::RCP< Thyra::ModelEvaluator< Scalar > > &)
 
virtual bool isInitialized ()
 True if stepper's member data is initialized. More...
 
virtual void checkInitialized ()
 Check initialization, and error out on failure. More...
 
void setStepperType (std::string s)
 
std::string getStepperType () const
 
void setUseFSAL (bool a)
 
bool getUseFSAL () const
 
virtual bool getUseFSALDefault () const
 
void setICConsistency (std::string s)
 
std::string getICConsistency () const
 
virtual std::string getICConsistencyDefault () const
 
void setICConsistencyCheck (bool c)
 
bool getICConsistencyCheck () const
 
virtual Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
getStepperX (Teuchos::RCP< SolutionState< Scalar > > state)
 Get x from SolutionState or Stepper storage. More...
 
virtual Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
getStepperXDot (Teuchos::RCP< SolutionState< Scalar > > state)
 Get xDot from SolutionState or Stepper storage. More...
 
virtual Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
getStepperXDotDot (Teuchos::RCP< SolutionState< Scalar > > state)
 Get xDotDot from SolutionState or Stepper storage. More...
 
virtual std::string description () const
 
virtual void createSubSteppers (std::vector< Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > >)
 
- Public Member Functions inherited from Tempus::StepperRKBase< Scalar >
virtual Teuchos::RCP< const
RKButcherTableau< Scalar > > 
getTableau () const
 
virtual Scalar getOrder () const
 
virtual Scalar getOrderMin () const
 
virtual Scalar getOrderMax () const
 
virtual int getNumberOfStages () const
 
virtual int getStageNumber () const
 
virtual void setStageNumber (int s)
 
virtual Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
getStageX ()
 
virtual Teuchos::RCP< const
Thyra::VectorBase< Scalar > > 
getStageX () const
 
virtual void setAppAction (Teuchos::RCP< StepperRKAppAction< Scalar > > appAction)
 
virtual Teuchos::RCP
< StepperRKAppAction< Scalar > > 
getAppAction () const
 

Protected Member Functions

void setupTableau ()
 
- Protected Member Functions inherited from Tempus::StepperDIRK< Scalar >
virtual void setupDefault ()
 Default setup for constructor. More...
 
virtual void setup (const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &wrapperModel, const Teuchos::RCP< StepperRKObserver< Scalar > > &obs, const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &solver, bool useFSAL, std::string ICConsistency, bool ICConsistencyCheck, bool useEmbedded, bool zeroInitialGuess)
 Setup for constructor. More...
 
virtual void setup (const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &wrapperModel, const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &solver, bool useFSAL, std::string ICConsistency, bool ICConsistencyCheck, bool useEmbedded, bool zeroInitialGuess, const Teuchos::RCP< StepperRKAppAction< Scalar > > &stepperRKAppAction)
 
- Protected Member Functions inherited from Tempus::Stepper< Scalar >
virtual void setStepperX (Teuchos::RCP< Thyra::VectorBase< Scalar > > x)
 Set x for Stepper storage. More...
 
virtual void setStepperXDot (Teuchos::RCP< Thyra::VectorBase< Scalar > > xDot)
 Set xDot for Stepper storage. More...
 
virtual void setStepperXDotDot (Teuchos::RCP< Thyra::VectorBase< Scalar > > xDotDot)
 Set x for Stepper storage. More...
 

Private Attributes

Scalar thetaDefault_
 
Scalar theta_
 

Additional Inherited Members

- Protected Attributes inherited from Tempus::StepperDIRK< Scalar >
std::vector< Teuchos::RCP
< Thyra::VectorBase< Scalar > > > 
stageXDot_
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
xTilde_
 
Teuchos::RCP
< StepperRKObserverComposite
< Scalar > > 
stepperObserver_
 
bool useEmbedded_
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
ee_
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
abs_u0
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
abs_u
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
sc
 
bool resetGuess_ = true
 
- Protected Attributes inherited from Tempus::StepperImplicit< Scalar >
Teuchos::RCP
< WrapperModelEvaluator
< Scalar > > 
wrapperModel_
 
Teuchos::RCP
< Thyra::NonlinearSolverBase
< Scalar > > 
solver_
 
Teuchos::RCP< const
Thyra::VectorBase< Scalar > > 
initialGuess_
 
bool zeroInitialGuess_
 
Teuchos::RCP< StepperObserver
< Scalar > > 
stepperObserver_
 
- Protected Attributes inherited from Tempus::Stepper< Scalar >
bool isInitialized_ = false
 True if stepper's member data is initialized. More...
 
- Protected Attributes inherited from Tempus::StepperRKBase< Scalar >
Teuchos::RCP< RKButcherTableau
< Scalar > > 
tableau_
 
int stageNumber_
 The current Runge-Kutta stage number, {0,...,s-1}. -1 indicates outside stage loop. More...
 
Teuchos::RCP
< Thyra::VectorBase< Scalar > > 
stageX_
 
Teuchos::RCP
< StepperRKAppAction< Scalar > > 
stepperRKAppAction_
 

Detailed Description

template<class Scalar>
class Tempus::StepperDIRK_1StageTheta< Scalar >

SDIRK 1 Stage Theta.

The tableau (order = 1 or 2) is

\[ \begin{array}{c|c} c & A \\ \hline & b^T \end{array} \;\;\;\;\mbox{ where }\;\;\;\; \begin{array}{c|c} \theta & \theta \\ & 1 \end{array} \]

Valid values are $0 < \theta \leq 1$, where $\theta$ = 0 implies Forward Euler (not avialble with this stepepr as it makes it explicit), $theta$ = 1/2 implies implicit midpoint method (default), and $theta$ = 1 implies Backward Euler. For $theta$ != 1/2, this method is first-order accurate, and with $theta$ = 1/2, it is second-order accurate. This method is A-stable, but becomes L-stable with $theta$ = 1. (A.K.A. Generalized Implicit Midpoint Method.)

Reference: Non-standard finite-difference methods in dynamical systems, P. Kama, Dissertation, University of Pretoria, pg. 49.

Definition at line 2654 of file Tempus_StepperRKButcherTableau.hpp.

Constructor & Destructor Documentation

template<class Scalar >
Tempus::StepperDIRK_1StageTheta< Scalar >::StepperDIRK_1StageTheta ( )
inline

Default constructor.

Requires subsequent setModel() and initialize() calls before calling takestep().

Definition at line 2663 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
Tempus::StepperDIRK_1StageTheta< Scalar >::StepperDIRK_1StageTheta ( const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &  appModel,
const Teuchos::RCP< StepperRKObserverComposite< Scalar > > &  obs,
const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &  solver,
bool  useFSAL,
std::string  ICConsistency,
bool  ICConsistencyCheck,
bool  useEmbedded,
bool  zeroInitialGuess,
Scalar  theta = Scalar(0.5) 
)
inline

Definition at line 2675 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
Tempus::StepperDIRK_1StageTheta< Scalar >::StepperDIRK_1StageTheta ( const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &  appModel,
const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &  solver,
bool  useFSAL,
std::string  ICConsistency,
bool  ICConsistencyCheck,
bool  useEmbedded,
bool  zeroInitialGuess,
const Teuchos::RCP< StepperRKAppAction< Scalar > > &  stepperRKAppAction,
Scalar  theta = Scalar(0.5) 
)
inline

Definition at line 2696 of file Tempus_StepperRKButcherTableau.hpp.

Member Function Documentation

template<class Scalar >
std::string Tempus::StepperDIRK_1StageTheta< Scalar >::getDescription ( ) const
inlinevirtual
template<class Scalar >
virtual bool Tempus::StepperDIRK_1StageTheta< Scalar >::getICConsistencyCheckDefault ( ) const
inlinevirtual

Reimplemented from Tempus::Stepper< Scalar >.

Definition at line 2744 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
Scalar Tempus::StepperDIRK_1StageTheta< Scalar >::getTheta ( )
inline

Definition at line 2728 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
Teuchos::RCP<const Teuchos::ParameterList> Tempus::StepperDIRK_1StageTheta< Scalar >::getValidParameters ( ) const
inlinevirtual

Reimplemented from Tempus::StepperDIRK< Scalar >.

Definition at line 2747 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
void Tempus::StepperDIRK_1StageTheta< Scalar >::setTheta ( Scalar  theta)
inline

Definition at line 2717 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
void Tempus::StepperDIRK_1StageTheta< Scalar >::setupTableau ( )
inlineprotectedvirtual

Member Data Documentation

template<class Scalar >
Scalar Tempus::StepperDIRK_1StageTheta< Scalar >::theta_
private

Definition at line 2787 of file Tempus_StepperRKButcherTableau.hpp.

template<class Scalar >
Scalar Tempus::StepperDIRK_1StageTheta< Scalar >::thetaDefault_
private

Definition at line 2786 of file Tempus_StepperRKButcherTableau.hpp.


The documentation for this class was generated from the following file: