source: src/comm/mpi/datatype.hpp@ a40eea

#ifndef DATATYPE_HPP_
#define DATATYPE_HPP_

#include <vector>

#include "base/index.hpp"
#include "grid/grid.hpp"

namespace VMG
{

namespace MPI
{

class Datatype
{
public:
  Datatype() :
    _sizes(0),
    _subsizes(0),
    _starts(0),
    _rank(-1),
    _type(MPI_DATATYPE_NULL)
  {}

  Datatype(Index sizes, Index subsizes, Index starts, const int& rank,
           const int& tag_send, const int& tag_receive) :
    _sizes(sizes),
    _subsizes(subsizes),
    _starts(starts),
    _rank(rank),
    _tag_send(tag_send),
    _tag_recv(tag_receive),
    _buffer(subsizes.Product())
  {
    InitDatatype();
  }

  Datatype(const GridIteratorSet& bounds, const Grid& grid, const int& rank,
           const int& tag_send, const int& tag_receive) :
    _sizes(grid.Local().SizeTotal()),
    _subsizes(bounds.Begin().GetEnd() - bounds.Begin().GetBegin()),
    _starts(bounds.Begin().GetBegin()),
    _rank(rank),
    _tag_send(tag_send),
    _tag_recv(tag_receive),
    _buffer(_subsizes.Product())
  {
    InitDatatype();
  }

  Datatype(const Datatype& other) :
    _sizes(other._sizes),
    _subsizes(other._subsizes),
    _starts(other._starts),
    _rank(other._rank),
    _tag_send(other._tag_send),
    _tag_recv(other._tag_recv),
    _buffer(_subsizes.Product())
  {
    InitDatatype();
  }

  ~Datatype()
  {
    if (_type != MPI_DATATYPE_NULL)
      MPI_Type_free(&_type);
  }

  void Set(const GridIteratorSet& bounds, const Grid& grid, const int& rank,
           const int& tag_send, const int& tag_receive);

  void Set(const Index& sizes, const Index& subsizes, const Index& starts, const int& rank,
           const int& tag_send, const int& tag_receive);

  const Index& Sizes() const {return _sizes;}
  const Index& Subsizes() const {return _subsizes;}
  const Index& Starts() const {return _starts;}

  const int& Rank() const {return _rank;}
  const int& TagSend() const {return _tag_send;}
  const int& TagReceive() const {return _tag_recv;}

  const MPI_Datatype& Type() const {return _type;}

  const std::vector<vmg_float>& Buffer() const {return _buffer;}

  void Send(Grid& grid, const int& tag, const MPI_Comm& comm) const;
  void Isend(Grid& grid, const int& tag, const MPI_Comm& comm, MPI_Request& request) const;
  void Recv(Grid& grid, const int& tag, const MPI_Comm& comm) const;
  void Irecv(Grid& grid, const int& tag, const MPI_Comm& comm, MPI_Request& request) const;

  void SendBuffered(const Grid& grid, const int& tag, const MPI_Comm& comm);
  void IsendBuffered(const Grid& grid, const int& tag, const MPI_Comm& comm, MPI_Request& request);
  void RecvBuffered(const int& tag, const MPI_Comm& comm);
  void IrecvBuffered(const int& tag, const MPI_Comm& comm, MPI_Request& request);

  void GridReplace(Grid& grid) const;
  void GridSum(Grid& grid) const;

  bool Feasible() const
  {
    return _sizes.Product() > 0 && _subsizes.Product() > 0;
  }

private:
  void InitDatatype();

  Index _sizes, _subsizes, _starts;
  int _rank, _tag_send, _tag_recv;
  MPI_Datatype _type;
  std::vector<vmg_float> _buffer;
};

}

}

#endif /* DATATYPE_HPP_ */