Tempus Version of the Day
Time Integration
Loading...
Searching...
No Matches
VanDerPol_IMEX_ExplicitModel_decl.hpp
Go to the documentation of this file.
1// @HEADER
2// ****************************************************************************
3// Tempus: Copyright (2017) Sandia Corporation
4//
5// Distributed under BSD 3-clause license (See accompanying file Copyright.txt)
6// ****************************************************************************
7// @HEADER
8
9#ifndef TEMPUS_TEST_VANDERPOL_IMEX_EXPLICITMODEL_DECL_HPP
10#define TEMPUS_TEST_VANDERPOL_IMEX_EXPLICITMODEL_DECL_HPP
11
12#include "Thyra_ModelEvaluator.hpp" // Interface
13#include "Thyra_StateFuncModelEvaluatorBase.hpp" // Implementation
14
15#include "Teuchos_ParameterListAcceptorDefaultBase.hpp"
16#include "Teuchos_ParameterList.hpp"
17
18namespace Tempus_Test {
19
106template<class Scalar>
108 : public Thyra::StateFuncModelEvaluatorBase<Scalar>,
109 public Teuchos::ParameterListAcceptorDefaultBase
110{
111 public:
112
113 // Constructor
115 Teuchos::RCP<Teuchos::ParameterList> pList = Teuchos::null,
116 bool useProductVector = false);
117
120 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_x_space() const;
121 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_f_space() const;
122 Thyra::ModelEvaluatorBase::InArgs<Scalar> getNominalValues() const;
123 Teuchos::RCP<Thyra::LinearOpBase<Scalar> > create_W_op() const;
124 Thyra::ModelEvaluatorBase::InArgs<Scalar> createInArgs() const;
125
126 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_p_space(int l) const;
127 Teuchos::RCP<const Teuchos::Array<std::string> > get_p_names(int l) const;
128 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_g_space(int j) const;
130
133 void setParameterList(Teuchos::RCP<Teuchos::ParameterList> const& paramList);
134 Teuchos::RCP<const Teuchos::ParameterList> getValidParameters() const;
136
137private:
138
139 void setupInOutArgs_() const;
140
143 Thyra::ModelEvaluatorBase::OutArgs<Scalar> createOutArgsImpl() const;
144 void evalModelImpl(
145 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs_bar,
146 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs_bar
147 ) const;
149
151 int dim_;
152 int Np_;
153 int np_;
154 int Ng_;
155 int ng_;
156 bool haveIC_;
159 mutable bool isInitialized_;
160 mutable Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs_;
161 mutable Thyra::ModelEvaluatorBase::OutArgs<Scalar> outArgs_;
162 mutable Thyra::ModelEvaluatorBase::InArgs<Scalar> nominalValues_;
163 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > xSpace_;
164 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > x_space_;
165 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > f_space_;
166 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > p_space_;
167 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > dxdp_space_;
168 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > g_space_;
169
170 // Parameters for the model:
171 Scalar epsilon_;
172 Scalar t0_ic_;
173 Scalar x0_ic_;
174 Scalar x1_ic_;
175};
176
177
178} // namespace Tempus_Test
179#endif // TEMPUS_TEST_VANDERPOL_IMEX_EXPLICITMODEL_DECL_HPP
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > xSpace_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > g_space_
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > inArgs_
bool haveIC_
false => no nominal values are provided (default=true)
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > f_space_
void setParameterList(Teuchos::RCP< Teuchos::ParameterList > const &paramList)
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs_bar, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs_bar) const
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int l) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int j) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > dxdp_space_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
int np_
Number of parameters in this vector (1)
Thyra::ModelEvaluatorBase::InArgs< Scalar > nominalValues_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
bool acceptModelParams_
Changes inArgs to require parameters.
Thyra::ModelEvaluatorBase::OutArgs< Scalar > outArgs_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > x_space_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int l) const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
int ng_
Number of elements in this observation function (0)
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > p_space_
bool useDfDpAsTangent_
Treat DfDp OutArg as tangent (df/dx*dx/dp+df/dp)