MPI_Gather发送接收数组错误

Question

    #include<mpi.h>
    #include<stdio.h>
    int main(int argc,char * argv[])
    {
      int rank,size,m;
      int arr[1000];
      int b[100];
      MPI_Init(&argc,&argv);
      MPI_Comm_rank(MPI_COMM_WORLD,&rank);
      MPI_Comm_size(MPI_COMM_WORLD,&size);

      if(rank == 0)
      {
        scanf("%d",&m);
        for(int i=0;i<size*m;i++)
        {
          scanf("%d",&arr[i]);
        }

      }
      MPI_Barrier(MPI_COMM_WORLD);
      MPI_Scatter(arr,m,MPI_INT,b,m,MPI_INT,0,MPI_COMM_WORLD);

      printf("in process %d \n",rank);
      for(int i=0;i<m;i++)
      {
        printf("%d ",b[i]);
      }
      printf("\n");

      MPI_Finalize();
      return 0;
    }

输入

mpiexec -n 4 ./three

实际输出

in process 2 
in process 1 
in process 3 
2
1 2 3 4 5 6 7 8 
in process 0 
1 2

预期输出

2
1 2 3 4 5 6 7 8 
in process 0
1 2
in process 1
3 4 
in process 2
5 6 
in process 3
7 8

Answer 1

您没有广播 m 的值，因此只有等级 0 具有正确的值 - 看起来它在其他等级上默认为 m=0，因此它们缺少输出。

不需要屏障 - MPI 集体自动执行所有必要的同步。

如果您将障碍调用替换为：

MPI_Bcast(&m, 1, MPI_INT, 0, MPI_COMM_WORLD);

那么您的代码将按预期工作（输出可能不像上面 Gilles 所指出的那样干净）。

MPI_Gather发送接收数组错误

MPI_Gather send and recieving array error

c

mpi