Fortran：计时代码部分的最佳方法？

有时在优化代码时需要对代码的某些部分进行计时，但我多年来一直在使用以下代码，但想知道是否有更简单/更好的方法来做到这一点？

call system_clock(count_rate=clock_rate) !Find the time rate
call system_clock(count=clock_start)     !Start Timer

call do_something_subroutine             !This is what gets timed

call system_clock(count=clock_stop)      ! Stop Timer

e_time = real(clock_stop-clock_start)/real(clock_rate)

还有其他几种方法，各有利弊：

• MPI_WTIME：这是高分辨率壁钟。这可能是最“值得信赖”的选择。它只是工作。缺点是，如果您的程序尚未使用MPI，则必须将MPI包装起来（这并不困难）。
• 使用fortran内部函数（如您所愿）：这可能是最简单的方法，通常就足够了，但是在奇怪的体系结构或并行作业上可能无法很好地工作。关于此堆栈溢出有一些讨论
• 包装C调用：Fortran和C是对象兼容的，因此它很容易为C调用编写包装器。我使用的代码使用getrusage，这可能是一个奇怪的选择。关于堆栈溢出有很多 讨论。

我个人的建议是MPI_WTIME，因为您知道只要有MPI，它都会很好地工作。这是快速搜索的示例：

  include 'mpif.h'
  DOUBLE PRECISION :: start, end
  start = MPI_Wtime()

  ! code to be timed

  end   = MPI_Wtime()
  write(*,*) 'That took ',end-start,' seconds'

如果使用GNU编译器，请签出gprof。

简而言之，您将-g标志添加到编译器中，如下所示：

g77 -g -pg -0 myprogram myprogram.F

然后，运行输出，一个名为gmon.out的文件将显示在您的目录中。然后打电话

gprof --line myprogram gmon.out

这将提供逐行的CPU时间配置文件。

如icurays1所述，分析是最好的。您也可以略微简化以上内容...

use utils
...
call tic()
   ! Section to be timed
call toc()
...
call tic()
   ! Section to be timed
call toc()
...

utils模块包含...

real(8) :: t1,t2
...
subroutine tic()
  implicit none
  call cpu_time(t1)
end subroutine tic

subroutine toc()
  implicit none
  call cpu_time(t2)
  ! if (rank==0) print*,"Time Taken -->", real(t2-t1)
  print*,"Time Taken -->", real(t2-t1)
end subroutine toc

如果您有很多这样的部分，则在toc中传递一个字符串，例如“ section_id”，以便它显示ID /名称以及时间。