hdiff output

r29123/quench.F 2015-11-17 23:33:19.916771853 +0000 r29122/quench.F 2015-11-17 23:33:20.120774589 +0000
318: ! Added timing for LBFGS call318: ! Added timing for LBFGS call
319:             CALL CPU_TIME(T_LBFGS_START)319:             CALL CPU_TIME(T_LBFGS_START)
320:             CALL MYLBFGS(NOPT,MUPDATE,P,.FALSE.,GMAX,CFLAG,EREAL,MAXIT,ITER,.TRUE.,NP)320:             CALL MYLBFGS(NOPT,MUPDATE,P,.FALSE.,GMAX,CFLAG,EREAL,MAXIT,ITER,.TRUE.,NP)
321:             CALL CPU_TIME(T_LBFGS_END)321:             CALL CPU_TIME(T_LBFGS_END)
322:             IF (FEBHT) THEN322:             IF (FEBHT) THEN
323:                WRITE(MYUNIT, '(A, F10.3)') 'Time to minimise (s):', T_LBFGS_END - T_LBFGS_START323:                WRITE(MYUNIT, '(A, F10.3)') 'Time to minimise (s):', T_LBFGS_END - T_LBFGS_START
324:             ENDIF324:             ENDIF
325:          ENDIF325:          ENDIF
326:          IF (EVAPREJECT) RETURN326:          IF (EVAPREJECT) RETURN
327:       ENDIF327:       ENDIF
328:       !328: 
 329: ! ds656> If (relative) chemical potentials are present, then we tack
 330: !        them on to the energy for semi-grand canonical BH.
 331:       IF (SEMIGRAND_MUT) THEN     
 332:          DO J1=2, NSPECIES(0)
 333:             EREAL = EREAL - NSPECIES(J1)*SEMIGRAND_MU(J1)
 334:          ENDDO
 335:       ENDIF
 336: 
329:       IF (FEBHT .AND. CFLAG) THEN337:       IF (FEBHT .AND. CFLAG) THEN
330:       ! Calculate the free energy338:       ! Calculate the free energy
331:          NUM_ZERO_EVS=6339:          NUM_ZERO_EVS=6
332:          IF (NATOMS.EQ.2) NUM_ZERO_EVS=5340:          IF (NATOMS.EQ.2) NUM_ZERO_EVS=5
333:          IF (ALLOCATED(HESS)) DEALLOCATE(HESS)341:          IF (ALLOCATED(HESS)) DEALLOCATE(HESS)
334:          ALLOCATE(HESS(3*NATOMS,3*NATOMS))342:          ALLOCATE(HESS(3*NATOMS,3*NATOMS))
335:          CALL POTENTIAL(P,GRAD,EREAL,.TRUE.,.TRUE.)343:          CALL POTENTIAL(P,GRAD,EREAL,.TRUE.,.TRUE.)
336:          CALL MASSWT()344:          CALL MASSWT()
337: #ifdef __SPARSE345: #ifdef __SPARSE
338:          IF (SPARSET) THEN346:          IF (SPARSET) THEN
478:          486:          
479:          IF (FETEMP .GT. 0.D0) THEN487:          IF (FETEMP .GT. 0.D0) THEN
480:             IF (DEBUG) WRITE(MYUNIT, '(A,F20.12)') 'Potential energy=', FEBH_POT_ENE488:             IF (DEBUG) WRITE(MYUNIT, '(A,F20.12)') 'Potential energy=', FEBH_POT_ENE
481:             EREAL=-FETEMP*(SYMFCTR+LOG(1.0D0/HORDER) - EREAL/FETEMP + (3*NATOMS-NUM_ZERO_EVS)*DLOG(FETEMP) - LOG_PROD_EV + ITDET)489:             EREAL=-FETEMP*(SYMFCTR+LOG(1.0D0/HORDER) - EREAL/FETEMP + (3*NATOMS-NUM_ZERO_EVS)*DLOG(FETEMP) - LOG_PROD_EV + ITDET)
482:             IF (DEBUG) WRITE(MYUNIT, '(A,F20.12)') 'Harmonic superposition contribution=', EREAL - FEBH_POT_ENE490:             IF (DEBUG) WRITE(MYUNIT, '(A,F20.12)') 'Harmonic superposition contribution=', EREAL - FEBH_POT_ENE
483:          END IF491:          END IF
484:       ELSE IF (FEBHT .AND. (.NOT. CFLAG)) THEN492:       ELSE IF (FEBHT .AND. (.NOT. CFLAG)) THEN
485:          WRITE(MYUNIT, '(A)') 'Quench did not converge, not calculating free energy and adding 1E10 to energy.'493:          WRITE(MYUNIT, '(A)') 'Quench did not converge, not calculating free energy and adding 1E10 to energy.'
486:          EREAL=EREAL + 1.0D10494:          EREAL=EREAL + 1.0D10
487:       ENDIF495:       ENDIF
488:  
489: ! ds656> If (relative) chemical potentials are present, then we tack 
490: !        them on to the (free) energy for semi-grand canonical (FE)BH. 
491:       IF (SEMIGRAND_MUT) THEN 
492:          DUMMY=0.0D0 
493:          DO J1=2, NSPECIES(0) 
494:             DUMMY = DUMMY + NSPECIES(J1)*SEMIGRAND_MU(J1) 
495:          ENDDO 
496:          EREAL = EREAL - DUMMY 
497:          IF(FEBHT) FEBH_POT_ENE = FEBH_POT_ENE - DUMMY  
498:       ENDIF 
499:  
500:       POTEL=EREAL496:       POTEL=EREAL
501: !fh301>{{{ 497: !fh301>{{{ 
502: ! for more effective restraints just compare the CamShift part of the energies498: ! for more effective restraints just compare the CamShift part of the energies
503: !      IF ((CHEMSHIFT2) .AND. (DABS(ENERGYCAMSHIFT).GE.0.0001)) THEN499: !      IF ((CHEMSHIFT2) .AND. (DABS(ENERGYCAMSHIFT).GE.0.0001)) THEN
504: !        FULLENERGY=POTEL500: !        FULLENERGY=POTEL
505: !        POTEL=ENERGYCAMSHIFT501: !        POTEL=ENERGYCAMSHIFT
506: !      ENDIF502: !      ENDIF
507: !fh301>}}}503: !fh301>}}}
508: 504: 
509:       IF (CFLAG) QDONE=1505:       IF (CFLAG) QDONE=1


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0