Tempus Version of the Day
Time Integration
Loading...
Searching...
No Matches
HarmonicOscillatorModel_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_HARMONIC_OSCILLATOR_MODEL_DECL_HPP
10#define TEMPUS_TEST_HARMONIC_OSCILLATOR_MODEL_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
51template<class Scalar>
53 : public Thyra::StateFuncModelEvaluatorBase<Scalar>,
54 public Teuchos::ParameterListAcceptorDefaultBase
55{
56 public:
57
58 // Constructor
59 HarmonicOscillatorModel(Teuchos::RCP<Teuchos::ParameterList> pList = Teuchos::null, const bool use_accel_IC = false);
60
61 // Exact solution
62 Thyra::ModelEvaluatorBase::InArgs<Scalar> getExactSolution(double t) const;
63
66
67 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_x_space() const;
68 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_f_space() const;
69 Thyra::ModelEvaluatorBase::InArgs<Scalar> getNominalValues() const;
70 Teuchos::RCP<Thyra::LinearOpWithSolveBase<Scalar> > create_W() const;
71 Teuchos::RCP<Thyra::LinearOpBase<Scalar> > create_W_op() const;
72 Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> > get_W_factory() const;
73 Thyra::ModelEvaluatorBase::InArgs<Scalar> createInArgs() const;
74
75 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_p_space(int l) const;
76 Teuchos::RCP<const Teuchos::Array<std::string> > get_p_names(int l) const;
77 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > get_g_space(int j) const;
78
80
83 void setParameterList(Teuchos::RCP<Teuchos::ParameterList> const& paramList);
84 Teuchos::RCP<const Teuchos::ParameterList> getValidParameters() const;
86
87private:
88
89 void setupInOutArgs_() const;
90
93 Thyra::ModelEvaluatorBase::OutArgs<Scalar> createOutArgsImpl() const;
94 void evalModelImpl(
95 const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs_bar,
96 const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs_bar
97 ) const;
99
100private:
101 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > x_space_;
102 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > p_space_;
103 Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > g_space_;
104 Teuchos::RCP<Thyra::VectorBase<Scalar> > x_vec_;
105 Teuchos::RCP<Thyra::VectorBase<Scalar> > x_dot_vec_;
106 Teuchos::RCP<Thyra::VectorBase<Scalar> > x_dot_dot_vec_;
107 Teuchos::RCP<Thyra::VectorBase<Scalar> > p_init_;
108 int vecLength_; //Number of state unknowns (1)
109 int numResponses_; //Number of responses (1)
110 mutable Thyra::ModelEvaluatorBase::InArgs<Scalar> inArgs_;
111 mutable Thyra::ModelEvaluatorBase::OutArgs<Scalar> outArgs_;
112 mutable Thyra::ModelEvaluatorBase::InArgs<Scalar> nominalValues_;
113 mutable bool isInitialized_;
114 double c_; //damping coefficient
115 double f_; //forcing coefficient
116 double k_; //coefficient multiplying x in ODE
117 double m_; //mass coefficient
118 Teuchos::RCP<Teuchos::FancyOStream> out_;
119};
120
121
122
123} // namespace Tempus_Test
124#endif // TEMPUS_TEST_HARMONIC_OSCILLATOR_MODEL_DECL_HPP
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > inArgs_
Teuchos::RCP< Thyra::VectorBase< Scalar > > p_init_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int l) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > g_space_
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int l) const
Teuchos::RCP< Thyra::VectorBase< Scalar > > x_dot_vec_
Teuchos::RCP< Thyra::VectorBase< Scalar > > x_vec_
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Thyra::ModelEvaluatorBase::InArgs< Scalar > getExactSolution(double t) const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > outArgs_
void setParameterList(Teuchos::RCP< Teuchos::ParameterList > const &paramList)
Teuchos::RCP< Thyra::LinearOpWithSolveBase< Scalar > > create_W() const
Teuchos::RCP< Teuchos::FancyOStream > out_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int j) const
void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs_bar, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs_bar) const
Teuchos::RCP< Thyra::VectorBase< Scalar > > x_dot_dot_vec_
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
HarmonicOscillatorModel(Teuchos::RCP< Teuchos::ParameterList > pList=Teuchos::null, const bool use_accel_IC=false)
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > x_space_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > p_space_
Thyra::ModelEvaluatorBase::InArgs< Scalar > nominalValues_