hdiff output

r30962/geopt.f 2016-08-04 16:30:08.187584333 +0100 r30961/geopt.f 2016-08-04 16:30:08.439587680 +0100
 24:       USE modcharmm 24:       USE modcharmm
 25:       use porfuncs 25:       use porfuncs
 26:       use binaryio 26:       use binaryio
 27:       USE MODEFOL 27:       USE MODEFOL
 28:       USE MODNEB,ONLY : NEWCONNECTT 28:       USE MODNEB,ONLY : NEWCONNECTT
 29:       USE AMHGLOBALS, ONLY : NMRES 29:       USE AMHGLOBALS, ONLY : NMRES
 30: ! hk286 30: ! hk286
 31:       USE GENRIGID 31:       USE GENRIGID
 32:       use internals_wrapper 32:       use internals_wrapper
 33:       USE MODCUDALBFGS, ONLY : CUDA_LBFGS_WRAPPER 33:       USE MODCUDALBFGS, ONLY : CUDA_LBFGS_WRAPPER
 34:       USE MODCUDABFGSTS, ONLY : CUDA_BFGSTS_WRAPPER 
 35:       USE CHIRALITY, ONLY: CIS_TRANS_CHECK, CHIRALITY_CHECK 
 36:  34: 
 37:       IMPLICIT NONE 35:       IMPLICIT NONE
 38:       INTEGER FCALL, INEG, ECALL, NPCALL, SCALL, ITDONE, INFO, NEXMODES, I1, NUMHB, J1, HORDER, J3, NSTEP, J2, K1, K2, 36:       INTEGER FCALL, INEG, ECALL, NPCALL, SCALL, ITDONE, INFO, NEXMODES, I1, NUMHB, J1, HORDER, J3, NSTEP, J2, K1, K2,
 39:      &        IMAX, IDONE, ITS  37:      &        IMAX, IDONE, ITS 
 40:       DOUBLE PRECISION  DIAG(3*NATOMS), ENERGY, EREAL, EVALMAX, EVALMIN, ETIME, FTIME, STIME, Q(3*NATOMS), OMAX, 38:       DOUBLE PRECISION  DIAG(3*NATOMS), ENERGY, EREAL, EVALMAX, EVALMIN, ETIME, FTIME, STIME, Q(3*NATOMS), OMAX,
 41:      1                  RMS, RMS2, VECS(3*NATOMS), VNEW(3*NATOMS), IT(3,3), IV(3,3), AMASS, ATMASSSAVE(NATOMS), DUMMY, 39:      1                  RMS, RMS2, VECS(3*NATOMS), VNEW(3*NATOMS), IT(3,3), IV(3,3), AMASS, ATMASSSAVE(NATOMS), DUMMY,
 42:      2                  EVALUES(3*NATOMS),TEMPA(9*NATOMS), DUMQ(3*NATOMS), ITX, ITY, ITZ, RMSD, RGYR, PROD, DPRAND, TSDISP, 40:      2                  EVALUES(3*NATOMS),TEMPA(9*NATOMS), DUMQ(3*NATOMS), ITX, ITY, ITZ, RMSD, RGYR, PROD, DPRAND, TSDISP,
 43:      &                  EVSAVE(3*NATOMS), ESAVE, TSDISPSAVE, MINFRQ2, MINCURVE, GRAD(3*NATOMS), RPBN, QFAC, TMPC(3*NATOMS) 41:      &                  EVSAVE(3*NATOMS), ESAVE, TSDISPSAVE, MINFRQ2, MINCURVE, GRAD(3*NATOMS), RPBN, QFAC, TMPC(3*NATOMS)
 44:       DOUBLE PRECISION DIFF, ORDERPLUS, ORDERMINUS, DIST, DIST2, RMAT(3,3), QPATH(3*NATOMS) 42:       DOUBLE PRECISION DIFF, ORDERPLUS, ORDERMINUS, DIST, DIST2, RMAT(3,3), QPATH(3*NATOMS)
 45:       DOUBLE PRECISION, ALLOCATABLE :: ORDERDERIV(:),ORDER(:),ORDERSAVE(:) 43:       DOUBLE PRECISION, ALLOCATABLE :: ORDERDERIV(:),ORDER(:),ORDERSAVE(:)
 46:       DOUBLE PRECISION QSAVE(3*NATOMS) 44:       DOUBLE PRECISION QSAVE(3*NATOMS)
 47:       CHARACTER(LEN=80) FNAMEF 45:       CHARACTER(LEN=80) FNAMEF
 48:       CHARACTER(LEN=20) EFNAMEF 46:       CHARACTER(LEN=20) EFNAMEF
 49:       CHARACTER(LEN=22) ITSTRING 47:       CHARACTER(LEN=22) ITSTRING
 50:       CHARACTER(LEN=5)  ADUMMY 48:       CHARACTER(LEN=5)  ADUMMY
 51:       LOGICAL MFLAG, ZT(3*NATOMS), PERMDISTLOCAL, EWARN 49:       LOGICAL MFLAG, ZT(3*NATOMS), PERMDISTLOCAL, EWARN
 52:       LOGICAL :: BTEST, LSELECT, LNATIVE, INERTIAT=.TRUE., GOODSTRUC1, GOODSTRUC2 50:       LOGICAL :: BTEST, LSELECT, LNATIVE, INERTIAT=.TRUE.
 53: !     AMH LOCAL VARIABLES 51: !     AMH LOCAL VARIABLES
 54:       INTEGER :: NRES,I_RES, GLY_COUNT, IPOT 52:       INTEGER :: NRES,I_RES, GLY_COUNT, IPOT
 55: !      CHARACTER(LEN=5) :: TARFL 53: !      CHARACTER(LEN=5) :: TARFL
 56:       CHARACTER(LEN=2) :: SDUMMY  54:       CHARACTER(LEN=2) :: SDUMMY 
 57:       REAL(8) QTEMP(3,NATOMS), Q4, Q6 55:       REAL(8) QTEMP(3,NATOMS), Q4, Q6
 58:  56: 
 59:       COMMON /PCALL/ NPCALL, ECALL, FCALL, SCALL, ETIME, FTIME, STIME 57:       COMMON /PCALL/ NPCALL, ECALL, FCALL, SCALL, ETIME, FTIME, STIME
 60:       LOGICAL KNOWE, KNOWG, KNOWH 58:       LOGICAL KNOWE, KNOWG, KNOWH
 61:       COMMON /KNOWN/ KNOWE, KNOWG, KNOWH 59:       COMMON /KNOWN/ KNOWE, KNOWG, KNOWH
 62:       LOGICAL CONNECTT, DUMPPATH, READPATH, CALCRATES, STOPFIRST, PROJGRAD 60:       LOGICAL CONNECTT, DUMPPATH, READPATH, CALCRATES, STOPFIRST, PROJGRAD
