MPI_Send中的细分错误，用于派生数据类型

Question

在下面的代码中，如果MPI_Get_address(&member, &offset[0]); 替换为offset[0] = 0; 该代码将按预期工作，否则将显示以下输出。 据我所知，要使用MPI_BOTTOM ，需要绝对内存地址，这就是为什么要使用MPI_Get_address()原因。 尽管struct Member MPI_Get_address()没有问题，但struct Family不能使用它。 问题是什么？

命令：

mpirun -n 2 ./out

输出：

Signal: Segmentation fault (11)
Signal code:  (128)
Failing at address: (nil)
mpirun noticed that process rank 0...

码：

#include <mpi.h>

struct Member
{
    double height;
    MPI_Datatype mpi_dtype;
    Member() { make_layout(); }
    void make_layout()
    {
        int nblock = 1;
        int block_count[nblock] = {1};
        MPI_Aint offset[nblock];
        MPI_Get_address(&height, &offset[0]);
        MPI_Datatype block_type[nblock] = {MPI_DOUBLE};
        MPI_Type_create_struct(nblock, block_count, offset, block_type, &mpi_dtype);
        MPI_Type_commit(&mpi_dtype);
    }
};

struct Family
{
    Member member;
    MPI_Datatype mpi_dtype;
    Family() { make_layout(); }
    void make_layout()
    {
        int nblock = 1;
        int block_count[nblock] = {1};
        MPI_Aint offset[nblock]; 
        MPI_Get_address(&member, &offset[0]);
        //offset[0] = 0; <-- HERE!!!!!!!
        MPI_Datatype block_type[nblock] = {member.mpi_dtype};
        MPI_Type_create_struct(nblock, block_count, offset, block_type, &mpi_dtype);
        MPI_Type_commit(&mpi_dtype);
    }
};

int main()
{
    int rank;

    MPI_Init(NULL, NULL);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);

    if (rank == 0)
    {
        Family f;
        MPI_Send(MPI_BOTTOM, 1, f.mpi_dtype, 1, 0, MPI_COMM_WORLD);
    }
    else if (rank == 1)
    {
        Family f;
        MPI_Recv(MPI_BOTTOM, 1, f.mpi_dtype, 0, 0, MPI_COMM_WORLD, NULL);
    }

    MPI_Finalize();

    return 0;
}

Answer 1

member.mpi_dtype已经在其typemap中带有Member::height的绝对地址作为偏移量。 在Family::make_layout()您指定的offset[0]等于member的地址，两个偏移量都加起来，这将导致错误的偏移量。 因此，在构造其他数据类型时不应使用具有绝对地址的MPI数据类型。

在这种情况下，绝对没有理由使用MPI_BOTTOM您的结构没有动态分配的字段，因此具有相对偏移的MPI数据类型就足够了。