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