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());