1

我正在尝试运行此 MPI Fortran 代码。有几个问题：

1）当我运行这段代码时，我希望程序在屏幕上写“输入间隔数：（0 退出）”然后问我 n。相反，它首先问我！为什么？

2）如果我不注释掉'goto 10'行，程序会一直问我n并且不会显示其他任何东西！

3）如果我注释掉'goto 10'，程序会问我n，然后写结果。但是，问题是每次程序写入部分结果而不是完整结果。它截断输出！以下是我运行程序连续三次的输出：

> mpiexec -n 40 ./a.out
10000000
 Enter the number of intervals: (0 quits)
 pi is    3.14159265358978      Error is   1.287858708565182E-014
 time is   1.687502861022949E-002  seconds

> mpiexec -n 40 ./a.out
10000000
 Enter the number of intervals: (0 quits)
 pi is    3.14159265358978      Error is   1.287858708565182E-014
 time is   1.68750286102

> mpiexec -n 40 ./a.out
10000000
 Enter the number of intervals: (0 quits)
 pi is    3.14159265358978      Error is   1.287858708565182E-014
 time is   1.687502861022949E-002  se

有人知道发生了什么吗？我提前感谢您的帮助。

      program main
      use mpi
      double precision starttime, endtime
      double precision  PI25DT
      parameter        (PI25DT = 3.141592653589793238462643d0)
      double precision  mypi, pi, h, sum, x, f, a
      double precision starttime, endtime
      integer n, myid, numprocs, i, ierr

      f(a) = 4.d0 / (1.d0 + a*a)   !   function to integrate

      call MPI_INIT(ierr)
      call MPI_COMM_RANK(MPI_COMM_WORLD, myid, ierr)
      call MPI_COMM_SIZE(MPI_COMM_WORLD, numprocs, ierr)

 10   if ( myid .eq. 0 ) then
         print *, 'Enter the number of intervals: (0 quits) '
         read(*,*) n
      endif
      starttime = MPI_WTIME()
!                                 broadcast n
      call MPI_BCAST(n,1,MPI_INTEGER,0,MPI_COMM_WORLD,ierr)
!                                 check for quit signal
      if ( n .le. 0 ) goto 30
!                                 calculate the interval size
      h   = 1.0d0/n
      sum = 0.0d0
      do 20 i = myid+1, n, numprocs
         x   = h * (dble(i) - 0.5d0)
         sum = sum + f(x)
 20   continue
      mypi = h * sum
!                                 collect all the partial sums
      call MPI_REDUCE(mypi,pi,1,MPI_DOUBLE_PRECISION,MPI_SUM,0, &
                      MPI_COMM_WORLD,ierr)
!                                 node 0 prints the answer.
      endtime = MPI_WTIME()
      if (myid .eq. 0) then
         print *, 'pi is ', pi, 'Error is ', abs(pi - PI25DT)
         print *, 'time is ', endtime-starttime, ' seconds'
      endif
      go to 10
 30   call MPI_FINALIZE(ierr)
      stop
      end
4

2 回答 2

3

该程序旨在通过最后的“goto 10”循环。打破这种情况的唯一方法是让 n 的值 <= 0，这将激活“goto 30”并分支通过“goto 10”。这是意图的其他线索是注释“检查退出信号”并且输入 n 的提示包括“（0退出）”。所以尝试输入0！

这不是现代 Fortran 的一个很好的例子。尽管清楚地使用了 Fortran 90 或更新的编译器（“use”语句显示了这一点），但它是以 FORTRAN 77 或更早的风格编写的。固定源代码布局，源代码行显然从第 7 列开始。第一列中的注释字符（旧的 Fortran 要求第一列中的 C）。“双精度”。程序逻辑大量使用 goto。（在我看来，其他人可能不同意，goto 语句有一个位置，但程序的基本控制流却没有。） 过时的 Fortran（在我看来）。

表达基本流程的现代 Fortran 方式：

MainLoop: do
   .....
   if (n .le. 0) exit MainLoop
   ....
end do MainLoop
于 2010-02-15T03:10:18.753 回答
1

你必须明确地刷新你的输出。我不记得fortran是否有标准的flush功能，如果flush不起作用，试试flush_。

基本上会发生什么，你的进程零缓冲输出，除非明确告诉它显示，否则你会得到一个有趣的东西

于 2010-02-14T18:09:01.890 回答