hdiff output

r32666/lbfgs.f90 2017-06-01 11:30:10.414445811 +0100 r32665/lbfgs.f90 2017-06-01 11:30:11.762463717 +0100
 83:      INTPTEST = .FALSE. 83:      INTPTEST = .FALSE.
 84:      IF (DESMDEBUG) MOREPRINTING = .TRUE. 84:      IF (DESMDEBUG) MOREPRINTING = .TRUE.
 85:      ADDATOM=.FALSE. 85:      ADDATOM=.FALSE.
 86:      NFAIL=0 86:      NFAIL=0
 87:      IF (FREEZENODEST) IMGFREEZE(1:NIMAGE)=.FALSE. 87:      IF (FREEZENODEST) IMGFREEZE(1:NIMAGE)=.FALSE.
 88:      BESTE=HUGE(1.0D0) 88:      BESTE=HUGE(1.0D0)
 89:  89: 
 90:      CALL CHECKINPUT 90:      CALL CHECKINPUT
 91:      IF (DEBUG) CALL DUMPFILES("b") 91:      IF (DEBUG) CALL DUMPFILES("b")
 92:      IF (DUMPNEBXYZ) CALL RWG("w",.False.,0) 92:      IF (DUMPNEBXYZ) CALL RWG("w",.False.,0)
 93: !     IF (DUMPNEBEOS) CALL WRITEPROFILE(-1)  ! sn402: moved this because we 93:      IF (DUMPNEBEOS) DVEC(1:NIMAGE+1)=0.0D0 ! otherwise uninitialised
 94: !     currently don't know the energies of the images 94:      IF (DUMPNEBEOS) CALL WRITEPROFILE(0)
 95:      IF (MOREPRINTING) THEN 95:      IF (MOREPRINTING) THEN
 96:           WRITE(*,'(A)') '                  Iter   Energy per image      RMS Force       Av.Dev    Path     Step per image' 96:           WRITE(*,'(A)') '                  Iter   Energy per image      RMS Force       Av.Dev    Path     Step per image'
 97:      ENDIF 97:      ENDIF
 98:      IF (PREVGRAD.LT.DNEBSWITCH) THEN 98:      IF (PREVGRAD.LT.DNEBSWITCH) THEN
 99:         CALL OLDNEBGRADIENT 99:         CALL OLDNEBGRADIENT
100:      ELSE100:      ELSE
101:         CALL NEBGRADIENT101:         CALL NEBGRADIENT
102:      ENDIF102:      ENDIF
103:      
104:      IF (DUMPNEBEOS) THEN 
105:         CALL DISTANCES  ! sn402: We need to know the distances between the images to write the profile 
106:         CALL WRITEPROFILE(-1) 
107:      ENDIF 
108: 103: 
109:      ENOLD=ETOTAL104:      ENOLD=ETOTAL
110:      IF (ETOTAL/NIMAGE.LT.COLDFUSIONLIMIT) THEN105:      IF (ETOTAL/NIMAGE.LT.COLDFUSIONLIMIT) THEN
111:         WRITE(*,'(A,2G20.10)') ' lbfgs> Cold fusion diagnosed - step discarded, energy, limit=',ETOTAL/NIMAGE,COLDFUSIONLIMIT106:         WRITE(*,'(A,2G20.10)') ' lbfgs> Cold fusion diagnosed - step discarded, energy, limit=',ETOTAL/NIMAGE,COLDFUSIONLIMIT
112:         IF (DEBUG) CALL DUMPFILES("e")107:         IF (DEBUG) CALL DUMPFILES("e")
113:         IF (KADJUSTFRQ.GT.0) PRINT '(A,G20.10)',' lbfgs> Final DNEB force constant ',NEWNEBK(1)108:         IF (KADJUSTFRQ.GT.0) PRINT '(A,G20.10)',' lbfgs> Final DNEB force constant ',NEWNEBK(1)
114:         RETURN109:         RETURN
115:      ENDIF110:      ENDIF
116: 111: 
117:      NITERDONE=1112:      NITERDONE=1


