12 #include "Teuchos_DefaultComm.hpp"
14 #include "Tempus_config.hpp"
15 #include "Tempus_IntegratorBasic.hpp"
16 #include "Tempus_StepperBackwardEuler.hpp"
18 #include "../TestModels/SinCosModel.hpp"
19 #include "../TestUtils/Tempus_ConvergenceTestUtils.hpp"
26 namespace Tempus_Test {
28 using Teuchos::getParametersFromXmlFile;
32 using Teuchos::rcp_const_cast;
33 using Teuchos::sublist;
44 std::vector<RCP<Thyra::VectorBase<double>>> solutions;
45 std::vector<RCP<Thyra::VectorBase<double>>> solutionsDot;
46 std::vector<double> StepSize;
47 std::vector<double> xErrorNorm;
48 std::vector<double> xDotErrorNorm;
49 const int nTimeStepSizes = 7;
52 for (
int n = 0; n < nTimeStepSizes; n++) {
55 getParametersFromXmlFile(
"Tempus_BackwardEuler_SinCos.xml");
70 pl->
sublist(
"Default Integrator")
72 .
set(
"Initial Time Step", dt);
73 integrator = Tempus::createIntegratorBasic<double>(pl, model);
80 model->getNominalValues().get_x()->clone_v();
81 integrator->initializeSolutionHistory(0.0, x0);
84 bool integratorStatus = integrator->advanceTime();
88 time = integrator->getTime();
89 double timeFinal = pl->sublist(
"Default Integrator")
90 .sublist(
"Time Step Control")
91 .
get<
double>(
"Final Time");
97 integrator->getSolutionHistory();
98 writeSolution(
"Tempus_BackwardEuler_SinCos.dat", solutionHistory);
101 for (
int i = 0; i < solutionHistory->getNumStates(); i++) {
102 double time_i = (*solutionHistory)[i]->getTime();
105 model->getExactSolution(time_i).get_x()),
107 model->getExactSolution(time_i).get_x_dot()));
108 state->setTime((*solutionHistory)[i]->getTime());
109 solnHistExact->addState(state);
111 writeSolution(
"Tempus_BackwardEuler_SinCos-Ref.dat", solnHistExact);
115 StepSize.push_back(dt);
116 auto solution = Thyra::createMember(model->get_x_space());
117 Thyra::copy(*(integrator->getX()), solution.ptr());
118 solutions.push_back(solution);
119 auto solutionDot = Thyra::createMember(model->get_x_space());
120 Thyra::copy(*(integrator->getXDot()), solutionDot.ptr());
121 solutionsDot.push_back(solutionDot);
122 if (n == nTimeStepSizes - 1) {
123 StepSize.push_back(0.0);
124 auto solutionExact = Thyra::createMember(model->get_x_space());
125 Thyra::copy(*(model->getExactSolution(time).get_x()),
126 solutionExact.ptr());
127 solutions.push_back(solutionExact);
128 auto solutionDotExact = Thyra::createMember(model->get_x_space());
129 Thyra::copy(*(model->getExactSolution(time).get_x_dot()),
130 solutionDotExact.ptr());
131 solutionsDot.push_back(solutionDotExact);
137 double xDotSlope = 0.0;
139 double order = stepper->getOrder();
140 writeOrderError(
"Tempus_BackwardEuler_SinCos-Error.dat", stepper, StepSize,
141 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.
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)
Sine-Cosine model problem from Rythmos. This is a canonical Sine-Cosine differential equation with a...
void writeSolution(const std::string filename, Teuchos::RCP< const Tempus::SolutionHistory< Scalar >> solutionHistory)
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)
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
ParameterList & sublist(const std::string &name, bool mustAlreadyExist=false, const std::string &docString="")
Solution state for integrators and steppers.