13 #include "Tempus_config.hpp"
14 #include "Tempus_IntegratorPseudoTransientForwardSensitivity.hpp"
15 #include "Tempus_IntegratorPseudoTransientAdjointSensitivity.hpp"
17 #include "Thyra_VectorStdOps.hpp"
18 #include "Thyra_MultiVectorStdOps.hpp"
20 #include "../TestModels/SteadyQuadraticModel.hpp"
22 #include "Thyra_DefaultMultiVectorProductVector.hpp"
23 #include "Thyra_DefaultProductVector.hpp"
28 namespace Tempus_Test {
30 using Teuchos::getParametersFromXmlFile;
33 using Teuchos::sublist;
41 RCP<ParameterList> pList =
42 getParametersFromXmlFile(
"Tempus_ExplicitRK_SteadyQuadratic.xml");
45 RCP<ParameterList> scm_pl = sublist(pList,
"SteadyQuadraticModel",
true);
46 scm_pl->set(
"Use DfDp as Tangent", use_dfdp_as_tangent);
47 RCP<SteadyQuadraticModel<double> > model =
51 RCP<ParameterList> pl = sublist(pList,
"Tempus",
true);
52 ParameterList &sens_pl = pl->sublist(
"Sensitivities");
53 sens_pl.set(
"Use DfDp as Tangent", use_dfdp_as_tangent);
56 RCP<Tempus::IntegratorPseudoTransientForwardSensitivity<double> > integrator =
57 Tempus::createIntegratorPseudoTransientForwardSensitivity<double>(pl,
61 bool integratorStatus = integrator->advanceTime();
65 double time = integrator->getTime();
66 double timeFinal = pl->sublist(
"Demo Integrator")
67 .sublist(
"Time Step Control")
68 .get<
double>(
"Final Time");
72 RCP<const Thyra::VectorBase<double> > x_vec = integrator->getX();
73 RCP<const Thyra::MultiVectorBase<double> > DxDp_vec = integrator->getDxDp();
74 const double x = Thyra::get_ele(*x_vec, 0);
75 const double dxdb = Thyra::get_ele(*(DxDp_vec->col(0)), 0);
76 const double x_exact = model->getSteadyStateSolution();
77 const double dxdb_exact = model->getSteadyStateSolutionSensitivity();
99 getParametersFromXmlFile(
"Tempus_ExplicitRK_SteadyQuadratic.xml");
112 Tempus::integratorPseudoTransientAdjointSensitivity<double>(pl, model);
115 bool integratorStatus = integrator->advanceTime();
119 double time = integrator->getTime();
120 double timeFinal = pl->sublist(
"Demo Integrator")
121 .sublist(
"Time Step Control")
122 .
get<
double>(
"Final Time");
128 const double x = Thyra::get_ele(*x_vec, 0);
129 const double dxdb = Thyra::get_ele(*(DxDp_vec->col(0)), 0);
130 const double x_exact = model->getSteadyStateSolution();
131 const double dxdb_exact = model->getSteadyStateSolutionSensitivity();
133 TEST_FLOATING_EQUALITY(x, x_exact, 1.0e-6);
134 TEST_FLOATING_EQUALITY(dxdb, dxdb_exact, 1.0e-6);
void test_pseudotransient_fsa(const bool use_dfdp_as_tangent, Teuchos::FancyOStream &out, bool &success)
#define TEST_FLOATING_EQUALITY(v1, v2, tol)
TEUCHOS_UNIT_TEST(BackwardEuler, SinCos_ASA)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Simple quadratic equation with a stable steady-state.