6

我想知道是否有人可以为我了解 MPI_Waital 函数。我有一个使用 MPI_Isend 和 MPI_Irecv 传递信息的程序。在所有发送和接收完成后,程序中的一个进程(在本例中为进程 0)将打印一条消息。我的 Isend/Irecv 正在工作,但消息在程序中的某个随机点打印出来;所以我试图使用 MPI_Waital 等到所有请求都完成后再打印消息。我收到以下错误消息:

Fatal error in PMPI_Waitall: Invalid MPI_Request, error stack:
PMPI_Waitall(311): MPI_Waitall(count=16, req_array=0x16f70d0, status_array=0x16f7260) failed
PMPI_Waitall(288): The supplied request in array element 1 was invalid (kind=0)

以下是一些相关代码:

MPI_Status *status;
MPI_Request *request;

MPI_Init(&argc, &argv);
MPI_Comm_rank(MPI_COMM_WORLD, &taskid);
MPI_Comm_size(MPI_COMM_WORLD, &numtasks);
status = (MPI_Status *) malloc(numtasks * sizeof(MPI_Status));
request = (MPI_Request *) malloc(numtasks * sizeof(MPI_Request));

/* Generate Data to send */
//Isend/Irecvs look like this:
MPI_Isend(&data, count, MPI_INT, dest, tag, MPI_COMM_WORLD, &request[taskid]);
MPI_Irecv(&data, count, MPI_INT, source, tag, MPI_COMM_WORLD, &request[taskid]);
MPI_Wait(&request[taskid], &status[taskid]

/* Calculations and such */

if (taskid == 0) {
        MPI_Waitall (numtasks, request, status);
        printf ("All done!\n");
}
MPI_Finalize();

如果不调用 MPI_Waitall,程序运行干净,但只要进程 0 的 Isend/Irecv 消息完成,就会打印“All done”消息,而不是在所有 Isend/Irecvs 完成后打印。

感谢您提供任何帮助。

4

1 回答 1

11

您只设置了request数组的一个元素,即request[taskid](顺便说一下,您用接收元素覆盖发送请求句柄,不可撤销地丢失前者)。请记住,MPI 用于对分布式内存机器进行编程,每个 MPI 进程都有自己的request数组副本。在等级中设置一个元素taskid不会神奇地将值传播到其他等级,即使这样做,请求也只有本地有效性。正确的实现是:

MPI_Status status[2];
MPI_Request request[2];

MPI_Init(&argc, &argv);
MPI_Comm_rank (MPI_COMM_WORLD, &taskid);
MPI_Comm_size (MPI_COMM_WORLD, &numtasks);

/* Generate Data to send */
//Isend/Irecvs look like this:
MPI_Isend (&data, count, MPI_INT, dest, tag, MPI_COMM_WORLD, &request[0]);
//          ^^^^
//           ||
//      data race !!
//           ||
//          vvvv
MPI_Irecv (&data, count, MPI_INT, source, tag, MPI_COMM_WORLD, &request[1]);
// Wait for both operations to complete
MPI_Waitall(2, request, status);

/* Calculations and such */

// Wait for all processes to reach this line in the code
MPI_Barrier(MPI_COMM_WORLD);

if (taskid == 0) {
  printf ("All done!\n");
}
MPI_Finalize();

顺便说一句,您的代码中存在数据竞争。两者MPI_IsendMPI_Irecv使用相同的数据缓冲区,这是不正确的。如果您只是尝试发送 to 的内容data然后dest从 接收到它source,那么MPI_Sendrecv_replace请改用并忘记非阻塞操作:

MPI_Status status;

MPI_Init(&argc, &argv);
MPI_Comm_rank (MPI_COMM_WORLD, &taskid);
MPI_Comm_size (MPI_COMM_WORLD, &numtasks);

/* Generate Data to send */
MPI_Sendrecv_replace (&data, count, MPI_INT, dest, tag, source, tag,
                      MPI_COMM_WORLD, &status);

/* Calculations and such */

// Wait for all processes to reach this line in the code
MPI_Barrier(MPI_COMM_WORLD);

if (taskid == 0) {
  printf ("All done!\n");
}
MPI_Finalize();
于 2013-03-15T20:38:31.717 回答