times.F90

Go to the documentation of this file.

module times
  integer,private,save :: last_count=0
  real, private,save :: last_cpucount=0
  logical, private,save :: alltimer_isactive=.false.
  
  integer, parameter :: nb_timer = 4
  integer, parameter :: timer_caldyn  = 1
  integer, parameter :: timer_vanleer = 2
  integer, parameter :: timer_dissip = 3
  integer, parameter :: timer_physic = 4
  integer, parameter :: stopped = 1
  integer, parameter :: running = 2
  integer, parameter :: suspended = 3 
  
  integer :: max_size
  real,    allocatable, dimension(:,:,:) :: timer_table
  real,    allocatable, dimension(:,:,:) :: timer_table_sqr 
  integer, allocatable, dimension(:,:,:) :: timer_iteration
  real,    allocatable, dimension(:,:,:) :: timer_average
  real,    allocatable, dimension(:,:,:) :: timer_delta
  real,    allocatable,dimension(:) :: timer_running, last_time
  integer, allocatable,dimension(:) :: timer_state
  
  contains
  
  subroutine init_timer
    use parallel
    implicit none
#include "dimensions.h"
#include "paramet.h"
    
    max_size=jjm+1
    allocate(timer_table(max_size,nb_timer,0:mpi_size-1))
    allocate(timer_table_sqr(max_size,nb_timer,0:mpi_size-1))
    allocate(timer_iteration(max_size,nb_timer,0:mpi_size-1))
    allocate(timer_average(max_size,nb_timer,0:mpi_size-1))
    allocate(timer_delta(max_size,nb_timer,0:mpi_size-1))
    allocate(timer_running(nb_timer))
    allocate(timer_state(nb_timer))
    allocate(last_time(nb_timer))
    
    timer_table(:,:,:)=0
    timer_table_sqr(:,:,:)=0
    timer_iteration(:,:,:)=0
    timer_average(:,:,:)=0
    timer_delta(:,:,:)=0
    timer_state(:)=stopped      
  end subroutine init_timer
  
  subroutine start_timer(no_timer)
    implicit none
    integer :: no_timer
    
    if (alltimer_isactive) then
    
      if (timer_state(no_timer)/=stopped) then
        stop 'start_timer :: timer is already running or suspended'
      else
        timer_state(no_timer)=running
      endif
      
      timer_running(no_timer)=0
      call cpu_time(last_time(no_timer))
    
    endif
    
  end subroutine start_timer
  
  subroutine suspend_timer(no_timer)
    implicit none
    integer :: no_timer
     
    if (alltimer_isactive) then   
      if (timer_state(no_timer)/=running) then
        stop 'suspend_timer :: timer is not running'
      else
        timer_state(no_timer)=suspended
      endif
    
      timer_running(no_timer)=timer_running(no_timer)-last_time(no_timer)
      call cpu_time(last_time(no_timer))
      timer_running(no_timer)=timer_running(no_timer)+last_time(no_timer)
    endif
  end subroutine suspend_timer
  
  subroutine resume_timer(no_timer)
    implicit none
    integer :: no_timer
     
    if (alltimer_isactive) then   
      if (timer_state(no_timer)/=suspended) then
        stop 'resume_timer :: timer is not suspended'
      else
        timer_state(no_timer)=running
      endif
      
      call cpu_time(last_time(no_timer))
    endif
    
  end subroutine resume_timer

  subroutine stop_timer(no_timer)
    use parallel
    implicit none
    integer :: no_timer
    integer :: n
    real :: v,v2
    
    if (alltimer_isactive) then
       
      if (timer_state(no_timer)/=running) then
        stop 'stop_timer :: timer is not running'
      else
        timer_state(no_timer)=stopped
      endif
    
      timer_running(no_timer)=timer_running(no_timer)-last_time(no_timer)
      call cpu_time(last_time(no_timer))
      timer_running(no_timer)=timer_running(no_timer)+last_time(no_timer)
    
      timer_table(jj_nb,no_timer,mpi_rank)=timer_table(jj_nb,no_timer,mpi_rank)+timer_running(no_timer)
      timer_table_sqr(jj_nb,no_timer,mpi_rank)=timer_table_sqr(jj_nb,no_timer,mpi_rank)+timer_running(no_timer)**2
      timer_iteration(jj_nb,no_timer,mpi_rank)=timer_iteration(jj_nb,no_timer,mpi_rank)+1
      timer_average(jj_nb,no_timer,mpi_rank)=timer_table(jj_nb,no_timer,mpi_rank)/timer_iteration(jj_nb,no_timer,mpi_rank)
      if (timer_iteration(jj_nb,no_timer,mpi_rank)>=2) then
        n=timer_iteration(jj_nb,no_timer,mpi_rank)
                v2=timer_table_sqr(jj_nb,no_timer,mpi_rank)
                v=timer_table(jj_nb,no_timer,mpi_rank)
                timer_delta(jj_nb,no_timer,mpi_rank)=sqrt(abs(v2-v*v/n)/(n-1)) 
      else
        timer_delta(jj_nb,no_timer,mpi_rank)=0
      endif
    endif
    
  end subroutine stop_timer
   
  subroutine allgather_timer
    use parallel
    implicit none
