12 #include "Teuchos_DefaultComm.hpp"
14 #include "Tempus_config.hpp"
15 #include "Tempus_IntegratorBasic.hpp"
16 #include "Tempus_StepperBackwardEuler.hpp"
18 #include "../TestModels/VanDerPolModel.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 = 4;
51 for (
int n = 0; n < nTimeStepSizes; n++) {
54 getParametersFromXmlFile(
"Tempus_BackwardEuler_VanDerPol.xml");
62 if (n == nTimeStepSizes - 1) dt /= 10.0;
68 .
set(
"Initial Time Step", dt);
69 integrator = Tempus::createIntegratorBasic<double>(pl, model);
72 bool integratorStatus = integrator->advanceTime();
76 double time = integrator->getTime();
77 double timeFinal = pl->sublist(
"Demo Integrator")
78 .sublist(
"Time Step Control")
79 .
get<
double>(
"Final Time");
80 double tol = 100.0 * std::numeric_limits<double>::epsilon();
84 StepSize.push_back(dt);
85 auto solution = Thyra::createMember(model->get_x_space());
86 Thyra::copy(*(integrator->getX()), solution.ptr());
87 solutions.push_back(solution);
88 auto solutionDot = Thyra::createMember(model->get_x_space());
89 Thyra::copy(*(integrator->getXDot()), solutionDot.ptr());
90 solutionsDot.push_back(solutionDot);
94 if ((n == 0) || (n == nTimeStepSizes - 1)) {
95 std::string fname =
"Tempus_BackwardEuler_VanDerPol-Ref.dat";
96 if (n == 0) fname =
"Tempus_BackwardEuler_VanDerPol.dat";
98 integrator->getSolutionHistory();
105 double xDotSlope = 0.0;
107 double order = stepper->getOrder();
108 writeOrderError(
"Tempus_BackwardEuler_VanDerPol-Error.dat", stepper, StepSize,
109 solutions, xErrorNorm, xSlope, solutionsDot, xDotErrorNorm,
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)
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...
van der Pol model problem for nonlinear electrical circuit.
ParameterList & sublist(const std::string &name, bool mustAlreadyExist=false, const std::string &docString="")
Solution state for integrators and steppers.