MueLu_IndexManager_decl.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//        MueLu: A package for multigrid based preconditioning
//
// Copyright 2012 NTESS and the MueLu contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef MUELU_INDEXMANAGER_DECL_HPP
#define MUELU_INDEXMANAGER_DECL_HPP

// use for Teuchos:Comm<T>
#include "Teuchos_CommHelpers.hpp"
#include "Teuchos_OrdinalTraits.hpp"

#include <Xpetra_Map_fwd.hpp>

#include "MueLu_ConfigDefs.hpp"
#include "MueLu_BaseClass.hpp"
#include "MueLu_IndexManager_fwd.hpp"

/*****************************************************************************

****************************************************************************/

namespace MueLu {

template <class LocalOrdinal, class GlobalOrdinal, class Node>
class IndexManager : public BaseClass {
#undef MUELU_INDEXMANAGER_SHORT
#include "MueLu_UseShortNamesOrdinal.hpp"

 private:
 protected:
  const RCP<const Teuchos::Comm<int> > comm_;  
  const bool coupled_;                         
  const bool singleCoarsePoint_;               
  const int numDimensions;                     
  const int interpolationOrder_;               

  Array<int> coarseRate;  
  Array<int> endRate;     

  GO gNumFineNodes;                  
  GO gNumFineNodes10;                
  const Array<GO> gFineNodesPerDir;  

  LO lNumFineNodes;                  
  LO lNumFineNodes10;                
  const Array<LO> lFineNodesPerDir;  

  GO gNumCoarseNodes;            
  GO gNumCoarseNodes10;          
  Array<GO> gCoarseNodesPerDir;  

  LO lNumCoarseNodes;            
  LO lNumCoarseNodes10;          
  Array<LO> lCoarseNodesPerDir;  

  LO numGhostNodes;              
  LO numGhostedNodes;            
  LO numGhostedNodes10;          
  Array<LO> ghostedNodesPerDir;  

  GO minGlobalIndex;                 
  Array<LO> offsets;                 
  Array<LO> coarseNodeOffsets;       
  Array<GO> startIndices;            
  Array<GO> startGhostedCoarseNode;  

  bool meshEdge[6]       = {false};  
  bool ghostInterface[6] = {false};  
  bool ghostedDir[6]     = {false};  

 public:
  IndexManager() = default;

  IndexManager(const RCP<const Teuchos::Comm<int> > comm, const bool coupled,
               const bool singleCoarsePoint, const int NumDimensions,
               const int interpolationOrder, const Array<GO> GFineNodesPerDir,
               const Array<LO> LFineNodesPerDir);

  virtual ~IndexManager() {}

  void computeMeshParameters();

  virtual void computeGlobalCoarseParameters() = 0;

  virtual void getGhostedNodesData(const RCP<const Map> fineMap,
                                   Array<LO>& ghostedNodeCoarseLIDs,
                                   Array<int>& ghostedNodeCoarsePIDs,
                                   Array<GO>& ghostedNodeCoarseGIDs) const = 0;

  virtual void getCoarseNodesData(const RCP<const Map> fineCoordinatesMap,
                                  Array<GO>& coarseNodeCoarseGIDs,
                                  Array<GO>& coarseNodeFineGIDs) const = 0;

  bool isAggregationCoupled() const { return coupled_; }

  bool isSingleCoarsePoint() const { return singleCoarsePoint_; }

  int getNumDimensions() const { return numDimensions; }

  int getInterpolationOrder() const { return interpolationOrder_; }

  GO getNumGlobalFineNodes() const { return gNumFineNodes; }

  GO getNumGlobalCoarseNodes() const { return gNumCoarseNodes; }

  LO getNumLocalFineNodes() const { return lNumFineNodes; }

  LO getNumLocalCoarseNodes() const { return lNumCoarseNodes; }

  LO getNumLocalGhostedNodes() const { return numGhostedNodes; }

