9 #ifndef Tempus_WrapperModelEvaluatorSecondOrder_decl_hpp
10 #define Tempus_WrapperModelEvaluatorSecondOrder_decl_hpp
13 #include "Tempus_config.hpp"
15 #include "Thyra_VectorBase.hpp"
32 template <
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") {
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";
75 #ifdef VERBOSE_DEBUG_OUTPUT
76 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
84 Scalar t, Scalar beta, Scalar gamma)
86 #ifdef VERBOSE_DEBUG_OUTPUT
87 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
97 #ifdef VERBOSE_DEBUG_OUTPUT
98 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
106 #ifdef VERBOSE_DEBUG_OUTPUT
107 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
114 #ifdef VERBOSE_DEBUG_OUTPUT
115 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
122 #ifdef VERBOSE_DEBUG_OUTPUT
123 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
130 #ifdef VERBOSE_DEBUG_OUTPUT
131 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
138 #ifdef VERBOSE_DEBUG_OUTPUT
139 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
149 #ifdef VERBOSE_DEBUG_OUTPUT
150 *
out_ <<
"DEBUG: " << __PRETTY_FUNCTION__ <<
"\n";
217 #endif // Tempus_WrapperModelEvaluatorSecondOrder_hpp
EVALUATION_TYPE
EVALUATION_TYPE indicates the evaluation to apply to the implicit ODE.
virtual Thyra::ModelEvaluatorBase::OutArgs< Scalar > getOutArgs()
Get OutArgs the wrapper ModelEvalutor.
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
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.
void setAppModel(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &me)
Set the underlying application ModelEvaluator.
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getAppModel() const
Get the underlying application model 'f'.
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int p) const
Get the p space.
Thyra::VectorBase< Scalar > Vector
Teuchos::RCP< const Vector > d_pred_
Teuchos::RCP< TimeDerivative< Scalar > > timeDer_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Get the x-solution space.
A ModelEvaluator for residual evaluations given a state. This ModelEvaluator takes a state...
Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > appModel_
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > wrapperOutArgs_
WrapperModelEvaluatorSecondOrder(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &appModel, const std::string schemeName)
Constructor.
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
A ModelEvaluator which wraps the application ModelEvaluator.
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Teuchos::RCP< const Vector > v_pred_
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > wrapperInArgs_
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int p) const
This interface defines the time derivative connection between an implicit Stepper and WrapperModelEva...
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::InArgs< Scalar > getInArgs()
Get InArgs the wrapper ModelEvalutor.
WrapperModelEvaluatorSecondOrder()
Default constructor - not allowed.
virtual void setInArgs(Thyra::ModelEvaluatorBase::InArgs< Scalar > inArgs)
Set InArgs the wrapper ModelEvalutor.
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Teuchos::RCP< Teuchos::FancyOStream > out_
Solve for x and determine xDot from x.
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int i) const
Get the g space.
virtual void setOutArgs(Thyra::ModelEvaluatorBase::OutArgs< Scalar > outArgs)
Set OutArgs the wrapper ModelEvalutor.