13 #include "Tempus_config.hpp"
14 #include "Tempus_IntegratorBasic.hpp"
16 #include "Tempus_StepperFactory.hpp"
17 #include "Tempus_StepperNewmarkImplicitAForm.hpp"
18 #include "Tempus_StepperNewmarkImplicitDForm.hpp"
20 #include "../TestModels/HarmonicOscillatorModel.hpp"
21 #include "../TestUtils/Tempus_ConvergenceTestUtils.hpp"
28 namespace Tempus_Test {
30 using Teuchos::getParametersFromXmlFile;
33 using Teuchos::rcp_const_cast;
34 using Teuchos::rcp_dynamic_cast;
35 using Teuchos::sublist;
45 std::vector<RCP<Thyra::VectorBase<double>>> solutions;
46 std::vector<RCP<Thyra::VectorBase<double>>> solutionsDot;
47 std::vector<double> StepSize;
48 std::vector<double> xErrorNorm;
49 std::vector<double> xDotErrorNorm;
50 const int nTimeStepSizes = 10;
55 "Tempus_Test_NewmarkImplicitDForm_HarmonicOscillator_Damped_SecondOrder."
69 double dt = pl->
sublist(
"Default Integrator")
71 .
get<
double>(
"Initial Time Step");
74 for (
int n = 0; n < nTimeStepSizes; n++) {
77 out <<
"\n \n time step #" << n <<
" (out of " << nTimeStepSizes - 1
78 <<
"), dt = " << dt <<
"\n";
79 pl->
sublist(
"Default Integrator")
81 .
set(
"Initial Time Step", dt);
82 integrator = Tempus::createIntegratorBasic<double>(pl, model);
85 bool integratorStatus = integrator->advanceTime();
89 time = integrator->getTime();
90 double timeFinal = pl->sublist(
"Default Integrator")
91 .sublist(
"Time Step Control")
92 .
get<
double>(
"Final Time");
98 integrator->getSolutionHistory();
100 "Tempus_Test_NewmarkImplicitDForm_HarmonicOscillator_Damped_"
106 for (
int i = 0; i < solutionHistory->getNumStates(); i++) {
107 double time_i = (*solutionHistory)[i]->getTime();
110 model->getExactSolution(time_i).get_x()),
112 model->getExactSolution(time_i).get_x_dot()));
113 state->setTime((*solutionHistory)[i]->getTime());
114 solnHistExact->addState(state);
117 "Tempus_Test_NewmarkImplicitDForm_HarmonicOscillator_Damped_"
118 "SecondOrder-Ref.dat",
123 StepSize.push_back(dt);
124 auto solution = Thyra::createMember(model->get_x_space());
125 Thyra::copy(*(integrator->getX()), solution.ptr());
126 solutions.push_back(solution);
127 auto solutionDot = Thyra::createMember(model->get_x_space());
128 Thyra::copy(*(integrator->getXDot()), solutionDot.ptr());
129 solutionsDot.push_back(solutionDot);
130 if (n == nTimeStepSizes - 1) {
131 StepSize.push_back(0.0);
132 auto solutionExact = Thyra::createMember(model->get_x_space());
133 Thyra::copy(*(model->getExactSolution(time).get_x()),
134 solutionExact.ptr());
135 solutions.push_back(solutionExact);
136 auto solutionDotExact = Thyra::createMember(model->get_x_space());
137 Thyra::copy(*(model->getExactSolution(time).get_x_dot()),
138 solutionDotExact.ptr());
139 solutionsDot.push_back(solutionDotExact);
145 double xDotSlope = 0.0;
147 double order = stepper->getOrder();
149 "Tempus_Test_NewmarkImplicitDForm_HarmonicOscillator_Damped_SecondOrder_"
151 stepper, StepSize, solutions, xErrorNorm, xSlope, solutionsDot,
152 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)
void writeSolution(const std::string filename, Teuchos::RCP< const Tempus::SolutionHistory< Scalar >> solutionHistory)
Consider the ODE: where is a constant, is a constant damping parameter, is a constant forcing par...
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.