Tempus  Version of the Day
Time Integration
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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 
17 namespace Tempus {
18 
25 template <typename Scalar>
27  : public SensitivityModelEvaluatorBase<Scalar>,
29  public:
33  forwardModel,
34  const bool is_pseudotransient,
35  const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null)
36  {
37  forwardModel_ = forwardModel;
38  appExplicitModel_ = forwardModel_->getExplicitModel();
39  appImplicitModel_ = forwardModel_->getImplicitModel();
42  is_pseudotransient, pList));
45  is_pseudotransient, pList));
47  }
48 
51 
53 
54 
57  const
58  {
59  return forwardModel_;
60  }
61 
65  {
66  fsaExplicitModel_->setForwardSolutionHistory(sh);
67  fsaImplicitModel_->setForwardSolutionHistory(sh);
68  }
69 
73  {
74  fsaExplicitModel_->setForwardSolutionState(s);
75  fsaImplicitModel_->setForwardSolutionState(s);
76  }
77 
79  virtual void setSolver(
81  const bool force_W_update)
82  {
83  fsaImplicitModel_->setSolver(solver, force_W_update);
84  }
85 
87 
88  private:
91 
92  protected:
95 
102 };
103 
104 } // namespace Tempus
105 
106 #endif // Tempus_ModelEvaluatorPairIMEX_StaggeredFSA_hpp
ModelEvaluator pair for implicit and explicit (IMEX) evaulations.
Specialization of IMEX ME for &quot;staggered&quot; FSA method.
void setup(const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &explicitModel, const Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > &implicitModel)
Setup ME when using default constructor – for derived classes.
virtual void setForwardSolutionState(const Teuchos::RCP< const Tempus::SolutionState< Scalar > > &s)
Set solution state from forward state evaluation (for frozen state)
virtual void setForwardSolutionHistory(const Teuchos::RCP< const Tempus::SolutionHistory< Scalar > > &sh)
Set solution history from forward state evaluation (for interpolation)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
WrapperModelEvaluatorPairIMEX_StaggeredFSA(const Teuchos::RCP< const WrapperModelEvaluatorPairIMEX_Basic< Scalar > > &forwardModel, const bool is_pseudotransient, const Teuchos::RCP< const Teuchos::ParameterList > &pList=Teuchos::null)
Constructor.
SolutionHistory is basically a container of SolutionStates. SolutionHistory maintains a collection of...
virtual Teuchos::RCP< const Thyra::ModelEvaluator< Scalar > > getForwardModel() const
Get the underlying forward model.
Teuchos::RCP< const WrapperModelEvaluatorPairIMEX_Basic< Scalar > > forwardModel_
Transform a ModelEvaluator&#39;s sensitivity equations to its residual.
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.
A ModelEvaluator decorator for sensitivity analysis.
Solution state for integrators and steppers.