实际上MPI_Waitall是如何工作的

Question

这个例子是否与手册相矛盾？该手册指出请求数组和状态数组必须具有相同的大小.更确切地说,两个数组应该至少与count参数指示的一样长.然而在下面的示例中,状态数组大小为2,而不是4.此外,该示例与手册中的此语句相矛盾

对于i = 0,...,count-1,MPI_Waitall(count,array_of_requests,array_of_statuses)的无错执行与执行MPI_Wait(&array_of_request [i],&array_of_statuses [i])具有相同的效果任意顺序.

#include "mpi.h"
#include <stdio.h>

int main(argc,argv)
int argc;
char *argv[];  {
int numtasks, rank, next, prev, buf[2], tag1=1, tag2=2;
MPI_Request reqs[4];
MPI_Status stats[2];

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

prev = rank-1;
next = rank+1;
if (rank == 0)  prev = numtasks - 1;
if (rank == (numtasks - 1))  next = 0;

MPI_Irecv(&buf[0], 1, MPI_INT, prev, tag1, MPI_COMM_WORLD, &reqs[0]);
MPI_Irecv(&buf[1], 1, MPI_INT, next, tag2, MPI_COMM_WORLD, &reqs[1]);

MPI_Isend(&rank, 1, MPI_INT, prev, tag2, MPI_COMM_WORLD, &reqs[2]);
MPI_Isend(&rank, 1, MPI_INT, next, tag1, MPI_COMM_WORLD, &reqs[3]);

      {  do some work  }

MPI_Waitall(4, reqs, stats);

MPI_Finalize();
}

PS定义主要看起来很奇怪.返回值丢失.是史前的C还是错字？

Answer 1

是的,这个例子与手册相矛盾.如果你将这个例子与Fortran版本进行比较,你会发现Fortran版本是正确的,因为状态数组足够大(奇怪的是,它是一个2D数组,但由于隐式接口和存储关联,它可以看作是大小为MPI_STATUS_SIZE*2的1D数组大于4,前提是MPI_STATUS_SIZE大于1(在我的系统上为5).

是的,缺少的return语句是一个错误; 但是有些编译器只是在main()中发出了忽略return语句的警告.此外,代码的史前性可以在参数的K&R样式声明中看到.