Fortran: melhor maneira de cronometrar seções do seu código?

Às vezes, ao otimizar o código, é necessário cronometrar determinadas partes do código, eu venho usando o seguinte há anos, mas fiquei imaginando se existe uma maneira melhor / mais simples de fazer isso?

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)

optimization fortran profiling Oscilação isópcica
fonte

Respostas:

Existem algumas outras maneiras de fazer isso, com vantagens e desvantagens:

MPI_WTIME : Este é um relógio de parede de alta resolução. É provavelmente a opção mais `confiável '; simplesmente funciona. A desvantagem é que, se o seu programa ainda não usa MPI, você terá que envolvê-lo (o que não é difícil).
Use um fortran intrínseco (como você tem): provavelmente é o mais fácil e geralmente suficiente, mas pode não funcionar tão bem em uma arquitetura estranha ou em trabalhos paralelos. Há um pouco de discussão sobre esse estouro de pilha
Finalizar uma chamada em C: Fortran e C são compatíveis com objetos, portanto é fácil escrever um wrapper em torno de chamadas em C. Um código com o qual trabalho usa getrusage, que pode ser uma escolha estranha. Há muitas discussões sobre isso no Stack Overflow.

Minha recomendação pessoal seria MPI_WTIME, como você sabe, funcionará bem onde houver MPI. Aqui está um exemplo de uma pesquisa rápida:

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

  ! code to be timed

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

Max Hutchinson
fonte

Se você usa o compilador GNU, consulte o gprof .

Em resumo, você adicionará o sinalizador -g ao seu compilador, assim:

g77 -g -pg -0 myprogram myprogram.F

Em seguida, execute a saída e um arquivo chamado gmon.out aparecerá no seu diretório. Em seguida, ligue

gprof --line myprogram gmon.out

Isso fornecerá um perfil de tempo de CPU linha por linha.

icurays1
fonte

Obrigado pela resposta, só preciso esclarecer que estava pedindo uma solução programática. Um criador de perfil é ótimo, mas é mais do que eu estava pedindo.

Oscilação isópica

a bandeira é -pg, -gé para símbolos de depuração (também interessante, mas não é obrigatório)

RSFalcon7

Ouvi em vários lugares que os horários fornecidos pelo gprof não são necessariamente precisos, como yosefk.com/blog/… , stackoverflow.com/questions/1777556/alternatives-to-gprof/… (e várias outras respostas de Mike Dunlavey estouro de pilha). Ferramentas como gprof e kcachegrind ainda são úteis, pois o número de chamadas de função ainda está correto e fornecem alguns dados de tempo, mas eu não os trataria como evangelho. O DOE tem algumas ferramentas para isso, mas não sei se elas são melhores do que inserir temporizadores.

precisa saber é o seguinte

Sério, @IsopycnalOscillation tenta usar o criador de perfil. É algo novo para aprender, mas irá ajudá-lo tremendamente (e limpar seu código!) A longo prazo.

tmarthal

thanks @tmarthal Eu já usei profilers antes e definitivamente usarei um para o meu próximo projeto - concordo totalmente com o que você disse.

Oscilação isópcica

Conforme mencionado pelo icurays1, a criação de perfil é a melhor. Você também pode simplificar um pouco o acima ...

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

onde o módulo utils contém ...

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

Se você tiver muitas seções, passe uma string, por exemplo, "section_id" em toc para que ela imprima o ID / nome junto com o tempo.

stali
fonte

Eu sugeriria não fazer t1e t2global, mas passar t1como parâmetro para ambas as funções, para permitir vários temporizadores. Você também pode simplesmente retornar a hora, não imprimir nada.

Pedro