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;
50 getParametersFromXmlFile(
"Tempus_BDF2_SinCos_AdaptDt.xml");
53 double dt = pList->
sublist(
"Tempus")
56 .
get<
double>(
"Initial Time Step");
61 const int nTimeStepSizes = scm_pl->
get<
int>(
"Number of Time Step Sizes", 7);
62 std::string output_file_string =
63 scm_pl->
get<std::string>(
"Output File Name",
"Tempus_BDF2_SinCos");
64 std::string output_file_name = output_file_string +
".dat";
65 std::string err_out_file_name = output_file_string +
"-Error.dat";
67 for (
int n = 0; n < nTimeStepSizes; n++) {
73 getParametersFromXmlFile(
"Tempus_BDF2_SinCos_AdaptDt.xml");
77 pl->
sublist(
"Default Integrator")
79 .
set(
"Initial Time Step", dt / 4.0);
82 pl->
sublist(
"Default Integrator")
84 .
set(
"Maximum Time Step", dt);
86 pl->
sublist(
"Default Integrator")
88 .
set(
"Minimum Time Step", dt / 4.0);
90 pl->
sublist(
"Default Integrator")
92 .
sublist(
"Time Step Control Strategy")
93 .
set(
"Minimum Value Monitoring Function", dt * 0.99);
94 integrator = Tempus::createIntegratorBasic<double>(pl, model);
101 model->getNominalValues().get_x()->clone_v();
102 integrator->initializeSolutionHistory(0.0, x0);
105 bool integratorStatus = integrator->advanceTime();
109 time = integrator->getTime();
110 double timeFinal = pl->sublist(
"Default Integrator")
111 .sublist(
"Time Step Control")
112 .
get<
double>(
"Final Time");
118 model->getExactSolution(time).get_x();
122 std::ofstream ftmp(output_file_name);
124 FILE *gold_file = fopen(
"Tempus_BDF2_SinCos_AdaptDt_gold.dat",
"r");
126 integrator->getSolutionHistory();
128 for (
int i = 0; i < solutionHistory->getNumStates(); i++) {
129 char time_gold_char[100];
130 fgets(time_gold_char, 100, gold_file);
132 sscanf(time_gold_char,
"%lf", &time_gold);
134 double time_i = solutionState->getTime();
137 TEST_FLOATING_EQUALITY(time_i, time_gold, 1.0e-5);
139 x_exact_plot = model->getExactSolution(time_i).get_x();
140 ftmp << time_i <<
" " << get_ele(*(x_plot), 0) <<
" "
141 << get_ele(*(x_plot), 1) <<
" " << get_ele(*(x_exact_plot), 0)
142 <<
" " << get_ele(*(x_exact_plot), 1) << std::endl;
148 StepSize.push_back(dt);
149 auto solution = Thyra::createMember(model->get_x_space());
150 Thyra::copy(*(integrator->getX()), solution.ptr());
151 solutions.push_back(solution);
152 auto solutionDot = Thyra::createMember(model->get_x_space());
153 Thyra::copy(*(integrator->getXDot()), solutionDot.ptr());
154 solutionsDot.push_back(solutionDot);
155 if (n == nTimeStepSizes - 1) {
156 StepSize.push_back(0.0);
157 auto solutionExact = Thyra::createMember(model->get_x_space());
158 Thyra::copy(*(model->getExactSolution(time).get_x()),
159 solutionExact.ptr());
160 solutions.push_back(solutionExact);
161 auto solutionDotExact = Thyra::createMember(model->get_x_space());
162 Thyra::copy(*(model->getExactSolution(time).get_x_dot()),
163 solutionDotExact.ptr());
164 solutionsDot.push_back(solutionDotExact);
169 if (nTimeStepSizes > 1) {
171 double xDotSlope = 0.0;
172 std::vector<double> xErrorNorm;
173 std::vector<double> xDotErrorNorm;
177 solutions, xErrorNorm, xSlope, solutionsDot, xDotErrorNorm,
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...
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.