r32666/minpermdist.f90 2017-06-01 11:30:10.866451815 +0100 r32665/minpermdist.f90 2017-06-01 11:30:12.202469562 +0100
472: IF (NFREEZE.LE.0 .AND. .NOT.MIEFT .AND. .NOT. NOTRANSROTT) THEN472: IF (NFREEZE.LE.0 .AND. .NOT.MIEFT .AND. .NOT. NOTRANSROTT) THEN
473:    IF (BULKT) THEN  ! we will not be doing any rotations.473:    IF (BULKT) THEN  ! we will not be doing any rotations.
474:       IF (BMTEST) BMCOORDSSV(1:3*NATOMS)=BMCOORDS(1:3*NATOMS)474:       IF (BMTEST) BMCOORDSSV(1:3*NATOMS)=BMCOORDS(1:3*NATOMS)
475:       NORBIT1=1; NORBIT2=1; NORBITB1=1; NORBITB2=1;475:       NORBIT1=1; NORBIT2=1; NORBITB1=1; NORBITB2=1;
476:       ! BULKMINDIST is used for a quick deterministic check for permutational isomers for periodic systems. 476:       ! BULKMINDIST is used for a quick deterministic check for permutational isomers for periodic systems. 
477:       ! If BMTEST is TRUE (i.e. ATOMMATCHDIST or ATOMMATCHFULL) then this routine also attempts to find the translational-permutational 477:       ! If BMTEST is TRUE (i.e. ATOMMATCHDIST or ATOMMATCHFULL) then this routine also attempts to find the translational-permutational 
478:       ! alignment that maximises the number of exact atom matches between the two structures. This is not necessarily the alignment with 478:       ! alignment that maximises the number of exact atom matches between the two structures. This is not necessarily the alignment with 
479:       ! the minimum cartesian distance.479:       ! the minimum cartesian distance.
480:       ! BULKMINDIST returns DUMMYB unchanged and BMCOORDS as the best trans+perm alignment of structure A. If a permutational isomer is480:       ! BULKMINDIST returns DUMMYB unchanged and BMCOORDS as the best trans+perm alignment of structure A. If a permutational isomer is
481:       ! detected, DUMMYA contains its coordinates and PITEST is TRUE. Otherwise, PITEST is FALSE and DUMMYA should be unchanged.481:       ! detected, DUMMYA contains its coordinates and PITEST is TRUE. Otherwise, PITEST is FALSE and DUMMYA should be unchanged.
 482:       write(*,*) "Before call to BULKMINDIST"
