hdiff output

r28822/geopt.f 2015-11-17 23:34:38.581826802 +0000 r28821/geopt.f 2015-11-17 23:34:38.997832382 +0000
360:                NEXMODES=6360:                NEXMODES=6
361:             ELSEIF (BFGSTST) THEN361:             ELSEIF (BFGSTST) THEN
362:                NEXMODES=7362:                NEXMODES=7
363:             ELSEIF (INR.EQ.2) THEN ! this will detect higher order saddles363:             ELSEIF (INR.EQ.2) THEN ! this will detect higher order saddles
364:                NEXMODES=7364:                NEXMODES=7
365:             ENDIF365:             ENDIF
366:             IF (BULKT) NEXMODES=NEXMODES-3366:             IF (BULKT) NEXMODES=NEXMODES-3
367:             IF (BULKT.AND.TWOD) NEXMODES=NATOMS+2367:             IF (BULKT.AND.TWOD) NEXMODES=NATOMS+2
368:             IF (PULLT.OR.EFIELDT) NEXMODES=4368:             IF (PULLT.OR.EFIELDT) NEXMODES=4
369: ! hk286369: ! hk286
370:             IF (GTHOMSONT .AND. (GTHOMMET .EQ. 5) ) NEXMODES = 3370:             IF (GTHOMSONT .AND. (GTHOMMET .EQ. 5) ) NEXMODES = 3 + NATOMS
371:             IF (GTHOMSONT .AND. (GTHOMMET < 5) ) NEXMODES = 1371:             IF (GTHOMSONT .AND. (GTHOMMET < 5) ) NEXMODES = 1
372:             IF (TWOD) NEXMODES=NEXMODES+NATOMS372:             IF (TWOD) NEXMODES=NEXMODES+NATOMS
373:             IF (FREEZE) THEN373:             IF (FREEZE) THEN
374:                NEXMODES=3*NFREEZE374:                NEXMODES=3*NFREEZE
375:             ENDIF375:             ENDIF
376:             IF (RBAAT) THEN376:             IF (RBAAT) THEN
377:                IF (EFIELDT) THEN377:                IF (EFIELDT) THEN
378:                   NEXMODES = 4378:                   NEXMODES = 4
379:                ELSE379:                ELSE
380:                   NEXMODES = 6380:                   NEXMODES = 6
1055:                   ! jmc put in this test mainly for pathsample purposes...1055:                   ! jmc put in this test mainly for pathsample purposes...
1056:                ENDIF 1056:                ENDIF 
1057:             ENDDO1057:             ENDDO
1058:          ELSE1058:          ELSE
1059:             IF (NENDHESS-NEXMODES.GT.0) THEN1059:             IF (NENDHESS-NEXMODES.GT.0) THEN
1060:                MINFRQ2=LOG(EVALUES(NENDHESS-NEXMODES))1060:                MINFRQ2=LOG(EVALUES(NENDHESS-NEXMODES))
1061:             ELSE1061:             ELSE
1062:                MINFRQ2=1.0D01062:                MINFRQ2=1.0D0
1063:             ENDIF1063:             ENDIF
1064:             EWARN=.FALSE.1064:             EWARN=.FALSE.
1065:             DO I1=1,NENDHESS - NEXMODES1065: ! jwrm2> For GTHOMSON, the eigenvalue list is padded with NATOMS 10^10 values, which occur at the top of the list.
 1066: !        These need to be skipped.
 1067:             IF (GTHOMSONT) THEN
 1068:                J1 = NATOMS + 1
 1069:                J2 = NATOMS + NENDHESS - NEXMODES
 1070:             ELSE
 1071:                J1 = 1
 1072:                J2 = NENDHESS - NEXMODES
 1073:             ENDIF
 1074:             DO I1=J1,J2
1066:                IF (I1.GT.J1) THEN1075:                IF (I1.GT.J1) THEN
1067:                   IF (EVALUES(I1-1).NE.0.0D0) THEN1076:                   IF (EVALUES(I1-1).NE.0.0D0) THEN
1068:                      IF (ABS(EVALUES(I1)/EVALUES(I1-1)).LT.1.0D-2) THEN1077:                      IF (ABS(EVALUES(I1)/EVALUES(I1-1)).LT.1.0D-2) THEN
1069:                         PRINT '(A,G20.10,A,G20.10)',' geopt> WARNING - decrease in magnitude of eigenvalues from ',EVALUES(I1-1),1078:                         PRINT '(A,G20.10,A,G20.10)',' geopt> WARNING - decrease in magnitude of eigenvalues from ',EVALUES(I1-1),
1070:      &                       ' to ',EVALUES(I1)1079:      &                       ' to ',EVALUES(I1)
1071:                         PRINT '(A)',' geopt> WARNING - this could indicate a stationary point of the wrong index'1080:                         PRINT '(A)',' geopt> WARNING - this could indicate a stationary point of the wrong index'
1072:                         EWARN=.TRUE.1081:                         EWARN=.TRUE.
1073:                      ENDIF1082:                      ENDIF
1074:                   ENDIF1083:                   ENDIF
1075:                ENDIF1084:                ENDIF


r28822/main.F 2015-11-17 23:34:38.773829380 +0000 r28821/main.F 2015-11-17 23:34:39.201835127 +0000
163:       NGLY=0163:       NGLY=0
164:       IF (TWOD) THEN164:       IF (TWOD) THEN
165:          NFSTART=1165:          NFSTART=1
166:          NFFINISH=2*NATOMS-3166:          NFFINISH=2*NATOMS-3
167:          KAPPA=2*NATOMS-3167:          KAPPA=2*NATOMS-3
168:          IF (BULKT) THEN168:          IF (BULKT) THEN
169:             NFSTART=1169:             NFSTART=1
170:             NFFINISH=2*NATOMS-2170:             NFFINISH=2*NATOMS-2
171:             KAPPA=2*NATOMS-2171:             KAPPA=2*NATOMS-2
172:          ENDIF172:          ENDIF
173:       ELSE IF (GTHOMSONT) THEN 
174:         NFSTART = 1 
175:         NFFINISH = 3*NATOMS - 3 
176:         KAPPA = 3*NATOMS - 3 
177:       ELSE IF (PHI4MODT) THEN173:       ELSE IF (PHI4MODT) THEN
178:          NFSTART=1174:          NFSTART=1
179:          NFFINISH=NATOMS175:          NFFINISH=NATOMS
180:          KAPPA=NATOMS176:          KAPPA=NATOMS
181:       ELSE IF (PULLT) THEN177:       ELSE IF (PULLT) THEN
182:          NFFINISH=3*NATOMS-4178:          NFFINISH=3*NATOMS-4
183:          KAPPA=3*NATOMS-4179:          KAPPA=3*NATOMS-4
184:       ELSE IF (BULKT) THEN180:       ELSE IF (BULKT) THEN
185:          NFFINISH=3*NATOMS-3181:          NFFINISH=3*NATOMS-3
186:          KAPPA=3*NATOMS-3182:          KAPPA=3*NATOMS-3


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0