GCC Code Coverage Report


Directory: ./
File: rad/gstats.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 0 147 0.0%
Branches: 0 122 0.0%

Line Branch Exec Source
1 SUBROUTINE GSTATS(KNUM,KSWITCH)
2
3 !**** *GSTATS* - Gather timing statistics
4
5 ! PURPOSE.
6 ! --------
7 ! To gather timings for subsequent output by routine STATS_OUTPUT
8
9
10 !** INTERFACE.
11 ! ----------
12 ! *CALL* *GSTATS(KNUM,KSWITCH)
13
14 ! EXPLICIT ARGUMENTS
15 ! --------------------
16 ! KNUM - timing event number (for list of already defined events
17 ! see routine STATS_OUTPUT)
18 ! KSWITCH - KSWITCH=0 - switch on timer
19 ! KSWITCH=1 - switch off timer
20 ! KSWITCH=2 - suspend timer
21 ! KSWITCH=3 - resume timer
22
23 ! IMPLICIT ARGUMENTS
24 ! --------------------
25 ! Module YOMSTATS
26
27 ! METHOD.
28 ! -------
29
30
31 ! EXTERNALS. USER_CLOCK - timing routine
32 ! ---------- MPL_BARRIER - syncronization of processors
33
34 ! REFERENCE.
35 ! ----------
36 ! ECMWF Research Department documentation of the IFS
37
38 ! AUTHOR.
39 ! -------
40 ! Mats Hamrud ECMWF
41
42 ! MODIFICATIONS.
43 ! --------------
44 ! ORIGINAL : 98-11-15
45 ! D.Salmond: 02-02-25 Return if not master thread when called from a
46 ! parallel region.
47 ! J.Hague: 03-06-11 Memory tracing (for NSTATS_MEM MPI tasks)
48 ! ------------------------------------------------------------------
49
50 USE PARKIND1 ,ONLY : JPIM ,JPRB ,JPIB
51 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK
52
53 USE YOMGSTATS
54 !USE MPL_MODULE ! MPL 28.11.08
55 USE YOMOML
56
57 IMPLICIT NONE
58
59 INTEGER(KIND=JPIM),INTENT(IN) :: KNUM
60 INTEGER(KIND=JPIM),INTENT(IN) :: KSWITCH
61
62 INTEGER(KIND=JPIM) :: IMOD,ICALL
63 INTEGER(KIND=JPIM) :: IIMEM, IIPAG, IIMEMC
64 INTEGER(KIND=JPIB) :: IMEM, IMEMH, IMEMS, IMEMC, IPAG, INUM
65 INTEGER(KIND=JPIB) :: GETRSS, GETHWM, GETSTK, GETCURHEAP, GETPAG
66 EXTERNAL GETRSS, GETHWM, GETSTK, GETCURHEAP, GETPAG
67 REAL(KIND=JPRB) :: ZTIMED,ZCLOCK,ZTIME,ZTCPU,ZVCPU
68 REAL(KIND=JPRB) :: ZLAST_PAR_TIME
69 LOGICAL :: LLFIRST=.TRUE.
70 LOGICAL :: LLMFIRST=.TRUE.
71 CHARACTER(LEN=32), SAVE :: CCDESC_DRHOOK(JPMAXSTAT)
72 CHARACTER(LEN=32), SAVE :: CCDESC_BARR(JPMAXSTAT)
73 SAVE ZLAST_PAR_TIME
74 SAVE IIMEM, IIPAG, IIMEMC
75
76 INTEGER(KIND=JPIM),SAVE :: NUM_THREADS
77 REAL(KIND=JPRB),ALLOCATABLE,SAVE :: ZHOOK_HANDLE(:)
78 REAL(KIND=JPRB),SAVE :: ZHOOK_HANDLE_COMMS, ZHOOK_HANDLE_COMMS1
79 REAL(KIND=JPRB),SAVE :: ZHOOK_HANDLE_TRANS
80 REAL(KIND=JPRB),SAVE :: ZHOOK_HANDLE_BARR
81 CHARACTER*4 CC
82
83 INTEGER(KIND=JPIM) :: NMAX_STATS, KULNAM
84
85 INTERFACE
86 SUBROUTINE USER_CLOCK(PELAPSED_TIME,PELAPSED_TIME_SINCE,PVECTOR_CP,PTOTAL_CP)
87
88 !**** *USER_CLOCK* - interface to system dependent timer routines
89
90 ! Purpose.
91 ! --------
92 ! Returns elapsed and CP from the start of execution.
93 ! Elapsed time is made relative to the first call to USER_CLOCK.
94
95 !** Interface.
96 ! ----------
97 ! ZTIME=USER_CLOCK(PELAPSED_TIME,PELAPSED_TIME_SINCE,
98 ! PVECTOR_CP,PTOTAL_CP)
99
100 ! Explicit arguments: (All are optional arguments)
101 ! PELAPSED_TIME=wall clock time (seconds)
102 ! PELAPSED_TIME_SINCE=wall clock time (seconds)
103 ! change from input value of this parameter
104 ! PVECTOR_CP=CP vector time (seconds)
105 ! PTOTAL_CP=total CP time (seconds)
106
107 ! Author.
108 ! -------
109 ! D.Dent *ECMWF*
110
111 ! External References:
112 ! -------------------
113
114 ! TIMEF,CPTIME
115
116 ! Modifications.
117 ! --------------
118 ! Original : 97-09-25
119 ! ----------------------------------------------------------
120
121
122 USE PARKIND1 ,ONLY : JPIM ,JPRB
123
124 IMPLICIT NONE
125
126 REAL(KIND=JPRB),INTENT(OUT) :: PELAPSED_TIME,PVECTOR_CP,PTOTAL_CP
127 REAL(KIND=JPRB),INTENT(INOUT) :: PELAPSED_TIME_SINCE
128 OPTIONAL PELAPSED_TIME,PELAPSED_TIME_SINCE
129 OPTIONAL PVECTOR_CP,PTOTAL_CP
130 REAL(KIND=JPRB) :: ZVECTOR_CP,ZTOTAL_CP,ZWALL
131 REAL(KIND=JPRB),EXTERNAL :: TIMEF
132
133
134 END SUBROUTINE USER_CLOCK
135 END INTERFACE
136
137 ! write(0,*) "GSTATS:LSTATS,JPMAXSTAT,LGSTATS_LABEL,KNUM=",LSTATS,JPMAXSTAT,LGSTATS_LABEL,KNUM
138
139 IF(LSTATS) THEN
140
141 IF(.NOT.ALLOCATED(ZHOOK_HANDLE))THEN
142 NUM_THREADS=OML_MAX_THREADS()
143 ALLOCATE(ZHOOK_HANDLE(NUM_THREADS))
144 ENDIF
145
146 IF(LGSTATS_LABEL)THEN
147 DO INUM=1,JPMAXSTAT
148 WRITE(CC,'(I4)')INUM
149 CCDESC_BARR(INUM)='>BAR-'//CCDESC(INUM)(1:21)//'('//CC//')'
150 ENDDO
151 DO INUM=1,JPMAXSTAT
152 WRITE(CC,'(I4)')INUM
153 ! write(6,*) inum,cctype(inum)
154 IF(CCTYPE(INUM).EQ."TRS".OR.CCTYPE(INUM).EQ.'MP-'.OR.CCTYPE(INUM).EQ.'MPL'&
155 & .OR.CCTYPE(INUM).EQ.'BAR'.OR.CCTYPE(INUM).EQ.'OMP') THEN
156 CCDESC_DRHOOK(INUM)='>'//CCTYPE(INUM)//'-'//CCDESC(INUM)(1:21)//'('//CC//')'
157 ENDIF
158 ENDDO
159 LGSTATS_LABEL=.FALSE.
160 ENDIF
161 ! ------------------------------------------------------------------
162
163 !J IF((KNUM > 1000 .AND.KNUM < 2001).AND.(.NOT.LSTATS_OMP))GOTO 99999
164 !J IF((KNUM > 500 .AND.KNUM < 1001).AND.(.NOT.LSTATS_COMMS))GOTO 99999
165 IF(OML_MY_THREAD() > 1)GOTO 99999
166
167 !J IF(KNUM/=0) THEN
168 !J IF(LSYNCSTATS .AND.(KSWITCH==0.OR. KSWITCH==2)) THEN
169 !J IF(.NOT.OML_IN_PARALLEL().AND. KNUM < 500 )THEN
170 !J IF(LHOOK)CALL DR_HOOK(CCDESC_BARR(KNUM),0,ZHOOK_HANDLE_BARR)
171 !J CALL MPL_BARRIER(CDSTRING='GSTATS:')
172 !J IF(LHOOK)CALL DR_HOOK(CCDESC_BARR(KNUM),1,ZHOOK_HANDLE_BARR)
173 !J ENDIF
174 !J ENDIF
175 !J ENDIF
176
177 IF (LHOOK .AND. (KSWITCH == 0 .OR. KSWITCH == 1)) THEN
178 ! write(0,*) "KNUM,SWITCH=",KNUM,KSWITCH
179 ! write(0,*) "CCTYPE=",CCTYPE(KNUM)
180 ! write(0,*) "CCDESC_DRHOOK=",CCDESC_DRHOOK(KNUM)
181 IF(CCTYPE(KNUM).EQ."TRS")THEN
182 CALL DR_HOOK(CCDESC_DRHOOK(KNUM),KSWITCH,ZHOOK_HANDLE_TRANS)
183 ELSEIF(CCTYPE(KNUM).EQ.'MP-')THEN
184 CALL DR_HOOK(CCDESC_DRHOOK(KNUM),KSWITCH,ZHOOK_HANDLE_COMMS)
185 ELSEIF(CCTYPE(KNUM).EQ.'MPL'.AND.KNUM.NE.682)THEN
186 CALL DR_HOOK(CCDESC_DRHOOK(KNUM),KSWITCH,ZHOOK_HANDLE_COMMS1)
187 ELSEIF(CCTYPE(KNUM).EQ.'OMP')THEN
188 ! CALL DR_HOOK(CCDESC_DRHOOK(KNUM),KSWITCH,ZHOOK_HANDLE)
189 ! MPL 1.12.08 : erreur a la compilation car zhook_handle est un tableau et non un
190 ! scalaire (PKEY)
191 CALL DR_HOOK(CCDESC_DRHOOK(KNUM),KSWITCH,ZTIME)
192 ELSEIF(CCTYPE(KNUM).EQ.'BAR')THEN
193 CALL DR_HOOK(CCDESC_DRHOOK(KNUM),KSWITCH,ZHOOK_HANDLE_BARR)
194 ENDIF
195 ENDIF
196
197 CALL USER_CLOCK(PELAPSED_TIME=ZCLOCK)
198 IF (LSTATSCPU.OR.KNUM==0) THEN
199 CALL USER_CLOCK(PTOTAL_CP=ZTCPU,PVECTOR_CP=ZVCPU)
200 ELSE
201 ZTCPU = 0.0_JPRB
202 ZVCPU = 0.0_JPRB
203 ENDIF
204
205 IF (LLFIRST) THEN
206
207 ! write(0,*) "JPMAXSTAT:2=",JPMAXSTAT
208
209 NCALLS(:) = 0
210 TIMESUM(:) = 0.0_JPRB
211 TIMESQSUM(:) = 0.0_JPRB
212 TIMEMAX(:) = 0.0_JPRB
213 TIMESUMB(:) = 0.0_JPRB
214 TTCPUSUM(:) = 0.0_JPRB
215 TVCPUSUM(:) = 0.0_JPRB
216 TIMELCALL(:) = ZCLOCK
217 CCDESC=""
218 CCTYPE=""
219 NTMEM = 0
220 NTMEM(:,5) = 99999999
221 IIMEM=0
222 IIPAG=0
223 IIMEMC=0
224 TIME_LAST_CALL = ZCLOCK
225 ZLAST_PAR_TIME=ZCLOCK
226 LLFIRST = .FALSE.
227 ENDIF
228
229 IF(KNUM < 0.OR. KNUM > JPMAXSTAT) CALL ABOR1('GSTATS')
230 IF(KSWITCH == 0.OR. KSWITCH == 1) THEN
231 NCALLS(KNUM) = NCALLS(KNUM)+1
232 ENDIF
233 IMOD = MOD(NCALLS(KNUM),2)
234 IF(.NOT.((KSWITCH == 0.AND. IMOD == 1) .OR.&
235 &(KSWITCH == 2.AND. IMOD == 1) .OR.&
236 &(KSWITCH == 3.AND. IMOD == 1) .OR.&
237 &(KSWITCH == 1.AND. IMOD == 0))) THEN
238 WRITE(JPERR,*) 'KNUM,KSWITCH,IMOD,NCALLS(KNUM)',&
239 &KNUM,KSWITCH,IMOD,NCALLS(KNUM)
240 CALL ABOR1('GSTATS')
241 ENDIF
242
243 IF( KSWITCH == 0 ) THEN
244 ! Start timing event
245 IF(KNUM < 500) THEN
246 ZTIMED = ZCLOCK-TIME_LAST_CALL
247 ELSE
248 ZTIMED = ZCLOCK - ZLAST_PAR_TIME
249 ENDIF
250 TIMESUMB(KNUM) = TIMESUMB(KNUM)+ZTIMED
251 THISTIME(KNUM) = 0.0_JPRB
252 TIMELCALL(KNUM) = ZCLOCK
253 TTCPULCALL(KNUM) = ZTCPU
254 TVCPULCALL(KNUM) = ZVCPU
255 THISTCPU(KNUM) = 0.0_JPRB
256 THISVCPU(KNUM) = 0.0_JPRB
257 IF(MYPROC_STATS.LE.NSTATS_MEM.AND.MYPROC_STATS.NE.0) THEN
258 ! CALL getrss(IMEM)
259 ! IMEM = getrss()/1024
260 IPAG = getpag()
261 IMEMH = gethwm()/1024
262 IMEMS = getstk()/1024
263 IMEMC = 0
264 IF(LSTATS_ALLOC) IMEMC = GETCURHEAP()/1024
265 IF(IMEM > IIMEM.OR.IPAG > IIPAG.OR.(LSTATS_ALLOC.AND.(IMEMC.NE.IIMEMC))) THEN
266 IF(LLMFIRST) THEN
267 WRITE(0,*) ".---------------------------------------------------------"
268 WRITE(0,*) "| Memory trace details"
269 WRITE(0,*) "| --------------------"
270 WRITE(0,*) "| Memory examined at each GSTATS call if NSTATS_MEM>0."
271 WRITE(0,*) "| Header for each trace line is:"
272 WRITE(0,*) "|"
273 WRITE(0,*) "| RSS_INC: Increase in RSS_MAX (KB)"
274 WRITE(0,*) "| RSS_MAX: Maximum real working set so far (KB)"
275 WRITE(0,*) "| HEAP_MX: High Water Mark for heap so far (KB)"
276 WRITE(0,*) "| STK: Current Stack usage (KB)"
277 WRITE(0,*) "| PGS: Page faults w I/O since last trace line"
278 WRITE(0,*) "| CALL: Number of gstats call"
279 WRITE(0,*) "| HEAP: Current malloc'd total (KB)"
280 WRITE(0,*) "|"
281 WRITE(0,*) "| Trace line written for NSTATS_MEM MPI tasks if RSS_MAX"
282 WRITE(0,*) "| RSS_MAX increases, PGS>0, or HEAP changed"
283 WRITE(0,*) "| (if LTATS_ALLOC=.TRUE.)"
284 WRITE(0,*) "`---------------------------------------------------------"
285 WRITE(0,*) ""
286 WRITE(0,'(A10,A5,21X,A7,2A8,A7,A5,A5,A8)') &
287 & "MEMORY "," KNUM","RSS_INC"," RSS_MAX"," HEAP_MX"," STK", &
288 & " PGS"," CALL"," HEAP"
289 LLMFIRST=.FALSE.
290 ENDIF
291 WRITE(0,'(A10,I5,1X,A20,1X,I6,2(1X,I7),1X,I6,1X,I4,1X,I4,1X,I7)') &
292 & "MEMORY bfr",KNUM,CCDESC(KNUM),IMEM-IIMEM,IMEM,IMEMH,IMEMS, &
293 & IPAG-IIPAG,(NCALLS(KNUM)+1)/2,IMEMC
294 ENDIF
295 NTMEM(KNUM,2)=IMEM
296 IIMEM=IMEM
297 IIPAG=IPAG
298 IIMEMC=IMEMC
299 ENDIF
300 ELSEIF( KSWITCH == 1 ) THEN
301 ! Finish timing event
302 ZTIME = THISTIME(KNUM)+(ZCLOCK-TIMELCALL(KNUM))
303 TIMESUM(KNUM) = TIMESUM(KNUM)+ZTIME
304 TIMESQSUM(KNUM) = TIMESQSUM(KNUM)+ZTIME**2
305 TIMEMAX(KNUM) = MAX(TIMEMAX(KNUM),ZTIME)
306 TTCPUSUM(KNUM) = TTCPUSUM(KNUM)+THISTCPU(KNUM)+ZTCPU-TTCPULCALL(KNUM)
307 TVCPUSUM(KNUM) = TVCPUSUM(KNUM)+THISVCPU(KNUM)+ZVCPU-TVCPULCALL(KNUM)
308 IF(MYPROC_STATS.LE.NSTATS_MEM.AND.MYPROC_STATS.NE.0) THEN
309 ! CALL getrss(IMEM)
310 IMEM = GETRSS()/1024
311 IPAG = GETPAG()
312 IMEMH = GETHWM()/1024
313 IMEMS = GETSTK()/1024
314 IMEMC = 0
315 IF(LSTATS_ALLOC) IMEMC = GETCURHEAP()/1024
316 IF(IMEM > IIMEM.OR.IPAG > IIPAG.OR.(LSTATS_ALLOC.AND.(IMEMC.NE.IIMEMC))) THEN
317 WRITE(0,'(A10,I5,1X,A20,1X,I6,2(1X,I7),1X,I6,1X,I4,1X,I4,1X,I7)') &
318 & "MEMORY aft ",KNUM,CCDESC(KNUM),IMEM-IIMEM,IMEM,IMEMH,IMEMS, &
319 & IPAG-IIPAG,NCALLS(KNUM)/2,IMEMC
320 ENDIF
321 IIMEM=IMEM
322 IIPAG=IPAG
323 IIMEMC=IMEMC
324 IMEM=IMEM-NTMEM(KNUM,2)
325 NTMEM(KNUM,4)=NTMEM(KNUM,4)+IMEM
326 IF(IMEM > NTMEM(KNUM,1)) THEN
327 NTMEM(KNUM,1)=IMEM
328 NTMEM(KNUM,3)=NCALLS(KNUM)
329 ENDIF
330 IF(IMEM < NTMEM(KNUM,5)) NTMEM(KNUM,5)=IMEM
331 ENDIF
332 ELSEIF( KSWITCH == 2 ) THEN
333 ! Suspend timing event
334 ZTIMED = ZCLOCK-TIMELCALL(KNUM)
335 THISTIME(KNUM) = THISTIME(KNUM)+ZTIMED
336 THISTCPU(KNUM) = THISTCPU(KNUM)+ZTCPU-TTCPULCALL(KNUM)
337 THISVCPU(KNUM) = THISVCPU(KNUM)+ZVCPU-TVCPULCALL(KNUM)
338 ELSEIF( KSWITCH == 3 ) THEN
339 ! Resume timing event
340 TIMELCALL(KNUM) = ZCLOCK
341 TTCPULCALL(KNUM) = ZTCPU
342 TVCPULCALL(KNUM) = ZVCPU
343 ENDIF
344 TIME_LAST_CALL = ZCLOCK
345 IF(KNUM > 500.OR.KNUM == 102.OR.KNUM == 103) ZLAST_PAR_TIME = ZCLOCK
346 ! Trace stats
347 NCALLS_TOTAL = NCALLS_TOTAL+1
348 IF (LTRACE_STATS .AND. NCALLS_TOTAL <= NTRACE_STATS) THEN
349 ICALL = NCALLS_TOTAL
350 TIME_TRACE(ICALL) = ZCLOCK
351 NCALL_TRACE(ICALL) = (JPMAXSTAT+1)*KSWITCH+KNUM
352 ENDIF
353
354 ENDIF
355
356 99999 CONTINUE
357 END SUBROUTINE GSTATS
358