Xpetra_TpetraMap_def.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//             Xpetra: A linear algebra interface package
//
// Copyright 2012 NTESS and the Xpetra contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef XPETRA_TPETRAMAP_DEF_HPP
#define XPETRA_TPETRAMAP_DEF_HPP

#include "Xpetra_TpetraMap_decl.hpp"

namespace Xpetra {



template <class LocalOrdinal, class GlobalOrdinal, class Node>
TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::
    TpetraMap(global_size_t numGlobalElements,
              GlobalOrdinal indexBase,
              const Teuchos::RCP<const Teuchos::Comm<int> > &comm,
              LocalGlobal lg)
  : map_(Teuchos::rcp(new Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node>(numGlobalElements,
                                                                         indexBase, comm,
                                                                         toTpetra(lg)))) {}


template <class LocalOrdinal, class GlobalOrdinal, class Node>
TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::
    TpetraMap(global_size_t numGlobalElements,
              size_t numLocalElements,
              GlobalOrdinal indexBase,
              const Teuchos::RCP<const Teuchos::Comm<int> > &comm)
  : map_(Teuchos::rcp(new Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node>(numGlobalElements,
                                                                         numLocalElements,
                                                                         indexBase, comm))) {}


template <class LocalOrdinal, class GlobalOrdinal, class Node>
TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::
    TpetraMap(global_size_t numGlobalElements,
              const Teuchos::ArrayView<const GlobalOrdinal> &elementList,
              GlobalOrdinal indexBase,
              const Teuchos::RCP<const Teuchos::Comm<int> > &comm)
  : map_(Teuchos::rcp(new Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node>(numGlobalElements,
                                                                         elementList,
                                                                         indexBase,
                                                                         comm))) {}

#ifdef HAVE_XPETRA_TPETRA

template <class LocalOrdinal, class GlobalOrdinal, class Node>
TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::
    TpetraMap(global_size_t numGlobalElements,
              const Kokkos::View<const GlobalOrdinal *, typename Node::device_type> &indexList,
              GlobalOrdinal indexBase,
              const Teuchos::RCP<const Teuchos::Comm<int> > &comm)
  : map_(Teuchos::rcp(new Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node>(numGlobalElements,
                                                                         indexList,
                                                                         indexBase,
                                                                         comm))) {}
#endif

template <class LocalOrdinal, class GlobalOrdinal, class Node>
TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::~TpetraMap() {}


template <class LocalOrdinal, class GlobalOrdinal, class Node>
global_size_t TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getGlobalNumElements() const {
  XPETRA_MONITOR("TpetraMap::getGlobalNumElements");
  return map_->getGlobalNumElements();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
size_t TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getLocalNumElements() const {
  XPETRA_MONITOR("TpetraMap::getLocalNumElements");
  return map_->getLocalNumElements();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
GlobalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getIndexBase() const {
  XPETRA_MONITOR("TpetraMap::getIndexBase");
  return map_->getIndexBase();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
LocalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMinLocalIndex() const {
  XPETRA_MONITOR("TpetraMap::getMinLocalIndex");
  return map_->getMinLocalIndex();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
LocalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMaxLocalIndex() const {
  XPETRA_MONITOR("TpetraMap::getMaxLocalIndex");
  return map_->getMaxLocalIndex();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
GlobalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMinGlobalIndex() const {
  XPETRA_MONITOR("TpetraMap::getMinGlobalIndex");
  return map_->getMinGlobalIndex();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
GlobalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMaxGlobalIndex() const {
  XPETRA_MONITOR("TpetraMap::getMaxGlobalIndex");
  return map_->getMaxGlobalIndex();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
GlobalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMinAllGlobalIndex() const {
  XPETRA_MONITOR("TpetraMap::getMinAllGlobalIndex");
  return map_->getMinAllGlobalIndex();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
GlobalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMaxAllGlobalIndex() const {
  XPETRA_MONITOR("TpetraMap::getMaxAllGlobalIndex");
  return map_->getMaxAllGlobalIndex();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
LocalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getLocalElement(GlobalOrdinal globalIndex) const {
  XPETRA_MONITOR("TpetraMap::getLocalElement");
  return map_->getLocalElement(globalIndex);
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
GlobalOrdinal TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getGlobalElement(LocalOrdinal localIndex) const {
  XPETRA_MONITOR("TpetraMap::getGlobalElement");
  return map_->getGlobalElement(localIndex);
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
LookupStatus TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getRemoteIndexList(const Teuchos::ArrayView<const GlobalOrdinal> &GIDList, const Teuchos::ArrayView<int> &nodeIDList, const Teuchos::ArrayView<LocalOrdinal> &LIDList) const {
  XPETRA_MONITOR("TpetraMap::getRemoteIndexList");
  return toXpetra(map_->getRemoteIndexList(GIDList, nodeIDList, LIDList));
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
LookupStatus TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getRemoteIndexList(const Teuchos::ArrayView<const GlobalOrdinal> &GIDList, const Teuchos::ArrayView<int> &nodeIDList) const {
  XPETRA_MONITOR("TpetraMap::getRemoteIndexList");
  return toXpetra(map_->getRemoteIndexList(GIDList, nodeIDList));
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
Teuchos::ArrayView<const GlobalOrdinal> TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getLocalElementList() const {
  XPETRA_MONITOR("TpetraMap::getLocalElementList");
  return map_->getLocalElementList();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
typename Map<LocalOrdinal, GlobalOrdinal, Node>::global_indices_array_device_type TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getMyGlobalIndicesDevice() const {
  XPETRA_MONITOR("TpetraMap::getMyGlobalIndicesDevice");
  return map_->getMyGlobalIndicesDevice();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::isNodeLocalElement(LocalOrdinal localIndex) const {
  XPETRA_MONITOR("TpetraMap::isNodeLocalElement");
  return map_->isNodeLocalElement(localIndex);
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::isNodeGlobalElement(GlobalOrdinal globalIndex) const {
  XPETRA_MONITOR("TpetraMap::isNodeGlobalElement");
  return map_->isNodeGlobalElement(globalIndex);
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::isContiguous() const {
  XPETRA_MONITOR("TpetraMap::isContiguous");
  return map_->isContiguous();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::isDistributed() const {
  XPETRA_MONITOR("TpetraMap::isDistributed");
  return map_->isDistributed();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::isCompatible(const Map<LocalOrdinal, GlobalOrdinal, Node> &map) const {
  XPETRA_MONITOR("TpetraMap::isCompatible");
  return map_->isCompatible(toTpetra(map));
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::isSameAs(const Map<LocalOrdinal, GlobalOrdinal, Node> &map) const {
  XPETRA_MONITOR("TpetraMap::isSameAs");
  return map_->isSameAs(toTpetra(map));
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
Teuchos::RCP<const Teuchos::Comm<int> > TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getComm() const {
  XPETRA_MONITOR("TpetraMap::getComm");
  return map_->getComm();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
std::string TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::description() const {
  XPETRA_MONITOR("TpetraMap::description");
  return map_->description();
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
void TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const {
  XPETRA_MONITOR("TpetraMap::describe");
  map_->describe(out, verbLevel);
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::removeEmptyProcesses() const {
  return toXpetra(map_->removeEmptyProcesses());
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::replaceCommWithSubset(const Teuchos::RCP<const Teuchos::Comm<int> > &newComm) const {
  return toXpetra(map_->replaceCommWithSubset(newComm));
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::TpetraMap(const Teuchos::RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > &map)
  : map_(map) {}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
UnderlyingLib TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::lib() const { return UseTpetra; }

template <class LocalOrdinal, class GlobalOrdinal, class Node>
RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getTpetra_Map() const { return map_; }

#ifdef HAVE_XPETRA_TPETRA

template <class LocalOrdinal, class GlobalOrdinal, class Node>
typename Map<LocalOrdinal, GlobalOrdinal, Node>::local_map_type TpetraMap<LocalOrdinal, GlobalOrdinal, Node>::getLocalMap() const {
  return map_->getLocalMap();
}
#endif

#ifdef HAVE_XPETRA_EPETRA

#if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \
     (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT))))

// specialization for Tpetra Map on EpetraNode and GO=int
template <>
class TpetraMap<int, int, EpetraNode>
  : public virtual Map<int, int, EpetraNode> {
 public:
  typedef int GlobalOrdinal;
  typedef int LocalOrdinal;
  typedef EpetraNode Node;




  TpetraMap(global_size_t numGlobalElements,
            GlobalOrdinal indexBase,
            const Teuchos::RCP<const Teuchos::Comm<int> > &comm,
            LocalGlobal lg = GloballyDistributed) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
  }


  TpetraMap(global_size_t numGlobalElements,
            size_t numLocalElements,
            GlobalOrdinal indexBase,
            const Teuchos::RCP<const Teuchos::Comm<int> > &comm) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
  }


  TpetraMap(global_size_t numGlobalElements,
            const Teuchos::ArrayView<const GlobalOrdinal> &elementList,
            GlobalOrdinal indexBase,
            const Teuchos::RCP<const Teuchos::Comm<int> > &comm) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
  }

  ~TpetraMap() {}




  global_size_t getGlobalNumElements() const { return 0; }

  size_t getLocalNumElements() const { return 0; }

  GlobalOrdinal getIndexBase() const { return 0; }

  LocalOrdinal getMinLocalIndex() const { return 0; }

  LocalOrdinal getMaxLocalIndex() const { return 0; }

  GlobalOrdinal getMinGlobalIndex() const { return 0; }

  GlobalOrdinal getMaxGlobalIndex() const { return 0; }

  GlobalOrdinal getMinAllGlobalIndex() const { return 0; }

  GlobalOrdinal getMaxAllGlobalIndex() const { return 0; }

  LocalOrdinal getLocalElement(GlobalOrdinal globalIndex) const { return 0; }

  GlobalOrdinal getGlobalElement(LocalOrdinal localIndex) const { return 0; }

  LookupStatus getRemoteIndexList(const Teuchos::ArrayView<const GlobalOrdinal> &GIDList, const Teuchos::ArrayView<int> &nodeIDList, const Teuchos::ArrayView<LocalOrdinal> &LIDList) const { return Xpetra::IDNotPresent; }

  LookupStatus getRemoteIndexList(const Teuchos::ArrayView<const GlobalOrdinal> &GIDList, const Teuchos::ArrayView<int> &nodeIDList) const { return Xpetra::IDNotPresent; }

  Teuchos::ArrayView<const GlobalOrdinal> getLocalElementList() const { return Teuchos::ArrayView<const GlobalOrdinal>(); }




  bool isNodeLocalElement(LocalOrdinal localIndex) const { return false; }

  bool isNodeGlobalElement(GlobalOrdinal globalIndex) const { return false; }

  bool isContiguous() const { return false; }

  bool isDistributed() const { return false; }

  bool isCompatible(const Map<LocalOrdinal, GlobalOrdinal, Node> &map) const { return false; }

  bool isSameAs(const Map<LocalOrdinal, GlobalOrdinal, Node> &map) const { return false; }




  Teuchos::RCP<const Teuchos::Comm<int> > getComm() const { return Teuchos::null; }




  std::string description() const { return std::string(""); }

  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const {}

  RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > removeEmptyProcesses() const { return Teuchos::null; }
  RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > replaceCommWithSubset(const Teuchos::RCP<const Teuchos::Comm<int> > &newComm) const { return Teuchos::null; }




  TpetraMap(const Teuchos::RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > &map) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
  }

  UnderlyingLib lib() const { return UseTpetra; }

  RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > getTpetra_Map() const { return Teuchos::null; }

#ifdef HAVE_XPETRA_TPETRA
  using local_map_type = typename Map<LocalOrdinal, GlobalOrdinal, Node>::local_map_type;
  local_map_type getLocalMap() const {
    return local_map_type();
  }
#endif


};  // TpetraMap class (specialization for GO=int and NO=EpetraNode)
#endif

#if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))) || \
     (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))))
// specialization for Tpetra Map on EpetraNode and GO=int
template <>
class TpetraMap<int, long long, EpetraNode>
  : public virtual Map<int, long long, EpetraNode> {
 public:
  typedef long long GlobalOrdinal;
  typedef int LocalOrdinal;
  typedef EpetraNode Node;



  TpetraMap(global_size_t numGlobalElements,
            GlobalOrdinal indexBase,
            const Teuchos::RCP<const Teuchos::Comm<int> > &comm,
            LocalGlobal lg = GloballyDistributed) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name());
  }

  TpetraMap(global_size_t numGlobalElements,
            size_t numLocalElements,
            GlobalOrdinal indexBase,
            const Teuchos::RCP<const Teuchos::Comm<int> > &comm) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name());
  }

  TpetraMap(global_size_t numGlobalElements,
            const Teuchos::ArrayView<const GlobalOrdinal> &elementList,
            GlobalOrdinal indexBase,
            const Teuchos::RCP<const Teuchos::Comm<int> > &comm) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name());
  }

  ~TpetraMap() {}




  global_size_t getGlobalNumElements() const { return 0; }

  size_t getLocalNumElements() const { return 0; }

  GlobalOrdinal getIndexBase() const { return 0; }

  LocalOrdinal getMinLocalIndex() const { return 0; }

  LocalOrdinal getMaxLocalIndex() const { return 0; }

  GlobalOrdinal getMinGlobalIndex() const { return 0; }

  GlobalOrdinal getMaxGlobalIndex() const { return 0; }

  GlobalOrdinal getMinAllGlobalIndex() const { return 0; }

  GlobalOrdinal getMaxAllGlobalIndex() const { return 0; }

  LocalOrdinal getLocalElement(GlobalOrdinal globalIndex) const { return 0; }

  GlobalOrdinal getGlobalElement(LocalOrdinal localIndex) const { return 0; }

  LookupStatus getRemoteIndexList(const Teuchos::ArrayView<const GlobalOrdinal> &GIDList, const Teuchos::ArrayView<int> &nodeIDList, const Teuchos::ArrayView<LocalOrdinal> &LIDList) const { return Xpetra::IDNotPresent; }

  LookupStatus getRemoteIndexList(const Teuchos::ArrayView<const GlobalOrdinal> &GIDList, const Teuchos::ArrayView<int> &nodeIDList) const { return Xpetra::IDNotPresent; }

  Teuchos::ArrayView<const GlobalOrdinal> getLocalElementList() const { return Teuchos::ArrayView<const GlobalOrdinal>(); }




  bool isNodeLocalElement(LocalOrdinal localIndex) const { return false; }

  bool isNodeGlobalElement(GlobalOrdinal globalIndex) const { return false; }

  bool isContiguous() const { return false; }

  bool isDistributed() const { return false; }

  bool isCompatible(const Map<LocalOrdinal, GlobalOrdinal, Node> &map) const { return false; }

  bool isSameAs(const Map<LocalOrdinal, GlobalOrdinal, Node> &map) const { return false; }




  Teuchos::RCP<const Teuchos::Comm<int> > getComm() const { return Teuchos::null; }




  std::string description() const { return std::string(""); }

  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const {}

  RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > removeEmptyProcesses() const { return Teuchos::null; }
  RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > replaceCommWithSubset(const Teuchos::RCP<const Teuchos::Comm<int> > &newComm) const { return Teuchos::null; }




  TpetraMap(const Teuchos::RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > &map) {
    XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMap<LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name());
  }

  UnderlyingLib lib() const { return UseTpetra; }

  RCP<const Tpetra::Map<LocalOrdinal, GlobalOrdinal, Node> > getTpetra_Map() const { return Teuchos::null; }

#ifdef HAVE_XPETRA_TPETRA
  using local_map_type = typename Map<LocalOrdinal, GlobalOrdinal, Node>::local_map_type;
  local_map_type getLocalMap() const {
    // We will never be here, this is a stub class
    return local_map_type();
  }
#endif

};  // TpetraMap class (specialization for GO=int and NO=EpetraNode)
#endif

#endif  // HAVE_XPETRA_EPETRA

}  // namespace Xpetra

// TODO: remove?
template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool operator==(const Xpetra::TpetraMap<LocalOrdinal, GlobalOrdinal, Node> &map1, const Xpetra::TpetraMap<LocalOrdinal, GlobalOrdinal, Node> &map2) {
  XPETRA_MONITOR("TpetraMap==TpetraMap");
  return map1.isSameAs(map2);
}

template <class LocalOrdinal, class GlobalOrdinal, class Node>
bool operator!=(const Xpetra::TpetraMap<LocalOrdinal, GlobalOrdinal, Node> &map1, const Xpetra::TpetraMap<LocalOrdinal, GlobalOrdinal, Node> &map2) {
  XPETRA_MONITOR("TpetraMap!=TpetraMap");
  return !map1.isSameAs(map2);
}

#endif  // XPETRA_TPETRAMAP_DEF_HPP