482:       CALL BULKMINDIST(DUMMYB,DUMMYA,BMCOORDS,NATOMS,DISTANCE,TWOD,DEBUG,BOXLX,BOXLY,BOXLZ,PITEST,.TRUE., TNMATCH, BMTEST)483:       CALL BULKMINDIST(DUMMYB,DUMMYA,BMCOORDS,NATOMS,DISTANCE,TWOD,DEBUG,BOXLX,BOXLY,BOXLZ,PITEST,.TRUE., TNMATCH, BMTEST)
483:       IF (PITEST) THEN484:       IF (PITEST) THEN
484:          ! The two structures are permutational isomers. The aligned structure A is contained in DUMMYA485:          ! The two structures are permutational isomers. The aligned structure A is contained in DUMMYA
485:          COORDSA(1:3*NATOMS)=DUMMYA(1:3*NATOMS)486:          COORDSA(1:3*NATOMS)=DUMMYA(1:3*NATOMS)
486:          ! No rotations (periodic system)487:          ! No rotations (periodic system)
487:          RMATBEST(1:3,1:3)=0.0D0488:          RMATBEST(1:3,1:3)=0.0D0
488:          RMATBEST(1,1)=1.0D0; RMATBEST(2,2)=1.0D0; RMATBEST(3,3)=1.0D0489:          RMATBEST(1,1)=1.0D0; RMATBEST(2,2)=1.0D0; RMATBEST(3,3)=1.0D0
489:          DISTANCE=SQRT(DISTANCE)490:          DISTANCE=SQRT(DISTANCE)
490:          IF (DEBUG) PRINT '(A,G20.10)',' minpermdist> after initial call to BULKMINDIST, distance=',DISTANCE491:          IF (DEBUG) PRINT '(A,G20.10)',' minpermdist> after initial call to BULKMINDIST, distance=',DISTANCE
491:          IF (DEBUG) PRINT '(A)',' minpermdist> permutation-inversion isomers identified by BULKMINDIST'492:          IF (DEBUG) PRINT '(A)',' minpermdist> permutation-inversion isomers identified by BULKMINDIST'
573:                 PRINT '(2(A,F25.15))',' final   energy for structure B=             ',ENERGY,' RMS=',RMS574:                 PRINT '(2(A,F25.15))',' final   energy for structure B=             ',ENERGY,' RMS=',RMS
574:                 IF (ABS(ENERGY-BINIT).GT.2*EDIFFTOL) THEN575:                 IF (ABS(ENERGY-BINIT).GT.2*EDIFFTOL) THEN
575:                    PRINT '(A)',' minpermdist> ERROR *** energy change for structure B is outside tolerance'576:                    PRINT '(A)',' minpermdist> ERROR *** energy change for structure B is outside tolerance'
576:                    STOP577:                    STOP
577:                 ENDIF578:                 ENDIF
578:              ENDIF579:              ENDIF
579: 580: 
580:              RETURN581:              RETURN
581:             ENDIF582:             ENDIF
582:          ELSE583:          ELSE
583:             WRITE(*,*) "minpermdist> Error in BMTEST (atom-matching) block. Stopping now"584:             WRITE(*,*) "minpermdist> Error in BMTEST (atom-matching) block."! Stopping now"
584:             STOP585:             !STOP
585:          ENDIF586:          ENDIF
586:       ELSE587:       ELSE
587:          CALL NEWMINDIST(DUMMYB,DUMMYA,NATOMS,DISTANCE,BULKT,TWOD,'AX    ',.FALSE.,RIGID,DEBUG,RMAT)588:          CALL NEWMINDIST(DUMMYB,DUMMYA,NATOMS,DISTANCE,BULKT,TWOD,'AX    ',.FALSE.,RIGID,DEBUG,RMAT)
588:          IF (DEBUG) PRINT '(A,G20.10)',' minpermdist> after initial call to BULK/NEWMINDIST distance=',DISTANCE589:          IF (DEBUG) PRINT '(A,G20.10)',' minpermdist> after initial call to BULK/NEWMINDIST distance=',DISTANCE
589:          DISTANCE=DISTANCE**2 ! minperdist returns the distance squared for historical reasons590:          DISTANCE=DISTANCE**2 ! minperdist returns the distance squared for historical reasons
590:       ENDIF591:       ENDIF
591:    ELSEIF (MKTRAPT) THEN592:    ELSEIF (MKTRAPT) THEN
592:       TMAT(1:3,1:3)=0.0D0593:       TMAT(1:3,1:3)=0.0D0
593:       TMAT(1,1)=INVERT*1.0D0; TMAT(2,2)=INVERT*1.0D0; TMAT(3,3)=INVERT*1.0D0594:       TMAT(1,1)=INVERT*1.0D0; TMAT(2,2)=INVERT*1.0D0; TMAT(3,3)=INVERT*1.0D0
594:       NORBIT1=1; NORBIT2=1; NORBITB1=1; NORBITB2=1;595:       NORBIT1=1; NORBIT2=1; NORBITB1=1; NORBITB2=1;


r32666/mylbfgs.f 2017-06-01 11:30:11.090454791 +0100 r32665/mylbfgs.f 2017-06-01 11:30:12.426472537 +0100
471:             G(J1)=VEC2(J1)471:             G(J1)=VEC2(J1)
472: C           G(J1)=VEC2(J1)/ENERGY472: C           G(J1)=VEC2(J1)/ENERGY
473:          ENDDO473:          ENDDO
474:       ELSE474:       ELSE
475:          EREAL=ENERGY475:          EREAL=ENERGY
476:          REALRMS=RMS476:          REALRMS=RMS
477:       ENDIF477:       ENDIF
478:       IF (PTEST) WRITE(*,'(A,2G20.10,A,I6,A)') 478:       IF (PTEST) WRITE(*,'(A,2G20.10,A,I6,A)') 
479:      1             ' mylbfgs> Energy and RMS force=',ENERGY,RMS,' after ',ITDONE,' steps'479:      1             ' mylbfgs> Energy and RMS force=',ENERGY,RMS,' after ',ITDONE,' steps'
480:       WRITE(ESTRING,16) ' mylbfgs> Energy for last cycle=',ENERGY,' '480:       WRITE(ESTRING,16) ' mylbfgs> Energy for last cycle=',ENERGY,' '
481: 16       FORMAT(A,26X,F20.10,A)481: 16       FORMAT(A,27X,F20.10,A)
482: 482: 
483: 10    CALL FLUSH(6)483: 10    CALL FLUSH(6)
484: 484: 
485: C     add additional check for gradient in internals. Otherwise rigid bonds might never converge.485: C     add additional check for gradient in internals. Otherwise rigid bonds might never converge.
486: C     GINT with bonds projected out is never transferred back to carthesians -> check in internals486: C     GINT with bonds projected out is never transferred back to carthesians -> check in internals
487:       MFLAG=.FALSE.487:       MFLAG=.FALSE.
488:       IF (RMS.LE.GMAX.OR.(INTMINT.AND.RMSINT.LE.GMAXINT)) THEN ! GMAX is in key module, so can be changed by changep call488:       IF (RMS.LE.GMAX.OR.(INTMINT.AND.RMSINT.LE.GMAXINT)) THEN ! GMAX is in key module, so can be changed by changep call
489:          IF (CHRMMT.AND.ACESOLV) THEN489:          IF (CHRMMT.AND.ACESOLV) THEN
490:             NCHENCALLS=ACEUPSTEP-1490:             NCHENCALLS=ACEUPSTEP-1
491:             CALL POTENTIAL(X,ENERGY,GSAVE,.TRUE.,.FALSE.,RMS,.FALSE.,.FALSE.)491:             CALL POTENTIAL(X,ENERGY,GSAVE,.TRUE.,.FALSE.,RMS,.FALSE.,.FALSE.)


