58 Teuchos::RCP<SolutionHistory<Scalar> > sh,
59 Teuchos::RCP<StepperRKBase<Scalar> > stepper,
65 const int stageNumber = stepper->getStageNumber();
66 Teuchos::SerialDenseVector<int,Scalar> c = stepper->getTableau()->c();
67 RCP<SolutionState<Scalar> > workingState = sh->getWorkingState();
68 const Scalar dt = workingState->getTimeStep();
69 Scalar time = sh->getCurrentState()->getTime();
70 if (stageNumber >= 0) time += c(stageNumber)*dt;
71 RCP<Thyra::VectorBase<Scalar> > x = workingState->getX();
107 time = workingState->getTime();
111 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
112 "Error - unknown action location.\n");
115 this->
modify(x, time, dt, stageNumber, modType);