0

我正在尝试使用MPI_Bcast. v1是一个普通的块数组。我也不确定是否需要将公共数组的计算值广播v1给所有其他进程，否则它们会因为通用而在每个进程中更改。以下是相关的一段代码：

  parameter (nprocz=48,nzro=1)

    do i=i101,i102
      dist  = 0.015*float(i-iv0)
      adamp = exp(-dist*dist)
      do j = je0, je1-1
      do k = ke0, ke1
        v1(k,j,i) = v1(k,j,i)*adamp
      end do
      end do
    end do

    nmpi01=floor((iv0-ie0-nzro)/(nprocz-1))
    if (mpirank .le. nprocz-2) then
       i101=ie0+(mpirank*nmpi01)
       i102=ie0+(mpirank+1)*nmpi01-1
    else
       i101=ie0+(mpirank*nmpi01)
       i102=iv0-1
    endif

   MPI_Bcast(v1(:,:,i101:i102),(ke1-ke0+1)*(je1-je0)*(i102-i101+1)
 & ,MPI_FLOAT,mpirank,MPI_COMM_WORLD,ierr01)

我收到错误消息：

PGFTN-S-0081-Matrix/vector v1 illegal as subprogram argument

传入的数组的大小是正确的。任何意见？

我更正了代码并遍历了排名并计算了每个排名中 rcount 和 displs 的所有元素：

integer :: myscount, myi101

do rank = 0, nprocz-1
  nmpi01=floor((iv0-ie0-nzro)/(nprocz-1))
  if (rank .le. nprocz-2) then
    i101=ie0+(rank*nmpi01)
    i102=ie0+(rank+1)*nmpi01-1
  else
    i101=ie0+(rank*nmpi01)
    i102=iv0-1
  endif
  scount=(i102-i101+1)*(je1-je0)*(ke1-ke0+1)
  rcount(rank+1)=scount
  displs(rank+1)=rank*scount+1
  if (rank .eq. mpirank) then
    myscount = scount
    myi101 = i101
 end if
end do

scount = myscount
i101 = myi101

call mpi_allgatherv(...)

但仍然是错误的结果。1-在我的情况下，每个部分的结果都用于下一部分，尤其是在之后。所以mpi_allgatherv我需要在每个部分之后添加mpi_barrier吗mpi_allgatherv？2-应该使用 mpi_in_place 吗？考虑我只有一个 3d 数组v1，每个子数组v1(1,1,i)都是通过某个过程计算的，我想将计算出的子数组放在同一数组的适当部分。3-我想我应该displs(i) = sum(rcount(1:i-1))+1为 i=>2 考虑到在 fortran77 中总是 displs(1)=1。所以我更正了这一点：循环之前，循环displs(1)=1内部和循环displs(rank+2)=rank*scount+1之后displs(nprocz+1)=0。我对吗？

4

1 回答 1

0

我记得，Fortran 77 对数组下标的限制比 Fortran 90 更严格，而 pgftn 是 Fortran 77 编译器。我会尝试传递v1(1,1,i101)给 mpi_bcast，而不是v1(:,:,i101:i102). （或者使用带有“-Mfixed”标志的 pgf95。）

如果每个进程都需要查看 v1，那么您确实需要使用 MPI 进行通信。MPI 任务之间不共享变量，即使是公共块中的任务也不共享。但是，如果每个进程都在计算 v1 的不同部分，那么每个进程都需要其他进程的一部分，那么您就不能mpi_bcast这样做了；改为使用mpi_allgather。

此外，如上所述，当您使用 MPI 程序时，您应该使用call它们，因为它们是子程序。

于 2013-08-18T10:17:54.213 回答