11 #include "Tempus_StepperFactory.hpp"
12 #include "Tempus_TimeStepControl.hpp"
15 #include "../TestModels/DahlquistTestModel.hpp"
18 namespace Tempus_Unit_Test {
22 using Teuchos::rcp_const_cast;
23 using Teuchos::rcp_dynamic_cast;
25 using Teuchos::sublist;
39 TEUCHOS_TEST_FOR_EXCEPT(tscs->getReductFactor() != 0.5);
40 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMinEta() != 0.0);
41 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMaxEta() != 1.0e+16);
44 tscs->setAmplFactor(1.33); tscs->initialize(); TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
45 tscs->setReductFactor(0.75); tscs->initialize(); TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
46 tscs->setMinEta(0.01); tscs->initialize(); TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
47 tscs->setMaxEta(0.05); tscs->initialize(); TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
49 TEUCHOS_TEST_FOR_EXCEPT(tscs->getAmplFactor() != 1.33);
50 TEUCHOS_TEST_FOR_EXCEPT(tscs->getReductFactor() != 0.75);
51 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMinEta() != 0.01);
52 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMaxEta() != 0.05);
61 1.33, 0.75, 0.01, 0.05));
67 TEUCHOS_TEST_FOR_EXCEPT(tscs->getReductFactor() != 0.75);
68 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMinEta() != 0.01);
69 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMaxEta() != 0.05);
77 auto pl = Tempus::getTimeStepControlStrategyBasicVS_PL<double>();
79 pl->set<
double>(
"Amplification Factor", 1.33);
80 pl->set<
double>(
"Reduction Factor" , 0.75);
81 pl->set<
double>(
"Minimum Value Monitoring Function", 0.01);
82 pl->set<
double>(
"Maximum Value Monitoring Function", 0.05);
84 auto tscs = Tempus::createTimeStepControlStrategyBasicVS<double>(pl);
89 TEUCHOS_TEST_FOR_EXCEPT(tscs->getReductFactor() != 0.75);
90 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMinEta() != 0.01);
91 TEUCHOS_TEST_FOR_EXCEPT(tscs->getMaxEta() != 0.05);
100 tscs->setAmplFactor(1.1);
101 tscs->setReductFactor(0.5);
102 tscs->setMinEta(0.01);
103 tscs->setMaxEta(0.05);
109 tsc->setTimeStepControlStrategy(tscs);
110 tsc->setInitTime(0.0);
111 tsc->setFinalTime(10.0);
112 tsc->setMinTimeStep (0.01);
113 tsc->setInitTimeStep(0.1);
114 tsc->setMaxTimeStep (1.0);
115 tsc->setFinalIndex(100);
122 auto inArgsIC = model->getNominalValues();
124 auto icState = Tempus::createSolutionStateX<double>(icSolution);
129 solutionHistory->addState(icState);
130 solutionHistory->getCurrentState()->setTimeStep(0.5);
131 solutionHistory->getCurrentState()->setTime(0.0);
132 solutionHistory->getCurrentState()->setIndex(0);
135 for (
int i = 0; i < 2; i++) {
136 solutionHistory->initWorkingState();
138 tsc->setNextTimeStep(solutionHistory, status);
141 auto currentState = solutionHistory->getCurrentState();
142 auto workingState = solutionHistory->getWorkingState();
143 auto xN = workingState->getX();
144 Thyra::Vp_S(xN.ptr(), 1.0);
146 workingState->computeNorms(currentState);
149 solutionHistory->promoteWorkingState();
153 auto currentState = solutionHistory->getCurrentState();
160 solutionHistory->clear();
161 solutionHistory->addState(icState);
162 solutionHistory->getCurrentState()->setTimeStep(0.5);
163 solutionHistory->getCurrentState()->setTime(0.0);
164 solutionHistory->getCurrentState()->setIndex(0);
167 for (
int i = 0; i < 2; i++) {
168 solutionHistory->initWorkingState();
170 tsc->setNextTimeStep(solutionHistory, status);
173 auto currentState = solutionHistory->getCurrentState();
174 auto workingState = solutionHistory->getWorkingState();
175 auto xN = workingState->getX();
176 Thyra::Vp_S(xN.ptr(), 0.0);
178 workingState->computeNorms(currentState);
181 solutionHistory->promoteWorkingState();
186 auto currentState = solutionHistory->getCurrentState();
199 auto pl = tscs->getValidParameters();
201 TEST_COMPARE ( pl->get<std::string>(
"Strategy Type"), ==,
"Basic VS");
208 std::ostringstream unusedParameters;
209 pl->unused(unusedParameters);
StepControlStrategy class for TimeStepControl.
The classic Dahlquist Test Problem.
#define TEST_COMPARE(v1, comp, v2)
#define TEST_FLOATING_EQUALITY(v1, v2, tol)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Status
Status for the Integrator, the Stepper and the SolutionState.
TEUCHOS_UNIT_TEST(BackwardEuler, Default_Construction)
TimeStepControl manages the time step size. There several mechanisms that effect the time step size a...
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)