c - 如何使用 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=en

Answer 1

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

与其创建显式 MPI 数据类型，还不如简单地发送结构本身，因为它的所有内容都在一 block 连续的内存中。诀窍是为 MPI_Send 操作提供正确的大小和数据类型。

使用你的结构，这是我过去所做的(假设变量 vertex_list list 已经被定义):

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

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

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

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

希望对您有所帮助!