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发送和接收数组错误

1 个答案: