Tpetra_Reindex_MultiVector_def.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//          Tpetra: Templated Linear Algebra Services Package
//
// Copyright 2008 NTESS and the Tpetra contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef TPETRA_REINDEX_MULTIVECTOR_DEF_HPP
#define TPETRA_REINDEX_MULTIVECTOR_DEF_HPP
 
#include <Tpetra_Reindex_MultiVector_decl.hpp>
 
 
namespace Tpetra {
 
template <class Scalar,
          class LocalOrdinal,
          class GlobalOrdinal,
          class Node>
Reindex_MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
    Reindex_MultiVector(Teuchos::RCP<Map<LocalOrdinal, GlobalOrdinal, Node> const> newRowMap)
  : ViewTransform<MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >()
  , newRowMap_(newRowMap) {
  // Nothing to do
}
 
template <class Scalar,
          class LocalOrdinal,
          class GlobalOrdinal,
          class Node>
Reindex_MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::~Reindex_MultiVector() {
  // Nothing to do
}
 
template <class Scalar,
          class LocalOrdinal,
          class GlobalOrdinal,
          class Node>
typename Reindex_MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::NewType
Reindex_MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::operator()(OriginalType const& origMultiVector) {
  this->origObj_ = origMultiVector;
  assert(origMultiVector->getMap()->getLocalNumElements() == this->newRowMap_->getLocalNumElements());
  assert(origMultiVector->isConstantStride() == true);  // So that it is valid to call origMultiVector->getStride()
 
  using mv_t    = MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>;
  this->newObj_ = Teuchos::RCP<mv_t>(new mv_t(*origMultiVector, Teuchos::DataAccess::Copy));
  this->newObj_->replaceMap(newRowMap_);
 
  return this->newObj_;
}
 
//
// Explicit instantiation macro
//
// Must be expanded from within the Tpetra namespace!
//
 
#define TPETRA_REINDEXMULTIVECTOR_INSTANT(SCALAR, LO, GO, NODE) \
  template class Reindex_MultiVector<SCALAR, LO, GO, NODE>;
 
}  // namespace Tpetra
 
#endif  // TPETRA_REINDEX_MULTIVECTOR_DEF_HPP