Zoltan2_Adapter.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//   Zoltan2: A package of combinatorial algorithms for scientific computing
//
// Copyright 2012 NTESS and the Zoltan2 contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef _ZOLTAN2_ADAPTER_HPP_
#define _ZOLTAN2_ADAPTER_HPP_
 
#include <Kokkos_Core.hpp>
#include <Zoltan2_Standards.hpp>
#include <Zoltan2_InputTraits.hpp>
#include <Zoltan2_PartitioningSolution.hpp>
#include <fstream>
 
namespace Zoltan2 {
 
//template<typename Adapter>
//class PartitioningSolution;
 
enum BaseAdapterType {
  InvalidAdapterType = 0, 
  IdentifierAdapterType,  
  VectorAdapterType,      
  MatrixAdapterType,      
  GraphAdapterType,       
  MeshAdapterType         
};
 
class BaseAdapterRoot {
public:
  virtual ~BaseAdapterRoot() = default; // required virtual declaration
 
  virtual size_t getLocalNumIDs() const = 0;
 
  virtual int getNumWeightsPerID() const { return 0; };
};
 
template <typename User>
  class BaseAdapter : public BaseAdapterRoot {
 
public:
  using lno_t = typename InputTraits<User>::lno_t;
  using gno_t = typename InputTraits<User>::gno_t;
  using scalar_t = typename InputTraits<User>::scalar_t;
  using node_t = typename InputTraits<User>::node_t;
  using part_t = typename InputTraits<User>::part_t;
  using offset_t = typename InputTraits<User>::offset_t;
  using host_t = typename Kokkos::HostSpace::memory_space;
  using device_t = typename node_t::device_type;
 
  using ConstIdsDeviceView = Kokkos::View<const gno_t *, device_t>;
  using ConstIdsHostView = typename ConstIdsDeviceView::host_mirror_type;
 
  using IdsDeviceView = Kokkos::View<gno_t *, device_t>;
  using IdsHostView = typename IdsDeviceView::host_mirror_type;
 
  using ConstOffsetsDeviceView = Kokkos::View<const offset_t *, device_t>;
  using ConstOffsetsHostView = typename ConstOffsetsDeviceView::host_mirror_type;
 
  using OffsetsDeviceView = Kokkos::View<offset_t *, device_t>;
  using OffsetsHostView = typename OffsetsDeviceView::host_mirror_type;
 
  using ConstScalarsDeviceView = Kokkos::View<const scalar_t *, device_t>;
  using ConstScalarsHostView = typename ConstScalarsDeviceView::host_mirror_type;
 
  using ScalarsDeviceView = Kokkos::View<scalar_t *, device_t>;
  using ScalarsHostView = typename ScalarsDeviceView::host_mirror_type;
 
  using ConstWeightsDeviceView1D = Kokkos::View<const scalar_t *, device_t>;
  using ConstWeightsHostView1D = typename ConstWeightsDeviceView1D::host_mirror_type;
 
  using WeightsDeviceView1D = Kokkos::View<scalar_t *, device_t>;
  using WeightsHostView1D = typename WeightsDeviceView1D::host_mirror_type;
 
  using ConstWeightsDeviceView = Kokkos::View<const scalar_t **, device_t>;
  using ConstWeightsHostView = typename ConstWeightsDeviceView::host_mirror_type;
 
  using WeightsDeviceView = Kokkos::View<scalar_t **, device_t>;
  using WeightsHostView = typename WeightsDeviceView::host_mirror_type;
 
  virtual enum BaseAdapterType adapterType() const = 0;
 
  virtual void getIDsView(const gno_t *&ids) const {
    // If adapter does not define getIDsView, getIDsKokkosView is called.
    // If adapter does not define getIDsKokkosView, getIDsView is called.
    // Allows forward and backwards compatibility.
    ConstIdsDeviceView kokkosIds;
    getIDsKokkosView(kokkosIds);
    //deep_copy(?)
    ids = kokkosIds.data();
  }
 
