13 #include "Thyra_VectorStdOps.hpp"
15 #include "Tempus_IntegratorBasic.hpp"
16 #include "Tempus_StepperFactory.hpp"
17 #include "Tempus_StepperOperatorSplit.hpp"
18 #include "Tempus_StepperForwardEuler.hpp"
19 #include "Tempus_StepperBackwardEuler.hpp"
21 #include "../TestModels/VanDerPol_IMEX_ExplicitModel.hpp"
22 #include "../TestModels/VanDerPol_IMEX_ImplicitModel.hpp"
23 #include "../TestUtils/Tempus_ConvergenceTestUtils.hpp"
28 namespace Tempus_Test {
30 using Teuchos::getParametersFromXmlFile;
34 using Teuchos::rcp_const_cast;
35 using Teuchos::sublist;
49 getParametersFromXmlFile(
"Tempus_OperatorSplit_VanDerPol.xml");
69 stepper->addStepper(subStepper1);
70 stepper->addStepper(subStepper2);
71 stepper->initialize();
77 timeStepControl->setInitIndex(tscPL.
get<
int>(
"Initial Time Index"));
78 timeStepControl->setInitTime(tscPL.
get<
double>(
"Initial Time"));
79 timeStepControl->setFinalTime(tscPL.
get<
double>(
"Final Time"));
80 timeStepControl->setInitTimeStep(dt);
81 timeStepControl->initialize();
84 auto inArgsIC = stepper->getModel()->getNominalValues();
88 icState->setTime(timeStepControl->getInitTime());
89 icState->setIndex(timeStepControl->getInitIndex());
90 icState->setTimeStep(0.0);
91 icState->setOrder(stepper->getOrder());
96 solutionHistory->setName(
"Forward States");
98 solutionHistory->setStorageLimit(2);
99 solutionHistory->addState(icState);
103 Tempus::createIntegratorBasic<double>();
104 integrator->setStepper(stepper);
105 integrator->setTimeStepControl(timeStepControl);
106 integrator->setSolutionHistory(solutionHistory);
108 integrator->initialize();
111 bool integratorStatus = integrator->advanceTime();
115 double time = integrator->getTime();
116 double timeFinal = pl->
sublist(
"Demo Integrator")
118 .
get<
double>(
"Final Time");
125 out <<
" Stepper = " << stepper->description() << std::endl;
126 out <<
" =========================" << std::endl;
127 out <<
" Computed solution: " << get_ele(*(x), 0) <<
" "
128 << get_ele(*(x), 1) << std::endl;
129 out <<
" =========================" << std::endl;
130 TEST_FLOATING_EQUALITY(get_ele(*(x), 0), -2.223910, 1.0e-4);
131 TEST_FLOATING_EQUALITY(get_ele(*(x), 1), 0.565441, 1.0e-4);
139 std::vector<RCP<Thyra::VectorBase<double>>> solutions;
140 std::vector<RCP<Thyra::VectorBase<double>>> solutionsDot;
141 std::vector<double> StepSize;
142 std::vector<double> xErrorNorm;
143 std::vector<double> xDotErrorNorm;
144 const int nTimeStepSizes = 4;
147 for (
int n = 0; n < nTimeStepSizes; n++) {
150 getParametersFromXmlFile(
"Tempus_OperatorSplit_VanDerPol.xml");
160 std::vector<RCP<const Thyra::ModelEvaluator<double>>> models;
161 models.push_back(explicitModel);
162 models.push_back(implicitModel);
166 if (n == nTimeStepSizes - 1) dt /= 10.0;
172 .
set(
"Initial Time Step", dt);
173 integrator = Tempus::createIntegratorBasic<double>(pl, models);
176 bool integratorStatus = integrator->advanceTime();
180 time = integrator->getTime();
181 double timeFinal = pl->sublist(
"Demo Integrator")
182 .sublist(
"Time Step Control")
183 .
get<
double>(
"Final Time");
184 double tol = 100.0 * std::numeric_limits<double>::epsilon();
188 StepSize.push_back(dt);
189 auto solution = Thyra::createMember(implicitModel->get_x_space());
190 Thyra::copy(*(integrator->getX()), solution.ptr());
191 solutions.push_back(solution);
192 auto solutionDot = Thyra::createMember(implicitModel->get_x_space());
193 Thyra::copy(*(integrator->getXDot()), solutionDot.ptr());
194 solutionsDot.push_back(solutionDot);
198 if ((n == 0) || (n == nTimeStepSizes - 1)) {
199 std::string fname =
"Tempus_OperatorSplit_VanDerPol-Ref.dat";
200 if (n == 0) fname =
"Tempus_OperatorSplit_VanDerPol.dat";
202 integrator->getSolutionHistory();
210 double xDotSlope = 0.0;
212 double order = stepper->getOrder();
213 writeOrderError(
"Tempus_OperatorSplit_VanDerPol-Error.dat", stepper, StepSize,
214 solutions, xErrorNorm, xSlope, solutionsDot, xDotErrorNorm,
Teuchos::RCP< SolutionState< Scalar > > createSolutionStateX(const Teuchos::RCP< Thyra::VectorBase< Scalar > > &x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &xdot=Teuchos::null, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &xdotdot=Teuchos::null)
Nonmember constructor from non-const solution vectors, x.
van der Pol model formulated for IMEX-RK.
T & get(const std::string &name, T def_value)
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
void writeOrderError(const std::string filename, Teuchos::RCP< Tempus::Stepper< Scalar >> stepper, std::vector< Scalar > &StepSize, std::vector< Teuchos::RCP< Thyra::VectorBase< Scalar >>> &solutions, std::vector< Scalar > &xErrorNorm, Scalar &xSlope, std::vector< Teuchos::RCP< Thyra::VectorBase< Scalar >>> &solutionsDot, std::vector< Scalar > &xDotErrorNorm, Scalar &xDotSlope, std::vector< Teuchos::RCP< Thyra::VectorBase< Scalar >>> &solutionsDotDot, std::vector< Scalar > &xDotDotErrorNorm, Scalar &xDotDotSlope, Teuchos::FancyOStream &out)
#define TEST_FLOATING_EQUALITY(v1, v2, tol)
void writeSolution(const std::string filename, Teuchos::RCP< const Tempus::SolutionHistory< Scalar >> solutionHistory)
OperatorSplit stepper loops through the Stepper list.
Teuchos::RCP< StepperForwardEuler< Scalar > > createStepperForwardEuler(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model, Teuchos::RCP< Teuchos::ParameterList > pl)
Nonmember constructor - ModelEvaluator and ParameterList.
TEUCHOS_UNIT_TEST(BackwardEuler, SinCos_ASA)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
static void summarize(Ptr< const Comm< int > > comm, std::ostream &out=std::cout, const bool alwaysWriteLocal=false, const bool writeGlobalStats=true, const bool writeZeroTimers=true, const ECounterSetOp setOp=Intersection, const std::string &filter="", const bool ignoreZeroTimers=false)
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...
Keep a fix number of states.
Teuchos::RCP< StepperBackwardEuler< Scalar > > createStepperBackwardEuler(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &model, Teuchos::RCP< Teuchos::ParameterList > pl)
Nonmember constructor - ModelEvaluator and ParameterList.
ParameterList & sublist(const std::string &name, bool mustAlreadyExist=false, const std::string &docString="")
van der Pol model formulated for IMEX.
Solution state for integrators and steppers.