Tempus_WrapperModelEvaluatorPairIMEX_CombinedFSA.hpp

Go to the documentation of this file.

//@HEADER
// *****************************************************************************
//          Tempus: Time Integration and Sensitivity Analysis Package
//
// Copyright 2017 NTESS and the Tempus contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
//@HEADER

#ifndef Tempus_ModelEvaluatorPairIMEX_CombinedFSA_hpp
#define Tempus_ModelEvaluatorPairIMEX_CombinedFSA_hpp

#include "Tempus_config.hpp"
#include "Tempus_SensitivityModelEvaluatorBase.hpp"
#include "Tempus_WrapperModelEvaluatorPairIMEX_Basic.hpp"
#include "Tempus_CombinedForwardSensitivityModelEvaluator.hpp"

namespace Tempus {

template <typename Scalar>
class WrapperModelEvaluatorPairIMEX_CombinedFSA
  : public SensitivityModelEvaluatorBase<Scalar>,
    public WrapperModelEvaluatorPairIMEX_Basic<Scalar> {
 public:
  WrapperModelEvaluatorPairIMEX_CombinedFSA(
      const Teuchos::RCP<const WrapperModelEvaluatorPairIMEX_Basic<Scalar> >&
          forwardModel,
      const Teuchos::RCP<const Teuchos::ParameterList>& pList = Teuchos::null)
  {
    forwardModel_     = forwardModel;
    appExplicitModel_ = forwardModel_->getExplicitModel();
    appImplicitModel_ = forwardModel_->getImplicitModel();
    fsaExplicitModel_ = rcp(new FSAME(appExplicitModel_, appExplicitModel_,
                                      appExplicitModel_, pList));
    fsaImplicitModel_ = rcp(new FSAME(appImplicitModel_, appImplicitModel_,
                                      appImplicitModel_, pList));
    Base::setup(fsaExplicitModel_, fsaImplicitModel_);
  }

  virtual ~WrapperModelEvaluatorPairIMEX_CombinedFSA() {}



  virtual Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > getForwardModel()
      const
  {
    return forwardModel_;
  }


 private:
  WrapperModelEvaluatorPairIMEX_CombinedFSA() {}

 protected:
  typedef WrapperModelEvaluatorPairIMEX_Basic<Scalar> Base;
  typedef CombinedForwardSensitivityModelEvaluator<Scalar> FSAME;

  Teuchos::RCP<const WrapperModelEvaluatorPairIMEX_Basic<Scalar> >
      forwardModel_;
  Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > appExplicitModel_;
  Teuchos::RCP<const Thyra::ModelEvaluator<Scalar> > appImplicitModel_;
  Teuchos::RCP<FSAME> fsaExplicitModel_;
  Teuchos::RCP<FSAME> fsaImplicitModel_;
};

}  // namespace Tempus

#endif  // Tempus_ModelEvaluatorPairIMEX_CombinedFSA_hpp