11#include "Tempus_StepperFactory.hpp"
12#include "Tempus_TimeStepControl.hpp"
15#include "../TestModels/DahlquistTestModel.hpp"
22using Teuchos::rcp_const_cast;
23using Teuchos::rcp_dynamic_cast;
24using Teuchos::ParameterList;
25using Teuchos::sublist;
33 TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
36 TEUCHOS_TEST_FOR_EXCEPT(tscs->getStrategyType() !=
"Basic VS");
37 TEUCHOS_TEST_FOR_EXCEPT(tscs->getStepType() !=
"Variable");
38 TEUCHOS_TEST_FOR_EXCEPT(tscs->getAmplFactor() != 1.75);
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));
62 TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
64 TEUCHOS_TEST_FOR_EXCEPT(tscs->getStrategyType() !=
"Basic VS");
65 TEUCHOS_TEST_FOR_EXCEPT(tscs->getStepType() !=
"Variable");
66 TEUCHOS_TEST_FOR_EXCEPT(tscs->getAmplFactor() != 1.33);
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);
86 TEUCHOS_TEST_FOR_EXCEPT(tscs->getStrategyType() !=
"Basic VS");
87 TEUCHOS_TEST_FOR_EXCEPT(tscs->getStepType() !=
"Variable");
88 TEUCHOS_TEST_FOR_EXCEPT(tscs->getAmplFactor() != 1.33);
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);
105 TEUCHOS_TEST_FOR_EXCEPT(!tscs->isInitialized());
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);
117 TEUCHOS_TEST_FOR_EXCEPT(!tsc->isInitialized());
122 auto inArgsIC = model->getNominalValues();
123 auto icSolution = rcp_const_cast<Thyra::VectorBase<double> >(inArgsIC.get_x());
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();
154 TEST_FLOATING_EQUALITY( currentState->getTimeStep(), 0.25, 1.0e-14);
155 TEST_FLOATING_EQUALITY( currentState->getTime(), 0.75, 1.0e-14);
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();
187 TEST_FLOATING_EQUALITY( currentState->getTimeStep(), 0.55, 1.0e-14);
188 TEST_FLOATING_EQUALITY( currentState->getTime(), 1.05, 1.0e-14);
199 auto pl = tscs->getValidParameters();
201 TEST_COMPARE ( pl->get<std::string>(
"Strategy Type"), ==,
"Basic VS");
202 TEST_FLOATING_EQUALITY( pl->get<
double>(
"Amplification Factor"), 1.75, 1.0e-14);
203 TEST_FLOATING_EQUALITY( pl->get<
double>(
"Reduction Factor"), 0.5, 1.0e-14);
204 TEST_FLOATING_EQUALITY( pl->get<
double>(
"Minimum Value Monitoring Function"), 0.0, 1.0e-14);
205 TEST_FLOATING_EQUALITY( pl->get<
double>(
"Maximum Value Monitoring Function"), 1.0e+16, 1.0e-14);
208 std::ostringstream unusedParameters;
209 pl->unused(unusedParameters);
210 TEST_COMPARE ( unusedParameters.str(), ==,
"");
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
StepControlStrategy class for TimeStepControl.
TimeStepControl manages the time step size. There several mechanisms that effect the time step size a...
The classic Dahlquist Test Problem.
TEUCHOS_UNIT_TEST(BackwardEuler, Default_Construction)
Status
Status for the Integrator, the Stepper and the SolutionState.