correction to the last commit - cleaning of homo restraints
[unres.git] / source / cluster / wham / src-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       double precision tcpu    ! function
9       include 'COMMON.TIME1'
10 C Diminish the assigned time limit a little so that there is some time to
11 C end a batch job
12 c     timlim=batime-150.0
13 C Calculate the initial time, if it is not zero (e.g. for the SUN).
14       stime=tcpu()
15       return 
16       end
17       logical function stopx(nf)
18 C
19 C     ..................................................................
20 C
21 C     *****PURPOSE...
22 C     THIS FUNCTION MAY SERVE AS THE STOPX (ASYNCHRONOUS INTERRUPTION)
23 C     FUNCTION FOR THE NL2SOL (NONLINEAR LEAST-SQUARES) PACKAGE AT
24 C     THOSE INSTALLATIONS WHICH DO NOT WISH TO IMPLEMENT A
25 C     DYNAMIC STOPX.
26 C
27 C     *****ALGORITHM NOTES...
28 C     AT INSTALLATIONS WHERE THE NL2SOL SYSTEM IS USED
29 C     INTERACTIVELY, THIS DUMMY STOPX SHOULD BE REPLACED BY A
30 C     FUNCTION THAT RETURNS .TRUE. IF AND ONLY IF THE INTERRUPT
31 C     (BREAK) KEY HAS BEEN PRESSED SINCE THE LAST CALL ON STOPX.
32 C
33 C     $$$ MODIFIED FOR USE AS  THE TIMER ROUTINE.
34 C     $$$                              WHEN THE TIME LIMIT HAS BEEN
35 C     $$$ REACHED     STOPX IS SET TO .TRUE  AND INITIATES (IN ITSUM)
36 C     $$$ AND ORDERLY EXIT OUT OF SUMSL.  IF ARRAYS IV AND V ARE
37 C     $$$ SAVED, THE SUMSL ROUTINES CAN BE RESTARTED AT THE SAME
38 C     $$$ POINT AT WHICH THEY WERE INTERRUPTED.
39 C
40 C     ..................................................................
41 C
42       include 'DIMENSIONS'
43       integer nf
44       logical ovrtim
45       include 'COMMON.IOUNITS'
46       include 'COMMON.TIME1'
47 #ifdef MPL
48       include 'COMMON.INFO'
49       integer Kwita
50
51 cd    print *,'Processor',MyID,' NF=',nf
52 #endif
53       if (ovrtim()) then
54 C Finish if time is up.
55          stopx = .true.
56 #ifdef MPL
57       else if (mod(nf,100).eq.0) then
58 C Other processors might have finished. Check this every 100th function 
59 C evaluation.
60 cd       print *,'Processor ',MyID,' is checking STOP: nf=',nf
61          call recv_stop_sig(Kwita)
62          if (Kwita.eq.-1) then
63            write (iout,'(a,i4,a,i5)') 'Processor',
64      &     MyID,' has received STOP signal in STOPX; NF=',nf
65            write (*,'(a,i4,a,i5)') 'Processor',
66      &     MyID,' has received STOP signal in STOPX; NF=',nf
67            stopx=.true.
68          else
69            stopx=.false.
70          endif
71 #endif
72       else
73          stopx = .false.
74       endif
75       return
76       end
77 C=========================================================================
78 C
79       logical function ovrtim()
80       double precision tcpu    ! function 
81       include 'COMMON.TIME1'
82 C Set a 100.0 secs. safety margin, so as to allow for the termination of 
83 C a batch job.
84 c      double  safety /150.0D0/
85       curtim= tcpu()
86 cd    print *,'curtim=',curtim,' timlim=',timlim
87 C  curtim is the current time in seconds.
88       ovrtim=(curtim .ge. timlim - safety )
89       return                                               
90       end
91 C=========================================================================
92 C
93       double precision function tcpu()
94       include 'COMMON.TIME1'
95 #ifdef ES9000 
96 ****************************
97 C Next definition for EAGLE (ibm-es9000)
98       real*8 micseconds
99       integer rcode
100       tcpu=cputime(micseconds,rcode)
101       tcpu=(micseconds/1.0E6) - stime
102 ****************************
103 #endif
104 #ifdef SUN
105 ****************************
106 C Next definitions for sun
107       integer seconds
108       call clock(seconds)
109       tcpu=seconds - stime
110 ****************************
111 #endif
112 #ifdef KSR
113 ****************************
114 C Next definitions for ksr
115 C this function uses the ksr timer ALL_SECONDS from the PMON library to
116 C return the elapsed time in seconds
117       tcpu= all_seconds() - stime
118 ****************************
119 #endif
120 #ifdef SGI
121 ****************************
122 C Next definitions for sgi
123       real timar(2), etime
124       seconds = etime(timar)
125 C      usrsec = timar(1)
126 C      syssec = timar(2)
127       tcpu=seconds - stime
128 ****************************
129 #endif
130 #ifdef CRAY
131 ****************************
132 C Next definitions for Cray
133 C     call date(curdat)
134 C     curdat=curdat(1:9)
135 C     call clock(curtim)
136 C     curtim=curtim(1:8)
137       cpusec = second()
138       tcpu=cpusec - stime
139 ****************************
140 #endif
141 #ifdef AIX
142 ****************************
143 C Next definitions for RS6000
144        integer*4 i1,mclock
145        i1 = mclock()
146        tcpu = (i1+0.0D0)/100.0D0
147 #endif
148 #ifdef LINUX
149 ****************************
150 C Next definitions for sgi
151       real timar(2), etime
152       seconds = etime(timar)
153 Cd    print *,'seconds=',seconds,' stime=',stime
154 C      usrsec = timar(1)
155 C      syssec = timar(2)
156       tcpu=seconds - stime
157 ****************************
158 #endif
159       return     
160       end  
161 *
162       subroutine dajczas(rntime,hrtime,mintime,sectime)
163       include 'COMMON.IOUNITS'
164       real*8 rntime,hrtime,mintime,sectime 
165       hrtime=rntime/3600.0D0 
166       hrtime=aint(hrtime)
167       mintime=aint((rntime-3600.0D0*hrtime)/60.0D0)
168       sectime=aint((rntime-3600.0D0*hrtime-60.0D0*mintime)+0.5D0)
169       if (sectime.eq.60.0D0) then
170         sectime=0.0D0
171         mintime=mintime+1.0D0
172       endif
173       ihr=hrtime
174       imn=mintime
175       isc=sectime
176       write (iout,328) ihr,imn,isc
177   328 FORMAT(//'***** Computation time: ',I4  ,' hours ',I2  ,
178      1         ' minutes ', I2  ,' seconds *****')       
179       return
180       end