  virtual void getIDsKokkosView(ConstIdsDeviceView &ids) const {
    // If adapter does not define getIDsView, getIDsKokkosView is called.
    // If adapter does not define getIDsKokkosView, getIDsView is called.
    // Allows forward and backwards compatibility.
    const gno_t * ptr_ids;
    getIDsView(ptr_ids);
 
    auto non_const_ids = IdsDeviceView("ptr_ids", getLocalNumIDs());
    auto host_ids = Kokkos::create_mirror_view(non_const_ids);
    for(size_t i = 0; i < this->getLocalNumIDs(); ++i) {
       host_ids(i) = ptr_ids[i];
    }
    Kokkos::deep_copy(non_const_ids, host_ids);
    ids = non_const_ids;
  }
 
  virtual void getIDsHostView(ConstIdsHostView& hostIds) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  virtual void getIDsDeviceView(ConstIdsDeviceView &deviceIds) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  // *    \param wgt on return a pointer to the weights for this idx
  // *    \param stride on return, the value such that
  // *       the \t nth weight should be found at <tt> wgt[n*stride] </tt>.
  // *    \param idx  the weight index, zero or greater
  // *   This function or getWeightsKokkosView must be implemented in
  // *   derived adapter if getNumWeightsPerID > 0.
  // *   This function should not be called if getNumWeightsPerID is zero.
  // */
  virtual void getWeightsView(const scalar_t *&wgt, int &stride,
                              int idx = 0) const {
    // If adapter does not define getWeightsView, getWeightsKokkosView is called.
    // If adapter does not define getWeightsKokkosView, getWeightsView is called.
    // Allows forward and backwards compatibility.
    Kokkos::View<scalar_t **, device_t> kokkos_wgts_2d;
    getWeightsKokkosView(kokkos_wgts_2d);
    Kokkos::View<scalar_t *, device_t> kokkos_wgts;
    wgt = Kokkos::subview(kokkos_wgts_2d, Kokkos::ALL, idx).data();
    stride = 1;
  }
 
  // *    \param wgt on return a Kokkos view of the weights for this idx
  // *   This function or getWeightsView must be implemented in
  // *   derived adapter if getNumWeightsPerID > 0.
  // *   This function should not be called if getNumWeightsPerID is zero.
  // */
  virtual void getWeightsKokkosView(Kokkos::View<scalar_t **,
    device_t> & wgt) const {
    // If adapter does not define getWeightsKokkosView, getWeightsView is called.
    // If adapter does not define getWeightsView, getWeightsKokkosView is called.
    // Allows forward and backwards compatibility.
    wgt = Kokkos::View<scalar_t **, device_t>(
      "wgts", getLocalNumIDs(), getNumWeightsPerID());
    auto host_wgt = Kokkos::create_mirror_view(wgt);
    for(int j = 0; j < this->getNumWeightsPerID(); ++j) {
      const scalar_t * ptr_wgts;
      int stride;
      getWeightsView(ptr_wgts, stride, j);
      size_t i = 0;
      for(size_t n = 0; n < this->getLocalNumIDs() * stride; n += stride) {
        host_wgt(i++,j) = ptr_wgts[n];
      }
    }
    Kokkos::deep_copy(wgt, host_wgt);
  }
 
  virtual void getWeightsHostView(WeightsHostView1D& hostWgts, int idx = 0) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  virtual void getWeightsHostView(WeightsHostView& hostWgts) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  virtual void getWeightsDeviceView(WeightsDeviceView1D& deviceWgts, int idx = 0) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  virtual void getWeightsDeviceView(WeightsDeviceView& deviceWgts) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  virtual void getPartsView(const part_t *&inputPart) const {
    // Default behavior:  return NULL for inputPart array;
    // assume input part == rank
    inputPart = NULL;
  }
 
  virtual void getPartsHostView(Kokkos::View<part_t*, host_t> &inputPart) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  virtual void getPartsDeviceView(Kokkos::View<part_t*, device_t> &inputPart) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
  template <typename Adapter>
    void applyPartitioningSolution(const User &in, User *&out,
      const PartitioningSolution<Adapter> &solution) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
 
protected:
 
  // Write Chaco-formatted graph and assign files echoing adapter input
  // This routine is serial and may be slow; use it only for debugging
  // This function does not write edge info to the graph file, as the
  // BaseAdapter does not know about edge info; it writes
  // only the Chaco header and vertex weights (if applicable).
  void generateWeightFileOnly(const char* fileprefix,
                              const Teuchos::Comm<int> &comm) const;
 
};
 
template <typename User>
class AdapterWithCoords : public BaseAdapter<User>
{
public:
  using scalar_t = typename BaseAdapter<User>::scalar_t;
  using device_t = typename BaseAdapter<User>::node_t::device_type;
  using host_t = typename Kokkos::HostSpace::memory_space;
 
