Tempus Version of the Day
Time Integration
Loading...
Searching...
No Matches
Tempus_WrapperModelEvaluatorPairIMEX_StaggeredFSA.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_ModelEvaluatorPairIMEX_StaggeredFSA_hpp
10#define Tempus_ModelEvaluatorPairIMEX_StaggeredFSA_hpp
11
12#include "Tempus_config.hpp"
14#include "Tempus_WrapperModelEvaluatorPairIMEX_Basic.hpp"
15#include "Tempus_StaggeredForwardSensitivityModelEvaluator.hpp"
16
17namespace Tempus {
18
25template <typename Scalar>
27 : public SensitivityModelEvaluatorBase<Scalar>,
29{
30public:
31
34 const Teuchos::RCP<const WrapperModelEvaluatorPairIMEX_Basic<Scalar> >& forwardModel,
35 const bool is_pseudotransient,
36 const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null)
37 {
38 forwardModel_ = forwardModel;
39 appExplicitModel_ = forwardModel_->getExplicitModel();
40 appImplicitModel_ = forwardModel_->getImplicitModel();
41 fsaExplicitModel_ = rcp(new FSAME(appExplicitModel_, appExplicitModel_, appExplicitModel_, is_pseudotransient, pList));
42 fsaImplicitModel_ = rcp(new FSAME(appImplicitModel_, appImplicitModel_, appImplicitModel_, is_pseudotransient, pList));
44 }
45
48
50
51
53 virtual Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> >
55 {
56 return forwardModel_;
57 }
58
61 const Teuchos::RCP<const Tempus::SolutionHistory<Scalar> >& sh)
62 {
63 fsaExplicitModel_->setForwardSolutionHistory(sh);
64 fsaImplicitModel_->setForwardSolutionHistory(sh);
65 }
66
69 const Teuchos::RCP<const Tempus::SolutionState<Scalar> >& s)
70 {
71 fsaExplicitModel_->setForwardSolutionState(s);
72 fsaImplicitModel_->setForwardSolutionState(s);
73 }
74
76 virtual void setSolver(
77 const Teuchos::RCP<Thyra::NonlinearSolverBase<Scalar> >& solver,
78 const bool force_W_update)
79 {
80 fsaImplicitModel_->setSolver(solver, force_W_update);
81 }
82
84
85private:
86
89
90protected:
91
92 typedef WrapperModelEvaluatorPairIMEX_Basic<Scalar> Base;
93 typedef StaggeredForwardSensitivityModelEvaluator<Scalar> FSAME;
94
95 Teuchos::RCP<const WrapperModelEvaluatorPairIMEX_Basic<Scalar> > forwardModel_;
96 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > appExplicitModel_;
97 Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > appImplicitModel_;
98 Teuchos::RCP<FSAME> fsaExplicitModel_;
99 Teuchos::RCP<FSAME> fsaImplicitModel_;
100};
101
102} // namespace Tempus
103
104#endif // Tempus_ModelEvaluatorPairIMEX_StaggeredFSA_hpp
A ModelEvaluator decorator for sensitivity analysis.
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
Solution state for integrators and steppers. SolutionState contains the metadata for solutions and th...
ModelEvaluator pair for implicit and explicit (IMEX) evaulations.
virtual void setForwardSolutionState(const Teuchos::RCP< const Tempus::SolutionState< Scalar > > &s)
Set solution state from forward state evaluation (for frozen state)
virtual Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getForwardModel() const
Get the underlying forward model.
virtual void setForwardSolutionHistory(const Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > &sh)
Set solution history from forward state evaluation (for interpolation)
virtual void setSolver(const Teuchos::RCP< Thyra::NonlinearSolverBase< Scalar > > &solver, const bool force_W_update)
Set the solver of the underlying model if you want to reuse it.
WrapperModelEvaluatorPairIMEX_StaggeredFSA(const Teuchos::RCP< const WrapperModelEvaluatorPairIMEX_Basic< Scalar > > &forwardModel, const bool is_pseudotransient, const Teuchos::RCP< const Teuchos::ParameterList > &pList=Teuchos::null)
Constructor.
Teuchos::RCP< const WrapperModelEvaluatorPairIMEX_Basic< Scalar > > forwardModel_