Teko_TpetraInverseFactoryOperator.hpp

// @HEADER
// *****************************************************************************
//      Teko: A package for block and physics based preconditioning
//
// Copyright 2010 NTESS and the Teko contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef __Teko_TpetraInverseFactoryOperator_hpp__
#define __Teko_TpetraInverseFactoryOperator_hpp__
 
#include "Teuchos_ConstNonconstObjectContainer.hpp"
 
#include "Teko_ConfigDefs.hpp"
#include "Teko_InverseFactory.hpp"
#include "Teko_TpetraInverseOpWrapper.hpp"
 
namespace Teko {
namespace TpetraHelpers {

class InverseFactoryOperator : public TpetraInverseOpWrapper {
 public:
  InverseFactoryOperator(const Teuchos::RCP<const InverseFactory> &bfp);

  virtual void initInverse(bool clearOld = false);

  virtual void buildInverseOperator(const Teuchos::RCP<const Tpetra::Operator<ST, LO, GO, NT> > &A,
                                    bool clear = true);

  virtual void buildInverseOperator(const Teuchos::RCP<Tpetra::Operator<ST, LO, GO, NT> > &A,
                                    bool clear = true);

  virtual void rebuildInverseOperator(
      const Teuchos::RCP<const Tpetra::Operator<ST, LO, GO, NT> > &A);

  virtual void rebuildInverseOperator(const Teuchos::RCP<Tpetra::Operator<ST, LO, GO, NT> > &A);

  Teuchos::RCP<const Tpetra::Operator<ST, LO, GO, NT> > getForwardOp() const {
    return fwdOp_.getConstObj();
  }

  Teuchos::RCP<Tpetra::Operator<ST, LO, GO, NT> > getNonconstForwardOp() const {
    return fwdOp_.getNonconstObj();
  }
 
 protected:
  Teuchos::RCP<const Thyra::LinearOpBase<ST> > extractLinearOp(
      const Teuchos::RCP<const Tpetra::Operator<ST, LO, GO, NT> > &A) const;
  Teuchos::RCP<const MappingStrategy> extractMappingStrategy(
      const Teuchos::RCP<const Tpetra::Operator<ST, LO, GO, NT> > &A) const;
 
  InverseFactoryOperator();
  InverseFactoryOperator(const InverseFactoryOperator &);
 
  Teuchos::RCP<const Teko::InverseFactory> inverseFactory_;
  Teko::ModifiableLinearOp invOperator_;
  bool firstBuildComplete_;
 
  Teuchos::ConstNonconstObjectContainer<Tpetra::Operator<ST, LO, GO, NT> > fwdOp_;
  bool setConstFwdOp_;
};
 
}  // end namespace TpetraHelpers
}  // end namespace Teko
 
#endif