12 #include "Teuchos_DefaultComm.hpp"
14 #include "Tempus_config.hpp"
15 #include "Tempus_IntegratorBasic.hpp"
16 #include "Tempus_StepperBDF2.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;
52 double dt = pList->
sublist(
"Tempus")
55 .
get<
double>(
"Initial Time Step");
60 const int nTimeStepSizes = scm_pl->
get<
int>(
"Number of Time Step Sizes", 7);
61 std::string output_file_string =
62 scm_pl->
get<std::string>(
"Output File Name",
"Tempus_BDF2_SinCos");
63 std::string output_file_name = output_file_string +
".dat";
64 std::string ref_out_file_name = output_file_string +
"-Ref.dat";
65 std::string err_out_file_name = output_file_string +
"-Error.dat";
67 for (
int n = 0; n < nTimeStepSizes; n++) {
74 getParametersFromXmlFile(
"Tempus_BDF2_SinCos.xml");
76 pl->
sublist(
"Default Integrator")
78 .
set(
"Initial Time Step", dt);
79 integrator = Tempus::createIntegratorBasic<double>(pl, model);
86 model->getNominalValues().get_x()->clone_v();
87 integrator->initializeSolutionHistory(0.0, x0);
90 bool integratorStatus = integrator->advanceTime();
94 time = integrator->getTime();
95 double timeFinal = pl->sublist(
"Default Integrator")
96 .sublist(
"Time Step Control")
97 .
get<
double>(
"Final Time");
103 integrator->getSolutionHistory();
107 for (
int i = 0; i < solutionHistory->getNumStates(); i++) {
108 double time_i = (*solutionHistory)[i]->getTime();
111 model->getExactSolution(time_i).get_x()),
113 model->getExactSolution(time_i).get_x_dot()));
114 state->setTime((*solutionHistory)[i]->getTime());
115 solnHistExact->addState(state);
121 StepSize.push_back(dt);
122 auto solution = Thyra::createMember(model->get_x_space());
123 Thyra::copy(*(integrator->getX()), solution.ptr());
124 solutions.push_back(solution);
125 auto solutionDot = Thyra::createMember(model->get_x_space());
126 Thyra::copy(*(integrator->getXDot()), solutionDot.ptr());
127 solutionsDot.push_back(solutionDot);
128 if (n == nTimeStepSizes - 1) {
129 StepSize.push_back(0.0);
130 auto solutionExact = Thyra::createMember(model->get_x_space());
131 Thyra::copy(*(model->getExactSolution(time).get_x()),
132 solutionExact.ptr());
133 solutions.push_back(solutionExact);
134 auto solutionDotExact = Thyra::createMember(model->get_x_space());
135 Thyra::copy(*(model->getExactSolution(time).get_x_dot()),
136 solutionDotExact.ptr());
137 solutionsDot.push_back(solutionDotExact);
142 if (nTimeStepSizes > 1) {
144 double xDotSlope = 0.0;
145 std::vector<double> xErrorNorm;
146 std::vector<double> xDotErrorNorm;
148 double order = stepper->getOrder();
149 writeOrderError(err_out_file_name, stepper, StepSize, solutions, xErrorNorm,
150 xSlope, solutionsDot, xDotErrorNorm, xDotSlope, out);
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.
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)
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.