MPI_Sendrecv使用2D矩阵的问题

时间:2015-11-03 15:54:44

标签: c mpi

据我所知,MPI_Sendrecv需要两个不同的缓冲区用于发送和接收。我的以下代码将(N / P)行的行发送到P-1个处理器但是它不起作用并且给我一个冻结的屏幕。我试图确保一切正确,但我没有看到问题出在哪里(我省略了变量声明以使其变短)

int **M, **FinalM, **M0;
M = malloc(N * sizeof (int *));
for (i = 0; i < N; i++) {
    M[i] = malloc(N * sizeof (int));
}
FinalM = malloc(N * sizeof (int *));
for (i = 0; i < N; i++) {
    FinalM[i] = malloc(n * sizeof (int));
}

M0 = malloc(N/P * sizeof (int *));
for (i = 0; i < N/P; i++) {
    M0[i] = malloc(N * sizeof (int));
}

c = N/P; // P is Number of Processors and N rows
if (rank == 0) {
    for (i = 0; i < P; i++) {
        k = i*c;
        k1 = (i + 1) * c;
        for (j = k; j < k1; j++) {
            MPI_Sendrecv(M[j], N, MPI_INT, i, TAG, FinalM[j], N, MPI_INT, i, TAG, MPI_COMM_WORLD, &status[d]);
        }
    }
} else {

    for (i = 0; i < (N / P); i++) {
        MPI_Recv(M0[i], N, MPI_INT, 0, TAG, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
    }
}

请说,有人可以给我一个解决这个问题的提示吗? 谢谢。

1 个答案:

答案 0 :(得分:1)

来自https://www.open-mpi.org/doc/v1.8/man3/MPI_Sendrecv.3.php#toc7的文档:

  

MPI_Sendrecv执行阻止发送和接收操作。

在您的情况下,主循环将矩阵的一部分发送给第一个工作人员并阻塞,等待工作人员做出响应。没有回应,所以它永远挂起。您需要在主循环中使用MPI_Send,然后在MPI_Recv使用另一个循环。工作人员必须使用MPI_Send发回一些内容。