如何使用 C 在 MPI 中使用 MPI_Datatype 发送嵌套结构

Question

我尝试使用 MPI_Datatype 发送以下结构，但 MPI_Send 在发送结构时崩溃。我想知道如何处理这种情况。以下是我编写的用于定义新 MPI 数据类型的代码：

   typedef struct
    {
       double x;
       double y;
    } vertex;

   typedef struct
   {
        int num_vertices;
        vertex vertex[2];
   } vertex_list;


       MPI_Datatype vertexType;
        MPI_Type_contiguous(2,MPI_DOUBLE,&vertexType);
        MPI_Type_commit(&vertexType);

      MPI_Datatype vertexListType;
        MPI_Datatype typev[3] = {MPI_INT, vertexType, MPI_UB};
        int blocklenv[3] = {1, 2, 1};
        MPI_Aint dispv[3];
        /* compute displacements of structure components */
        MPI_Address( vertexl, dispv);
        MPI_Address( vertexl[0].vertex, dispv+1);
        MPI_Address(  vertexl+1, dispv+2);
        base = dispv[0];

        for (i=0; i <3; i++)
         dispv[i] -= base;

        /* build datatype describing structure */
        MPI_Type_struct( 3, blocklenv, dispv, typev, &vertexListType);
        MPI_Type_commit(&vertexListType);

https： //docs.google.com/document/d/1OQFtx0ClkKQx7X91BlVgiizs5D9jShhtgsKafrgC7hk/edit?hl=zh-CN

Answer 1

我不确定您的结构是否会按照您的预期方式工作，但我可以分享我使用 MPI_Send 发送结构的经验。

无需创建显式 MPI 数据类型，只需发送结构本身即可，因为它的所有内容都位于连续的内存中。诀窍是为 MPI_Send 操作提供正确的大小和数据类型。

使用您的结构，这就是我过去所做的（假设变量 vertex_list list 已经定义）：

MPI_Send(&list, sizeof(vertex_list), MPI_BYTE, <rankToSendTo>, <tagInteger>, <comm>);

因此要发送的数据缓冲区是指向列表结构的指针， buffer 是 vertex_list 的大小（以字节为单位），MPI 数据类型只是字节。

在接收端，您只需要提供一个 vertex_list 引用作为接收缓冲区：

vertex_list receiveList;
MPI_Recv(&receiveList, sizeof(vertex_list), MPI_BYTE, <rankOfSender>, <tagInteger>, <comm>, <status>);

希望有帮助！