使用 MPI_Send 和 Recv 发送 2D Int 数组

Question

我正在尝试从从属进程向主进程发送一个任意长度的二维整数数组，但我不断收到分段错误。由于 MPI 很难调试，我不确定问题是否与发送/接收有关，但如果不是，那么它必须与我自己分配数组的方式有关。

我在这里关注了一个关于确保分配给数组的内存是连续的问题，但这仍然没有解决分段错误。

以下是我的部分代码：

创建数组：

int** create2DArray(int sizeX, int sizeY)
{
    int* data = (int *) malloc(sizeX * sizeY * sizeof(int));
    int** array= (int **) malloc(sizeX * sizeof(int*));

    int i;

    for (i=0; i<sizeX; i++)
    {
        array[i] = &(data[sizeY * i]);
    }

    return array;
}

初始化数组：

if(rank==0)
{
    display = x11setup(&win, &gc, width, height);
    pixels = create2DArray(X_RESN, Y_RESN);
}
else
{
    xStart = xPixels * (rank - 1);
    xFinish = xStart + xPixels;
    pixels = create2DArray(xPixels, Y_RESN);
} 

发送：

MPI_Send(&pixels[0][0], xPixels * Y_RESN, MPI_INT, 0, type, MPI_COMM_WORLD);

接收：

for(i = 1; i < processes; i++)
{
    int** pixelChunk = create2DArray(xPixels, Y_RESN);
    MPI_Recv(&pixelChunk[0][0], xPixels * Y_RESN, MPI_INT, i, type, MPI_COMM_WORLD, &status);

    int xStart = xPixels * (i - 1);
    int xFinish = xStart + xPixels;

    int k;

    for(j = xStart; j < xFinish; j++)
    {
        for(k = 0; k < Y_RESN; k++)
        {
            pixels[j][k] = pixelChunk[j - (xPixels * i - 1)][k];
        }               
    }
}

Answer 1

这条线看起来很可疑：

pixels[j][k] = pixelChunk[j - (xPixels * i - 1)][k];

例如，假设我们有np = 2，那么我们只剩下一个块，那么

i = 1;
xStart = 0;
j = 0;
xPixels = 600;
pixelChunk[0 - (600 * 1 - 1)[k] == pixelChunk[-599][k]

看起来不太对，是吗？

这个？

pixels[j][k] = pixelChunk[j - xPixels * (i - 1)][k];

发送/接收代码大概没问题。