source: LMDZ5/trunk/libf/dyn3dmem/times.F90 @ 1646

Last change on this file since 1646 was 1632, checked in by Laurent Fairhead, 12 years ago

Import initial du répertoire dyn3dmem

Attention! ceci n'est qu'une version préliminaire du code "basse mémoire":
le code contenu dans ce répertoire est basé sur la r1320 et a donc besoin
d'être mis à jour par rapport à la dynamique parallèle d'aujourd'hui.
Ce code est toutefois mis à disposition pour circonvenir à des problèmes
de mémoire que certaines configurations du modèle pourraient rencontrer.
Dans l'état, il compile et tourne sur vargas et au CCRT

Initial import of dyn3dmem

Warning! this is just a preliminary version of the memory light code:
it is based on r1320 of the code and thus needs to be updated before
it can replace the present dyn3dpar code. It is nevertheless put at your
disposal to circumvent some memory problems some LMDZ configurations may
encounter. In its present state, it will compile and run on vargas and CCRT

File size: 7.1 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#ifdef CPP_MPI   
141    include 'mpif.h'
142#endif
143    integer :: ierr
144    integer :: data_size
145    real, allocatable,dimension(:,:) :: tmp_table
146
147    IF (using_mpi) THEN   
148   
149      if (AllTimer_IsActive) then
150   
151   
152      allocate(tmp_table(max_size,nb_timer))
153   
154      data_size=max_size*nb_timer
155   
156      tmp_table(:,:)=timer_table(:,:,mpi_rank)
157#ifdef CPP_MPI
158      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)
159#endif
160      tmp_table(:,:)=timer_table_sqr(:,:,mpi_rank)
161#ifdef CPP_MPI
162      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)
163#endif   
164      deallocate(tmp_table)
165   
166      endif
167     
168    ENDIF ! using_mpi
169   
170  end subroutine allgather_timer
171 
172  subroutine allgather_timer_average
173    use parallel
174    implicit none
175#ifdef CPP_MPI
176    include 'mpif.h'
177#endif
178    integer :: ierr
179    integer :: data_size
180    real, allocatable,dimension(:,:),target :: tmp_table
181    integer, allocatable,dimension(:,:),target :: tmp_iter
182    integer :: istats
183
184    IF (using_mpi) THEN
185       
186      if (AllTimer_IsActive) then
187   
188      allocate(tmp_table(max_size,nb_timer))
189      allocate(tmp_iter(max_size,nb_timer))
190   
191      data_size=max_size*nb_timer
192
193      tmp_table(:,:)=timer_average(:,:,mpi_rank)
194#ifdef CPP_MPI
195      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)
196#endif
197      tmp_table(:,:)=timer_delta(:,:,mpi_rank)
198#ifdef CPP_MPI
199      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)
200#endif
201      tmp_iter(:,:)=timer_iteration(:,:,mpi_rank)
202#ifdef CPP_MPI
203      call mpi_allgather(tmp_iter(1,1),data_size,MPI_INTEGER,timer_iteration(1,1,0),data_size,MPI_INTEGER,COMM_LMDZ,ierr)
204#endif   
205      deallocate(tmp_table)
206   
207      endif
208     
209    ENDIF  ! using_mp�
210  end subroutine allgather_timer_average
211 
212  subroutine InitTime
213  implicit none
214    integer :: count,count_rate,count_max
215   
216    AllTimer_IsActive=.TRUE.
217    if (AllTimer_IsActive) then
218      call system_clock(count,count_rate,count_max)
219      call cpu_time(Last_cpuCount)
220      Last_Count=count
221    endif
222  end subroutine InitTime
223 
224  function DiffTime()
225  implicit none
226    double precision :: DiffTime
227    integer :: count,count_rate,count_max
228 
229    call system_clock(count,count_rate,count_max)
230    if (Count>=Last_Count) then
231      DiffTime=(1.*(Count-last_Count))/count_rate
232    else
233      DiffTime=(1.*(Count-last_Count+Count_max))/count_rate
234    endif
235    Last_Count=Count
236  end function DiffTime
237 
238  function DiffCpuTime()
239  implicit none
240    real :: DiffCpuTime
241    real :: Count
242   
243    call cpu_time(Count)
244    DiffCpuTime=Count-Last_cpuCount
245    Last_cpuCount=Count
246  end function DiffCpuTime
247
248end module times
Note: See TracBrowser for help on using the repository browser.