MPI自定义数据类型

自定义数据类型

1.数据类型辅助函数

• MPI_Type_commit

int MPI_Type_commit(
  MPI_Datatype *datatype
);

在通信中使用数据类型对象之前，必须提交数据类型对象。提交的数据类型仍可用作数据类型构造函数中的参数。无需提交基本数据类型。它们是"预先承诺的"。

• MPI_Type_free

int MPI_Type_free(
  MPI_Datatype *datatype
);

标记与处理数据的数据类型关联的数据类型对象，并将数据类型设置到MPI_DATATYPE_NULL。当前使用此数据类型的任何通信都将正常完成。从释放的数据类型定义的派生数据类型不受影响。

• MPI_Type_extent（已弃用，由替换MPI_Type_get_extent）

int MPI_Type_extent(
  MPI_Datatype datatype,
  MPI_Aint *extent
);

此例程既具有线程安全，也具有中断安全。这意味着此例程可能安全地由多个线程使用，也可以在信号处理程序中使用。

• MPI_Address

int MPI_Address(
  void *location,
  MPI_Aint *address
);

获取内存中的某位置的地址，返回指定变量内存中的“绝对”地址。

2.Contiguous数据类型

• MPI_Type_contiguous

int MPI_Type_contiguous(
  int count,
  MPI_Datatype old_type,
  MPI_Datatype *new_type_p
);

连续复制：将原数据类型oldtype按顺序依次连续复制后，得到一个新的数据类型。比如定义连续的2个实数为一个新的数据类型。

 MPI_Type_contiguous( 2, MPI_CHAR, &type );

3.Vector数据类型

• MPI_Type_vector

int MPI_Type_vector(
  int count,
  int blocklength,
  int stride,
  MPI_Datatype old_type,
  MPI_Datatype *newtype_p
);

创建向量（步幅）数据类型。先连续复制blocklength个oldtype类型的数据，形成一个数据块；再通过等间隔的复制count个该数据块形成新的数据类型；相邻两个数据的的起始位置的位移相差stride乘extent(oldtype)个字节。

例如有oldtype数据

1，2，3，4，5，6，7，8

现在要求新的数据类型3个，每连续2个作为一个新的数据类型，间隔为1.

MPI_Type_vector(3,2,1,old,new)

• MPI_Type_hvector

int MPI_Type_hvector(
  int count,
  int blocklen,
  MPI_Aint stride,
  MPI_Datatype old_type,
  MPI_Datatype *newtype_p
);

功能同MPI_Type_vector相同，唯一区别就是stride的单位是字节。

4.Index数据类型

• MPI_Type_indexed

int MPI_Type_indexed(
  int count,
  int blocklens[],
  int indices[],//位移，和首地址的位移
  MPI_Datatype old_type,
  MPI_Datatype *newtype
);

创建索引数据类型：该函数生成的新数据类型有count个数据块构成，第i个数据块包含blocklens[i]个连续存放的oldtype,第i个数据块与首地址的偏移量（字节数）为indices[i]乘extent(oldtype)。

可以看作是MPI_Type_vector的扩展，区别是买个数据块的长度可以不同，数据块之间的间隔也可以不同。

• MPI_Type_hindexed

int MPI_Type_hindexed(
  int count,
  int blocklens[],
  MPI_Aint indices[],
  MPI_Datatype old_type,
  MPI_Datatype *newtype
);

功能同MPI_Type_vector相同，唯一区别就是indices数组存放的数据单位是字节。

• MPI_Type_struct

int MPI_Type_struct(
  int count,
  int blocklens[],
  MPI_Aint indices[],
  MPI_Datatype old_types[],
  MPI_Datatype *newtype
);

创建结构数据类型：与MPI_Type_indexed的区别在于每个数据块的数据类型可以不同，这里的indices数组存访的单位是字节。

该函数时最一般的新数据类型的构造函数，也是最广泛的一个。

MPI的数据打包与拆包

打包(Pack)和解包(Unpack)操作是为了发送不连续的数据， 在发送前显示地把数据包装大一个连续的缓冲区， 在接收之后从连续的缓冲区中解包。

• MPI_Pack

int MPI_Pack(
  void *inbuf,
  int incount,
  MPI_Datatype datatype,
  void *outbuf,
  int outcount,
  int *position,
  MPI_Comm comm
);

MPI_PACK把由inbuf,incount, datatype指定的发送缓冲区中的incount个datatype类型的消息放到起始为 outbuf 的连续空间，该空间共有 outcount 个字节。 输入缓冲区可以是 MPI_Send 允许的任何通信缓冲区。入口参数 position 的值是输出缓冲区中用于打包的起始地址，打包后它的值根据打包消息的大小来增加，出口参数 position 的值是被打包的消息占用的输出缓冲区后面的第一个地址。通过连续几次对不同位置的消息调用打包操作， 就将不连续的消息放到了一个连续的空间。comm参数是将在后面用于发送打包的消息时用的通信域。 (NOTE:这里要连续多次调用打包操作)

• MPI_Unpack

int MPI_Unpack(
  void *inbuf,
  int insize,
  int *position,
  void *outbuf,
  int outcount,
  MPI_Datatype datatype,
  MPI_Comm comm
);

MPI_UNPACK 和 MPI_PACK 对应， 它从 inbuf 和 insize 指定的缓冲区空间将不连续的消息解开，放到 outbuf, outcount, datatype 指定的缓冲区中。 输出缓冲区可以是 MPI_RECV 允许的任何通信缓冲区。 输入缓冲区是一个连续的存储空间，大小为insize字节， 开始地址为 inbuf。入口参数 position的初始值是输出缓冲区中被打包消息占用的起始地址， 解包后它的值根据打包消息的大小来增加，因此出口参数 position的值是输出缓冲区中被解包的消息占用空间后面的第一个地址。 通过连续几次对已打包的消息调用与打包时相应的解包操作，就可以将连续的消息解开放到一个不连续的空间。comm参数是用于接收消息的通信域。

• MPI_Pack_size

int MPI_Pack_size(
  int incount,
  MPI_Datatype datatype,
  MPI_Comm comm,
  int *size
);

描述MPI_Pack和MPI_Unpack的最大值。