Fortran: Cách tốt nhất để các phần thời gian trong mã của bạn?

Đôi khi, trong khi tối ưu hóa mã, cần phải có thời gian cho một số phần nhất định của mã, tôi đã sử dụng những điều sau đây trong nhiều năm nhưng tự hỏi liệu có cách nào đơn giản hơn / tốt hơn để làm điều đó không?

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)

Có một vài cách khác để làm điều này, với những ưu điểm và nhược điểm:

• MPI_WTIME : Đây là đồng hồ treo tường có độ phân giải cao. Nó có lẽ là tùy chọn 'đáng tin cậy' nhất; nó chỉ hoạt động. Nhược điểm là nếu chương trình của bạn chưa sử dụng MPI, bạn sẽ phải bọc MPI xung quanh nó (điều này không khó).
• Sử dụng một fortranic (như bạn có): Đây có thể là dễ nhất và nói chung là đủ, nhưng có thể không hoạt động tốt trên một kiến ​​trúc lạ hoặc cho các công việc song song. Có một chút thảo luận về Stack Overflow này
• Kết thúc cuộc gọi C: Fortran và C tương thích với đối tượng, do đó, đủ dễ dàng để viết một trình bao bọc xung quanh các cuộc gọi C. Một mã tôi làm việc với việc sử dụng getrusage, có thể là một lựa chọn kỳ lạ. Có nhiều cuộc thảo luận về điều này trên Stack Overflow.

Đề xuất cá nhân của tôi sẽ là MPI_WTIME, vì bạn biết nó sẽ hoạt động tốt ở bất cứ nơi nào có MPI. Đây là một ví dụ từ một tìm kiếm nhanh:

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

  ! code to be timed

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

Nếu bạn sử dụng trình biên dịch GNU, hãy xem gprof .

Nói tóm lại, bạn sẽ thêm cờ -g vào trình biên dịch của mình, như vậy:

g77 -g -pg -0 myprogram myprogram.F

Sau đó, chạy đầu ra và một tệp có tên gmon.out sẽ hiển thị trong thư mục của bạn. Sau đó gọi

gprof --line myprogram gmon.out

Điều này sẽ cung cấp một hồ sơ thời gian CPU từng dòng.

Như được đề cập bởi icurays1 profiling là tốt nhất. Bạn cũng có thể đơn giản hóa một chút ở trên ...

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

nơi mô-đun utils chứa ...

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

Nếu bạn có nhiều phần như vậy thì hãy truyền một chuỗi, ví dụ: "part_id" trong toc để nó in id / name cùng với thời gian.