11 #include "../TestModels/DahlquistTestModel.hpp"
14 namespace Tempus_Unit_Test {
18 using Teuchos::rcp_const_cast;
19 using Teuchos::rcp_dynamic_cast;
31 const auto rk_fe = stepper->getTableau();
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);
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);
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.
void testExplicitRKAccessorsFullConstruction(const RCP< Tempus::StepperExplicitRK< double > > &stepper)
Unit test utility for ExplicitRK Stepper construction and accessors.
void testFactoryConstruction(std::string stepperType, const Teuchos::RCP< const Thyra::ModelEvaluator< double > > &model)
Unit test utility for Stepper construction through StepperFactory.
The classic Dahlquist Test Problem.
#define TEST_FLOATING_EQUALITY(v1, v2, tol)
Sine-Cosine model problem from Rythmos. This is a canonical Sine-Cosine differential equation with a...
Teuchos::RCP< SolutionHistory< Scalar > > createSolutionHistoryME(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model)
Nonmember contructor from a Thyra ModelEvaluator.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
TEUCHOS_UNIT_TEST(BackwardEuler, Default_Construction)
Forward Euler Runge-Kutta Butcher Tableau.
#define TEUCHOS_ASSERT(assertion_test)