r32666/nnutils.f90 2017-06-01 11:30:10.638448787 +0100 r32665/nnutils.f90 2017-06-01 11:30:11.982466639 +0100
1173:           DOUBLE PRECISION :: DUMMY1173:           DOUBLE PRECISION :: DUMMY
1174:           CHARACTER(LEN=20) :: FILENAME1174:           CHARACTER(LEN=20) :: FILENAME
1175:           1175:           
1176: !         IF (PRESENT(UNITIN)) THEN1176: !         IF (PRESENT(UNITIN)) THEN
1177: !              UNIT=UNITIN1177: !              UNIT=UNITIN
1178: !         ELSE1178: !         ELSE
1179:                UNIT=9921179:                UNIT=992
1180:                IF (NITER.GT.0) THEN1180:                IF (NITER.GT.0) THEN
1181:                   WRITE(FILENAME,'(I8)') NITER1181:                   WRITE(FILENAME,'(I8)') NITER
1182:                   FILENAME='neb.EofS.' // TRIM(ADJUSTL(FILENAME))1182:                   FILENAME='neb.EofS.' // TRIM(ADJUSTL(FILENAME))
1183:                ELSEIF (NITER.EQ.-1) THEN 
1184:                   FILENAME='neb.EofS.0' 
1185:                ELSE   1183:                ELSE   
1186:                   FILENAME='neb.EofS'1184:                   FILENAME='neb.EofS'
1187:                ENDIF1185:                ENDIF
1188:                IF (.NOT.FILTH==0) THEN1186:                IF (.NOT.FILTH==0) THEN
1189:                     FILENAME=TRIM(FILENAME)//'.'//TRIM(ADJUSTL(FILTHSTR))1187:                     FILENAME=TRIM(FILENAME)//'.'//TRIM(ADJUSTL(FILTHSTR))
1190:                ENDIF1188:                ENDIF
1191:                OPEN(UNIT=UNIT,FILE=FILENAME,STATUS='replace')1189:                OPEN(UNIT=UNIT,FILE=FILENAME,STATUS='replace')
1192: !         ENDIF1190: !         ENDIF
1193: 1191: 
1194:           DUMMY=0.0D01192:           DUMMY=0.0D0


