在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 * ，所述sizeof(float)在phie+iter*20*size*sizeof(float)是多余的，你会有效存储器的外部。

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

您还应该删除所有指针强制转换，它们都是多余的并且可以隐藏问题。 仅在知道需要时才进行投射。

关于你提到的另一个问题：该recvbuf的参数MPI_Gather只在根进程显著。 因此，仅根过程的iter重要。 但是，由于您的所有进程都经过运行MPI_Gather的相同循环，并且彼此之间无法相互超越，因此在匹配收集期间它们始终具有相同的iter 。