  // Coordinates in MJ are LayoutLeft since Tpetra Multivector gives LayoutLeft
  using CoordsDeviceView = Kokkos::View<scalar_t **, Kokkos::LayoutLeft, device_t>;
  using CoordsHostView = typename CoordsDeviceView::host_mirror_type;
 
public:
  virtual void getCoordinatesView(const scalar_t *&coords, int &stride,
                                  int coordDim) const = 0;
  virtual void getCoordinatesKokkosView(CoordsDeviceView &elements) const = 0;
 
  virtual void getCoordinatesHostView(CoordsHostView &) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
  virtual void getCoordinatesDeviceView(CoordsDeviceView &elements) const {
    Z2_THROW_NOT_IMPLEMENTED
  }
};
 
// Forward declare
template <typename User>
class VectorAdapter;
 
template <typename User, typename UserCoord=User>
class AdapterWithCoordsWrapper : public BaseAdapter<User>
{
public:
  virtual void setCoordinateInput(VectorAdapter<UserCoord> *coordData) = 0;
  virtual VectorAdapter<UserCoord> *getCoordinateInput() const = 0;
};
 
template <typename User>
void BaseAdapter<User>::generateWeightFileOnly(
  const char *fileprefix,
  const Teuchos::Comm<int> &comm
) const
{
  int np = comm.getSize();
  int me = comm.getRank();
 
  size_t nLocalIDs = this->getLocalNumIDs();
 
  // Write .graph file:  header and weights only (no edges)
  // Adapters with edges have to implement their own generateFiles function
  // to provide edge info.
  {
    // append suffix to filename
    std::string filenamestr = fileprefix;
    filenamestr = filenamestr + ".graph";
    const char *filename = filenamestr.c_str();
 
    size_t nGlobalIDs;
    Teuchos::reduceAll(comm, Teuchos::REDUCE_SUM, 1, &nLocalIDs, &nGlobalIDs);
 
    int nWgts = this->getNumWeightsPerID();
 
    for (int p = 0; p < np; p++) {
 
      // Is it this processor's turn to write to files?
      if (me == p) {
 
        std::ofstream fp;
 
        if (me == 0) {
          // open file for writing
          fp.open(filename, std::ios::out);
          // write Chaco header info
          // this function assumes no edges
          fp << nGlobalIDs << " " << 0 << " "
             << (nWgts ? "010" : "000") << " "
             << (nWgts > 1 ? std::to_string(nWgts) : " ") << std::endl;
        }
        else {
          // open file for appending
          fp.open(filename, std::ios::app);
        }
 
        if (nWgts) {
 
          // get weight data
          const scalar_t **wgts = new const scalar_t *[nWgts];
          int *strides = new int[nWgts];
          for (int n = 0; n < nWgts; n++)
            getWeightsView(wgts[n], strides[n], n);
 
          // write weights to file
          for (size_t i = 0; i < nLocalIDs; i++) {
            for (int n = 0; n < nWgts; n++)
              fp << wgts[n][i*strides[n]] << " ";
            fp << "\n";
          }
 
          delete [] strides;
          delete [] wgts;
        }
 
        fp.close();
      }
 
      comm.barrier();
    }
  }
 
  // write assignments file
  {
    std::string filenamestr = fileprefix;
    filenamestr = filenamestr + ".assign";
    const char *filename = filenamestr.c_str();
 
    for (int p = 0; p < np; p++) {
 
      // Is it this processor's turn to write to files?
      if (me == p) {
 
        std::ofstream fp;
 
        if (me == 0) {
          // open file for writing
          fp.open(filename, std::ios::out);
        }
        else {
          // open file for appending
          fp.open(filename, std::ios::app);
        }
 
        const part_t *parts;
        this->getPartsView(parts);
 
        for (size_t i = 0; i < nLocalIDs; i++) {
          fp << (parts != NULL ? parts[i] : me) << "\n";
        }
        fp.close();
      }
 
      comm.barrier();
    }
  }
}
 
}  //namespace Zoltan2
 
#endif