Question

尝试在MPI进程之间传递数组指针并将其接收到动态分配的内存。它不断给出分段错误，我们认为这是由于我们在进程之间发送它的方式。我们的代码如下：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include "mpi.h"
#include <math.h>

int main(int argc, char* argv[])
{
    int my_rank, i, j, p;
    MPI_Request     request, request2;
    MPI_Status      mpi_status;

MPI_Init(&argc, &argv);

MPI_Comm_rank(MPI_COMM_WORLD, &my_rank);
MPI_Comm_size(MPI_COMM_WORLD, &p);
MPI_Barrier(MPI_COMM_WORLD);

    if (my_rank == 0) 
    {
            int size = 5;
            int **slice;
             slice = (int **)malloc(sizeof(int*)* size);
             slice[0] = (int *)malloc(sizeof(int)* size*size);

             for(i = 0; i< size;i++)
             {
                    slice[i] = (*slice + size*i);
            }

            for (i = 0; i < size; i++)
            {
                    for (j = 0; j < size; j++)
                    {
                            slice[i][j] = i*j;
                    }
            }
           for (i = 0; i < size; i++)
            {
                    for (j = 0; j < size; j++) 
                    {
                            printf("slice[%d][%d]: %d\n", i, j, slice[i][j]);
                    }
            }

            MPI_Send(&slice, size * size, MPI_INT, 1, 1, MPI_COMM_WORLD);

    } else {
            int local_size=5;
            int **local_slice;
            local_slice = (int **)malloc(sizeof(int*)* local_size);
            local_slice[0] = (int *)malloc(sizeof(int)* local_size*local_size);

             for(i = 0; i< local_size;i++)
             {
                    local_slice[i] = (*local_slice + local_size*i);
            }

            MPI_Recv(&local_slice, local_size * local_size, MPI_INT, 0, 1, MPI_COMM_WORLD, &mpi_status);

           for (i = 0; i < local_size; i++)
            {
                    for (j = 0; j < local_size; j++) 
                    {
                            printf("local_slice[%d][%d]: %d\n", i, j, local_slice[i][j]); 
                    }
            }

    }

    MPI_Finalize();
    return 0;
}

有人可以解释如何在MPI进程之间正确传递这种类型的数组吗？

Answer 1

您似乎需要将MPI_Send的第一个参数从&slice更改为slice[0]并对MPI_Recv执行相同操作。

MPI_Send和MPI_Recv分段错误

1 个答案: