source: LMDZ4/trunk/libf/dyn3dpar/times.F90 @ 801

Last change on this file since 801 was 792, checked in by Laurent Fairhead, 17 years ago

Modifications suite a la transformation des fichiers include pour
qu'ils soient compatibles a la fois au format fixe et au format libre
Un bon nombre de fichier *.inc du coup disparaissent
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.8 KB
Line 
1module times
2  integer,private,save :: Last_Count=0
3  real, private,save :: Last_cpuCount=0
4  logical, private,save :: AllTimer_IsActive=.false.
5 
6  integer, parameter :: nb_timer = 4
7  integer, parameter :: timer_caldyn  = 1
8  integer, parameter :: timer_vanleer = 2
9  integer, parameter :: timer_dissip = 3
10  integer, parameter :: timer_physic = 4
11  integer, parameter :: stopped = 1
12  integer, parameter :: running = 2
13  integer, parameter :: suspended = 3
14 
15  integer :: max_size
16  real,    allocatable, dimension(:,:,:) :: timer_table
17  real,    allocatable, dimension(:,:,:) :: timer_table_sqr
18  integer, allocatable, dimension(:,:,:) :: timer_iteration
19  real,    allocatable, dimension(:,:,:) :: timer_average
20  real,    allocatable, dimension(:,:,:) :: timer_delta
21  real,    allocatable,dimension(:) :: timer_running, last_time
22  integer, allocatable,dimension(:) :: timer_state
23 
24  contains
25 
26  subroutine init_timer
27    use parallel
28    implicit none
29#include "dimensions.h"
30#include "paramet.h"
31   
32    max_size=jjm+1
33    allocate(timer_table(max_size,nb_timer,0:mpi_size-1))
34    allocate(timer_table_sqr(max_size,nb_timer,0:mpi_size-1))
35    allocate(timer_iteration(max_size,nb_timer,0:mpi_size-1))
36    allocate(timer_average(max_size,nb_timer,0:mpi_size-1))
37    allocate(timer_delta(max_size,nb_timer,0:mpi_size-1))
38    allocate(timer_running(nb_timer))
39    allocate(timer_state(nb_timer))
40    allocate(last_time(nb_timer))
41   
42    timer_table(:,:,:)=0
43    timer_table_sqr(:,:,:)=0
44    timer_iteration(:,:,:)=0
45    timer_average(:,:,:)=0
46    timer_delta(:,:,:)=0
47    timer_state(:)=stopped     
48  end subroutine init_timer
49 
50  subroutine start_timer(no_timer)
51    implicit none
52    integer :: no_timer
53   
54    if (AllTimer_IsActive) then
55   
56      if (timer_state(no_timer)/=stopped) then
57        stop 'start_timer :: timer is already running or suspended'
58      else
59        timer_state(no_timer)=running
60      endif
61     
62      timer_running(no_timer)=0
63      call cpu_time(last_time(no_timer))
64   
65    endif
66   
67  end subroutine start_timer
68 
69  subroutine suspend_timer(no_timer)
70    implicit none
71    integer :: no_timer
72     
73    if (AllTimer_IsActive) then   
74      if (timer_state(no_timer)/=running) then
75        stop 'suspend_timer :: timer is not running'
76      else
77        timer_state(no_timer)=suspended
78      endif
79   
80      timer_running(no_timer)=timer_running(no_timer)-last_time(no_timer)
81      call cpu_time(last_time(no_timer))
82      timer_running(no_timer)=timer_running(no_timer)+last_time(no_timer)
83    endif
84  end subroutine suspend_timer
85 
86  subroutine resume_timer(no_timer)
87    implicit none
88    integer :: no_timer
89     
90    if (AllTimer_IsActive) then   
91      if (timer_state(no_timer)/=suspended) then
92        stop 'resume_timer :: timer is not suspended'
93      else
94        timer_state(no_timer)=running
95      endif
96     
97      call cpu_time(last_time(no_timer))
98    endif
99   
100  end subroutine resume_timer
101
102  subroutine stop_timer(no_timer)
103    use parallel
104    implicit none
105    integer :: no_timer
106    integer :: N
107    real :: V,V2
108   
109    if (AllTimer_IsActive) then
110       
111      if (timer_state(no_timer)/=running) then
112        stop 'stop_timer :: timer is not running'
113      else
114        timer_state(no_timer)=stopped
115      endif
116   
117      timer_running(no_timer)=timer_running(no_timer)-last_time(no_timer)
118      call cpu_time(last_time(no_timer))
119      timer_running(no_timer)=timer_running(no_timer)+last_time(no_timer)
120   
121      timer_table(jj_nb,no_timer,mpi_rank)=timer_table(jj_nb,no_timer,mpi_rank)+timer_running(no_timer)
122      timer_table_sqr(jj_nb,no_timer,mpi_rank)=timer_table_sqr(jj_nb,no_timer,mpi_rank)+timer_running(no_timer)**2
123      timer_iteration(jj_nb,no_timer,mpi_rank)=timer_iteration(jj_nb,no_timer,mpi_rank)+1
124      timer_average(jj_nb,no_timer,mpi_rank)=timer_table(jj_nb,no_timer,mpi_rank)/timer_iteration(jj_nb,no_timer,mpi_rank)
125      if (timer_iteration(jj_nb,no_timer,mpi_rank)>=2) then
126        N=timer_iteration(jj_nb,no_timer,mpi_rank)
127        V2=timer_table_sqr(jj_nb,no_timer,mpi_rank)
128        V=timer_table(jj_nb,no_timer,mpi_rank)
129        timer_delta(jj_nb,no_timer,mpi_rank)=sqrt(ABS(V2-V*V/N)/(N-1))
130      else
131        timer_delta(jj_nb,no_timer,mpi_rank)=0
132      endif
133    endif
134   
135  end subroutine stop_timer
136   
137  subroutine allgather_timer
138    use parallel
139    implicit none
140    include 'mpif.h'
141    integer :: ierr
142    integer :: data_size
143    real, allocatable,dimension(:,:) :: tmp_table
144   
145    if (AllTimer_IsActive) then
146   
147    allocate(tmp_table(max_size,nb_timer))
148   
149    data_size=max_size*nb_timer
150   
151    tmp_table(:,:)=timer_table(:,:,mpi_rank)
152    call mpi_allgather(tmp_table(1,1),data_size,MPI_REAL8,timer_table(1,1,mpi_rank),data_size,MPI_REAL8,COMM_LMDZ,ierr)
153
154    tmp_table(:,:)=timer_table_sqr(:,:,mpi_rank)
155    call mpi_allgather(tmp_table(1,1),data_size,MPI_REAL8,timer_table_sqr(1,1,mpi_rank),data_size,MPI_REAL8,COMM_LMDZ,ierr)
156   
157    deallocate(tmp_table)
158   
159    endif
160   
161  end subroutine allgather_timer
162 
163  subroutine allgather_timer_average
164    use parallel
165    implicit none
166    include 'mpif.h'
167    integer :: ierr
168    integer :: data_size
169    real, allocatable,dimension(:,:),target :: tmp_table
170    integer, allocatable,dimension(:,:),target :: tmp_iter
171    integer :: istats
172   
173    if (AllTimer_IsActive) then
174   
175    allocate(tmp_table(max_size,nb_timer))
176    allocate(tmp_iter(max_size,nb_timer))
177   
178    data_size=max_size*nb_timer
179
180    tmp_table(:,:)=timer_average(:,:,mpi_rank)
181    call mpi_allgather(tmp_table(1,1),data_size,MPI_REAL8,timer_average(1,1,0),data_size,MPI_REAL8,COMM_LMDZ,ierr)
182
183    tmp_table(:,:)=timer_delta(:,:,mpi_rank)
184    call mpi_allgather(tmp_table(1,1),data_size,MPI_REAL8,timer_delta(1,1,0),data_size,MPI_REAL8,COMM_LMDZ,ierr)
185
186    tmp_iter(:,:)=timer_iteration(:,:,mpi_rank)
187    call mpi_allgather(tmp_iter(1,1),data_size,MPI_INTEGER,timer_iteration(1,1,0),data_size,MPI_INTEGER,COMM_LMDZ,ierr)
188   
189    deallocate(tmp_table)
190   
191    endif
192  end subroutine allgather_timer_average
193 
194  subroutine InitTime
195  implicit none
196    integer :: count,count_rate,count_max
197   
198    AllTimer_IsActive=.TRUE.
199    if (AllTimer_IsActive) then
200      call system_clock(count,count_rate,count_max)
201      call cpu_time(Last_cpuCount)
202      Last_Count=count
203    endif
204  end subroutine InitTime
205 
206  function DiffTime
207  implicit none
208    double precision :: DiffTime
209    integer :: count,count_rate,count_max
210 
211    call system_clock(count,count_rate,count_max)
212    if (Count>=Last_Count) then
213      DiffTime=(1.*(Count-last_Count))/count_rate
214    else
215      DiffTime=(1.*(Count-last_Count+Count_max))/count_rate
216    endif
217    Last_Count=Count
218  end function DiffTime
219 
220  function DiffCpuTime
221  implicit none
222    real :: DiffCpuTime
223    real :: Count
224   
225    call cpu_time(Count)
226    DiffCpuTime=Count-Last_cpuCount
227    Last_cpuCount=Count
228  end function DiffCpuTime
229
230end module times
Note: See TracBrowser for help on using the repository browser.