r32666/secdiag.f90 2017-06-01 11:30:11.310457712 +0100 r32665/secdiag.f90 2017-06-01 11:30:12.642475406 +0100
 48:       DOUBLE PRECISION, INTENT(IN) :: ENERGY ! the energy at point coords 48:       DOUBLE PRECISION, INTENT(IN) :: ENERGY ! the energy at point coords
 49:       DOUBLE PRECISION, INTENT(IN) :: GRAD(3*NATOMS) ! the gradient at coords (will be nonsense unless nsecdiag > 2) 49:       DOUBLE PRECISION, INTENT(IN) :: GRAD(3*NATOMS) ! the gradient at coords (will be nonsense unless nsecdiag > 2)
 50:       DOUBLE PRECISION, INTENT(OUT) :: DIAG ! the curvature 50:       DOUBLE PRECISION, INTENT(OUT) :: DIAG ! the curvature
 51:       DOUBLE PRECISION, INTENT(OUT) :: GL(3*NATOMS) ! the gradient of the curvature 51:       DOUBLE PRECISION, INTENT(OUT) :: GL(3*NATOMS) ! the gradient of the curvature
 52:       DOUBLE PRECISION, INTENT(OUT) :: XRMS ! the rms of GL 52:       DOUBLE PRECISION, INTENT(OUT) :: XRMS ! the rms of GL
 53:       DOUBLE PRECISION DUMMY3(3*NATOMS), DIFF, DIAG2, DIAG3, & 53:       DOUBLE PRECISION DUMMY3(3*NATOMS), DIFF, DIAG2, DIAG3, &
 54:                        EPLUS,EMINUS,GRAD1(3*NATOMS),LOCALV(3*NATOMS), & 54:                        EPLUS,EMINUS,GRAD1(3*NATOMS),LOCALV(3*NATOMS), &
 55:                        RMS,GRAD2(3*NATOMS), VECL, ZETA, PROJ 55:                        RMS,GRAD2(3*NATOMS), VECL, ZETA, PROJ
 56:       double precision diag4 ! the curvature computed with the first order forward finite differences method 56:       double precision diag4 ! the curvature computed with the first order forward finite differences method
 57:  57: 
 58:       DIAG = 0.0D0 
 59:       DIAG2 = 0.0D0 58:       DIAG2 = 0.0D0
 60:       DIAG3 = 0.0D0 59:       DIAG3 = 0.0D0
 61: !      EPLUS = 0.0D0 60: !      EPLUS = 0.0D0
 62:       EMINUS = 0.0D0 61:       EMINUS = 0.0D0
 63: !      XRMS = 0.0D0 62: !      XRMS = 0.0D0
 64:       DIFF=1.0D-3 63:       DIFF=1.0D-3
 65:       IF (NIH2LEPST .or. NIMET .or. NIHEAM7T .or. NIHLEPST .or. NIHPAIRONLYT) DIFF=1.0D-5 64:       IF (NIH2LEPST .or. NIMET .or. NIHEAM7T .or. NIHLEPST .or. NIHPAIRONLYT) DIFF=1.0D-5
 66:       IF (PYADD2T) DIFF=1.0D-5 65:       IF (PYADD2T) DIFF=1.0D-5
 67:       IF (CHARMMDFTBT) DIFF=1.0D-2 66:       IF (CHARMMDFTBT) DIFF=1.0D-2
 68:       IF (AMHT) DIFF=1.0D-2 67:       IF (AMHT) DIFF=1.0D-2
217:          DO J1=1,NOPT216:          DO J1=1,NOPT
218:             GL(J1)=GL(J1)-PROJ*LOCALV(J1)217:             GL(J1)=GL(J1)-PROJ*LOCALV(J1)
219:          ENDDO218:          ENDDO
220:          XRMS=0.0D0219:          XRMS=0.0D0
221:          DO J1=1,NOPT220:          DO J1=1,NOPT
222:             XRMS=XRMS+GL(J1)**2221:             XRMS=XRMS+GL(J1)**2
223:          ENDDO222:          ENDDO
224:          XRMS=DSQRT(XRMS/NOPT)223:          XRMS=DSQRT(XRMS/NOPT)
225: !        PRINT *,'secdiag> after proj stuff GL:'224: !        PRINT *,'secdiag> after proj stuff GL:'
226: !        PRINT '(3F20.10)',GL(1:NOPT)225: !        PRINT '(3F20.10)',GL(1:NOPT)
227:          IF (DEBUG) THEN226:          IF (DEBUG) WRITE(*,'(A,3G15.5,3G20.12,G10.3)') 'D,D2,D3,E+,E-,E,RMS=',DIAG,DIAG2,DIAG3,EPLUS,EMINUS,ENERGY,XRMS
228:             IF(NSECDIAG.NE.3) THEN227:          IF (DEBUG.AND.(NSECDIAG.NE.3)) WRITE(*,'(A,G20.10)') 'predicted gradient component=',(EPLUS-EMINUS)/(2*ZETA)
229:                WRITE(*,'(A,3G15.5,3G20.12,G10.3)') 'D,D2,D3,E+,E-,E,RMS=',DIAG,DIAG2,DIAG3,EPLUS,EMINUS,ENERGY,XRMS 
230:                WRITE(*,'(A,G20.10)') 'predicted gradient component=',(EPLUS-EMINUS)/(2*ZETA) 
231:             ELSE  ! forward difference method: don't need to print so much 
232:                WRITE(*,'(A,3G15.5,3G20.12,G10.3)') 'D4,E+,E,RMS=',DIAG4,EPLUS,ENERGY,XRMS 
233:             ENDIF 
234:          ENDIF 
235:       ENDIF228:       ENDIF
236: !     PRINT '(A)','LOCALV:'229: !     PRINT '(A)','LOCALV:'
237: !     PRINT '(3G20.10)',LOCALV(1:NOPT)230: !     PRINT '(3G20.10)',LOCALV(1:NOPT)
238:       IF (NSECDIAG.EQ.2) DIAG=DIAG2231:       IF (NSECDIAG.EQ.2) DIAG=DIAG2
239:       IF (NSECDIAG.GT.2) DIAG=DIAG4232:       IF (NSECDIAG.GT.2) DIAG=DIAG4
240: 233: 
241:       RETURN234:       RETURN
242:       END235:       END


