9#ifndef Tempus_PhysicsStateTest_StepperForwardEuler_hpp
10#define Tempus_PhysicsStateTest_StepperForwardEuler_hpp
12#include "Tempus_config.hpp"
13#include "Tempus_StepperForwardEuler.hpp"
44 {
return Teuchos::null; }
56 TEMPUS_FUNC_TIME_MONITOR(
"Tempus::StepperPhysicsStateTest::takeStep()");
58 RCP<Tempus::SolutionState<Scalar> > currentState =
59 solutionHistory->getCurrentState();
61 typedef Thyra::ModelEvaluatorBase MEB;
62 this->
inArgs_.set_x(currentState->getX());
63 if (this->
inArgs_.supports(MEB::IN_ARG_t))
64 this->
inArgs_.set_t(currentState->getTime());
71 if (this->
inArgs_.supports(MEB::IN_ARG_x_dot))
72 this->
inArgs_.set_x_dot(Teuchos::null);
73 this->
outArgs_.set_f(currentState->getXDot());
78 RCP<Tempus::SolutionState<Scalar> > workingState =
79 solutionHistory->getWorkingState();
80 const Scalar dt = workingState->getTimeStep();
81 Thyra::V_VpStV(Teuchos::outArg(*(workingState->getX())),
82 *(currentState->getX()),dt,*(currentState->getXDot()));
84 RCP<PhysicsStateCounter<Scalar> > pSC =
85 Teuchos::rcp_dynamic_cast<PhysicsStateCounter<Scalar> >
86 (workingState->getPhysicsState());
87 int counter = pSC->getCounter();
89 pSC->setCounter(counter);
92 workingState->setOrder(this->
getOrder());
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
Thyra Base interface for implicit time steppers.
Thyra::ModelEvaluatorBase::InArgs< Scalar > inArgs_
Thyra::ModelEvaluatorBase::OutArgs< Scalar > outArgs_
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > appModel_
Explicit ODE ModelEvaluator.
virtual void setModel(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel)
Set model.
void setICConsistencyCheck(bool c)
virtual std::string description() const
virtual void setUseFSAL(bool a)
void setStepperType(std::string s)
Set the stepper type.
void setICConsistency(std::string s)
This is a Forward Euler time stepper to test the PhysicsState.
Scalar getOrderMax() const
StepperPhysicsStateTest(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel)
Constructor.
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Scalar getOrderMin() const
Teuchos::RCP< Tempus::StepperState< Scalar > > getDefaultStepperState()
virtual void takeStep(const Teuchos::RCP< Tempus::SolutionHistory< Scalar > > &solutionHistory)
Take the specified timestep, dt, and return true if successful.
Tempus::OrderODE getOrderODE() const
void initialize()
Initialize after construction and changing input parameters.
@ FIRST_ORDER_ODE
Stepper integrates first-order ODEs.