9#ifndef Tempus_WrapperModelEvaluatorSecondOrder_decl_hpp
10#define Tempus_WrapperModelEvaluatorSecondOrder_decl_hpp
13#include "Tempus_config.hpp"
15#include "Thyra_VectorBase.hpp"
32template <
typename Scalar>
43 out_(
Teuchos::VerboseObjectBase::getDefaultOStream())
45#ifdef VERBOSE_DEBUG_OUTPUT
46 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
47 *
out_ <<
" schemeName = " << schemeName <<
"\n";
49 if (schemeName ==
"Newmark Implicit a-Form" || schemeName ==
"HHT-Alpha") {
52 else if (schemeName ==
"Newmark Implicit d-Form") {
56 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
57 "Error: WrapperModelEvaluatorSecondOrder called with unsopported schemeName = " << schemeName
58 <<
"! Supported schemeNames are: 'Newmark Implicit a-Form' and 'HHT-Alpha'.\n");
66#ifdef VERBOSE_DEBUG_OUTPUT
67 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
73 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
getAppModel()
const
75#ifdef VERBOSE_DEBUG_OUTPUT
76 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
83 Teuchos::RCP<const Vector> d_pred, Scalar delta_t,
84 Scalar t, Scalar beta, Scalar gamma)
86#ifdef VERBOSE_DEBUG_OUTPUT
87 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
95 Teuchos::RCP<Thyra::LinearOpBase<Scalar> >
create_W_op()
const
97#ifdef VERBOSE_DEBUG_OUTPUT
98 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
103 Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> >
106#ifdef VERBOSE_DEBUG_OUTPUT
107 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
112 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_f_space()
const
114#ifdef VERBOSE_DEBUG_OUTPUT
115 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
120 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_p_space(
int p)
const
122#ifdef VERBOSE_DEBUG_OUTPUT
123 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
128 Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int p)
const
130#ifdef VERBOSE_DEBUG_OUTPUT
131 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
136 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_x_space()
const
138#ifdef VERBOSE_DEBUG_OUTPUT
139 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
144 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> >
get_g_space(
int i)
const
149#ifdef VERBOSE_DEBUG_OUTPUT
150 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
156 virtual void setInArgs(Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs)
160 virtual Thyra::ModelEvaluatorBase::InArgs<Scalar>
getInArgs()
164 virtual void setOutArgs(Thyra::ModelEvaluatorBase::OutArgs<Scalar> outArgs)
168 virtual Thyra::ModelEvaluatorBase::OutArgs<Scalar>
getOutArgs()
172 virtual void setForSolve(Teuchos::RCP<TimeDerivative<Scalar> > timeDer,
173 Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs,
174 Thyra::ModelEvaluatorBase::OutArgs<Scalar> outArgs,
182 Thyra::ModelEvaluatorBase::InArgs<Scalar>
createInArgs()
const;
186 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
187 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs)
const;
199 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
appModel_;
206 Teuchos::RCP<Teuchos::FancyOStream>
out_;
A ModelEvaluator for residual evaluations given a state. This ModelEvaluator takes a state,...
void initializeNewmark(Teuchos::RCP< const Vector > v_pred, Teuchos::RCP< const Vector > d_pred, Scalar delta_t, Scalar t, Scalar beta, Scalar gamma)
Set values needed in evalModelImpl.
WrapperModelEvaluatorSecondOrder()
Default constructor - not allowed.
Thyra::ModelEvaluatorBase::InArgs< Scalar > wrapperInArgs_
virtual Thyra::ModelEvaluatorBase::InArgs< Scalar > getInArgs()
Get InArgs the wrapper ModelEvalutor.
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int i) const
Get the g space.
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int p) const
Get the p space.
Teuchos::RCP< const Vector > d_pred_
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
void setAppModel(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &me)
Set the underlying application ModelEvaluator.
Teuchos::RCP< TimeDerivative< Scalar > > timeDer_
Teuchos::RCP< const Vector > v_pred_
virtual void setForSolve(Teuchos::RCP< TimeDerivative< Scalar > > timeDer, Thyra::ModelEvaluatorBase::InArgs< Scalar > inArgs, Thyra::ModelEvaluatorBase::OutArgs< Scalar > outArgs, EVALUATION_TYPE=SOLVE_FOR_X)
Set parameters for application implicit ModelEvaluator solve.
virtual Thyra::ModelEvaluatorBase::OutArgs< Scalar > getOutArgs()
Get OutArgs the wrapper ModelEvalutor.
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Get the x-solution space.
virtual void setOutArgs(Thyra::ModelEvaluatorBase::OutArgs< Scalar > outArgs)
Set OutArgs the wrapper ModelEvalutor.
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int p) const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > appModel_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
virtual void setInArgs(Thyra::ModelEvaluatorBase::InArgs< Scalar > inArgs)
Set InArgs the wrapper ModelEvalutor.
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getAppModel() const
Get the underlying application model 'f'.
Thyra::ModelEvaluatorBase::OutArgs< Scalar > wrapperOutArgs_
WrapperModelEvaluatorSecondOrder(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel, const std::string schemeName)
Constructor.
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Teuchos::RCP< Teuchos::FancyOStream > out_
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Thyra::VectorBase< Scalar > Vector
A ModelEvaluator which wraps the application ModelEvaluator.
EVALUATION_TYPE
EVALUATION_TYPE indicates the evaluation to apply to the implicit ODE.
@ SOLVE_FOR_X
Solve for x and determine xDot from x.