r32666/xmylbfgs.f 2017-06-01 11:30:11.530460635 +0100 r32665/xmylbfgs.f 2017-06-01 11:30:12.862478328 +0100
350:       EVPERCENT=100.0D0*ABS(ENEW-ENERGY)/MAX(ABS(ENEW),1.0D-100)350:       EVPERCENT=100.0D0*ABS(ENEW-ENERGY)/MAX(ABS(ENEW),1.0D-100)
351: C     IF (((ENEW-ENERGY)/ABS(ENEW).LE.XMAXERISE).OR.(CASTEP.AND.((ENEW-ENERGY)/ABS(ENEW).LE.1.0D-1))) THEN351: C     IF (((ENEW-ENERGY)/ABS(ENEW).LE.XMAXERISE).OR.(CASTEP.AND.((ENEW-ENERGY)/ABS(ENEW).LE.1.0D-1))) THEN
352:       IF (((ENEW-ENERGY)/ABS(ENEW).LE.XMAXERISE)) THEN352:       IF (((ENEW-ENERGY)/ABS(ENEW).LE.XMAXERISE)) THEN
353: C     IF (.TRUE.) THEN353: C     IF (.TRUE.) THEN
354:          ITER=ITER+1354:          ITER=ITER+1
355:          ITDONE=ITDONE+1355:          ITDONE=ITDONE+1
356:          ENERGY=ENEW356:          ENERGY=ENEW
357:          DO J1=1,N357:          DO J1=1,N
358:             G(J1)=GNEW(J1)358:             G(J1)=GNEW(J1)
359:          ENDDO359:          ENDDO
360:          IF (PTEST) WRITE(*,'(A,3G19.10,A,I6,A,F13.10)') 360:          IF (PTEST) WRITE(*,'(A,3G20.10,A,I6,A,F13.10)') 
361:      1          'xmylbfgs> Eigenvalue, RMS, % change=',ENERGY,RMS,EVPERCENT,' after ',ITDONE,' steps, step:',STP*SLENGTH361:      1          'xmylbfgs> Eigenvalue, RMS, % change=',ENERGY,RMS,EVPERCENT,' after ',ITDONE,' steps, step:',STP*SLENGTH
362: !        PRINT '(A)','X for this eigenvalue:'362: !        PRINT '(A)','X for this eigenvalue:'
363: !        PRINT '(3G20.10)',X(1:15)363: !        PRINT '(3G20.10)',X(1:15)
364:       ELSE364:       ELSE
365: C365: C
366: C  Eigenvalue increased - try again with a smaller step size?366: C  Eigenvalue increased - try again with a smaller step size?
367: C367: C
368:          IF (NDECREASE.GT.3) THEN   ! decreasing the step size for castep doesn;t seem to help.368:          IF (NDECREASE.GT.3) THEN   ! decreasing the step size for castep doesn;t seem to help.
369: !        IF ((NDECREASE.GT.5).OR.(CASTEP.OR.ONETEP.OR.CP2K)) THEN   ! decreasing the step size for castep doesn;t seem to help.369: !        IF ((NDECREASE.GT.5).OR.(CASTEP.OR.ONETEP.OR.CP2K)) THEN   ! decreasing the step size for castep doesn;t seem to help.
370:                                                                                  ! could try this again one day!370:                                                                                  ! could try this again one day!


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0