hdiff output

r33343/intlbfgs.f90 2017-09-24 13:30:11.098869759 +0100 r33342/intlbfgs.f90 2017-09-24 13:30:11.354873122 +0100
307:    WRITE(*,'(12I8)') NREPI(1:NNREPULSIVE)307:    WRITE(*,'(12I8)') NREPI(1:NNREPULSIVE)
308:    READ(LUNIT,*) NREPJ(1:NNREPULSIVE)308:    READ(LUNIT,*) NREPJ(1:NNREPULSIVE)
309:    WRITE(*,'(A)') ' intlbfgs> read NREPJ:'309:    WRITE(*,'(A)') ' intlbfgs> read NREPJ:'
310:    WRITE(*,'(12I8)') NREPJ(1:NNREPULSIVE)310:    WRITE(*,'(12I8)') NREPJ(1:NNREPULSIVE)
311:    READ(LUNIT,*) REPCUT(1:NREPULSIVE)311:    READ(LUNIT,*) REPCUT(1:NREPULSIVE)
312:    WRITE(*,'(A)') ' intlbfgs> read REPCUT:'312:    WRITE(*,'(A)') ' intlbfgs> read REPCUT:'
313:    WRITE(*,'(6G20.10)') REPCUT(1:NREPULSIVE)313:    WRITE(*,'(6G20.10)') REPCUT(1:NREPULSIVE)
314:    READ(LUNIT,*) NREPCUT(1:NNREPULSIVE)314:    READ(LUNIT,*) NREPCUT(1:NNREPULSIVE)
315:    WRITE(*,'(A)') ' intlbfgs> read NREPCUT:'315:    WRITE(*,'(A)') ' intlbfgs> read NREPCUT:'
316:    WRITE(*,'(6G20.10)') NREPCUT(1:NNREPULSIVE)316:    WRITE(*,'(6G20.10)') NREPCUT(1:NNREPULSIVE)
317:    READ(LUNIT,*) INTFROZEN(1:NATOMS) 
318:    WRITE(*,'(A)') ' intlbfgs> read INTFROZEN' 
319:    WRITE(*,'(12L5)') INTFROZEN(1:NATOMS) 
320:    CLOSE(LUNIT)317:    CLOSE(LUNIT)
321: 318: 
322:    GLAST(1:D)=G(1:D)319:    GLAST(1:D)=G(1:D)
323:    XSAVE(1:D)=X(1:D)320:    XSAVE(1:D)=X(1:D)
324:    GOTO 986321:    GOTO 986
325: ENDIF322: ENDIF
326: IF (INTFREEZET) THEN323: IF (INTFREEZET) THEN
327:    DO J1=1,NATOMS324:    DO J1=1,NATOMS
328:       IF (INTFROZEN(J1)) THEN325:       IF (INTFROZEN(J1)) THEN
329: ! 326: ! 
1166:          JMAXRMS=J11163:          JMAXRMS=J1
1167:       ENDIF1164:       ENDIF
1168:    ENDDO1165:    ENDDO
1169:    MAXRMS=SQRT(MAXRMS/(3*NACTIVE))1166:    MAXRMS=SQRT(MAXRMS/(3*NACTIVE))
1170:    SUMEEE=SUMEEE/INTIMAGE1167:    SUMEEE=SUMEEE/INTIMAGE
1171:    SUMEEE2=SUMEEE2+EEE(J1)**21168:    SUMEEE2=SUMEEE2+EEE(J1)**2
1172:    SUMEEE2=SUMEEE2/INTIMAGE1169:    SUMEEE2=SUMEEE2/INTIMAGE
1173:    SIGMAEEE=SQRT(SUMEEE2-SUMEEE**2)1170:    SIGMAEEE=SQRT(SUMEEE2-SUMEEE**2)
1174:    REMOVEIMAGE=.FALSE.1171:    REMOVEIMAGE=.FALSE.
1175: !  IF (ABS(MAXEEE-SUMEEE).GT.3.0D0*SIGMAEEE) THEN1172: !  IF (ABS(MAXEEE-SUMEEE).GT.3.0D0*SIGMAEEE) THEN
1176: !  IF (MAXEEE.GT.1.0D2*MINEEE) THEN1173:    IF (MAXEEE.GT.1.0D2*MINEEE) THEN
1177:       WRITE(*,'(A,I8,A,G20.10,A,G20.10,A)') ' intlbfgs> Highest image ',JMAXEEE,' energy ',MAXEEE,' is ',ABS(MAXEEE-SUMEEE)/SIGMAEEE, &1174:       WRITE(*,'(A,I8,A,G20.10,A)') ' intlbfgs> warning *** highest image ',JMAXEEE,' energy ',MAXEEE,' is > 100 * higher than minimum '  
1178:   &                        ' sigma from the mean' 
1179: !     REMOVEIMAGE=.TRUE.1175: !     REMOVEIMAGE=.TRUE.
1180: !     IF (NITERDONE-NLASTGOODE.LT.40) REMOVEIMAGE=.FALSE.1176: !     IF (NITERDONE-NLASTGOODE.LT.40) REMOVEIMAGE=.FALSE.
1181: !  ENDIF1177:    ENDIF
1182: 1178: 
1183:    IF (DEBUG) THEN1179:    IF (DEBUG) THEN
1184:       WRITE(*,'(A,2G20.10)') ' intlbfgs> mean and sigma of image energies=',SUMEEE,SIGMAEEE1180:       WRITE(*,'(A,2G20.10)') ' intlbfgs> mean and sigma of image energies=',SUMEEE,SIGMAEEE
1185:       WRITE(*,'(A,I6,2G20.10,3(G20.10,I8))') ' intlbfgs> steps: ',NITERDONE,ETOTAL/INTIMAGE,RMS,STEPTOT,NACTIVE, &1181:       WRITE(*,'(A,I6,2G20.10,3(G20.10,I8))') ' intlbfgs> steps: ',NITERDONE,ETOTAL/INTIMAGE,RMS,STEPTOT,NACTIVE, &
1186:   &                                                        MAXEEE,JMAXEEE,MAXRMS,JMAXRMS1182:   &                                                        MAXEEE,JMAXEEE,MAXRMS,JMAXRMS
1187:       CALL FLUSH(6)1183:       CALL FLUSH(6)
1188:    ENDIF1184:    ENDIF
1189: 1185: 
1190:    IF (.NOT.SWITCHED) THEN1186:    IF (.NOT.SWITCHED) THEN
1191: !     IF ((NITERDONE-NLASTGOODE.GT.INTRELSTEPS).AND.((ETOTAL.GT.LASTGOODE).OR.(ETOTAL/INTIMAGE.GT.MAXCONE*1.0D8))) THEN1187: !     IF ((NITERDONE-NLASTGOODE.GT.INTRELSTEPS).AND.((ETOTAL.GT.LASTGOODE).OR.(ETOTAL/INTIMAGE.GT.MAXCONE*1.0D8))) THEN
1475:    ENDDO 1471:    ENDDO 
1476: 246 CONTINUE1472: 246 CONTINUE
1477: ENDDO1473: ENDDO
1478: IF (DEBUG) WRITE(*,'(A,2I8)') ' checkrep> number of active repulsions and total=',NNREPULSIVE,NREPULSIVE1474: IF (DEBUG) WRITE(*,'(A,2I8)') ' checkrep> number of active repulsions and total=',NNREPULSIVE,NREPULSIVE
1479: 1475: 
1480: END SUBROUTINE CHECKREP1476: END SUBROUTINE CHECKREP
1481: 1477: 
1482: SUBROUTINE INTRWG(NACTIVE,NITER,INTIMAGE,XYZ,TURNONORDER)1478: SUBROUTINE INTRWG(NACTIVE,NITER,INTIMAGE,XYZ,TURNONORDER)
1483: USE PORFUNCS1479: USE PORFUNCS
1484: USE KEY,ONLY: STOCKT,STOCKAAT, RBAAT, ATOMACTIVE, NCONSTRAINT, CONACTIVE, NREPULSIVE, NNREPULSIVE, REPI, REPJ, REPCUT, NREPCUT, &1480: USE KEY,ONLY: STOCKT,STOCKAAT, RBAAT, ATOMACTIVE, NCONSTRAINT, CONACTIVE, NREPULSIVE, NNREPULSIVE, REPI, REPJ, REPCUT, NREPCUT, &
1485:   &           NREPMAX, NREPI, NREPJ, INTFROZEN1481:   &           NREPMAX, NREPI, NREPJ
1486: USE COMMONS, ONLY: NATOMS1482: USE COMMONS, ONLY: NATOMS
1487: IMPLICIT NONE1483: IMPLICIT NONE
1488: CHARACTER(LEN=10) :: XYZFILE   = 'int.xyz   '1484: CHARACTER(LEN=10) :: XYZFILE   = 'int.xyz   '
1489: CHARACTER(LEN=10) :: QCIFILE   = 'QCIdump   '1485: CHARACTER(LEN=10) :: QCIFILE   = 'QCIdump   '
1490: INTEGER,INTENT(IN) :: NITER, TURNONORDER(NATOMS)1486: INTEGER,INTENT(IN) :: NITER, TURNONORDER(NATOMS)
1491: INTEGER :: J1,J2,INTIMAGE,J3,NACTIVE,LUNIT,GETUNIT1487: INTEGER :: J1,J2,INTIMAGE,J3,NACTIVE,LUNIT,GETUNIT
1492: CHARACTER(LEN=80) :: FILENAME,DUMMYS1488: CHARACTER(LEN=80) :: FILENAME,DUMMYS
1493: DOUBLE PRECISION XYZ((3*NATOMS)*(INTIMAGE+2))1489: DOUBLE PRECISION XYZ((3*NATOMS)*(INTIMAGE+2))
1494: 1490: 
1495: FILENAME=XYZFILE1491: FILENAME=XYZFILE
1542:    WRITE(LUNIT,'(12I8)') NREPI(1:NNREPULSIVE)1538:    WRITE(LUNIT,'(12I8)') NREPI(1:NNREPULSIVE)
1543:    WRITE(*,'(A)') ' intlbfgs> dumped NREPI:'1539:    WRITE(*,'(A)') ' intlbfgs> dumped NREPI:'
1544:    WRITE(LUNIT,'(12I8)') NREPJ(1:NNREPULSIVE)1540:    WRITE(LUNIT,'(12I8)') NREPJ(1:NNREPULSIVE)
1545:    WRITE(*,'(A)') ' intlbfgs> dumped NREPJ:'1541:    WRITE(*,'(A)') ' intlbfgs> dumped NREPJ:'
1546: 1542: 
1547:    WRITE(LUNIT,'(6G20.10)') REPCUT(1:NREPULSIVE)1543:    WRITE(LUNIT,'(6G20.10)') REPCUT(1:NREPULSIVE)
1548:    WRITE(*,'(A)') ' intlbfgs> dumped REPCUT:'1544:    WRITE(*,'(A)') ' intlbfgs> dumped REPCUT:'
1549:    WRITE(LUNIT,'(6G20.10)') NREPCUT(1:NNREPULSIVE)1545:    WRITE(LUNIT,'(6G20.10)') NREPCUT(1:NNREPULSIVE)
1550:    WRITE(*,'(A)') ' intlbfgs> dumped NREPCUT:'1546:    WRITE(*,'(A)') ' intlbfgs> dumped NREPCUT:'
1551: 1547: 
1552:    WRITE(LUNIT,'(12L5)') INTFROZEN(1:NATOMS) 
1553:    WRITE(*,'(A)') ' intlbfgs> dumped INTFROZEN' 
1554:  
1555: CLOSE(LUNIT)1548: CLOSE(LUNIT)
1556: 1549: 
1557: END SUBROUTINE INTRWG1550: END SUBROUTINE INTRWG
1558: 1551: 
1559: SUBROUTINE WRITEPROFILE(NITER,EEE,INTIMAGE)1552: SUBROUTINE WRITEPROFILE(NITER,EEE,INTIMAGE)
1560: IMPLICIT NONE 1553: IMPLICIT NONE 
1561: INTEGER,INTENT(IN) :: NITER, INTIMAGE1554: INTEGER,INTENT(IN) :: NITER, INTIMAGE
1562: INTEGER :: I,LUNIT,GETUNIT1555: INTEGER :: I,LUNIT,GETUNIT
1563: DOUBLE PRECISION :: EEE(INTIMAGE+2)1556: DOUBLE PRECISION :: EEE(INTIMAGE+2)
1564: CHARACTER(LEN=20) :: FILENAME1557: CHARACTER(LEN=20) :: FILENAME


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0