13 #include "../TestModels/DahlquistTestModel.hpp"
16 namespace Tempus_Unit_Test {
20 using Teuchos::rcp_const_cast;
21 using Teuchos::rcp_dynamic_cast;
23 using Teuchos::sublist;
60 class StepperRKModifierEDIRK_TrapezoidaTest
67 : out(Out), success(Success)
72 virtual ~StepperRKModifierEDIRK_TrapezoidaTest(){}
80 const double relTol = 1.0e-14;
81 auto stageNumber = stepper->getStageNumber();
84 auto currentState = sh->getCurrentState();
85 auto workingState = sh->getWorkingState();
86 const double dt = workingState->getTimeStep();
87 double time = currentState->getTime();
88 if (stageNumber >= 0) time += c(stageNumber)*dt;
90 auto x = workingState->getX();
91 auto xDot = workingState->getXDot();
92 if (xDot == Teuchos::null) xDot = stepper->getStepperXDot();
98 auto DME = Teuchos::rcp_dynamic_cast<
104 const double x_0 = get_ele(*(x), 0);
105 const double xDot_0 = get_ele(*(xDot), 0);
116 const double X_i = get_ele(*(x), 0);
117 const double f_i = get_ele(*(xDot), 0);
119 if (stageNumber == 0) {
123 }
else if (stageNumber == 1) {
137 const double X_i = get_ele(*(x), 0);
138 const double f_i = get_ele(*(xDot), 0);
140 if (stageNumber == 0) {
145 }
else if (stageNumber == 1) {
158 const double x_1 = get_ele(*(x), 0);
159 time = workingState->getTime();
165 if (stepper->getUseEmbedded() ==
true) {
167 TEST_ASSERT(std::abs(workingState->getTolAbs()) < relTol);
169 TEST_ASSERT(std::abs(workingState->getErrorRel()) < relTol);
194 auto modifier =
rcp(
new StepperRKModifierEDIRK_TrapezoidaTest(out, success));
196 stepper->setModel(model);
197 stepper->setAppAction(modifier);
198 stepper->setICConsistency(
"Consistent");
199 stepper->setUseFSAL(
false);
200 stepper->initialize();
206 stepper->setInitialConditions(solutionHistory);
207 solutionHistory->initWorkingState();
209 solutionHistory->getWorkingState()->setTimeStep(dt);
210 solutionHistory->getWorkingState()->setTime(dt);
211 stepper->takeStep(solutionHistory);
At the beginning of the stage.
At the beginning of the step.
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 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_COMPARE(v1, comp, v2)
#define TEST_FLOATING_EQUALITY(v1, v2, tol)
Before the implicit solve.
Sine-Cosine model problem from Rythmos. This is a canonical Sine-Cosine differential equation with a...
RK Trapezoidal Rule (A.K.A. RK Crank-Nicolson)
Teuchos::RCP< SolutionHistory< Scalar > > createSolutionHistoryME(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model)
Nonmember contructor from a Thyra ModelEvaluator.
Base class for Runge-Kutta methods, ExplicitRK, DIRK and IMEX.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
TEUCHOS_UNIT_TEST(BackwardEuler, Default_Construction)
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
ACTION_LOCATION
Indicates the location of application action (see algorithm).
void testDIRKAccessorsFullConstruction(const RCP< Tempus::StepperDIRK< double > > &stepper)
Unit test utility for ExplicitRK Stepper construction and accessors.
After the implicit solve.
Before the explicit evaluation.
#define TEUCHOS_ASSERT(assertion_test)
Base modifier for StepperRK.
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)