177:                 VECS(J1)=DPRAND()*2-1.0D0175:                 VECS(J1)=DPRAND()*2-1.0D0
178:              ENDDO176:              ENDDO
179:              CALL VECNORM(VECS,NINTS)177:              CALL VECNORM(VECS,NINTS)
180:              CALL INTBFGSTS(NSTEPS,Q,ENERGY,VNEW,MFLAG,RMS,EVALMIN,EVALMAX,VECS,ITDONE,.TRUE.,.TRUE.) 178:              CALL INTBFGSTS(NSTEPS,Q,ENERGY,VNEW,MFLAG,RMS,EVALMIN,EVALMAX,VECS,ITDONE,.TRUE.,.TRUE.) 
181:           ELSE179:           ELSE
182:              PRINT '(A)', ' geopt> setting random initial vector for eigenvector'180:              PRINT '(A)', ' geopt> setting random initial vector for eigenvector'
183:              DO J1=1,NOPT181:              DO J1=1,NOPT
184:                 VECS(J1)=DPRAND()*2-1.0D0182:                 VECS(J1)=DPRAND()*2-1.0D0
185:              ENDDO 183:              ENDDO 
186:              CALL VECNORM(VECS,NOPT)184:              CALL VECNORM(VECS,NOPT)
187: 185:              CALL BFGSTS(NSTEPS,Q,ENERGY,VNEW,MFLAG,RMS,EVALMIN,EVALMAX,VECS,ITDONE,.TRUE.,.TRUE.)
188:              IF (CUDAT) THEN 
189:                  CALL CUDA_BFGSTS_WRAPPER(NSTEPS,Q,ENERGY,MFLAG,RMS,EVALMIN,VECS,ITDONE) 
190:              ELSE 
191:                  CALL BFGSTS(NSTEPS,Q,ENERGY,VNEW,MFLAG,RMS,EVALMIN,EVALMAX,VECS,ITDONE,.TRUE.,.TRUE.) 
192:              END IF 
193:              IF (MFLAG) THEN 
194:                  WRITE(*,'(1X,A,I6)') 'Converged to TS (number of iterations):     ',ITDONE 
195:              ELSE 
196:                  WRITE(*,'(1X,A,I6)') 'DID NOT converge to TS (number of iterations):     ',ITDONE 
197:              END IF 
198:              WRITE(*,'(A,F20.10)') 'Energy of TS is:                    ',ENERGY 
199:  
200:           ENDIF186:           ENDIF
201: 187: 
202: C188: C
203: C  We may now have the energy and gradient, so the first energy and gradient call in189: C  We may now have the energy and gradient, so the first energy and gradient call in
204: C  mylbfgs for BFGSSTEP may be unnecessary190: C  mylbfgs for BFGSSTEP may be unnecessary
205: C191: C
206:           IF (BFGSSTEP) THEN192:           IF (BFGSSTEP) THEN
207: C193: C
208: C  Switch to appropriate minimization after pushoff.194: C  Switch to appropriate minimization after pushoff.
209: C195: C
1373: ! containing one frame per mode and scaled atomic and residue1359: ! containing one frame per mode and scaled atomic and residue
1374: ! displacements - coords are in Q here!1360: ! displacements - coords are in Q here!
1375:                 IF (CHRMMT.AND.MWVECTORS) CALL CHARMMDUMPMODES(Q,DIAG,ZT,NOPT,3*NATOMS)1361:                 IF (CHRMMT.AND.MWVECTORS) CALL CHARMMDUMPMODES(Q,DIAG,ZT,NOPT,3*NATOMS)
1376: ! AMBER call to routine in amberinterface.f1362: ! AMBER call to routine in amberinterface.f
1377:                 IF ((AMBERT.OR.NABT).AND.MWVECTORS) THEN1363:                 IF ((AMBERT.OR.NABT).AND.MWVECTORS) THEN
1378:                         CALL A9DUMPMODES(DIAG,ZT,NOPT,3*NATOMS)1364:                         CALL A9DUMPMODES(DIAG,ZT,NOPT,3*NATOMS)
1379:                 ENDIF1365:                 ENDIF
1380:               IF(NORMALMODET) CALL VISNORMODES(Q)1366:               IF(NORMALMODET) CALL VISNORMODES(Q)
1381:          ENDIF1367:          ENDIF
1382: 1368: 
1383:          GOODSTRUC1 = .TRUE. 
1384:          GOODSTRUC2 = .TRUE. 
1385:          IF (AMBER12T) THEN 
1386:             CALL CIS_TRANS_CHECK(Q, GOODSTRUC1) 
1387:             CALL CHIRALITY_CHECK(Q, GOODSTRUC2) 
1388:          END IF 
1389:          IF (.NOT. GOODSTRUC1) THEN 
1390:             WRITE(*,'(A)') ' geopt> cistrans problem'  
1391:          END IF 
1392:          IF (.NOT. GOODSTRUC2) THEN 
1393:             WRITE(*,'(A)') ' geopt> chirality problem' 
1394:          END IF 
1395:  
1396:          PRINT*1369:          PRINT*
1397:          WRITE(*,'(A)') ' geopt>                          **** CONVERGED ****'1370:          WRITE(*,'(A)') ' geopt>                          **** CONVERGED ****'
1398:          PRINT*1371:          PRINT*
1399:          CALL FLUSH(6)1372:          CALL FLUSH(6)
1400: 1373: 
1401: !1374: !
1402: ! Print TIP frequencies in cm-1 for CoM/Euler angle representation for debugging angle-axis.1375: ! Print TIP frequencies in cm-1 for CoM/Euler angle representation for debugging angle-axis.
1403: !1376: !
1404:          IF (DEBUG.AND.(ZSYM(1)(1:1).EQ.'W')) THEN1377:          IF (DEBUG.AND.(ZSYM(1)(1:1).EQ.'W')) THEN
1405:             IF (ZSYM(1)(1:2).EQ.'W4') IPOT=41378:             IF (ZSYM(1)(1:2).EQ.'W4') IPOT=4


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0