hdiff output

r29834/mymylbfgs.f 2016-01-25 15:30:07.172397890 +0000 r29833/mymylbfgs.f 2016-01-25 15:30:07.400400943 +0000
432:          STP=MIN(1.0D0/GNORM,GNORM)432:          STP=MIN(1.0D0/GNORM,GNORM)
433:       ELSE 433:       ELSE 
434:          BOUND=ITER434:          BOUND=ITER
435:          IF (ITER.GT.M) BOUND=M435:          IF (ITER.GT.M) BOUND=M
436:          YS= DDOT(N,W(IYPT+NPT+1),1,W(ISPT+NPT+1),1)436:          YS= DDOT(N,W(IYPT+NPT+1),1,W(ISPT+NPT+1),1)
437: !437: !
438: !  Update estimate of diagonal inverse Hessian elements438: !  Update estimate of diagonal inverse Hessian elements
439: !439: !
440:          IF (.NOT.DIAGCO) THEN440:          IF (.NOT.DIAGCO) THEN
441:             YY= DDOT(N,W(IYPT+NPT+1),1,W(IYPT+NPT+1),1)441:             YY= DDOT(N,W(IYPT+NPT+1),1,W(IYPT+NPT+1),1)
442:             IF (ABS(YY).LT.1.0D-30) THEN442:             IF (ABS(YY).LT.1.0D-20) THEN
443:                YY=SIGN(1.0D-30,YY)443:                YY=SIGN(1.0D-20,YY)
444:                IF (DEBUG) WRITE(MYUNIT,'(A,G20.10,A)') 'WARNING, resetting YY to ',YY,' in mymylbfgs'444:                IF (DEBUG) WRITE(MYUNIT,'(A,G20.10,A)') 'WARNING, resetting YY to ',YY,' in mymylbfgs'
445:             ENDIF445:             ENDIF
446:             IF (ABS(YS).LT.1.0D-30) THEN446:             IF (ABS(YS).LT.1.0D-20) THEN
447:                YS=SIGN(1.0D-30,YS)447:                YS=SIGN(1.0D-20,YS)
448:                IF (DEBUG) WRITE(MYUNIT,'(A,G20.10,A)') 'WARNING, resetting YS to ',YS,' in mymylbfgs'448:                IF (DEBUG) WRITE(MYUNIT,'(A,G20.10,A)') 'WARNING, resetting YS to ',YS,' in mymylbfgs'
449:             ENDIF449:             ENDIF
450: !           WRITE(MYUNIT,'(A,2G20.10)') 'YS/YY=',YS/YY450: !           WRITE(MYUNIT,'(A,2G20.10)') 'YS/YY=',YS/YY
451:             DO J1=1,N451:             DO J1=1,N
452: !              DIAG(J1)= ABS(YS/YY) ! messes up after step reversals!452: !              DIAG(J1)= ABS(YS/YY) ! messes up after step reversals!
453:                DIAG(J1)= YS/YY453:                DIAG(J1)= YS/YY
454:             ENDDO454:             ENDDO
455:          ELSE455:          ELSE
456:             WRITE(MYUNIT,'(A)') 'using estimate of the inverse diagonal elements'456:             WRITE(MYUNIT,'(A)') 'using estimate of the inverse diagonal elements'
457:             DO J1=1,N457:             DO J1=1,N


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0