#ifdef CPP_MPI    
    include 'mpif.h'
#endif
    integer :: ierr
    integer :: data_size
    real, allocatable,dimension(:,:) :: tmp_table

    IF (using_mpi) THEN    
   
      if (alltimer_isactive) then
    
    
      allocate(tmp_table(max_size,nb_timer))
    
      data_size=max_size*nb_timer
    
      tmp_table(:,:)=timer_table(:,:,mpi_rank)
#ifdef CPP_MPI 
      call mpi_allgather(tmp_table(1,1),data_size,mpi_real_lmdz,timer_table(1,1,0),data_size,mpi_real_lmdz,comm_lmdz,ierr)
#endif
      tmp_table(:,:)=timer_table_sqr(:,:,mpi_rank)
#ifdef CPP_MPI
      call mpi_allgather(tmp_table(1,1),data_size,mpi_real_lmdz,timer_table_sqr(1,1,0),data_size,mpi_real_lmdz,comm_lmdz,ierr)
#endif    
      deallocate(tmp_table)
    
      endif
      
    ENDIF ! using_mpi
    
  end subroutine allgather_timer
  
  subroutine allgather_timer_average
    use parallel
    implicit none
#ifdef CPP_MPI
    include 'mpif.h'
#endif
    integer :: ierr
    integer :: data_size
    real, allocatable,dimension(:,:),target :: tmp_table
    integer, allocatable,dimension(:,:),target :: tmp_iter
    integer :: istats

    IF (using_mpi) THEN
        
      if (alltimer_isactive) then
    
      allocate(tmp_table(max_size,nb_timer))
      allocate(tmp_iter(max_size,nb_timer))
   
      data_size=max_size*nb_timer

      tmp_table(:,:)=timer_average(:,:,mpi_rank)
#ifdef CPP_MPI
      call mpi_allgather(tmp_table(1,1),data_size,mpi_real_lmdz,timer_average(1,1,0),data_size,mpi_real_lmdz,comm_lmdz,ierr)
#endif
      tmp_table(:,:)=timer_delta(:,:,mpi_rank)
#ifdef CPP_MPI
      call mpi_allgather(tmp_table(1,1),data_size,mpi_real_lmdz,timer_delta(1,1,0),data_size,mpi_real_lmdz,comm_lmdz,ierr)
#endif
      tmp_iter(:,:)=timer_iteration(:,:,mpi_rank)
#ifdef CPP_MPI
      call mpi_allgather(tmp_iter(1,1),data_size,mpi_integer,timer_iteration(1,1,0),data_size,mpi_integer,comm_lmdz,ierr)
#endif    
      deallocate(tmp_table)
    
      endif
      
    ENDIF  ! using_mp�
  end subroutine allgather_timer_average
  
  subroutine inittime
  implicit none
    integer :: count,count_rate,count_max
    
    alltimer_isactive=.true.
    if (alltimer_isactive) then
      call system_clock(count,count_rate,count_max)
      call cpu_time(last_cpucount)
      last_count=count
    endif
  end subroutine inittime
  
  function difftime()
  implicit none
    double precision :: difftime
    integer :: count,count_rate,count_max
  
    call system_clock(count,count_rate,count_max)
    if (count>=last_count) then
      difftime=(1.*(count-last_count))/count_rate
    else
      difftime=(1.*(count-last_count+count_max))/count_rate
    endif
    last_count=count 
  end function difftime
  
  function diffcputime()
  implicit none
    real :: diffcputime
    real :: count
    
    call cpu_time(count)
    diffcputime=count-last_cpucount
    last_cpucount=count 
  end function diffcputime

end module times