- added UNRES Multichain examples
[unres.git] / source / unres / src_MD-M / timing.F
1 C $Date: 1994/10/05 16:41:52 $
2 C $Revision: 2.2 $
3 C
4 C
5 C
6       subroutine set_timers
7 c
8       implicit none
9       double precision tcpu
10       include 'COMMON.TIME1'
11 #ifdef MP
12       include 'mpif.h'
13 #endif
14 C Diminish the assigned time limit a little so that there is some time to
15 C end a batch job
16 c     timlim=batime-150.0
17 C Calculate the initial time, if it is not zero (e.g. for the SUN).
18       stime=tcpu()
19 #ifdef MPI
20       walltime=MPI_WTIME()
21       time_reduce=0.0d0
22       time_allreduce=0.0d0
23       time_bcast=0.0d0
24       time_gather=0.0d0
25       time_sendrecv=0.0d0
26       time_scatter=0.0d0
27       time_scatter_fmat=0.0d0
28       time_scatter_ginv=0.0d0
29       time_scatter_fmatmult=0.0d0
30       time_scatter_ginvmult=0.0d0
31       time_barrier_e=0.0d0
32       time_barrier_g=0.0d0
33       time_enecalc=0.0d0
34       time_sumene=0.0d0
35       time_lagrangian=0.0d0
36       time_sumgradient=0.0d0
37       time_intcartderiv=0.0d0
38       time_inttocart=0.0d0
39       time_ginvmult=0.0d0
40       time_fricmatmult=0.0d0
41       time_cartgrad=0.0d0
42       time_bcastc=0.0d0
43       time_bcast7=0.0d0
44       time_bcastw=0.0d0
45       time_intfcart=0.0d0
46       time_vec=0.0d0
47       time_mat=0.0d0
48       time_fric=0.0d0
49       time_stoch=0.0d0
50       time_fricmatmult=0.0d0
51       time_fsample=0.0d0
52 #endif
53 cd    print *,' in SET_TIMERS stime=',stime
54       return 
55       end
56 C------------------------------------------------------------------------------
57       logical function stopx(nf)
58 C This function returns .true. if one of the following reasons to exit SUMSL
59 C occurs. The "reason" code is stored in WHATSUP passed thru a COMMON block:
60 C
61 C... WHATSUP = 0 - go on, no reason to stop. Stopx will return .false.
62 C...           1 - Time up in current node;
63 C...           2 - STOP signal was received from another node because the
64 C...               node's task was accomplished (parallel only);
65 C...          -1 - STOP signal was received from another node because of error;
66 C...          -2 - STOP signal was received from another node, because 
67 C...               the node's time was up.
68       implicit real*8 (a-h,o-z)
69       include 'DIMENSIONS'
70       integer nf
71       logical ovrtim
72 #ifdef MP
73       include 'mpif.h'
74       include 'COMMON.INFO'
75 #endif
76       include 'COMMON.IOUNITS'
77       include 'COMMON.TIME1'
78       integer Kwita
79
80 cd    print *,'Processor',MyID,' NF=',nf
81 #ifndef MPI
82       if (ovrtim()) then
83 C Finish if time is up.
84          stopx = .true.
85          WhatsUp=1
86 #ifdef MPL
87       else if (mod(nf,100).eq.0) then
88 C Other processors might have finished. Check this every 100th function 
89 C evaluation.
90 C Master checks if any other processor has sent accepted conformation(s) to it. 
91          if (MyID.ne.MasterID) call receive_mcm_info
92          if (MyID.eq.MasterID) call receive_conf
93 cd       print *,'Processor ',MyID,' is checking STOP: nf=',nf
94          call recv_stop_sig(Kwita)
95          if (Kwita.eq.-1) then
96            write (iout,'(a,i4,a,i5)') 'Processor',
97      &     MyID,' has received STOP signal in STOPX; NF=',nf
98            write (*,'(a,i4,a,i5)') 'Processor',
99      &     MyID,' has received STOP signal in STOPX; NF=',nf
100            stopx=.true.
101            WhatsUp=2
102          elseif (Kwita.eq.-2) then
103            write (iout,*)
104      &    'Processor',MyID,' received TIMEUP-STOP signal in SUMSL.'
105            write (*,*)
106      &    'Processor',MyID,' received TIMEUP-STOP signal in SUMSL.'
107            WhatsUp=-2
108            stopx=.true.  
109          else if (Kwita.eq.-3) then
110            write (iout,*)
111      &    'Processor',MyID,' received ERROR-STOP signal in SUMSL.'
112            write (*,*)
113      &    'Processor',MyID,' received ERROR-STOP signal in SUMSL.'
114            WhatsUp=-1
115            stopx=.true.
116          else
117            stopx=.false.
118            WhatsUp=0
119          endif
120 #endif
121       else
122          stopx = .false.
123          WhatsUp=0
124       endif
125 #else
126       stopx=.false.
127 #endif
128
129 #ifdef OSF
130 c Check for FOUND_NAN flag
131       if (FOUND_NAN) then
132         write(iout,*)"   ***   stopx : Found a NaN"
133         stopx=.true.
134       endif
135 #endif
136
137       return
138       end
139 C--------------------------------------------------------------------------
140       logical function ovrtim() 
141       include 'DIMENSIONS'
142       include 'COMMON.IOUNITS'
143       include 'COMMON.TIME1'
144       real*8 tcpu
145 #ifdef MPI
146       include "mpif.h"
147       curtim = MPI_Wtime()-walltime
148 #else
149       curtim= tcpu()
150 #endif
151 C  curtim is the current time in seconds.
152 c      write (iout,*) "curtim",curtim," timlim",timlim," safety",safety
153       if (curtim .ge. timlim - safety) then
154         write (iout,'(a,f10.2,a,f10.2,a,f10.2,a)') 
155      &  "***************** Elapsed time (",curtim,
156      &  " s) is within the safety limit (",safety,
157      &  " s) of the allocated time (",timlim," s). Terminating."
158         ovrtim=.true.
159       else
160         ovrtim=.false.
161       endif
162       return                                               
163       end
164 **************************************************************************      
165       double precision function tcpu()
166       include 'COMMON.TIME1'
167 #ifdef ES9000 
168 ****************************
169 C Next definition for EAGLE (ibm-es9000)
170       real*8 micseconds
171       integer rcode
172       tcpu=cputime(micseconds,rcode)
173       tcpu=(micseconds/1.0E6) - stime
174 ****************************
175 #endif
176 #ifdef SUN
177 ****************************
178 C Next definitions for sun
179       REAL*8  ECPU,ETIME,ETCPU
180       dimension tarray(2)
181       tcpu=etime(tarray)
182       tcpu=tarray(1)
183 ****************************
184 #endif
185 #ifdef KSR
186 ****************************
187 C Next definitions for ksr
188 C this function uses the ksr timer ALL_SECONDS from the PMON library to
189 C return the elapsed time in seconds
190       tcpu= all_seconds() - stime
191 ****************************
192 #endif
193 #ifdef SGI
194 ****************************
195 C Next definitions for sgi
196       real timar(2), etime
197       seconds = etime(timar)
198 Cd    print *,'seconds=',seconds,' stime=',stime
199 C      usrsec = timar(1)
200 C      syssec = timar(2)
201       tcpu=seconds - stime
202 ****************************
203 #endif
204
205 #ifdef LINUX
206 ****************************
207 C Next definitions for sgi
208       real timar(2), etime
209       seconds = etime(timar)
210 Cd    print *,'seconds=',seconds,' stime=',stime
211 C      usrsec = timar(1)
212 C      syssec = timar(2)
213       tcpu=seconds - stime
214 ****************************
215 #endif
216
217
218 #ifdef CRAY
219 ****************************
220 C Next definitions for Cray
221 C     call date(curdat)
222 C     curdat=curdat(1:9)
223 C     call clock(curtim)
224 C     curtim=curtim(1:8)
225       cpusec = second()
226       tcpu=cpusec - stime
227 ****************************
228 #endif
229 #ifdef AIX
230 ****************************
231 C Next definitions for RS6000
232        integer*4 i1,mclock
233        i1 = mclock()
234        tcpu = (i1+0.0D0)/100.0D0
235 #endif
236 #ifdef WINPGI
237 ****************************
238 c next definitions for windows NT Digital fortran
239        real time_real
240        call cpu_time(time_real)
241        tcpu = time_real
242 #endif
243 #ifdef WINIFL
244 ****************************
245 c next definitions for windows NT Digital fortran
246        real time_real
247        call cpu_time(time_real)
248        tcpu = time_real
249 #endif
250
251       return     
252       end  
253 C---------------------------------------------------------------------------
254       subroutine dajczas(rntime,hrtime,mintime,sectime)
255       include 'COMMON.IOUNITS'
256       real*8 rntime,hrtime,mintime,sectime 
257       hrtime=rntime/3600.0D0 
258       hrtime=aint(hrtime)
259       mintime=aint((rntime-3600.0D0*hrtime)/60.0D0)
260       sectime=aint((rntime-3600.0D0*hrtime-60.0D0*mintime)+0.5D0)
261       if (sectime.eq.60.0D0) then
262         sectime=0.0D0
263         mintime=mintime+1.0D0
264       endif
265       ihr=hrtime
266       imn=mintime
267       isc=sectime
268       write (iout,328) ihr,imn,isc
269   328 FORMAT(//'***** Computation time: ',I4  ,' hours ',I2  ,
270      1         ' minutes ', I2  ,' seconds *****')       
271       return
272       end
273 C---------------------------------------------------------------------------
274       subroutine print_detailed_timing
275       implicit real*8 (a-h,o-z)
276       include 'DIMENSIONS'
277 #ifdef MPI
278       include 'mpif.h'
279 #endif
280       include 'COMMON.IOUNITS'
281       include 'COMMON.TIME1'
282       include 'COMMON.SETUP'
283       time1=MPI_WTIME()
284          write (iout,'(80(1h=)/a/(80(1h=)))') 
285      &    "Details of FG communication time"
286          write (*,'(7(a40,1pe15.5/),40(1h-)/a40,1pe15.5/80(1h=))') 
287      &    "BROADCAST:",time_bcast,"REDUCE:",time_reduce,
288      &    "GATHER:",time_gather,
289      &    "SCATTER:",time_scatter,"SENDRECV:",time_sendrecv,
290      &    "BARRIER ene",time_barrier_e,
291      &    "BARRIER grad",time_barrier_g,
292      &    "TOTAL:",
293      &    time_bcast+time_reduce+time_gather+time_scatter+time_sendrecv
294          write (*,*) fg_rank,myrank,
295      &     ': Total wall clock time',time1-walltime,' sec'
296          write (*,*) "Processor",fg_rank,myrank,
297      &     ": BROADCAST time",time_bcast," REDUCE time",
298      &      time_reduce," GATHER time",time_gather," SCATTER time",
299      &      time_scatter,
300      &     " SCATTER fmatmult",time_scatter_fmatmult,
301      &     " SCATTER ginvmult",time_scatter_ginvmult,
302      &     " SCATTER fmat",time_scatter_fmat,
303      &     " SCATTER ginv",time_scatter_ginv,
304      &      " SENDRECV",time_sendrecv,
305      &      " BARRIER ene",time_barrier_e,
306      &      " BARRIER GRAD",time_barrier_g,
307      &      " BCAST7",time_bcast7," BCASTC",time_bcastc,
308      &      " BCASTW",time_bcastw," ALLREDUCE",time_allreduce,
309      &      " TOTAL",
310      &      time_bcast+time_reduce+time_gather+time_scatter+
311      &      time_sendrecv+time_barrier+time_bcastc
312          write (*,*) "Processor",fg_rank,myrank," enecalc",time_enecalc
313          write (*,*) "Processor",fg_rank,myrank," sumene",time_sumene
314          write (*,*) "Processor",fg_rank,myrank," intfromcart",
315      &     time_intfcart
316          write (*,*) "Processor",fg_rank,myrank," vecandderiv",
317      &     time_vec
318          write (*,*) "Processor",fg_rank,myrank," setmatrices",
319      &     time_mat
320          write (*,*) "Processor",fg_rank,myrank," ginvmult",
321      &     time_ginvmult
322          write (*,*) "Processor",fg_rank,myrank," fricmatmult",
323      &     time_fricmatmult
324          write (*,*) "Processor",fg_rank,myrank," inttocart",
325      &     time_inttocart
326          write (*,*) "Processor",fg_rank,myrank," sumgradient",
327      &     time_sumgradient
328          write (*,*) "Processor",fg_rank,myrank," intcartderiv",
329      &     time_intcartderiv
330          if (fg_rank.eq.0) then
331            write (*,*) "Processor",fg_rank,myrank," lagrangian",
332      &       time_lagrangian
333            write (*,*) "Processor",fg_rank,myrank," cartgrad",
334      &       time_cartgrad
335          endif
336       return
337       end