在 MPI_Gather C 中寻址内存

Question

我正在尝试将数据传递给MPI_Gather。我分配内存如下：

float *phie, *phitemp;
MPI_Comm_size(MPI_COMM_WORLD, &size);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);

phitemp=(float *) malloc(20*sizeof(float));
if (rank==1) phie=(float *) malloc(itermax*20*size*sizeof(float));

然后让所有进程使用MPI_Gather() 将数据发送到排名 1，如下所示：

for (iter=0;iter<itermax;iter++) {    
   MPI_Gather((float *) phitemp, 20, MPI_FLOAT, (float *) (phie+iter*20*size*sizeof(float)), 20, MPI_FLOAT, 1, MPI_COMM_WORLD);
   iter=0;

}

我收到错误消息，提示我没有正确分配内存。

Answer 1

指针运算是根据指针指向的字的大小进行的。由于phie 是float *，因此phie+iter*20*size*sizeof(float) 中的sizeof(float) 是多余的，您将超出有效内存。

删除sizeof(float)，或者更好地将其更改为更清晰的数组索引：&(phie[iter * 20 * size])

您还应该删除所有指针强制转换，它们都是多余的并且可以隐藏问题。只有在知道自己需要时才应该施放。

关于您的附加问题：MPI_Gather 的 recvbuf 参数仅在根进程中有意义。因此，只有根进程的iter 很重要。尽管如此，由于您的所有进程都通过相同的循环运行MPI_Gather，并且它们不能相互超越，因此在匹配聚集期间它们将始终具有相同的iter。