62 Teuchos::RCP<const Thyra::ModelEvaluator<double> >
63 model = rcp(
new Tempus_Test::DahlquistTestModel<double>(-1.0,
true));
65 stepper->setModel(model);
66 stepper->setICConsistency(
"Consistent");
67 stepper->setUseFSAL(
true);
68 stepper->initialize();
74 stepper->setInitialConditions(solutionHistory);
75 solutionHistory->initWorkingState();
77 solutionHistory->getWorkingState()->setTimeStep(dt);
78 solutionHistory->getWorkingState()->setTime(dt);
79 stepper->takeStep(solutionHistory);
82 const double relTol = 1.0e-14;
85 auto currentState = solutionHistory->getCurrentState();
86 const double x_0 = get_ele(*(currentState->getX()), 0);
87 const double xDot_0 = get_ele(*(currentState->getXDot()), 0);
88 TEST_FLOATING_EQUALITY(x_0, 1.0, relTol);
89 TEST_FLOATING_EQUALITY(xDot_0, -1.0, relTol);
90 TEST_ASSERT(std::abs(currentState->getTime()) < relTol);
93 auto workingState = solutionHistory->getWorkingState();
94 const double x_1 = get_ele(*(workingState->getX()), 0);
95 const double xDot_1 = get_ele(*(workingState->getXDot()), 0);
97 TEST_ASSERT(std::abs(x_1 ) < relTol);
98 TEST_ASSERT(std::abs(xDot_1) < relTol);
99 TEST_FLOATING_EQUALITY(workingState->getTime(), 1.0, relTol);
void testRKAppAction(const Teuchos::RCP< Tempus::StepperRKBase< double > > &stepper, const Teuchos::RCP< const Thyra::ModelEvaluator< double > > &model, Teuchos::FancyOStream &out, bool &success)
Unit test utility for Stepper RK AppAction.