Teko_TpetraBlockedMappingStrategy.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_TpetraBlockedMappingStrategy_hpp__
#define __Teko_TpetraBlockedMappingStrategy_hpp__

// stl includes
#include <vector>

// Teuchos includes
#include "Teuchos_RCP.hpp"

// Thyra includes
#include "Thyra_MultiVectorBase.hpp"
#include "Thyra_LinearOpBase.hpp"
#include "Thyra_BlockedLinearOpBase.hpp"

// Tpetra includes
#include "Tpetra_Map.hpp"
#include "Tpetra_MultiVector.hpp"
#include "Tpetra_Import.hpp"
#include "Tpetra_Export.hpp"
#include "Tpetra_CrsMatrix.hpp"

// Teko includes
#include "Teko_TpetraOperatorWrapper.hpp"
#include "Teko_BlockingTpetra.hpp"
#include "Teko_ConfigDefs.hpp"

namespace Teko {
namespace TpetraHelpers {

class TpetraBlockedMappingStrategy : public MappingStrategy {
 public:


  TpetraBlockedMappingStrategy(const std::vector<std::vector<GO>>& vars,
                               const Teuchos::RCP<const Tpetra::Map<LO, GO, NT>>& baseMap,
                               const Teuchos::Comm<int>& comm);



  virtual void copyTpetraIntoThyra(const Tpetra::MultiVector<ST, LO, GO, NT>& tpetra_X,
                                   const Teuchos::Ptr<Thyra::MultiVectorBase<ST>>& thyra_X) const;

  virtual void copyThyraIntoTpetra(const Teuchos::RCP<const Thyra::MultiVectorBase<ST>>& thyra_Y,
                                   Tpetra::MultiVector<ST, LO, GO, NT>& tpetra_Y) const;

  virtual const Teuchos::RCP<const Tpetra::Map<LO, GO, NT>> domainMap() const { return domainMap_; }

  virtual const Teuchos::RCP<const Tpetra::Map<LO, GO, NT>> rangeMap() const { return rangeMap_; }

  virtual std::string toString() const { return std::string("TpetraBlockedMappingStrategy"); }




  void buildBlockTransferData(const std::vector<std::vector<GO>>& vars,
                              const Teuchos::RCP<const Tpetra::Map<LO, GO, NT>>& baseMap,
                              const Teuchos::Comm<int>& comm);

  const std::vector<Blocking::MapPair>& getMaps() const { return blockMaps_; }

  const Teuchos::RCP<Thyra::BlockedLinearOpBase<ST>> buildBlockedThyraOp(
      const Teuchos::RCP<const Tpetra::CrsMatrix<ST, LO, GO, NT>>& mat,
      const std::string& label = "<ANYM>") const;

  void rebuildBlockedThyraOp(const RCP<const Tpetra::CrsMatrix<ST, LO, GO, NT>>& mat,
                             const RCP<Thyra::BlockedLinearOpBase<ST>>& A) const;


 protected:
  // member variables


  Teuchos::RCP<const Tpetra::Map<LO, GO, NT>> domainMap_;
  Teuchos::RCP<const Tpetra::Map<LO, GO, NT>> rangeMap_;


  std::vector<Blocking::MapPair> blockMaps_;
  std::vector<Teuchos::RCP<Tpetra::Import<LO, GO, NT>>> blockImport_;
  std::vector<Teuchos::RCP<Tpetra::Export<LO, GO, NT>>> blockExport_;

  mutable std::vector<RCP<Tpetra::Vector<GO, LO, GO, NT>>> plocal2ContigGIDs;
};

}  // end namespace TpetraHelpers
}  // end namespace Teko

#endif