  Array<int> getCoarseningRates() const { return coarseRate; }

  int getCoarseningRate(const int dim) const { return coarseRate[dim]; }

  Array<int> getCoarseningEndRates() const { return endRate; }

  int getCoarseningEndRate(const int dim) const { return endRate[dim]; }

  bool getMeshEdge(const int dir) const { return meshEdge[dir]; }

  bool getGhostInterface(const int dir) const { return ghostInterface[dir]; }

  Array<LO> getOffsets() const { return offsets; }

  LO getOffset(int const dim) const { return offsets[dim]; }

  Array<LO> getCoarseNodeOffsets() const { return coarseNodeOffsets; }

  LO getCoarseNodeOffset(int const dim) const { return coarseNodeOffsets[dim]; }

  Array<GO> getStartIndices() const { return startIndices; }

  GO getStartIndex(int const dim) const { return startIndices[dim]; }

  Array<GO> getStartGhostedCoarseNodes() const { return startGhostedCoarseNode; }

  GO getStartGhostedCoarseNode(int const dim) const { return startGhostedCoarseNode[dim]; }

  Array<LO> getLocalFineNodesPerDir() const { return lFineNodesPerDir; }

  LO getLocalFineNodesInDir(const int dim) const { return lFineNodesPerDir[dim]; }

  Array<GO> getGlobalFineNodesPerDir() const { return gFineNodesPerDir; }

  GO getGlobalFineNodesInDir(const int dim) const { return gFineNodesPerDir[dim]; }

  Array<LO> getLocalCoarseNodesPerDir() const { return lCoarseNodesPerDir; }

  LO getLocalCoarseNodesInDir(const int dim) const { return lCoarseNodesPerDir[dim]; }

  Array<GO> getGlobalCoarseNodesPerDir() const { return gCoarseNodesPerDir; }

  GO getGlobalCoarseNodesInDir(const int dim) const { return gCoarseNodesPerDir[dim]; }

  Array<LO> getGhostedNodesPerDir() const { return ghostedNodesPerDir; }

  LO getGhostedNodesInDir(const int dim) const { return ghostedNodesPerDir[dim]; }

  virtual std::vector<std::vector<GO> > getCoarseMeshData() const = 0;

  virtual void getFineNodeGlobalTuple(const GO myGID, GO& i, GO& j, GO& k) const = 0;

  virtual void getFineNodeLocalTuple(const LO myLID, LO& i, LO& j, LO& k) const = 0;

  virtual void getFineNodeGhostedTuple(const LO myLID, LO& i, LO& j, LO& k) const = 0;

  virtual void getFineNodeGID(const GO i, const GO j, const GO k, GO& myGID) const = 0;

  virtual void getFineNodeLID(const LO i, const LO j, const LO k, LO& myLID) const = 0;

  virtual void getCoarseNodeGlobalTuple(const GO myGID, GO& i, GO& j, GO& k) const = 0;

  virtual void getCoarseNodeLocalTuple(const LO myLID, LO& i, LO& j, LO& k) const = 0;

  virtual void getCoarseNodeGID(const GO i, const GO j, const GO k, GO& myGID) const = 0;

  virtual void getCoarseNodeLID(const LO i, const LO j, const LO k, LO& myLID) const = 0;

  virtual void getCoarseNodeGhostedLID(const LO i, const LO j, const LO k, LO& myLID) const = 0;

  virtual void getCoarseNodeFineLID(const LO i, const LO j, const LO k, LO& myLID) const = 0;

  virtual void getGhostedNodeFineLID(const LO i, const LO j, const LO k, LO& myLID) const = 0;

  virtual void getGhostedNodeCoarseLID(const LO i, const LO j, const LO k, LO& myLID) const = 0;
};

}  // namespace MueLu

#define MUELU_INDEXMANAGER_SHORT
#endif  // MUELU_INDEXMANAGER_DECL_HPP