MPI状态FORTRAN

3

我有这段代码,但是出现了错误。

我尝试将status声明为

INTEGER :: status

但是这会改变MPI_SENDRECV之后我的进程等级的值(即所有处理器的rank=0)。

PROGRAM testsendrecv
  IMPLICIT NONE

  INTEGER :: i, k, nx, nz
  INTEGER :: ierror, comm, p, rank, npr, prev
  INTEGER :: status(MPI_STATUS_SIZE)
  REAL(KIND = 8), ALLOCATABLE :: A(:,:), B(:), C(:)

  include 'mpif.h'

  nx = 5
  nz = 5

  ALLOCATE(A(nx,nz), B(nx))

  CALL MPI_INIT(ierror)
  comm = MPI_COMM_WORLD
  !Get rank
  CALL MPI_COMM_RANK(comm, rank, ierror)
  !Get number of processors
  CALL MPI_COMM_SIZE(comm, p, ierror)

  A(:,:) = rank

  IF(rank==0) THEN
     prev = p-1
  ELSE
     prev = rank-1
  END IF


  CALL MPI_SENDRECV(A(:,1), nx, MPI_DOUBLE_PRECISION, MOD(rank+1,p), 1,  &
       B(:), nx, MPI_DOUBLE_PRECISION, prev, 1, comm, status, ierror)

  WRITE(*,*) rank
  WRITE(*,*) B(1)

  CALL MPI_FINALIZE(ierror)

END PROGRAM testsendrecv

上述代码给我返回了以下错误信息。
bash-4.1$ mpif90 testsendr.f90
mpif.h:79.35:
    Included at testsendr.f90:9:

       PARAMETER (MPI_STATUS_SIZE=5)
                                   1
Error: VARIABLE attribute of 'mpi_status_size' conflicts with PARAMETER attribute at (1)
mpif.h:80.33:
    Included at testsendr.f90:9:

       INTEGER MPI_STATUS_IGNORE(MPI_STATUS_SIZE)
                                 1
Error: Variable 'mpi_status_size' cannot appear in the expression at (1)
mpif.h:80.49:
    Included at testsendr.f90:9:

       INTEGER MPI_STATUS_IGNORE(MPI_STATUS_SIZE)
                                                 1
Error: The module or main program array 'mpi_status_ignore' at (1) must have constant shape
mpif.h:81.35:
    Included at testsendr.f90:9:

       INTEGER MPI_STATUSES_IGNORE(MPI_STATUS_SIZE,1)
                                   1
Error: Variable 'mpi_status_size' cannot appear in the expression at (1)
mpif.h:81.53:
    Included at testsendr.f90:9:

       INTEGER MPI_STATUSES_IGNORE(MPI_STATUS_SIZE,1)
                                                     1
Error: The module or main program array 'mpi_statuses_ignore' at (1) must have constant shape
testsendr.f90:6.20:

  INTEGER :: status(MPI_STATUS_SIZE)
                    1
Error: Variable 'mpi_status_size' cannot appear in the expression at (1)
testsendr.f90:6.36:

  INTEGER :: status(MPI_STATUS_SIZE)
                                    1
Error: The module or main program array 'status' at (1) must have constant shape

任何想法。这是一个非常简单的程序。
谢谢。
1个回答

4
我认为你的问题来自于程序语句顺序错误。在此行之前:

include mpif.h

您在代码行中声明了一个变量,使用了该文件中定义的常量之一。

INTEGER :: status(MPI_STATUS_SIZE)

要么将include语句移动到IMPLICIT NONE之后,或者更好的方法是完全删除该语句并在隐式声明语句之前插入USE MPI,并解决修订代码的链接问题。

网页内容由stack overflow 提供, 点击上面的
可以查看英文原文,
原文链接