hdiff output

r28887/persistence.f90 2015-11-17 23:34:35.801789524 +0000 r28886/persistence.f90 2015-11-17 23:34:36.029792586 +0000
  5: !  superbasin with a lower minimum in it. The corresponding downhill   5: !  superbasin with a lower minimum in it. The corresponding downhill 
  6: !  barrier height is called the culminance of the minimum in question.  6: !  barrier height is called the culminance of the minimum in question.
  7: !  7: !
  8: !  Usual superbasin analysis at fine grained total energies.  8: !  Usual superbasin analysis at fine grained total energies.
  9: !  Could use the sorted list of ts energies?  9: !  Could use the sorted list of ts energies?
 10: !  Check if minimum has a threshold assigned. 10: !  Check if minimum has a threshold assigned.
 11: !  If not, check if it is the lowest minimum in the superbasin. 11: !  If not, check if it is the lowest minimum in the superbasin.
 12: !  If not, assign that threshold as its culminance, and mark it done. 12: !  If not, assign that threshold as its culminance, and mark it done.
 13: ! 13: !
 14: SUBROUTINE PERSISTENCE 14: SUBROUTINE PERSISTENCE
 15: USE COMMONS,ONLY : NMIN, NTS, ETS, EMIN, PLUS, MINUS, EINC, DEBUG, TSTHRESH, FVIBMIN, HORDERMIN, TEMPERATURE, TSTHRESH, & 15: USE COMMONS,ONLY : NMIN, NTS, ETS, EMIN, PLUS, MINUS, EINC, DEBUG, TSTHRESH, PFMIN, FVIBMIN, HORDERMIN, TEMPERATURE, TSTHRESH, &
 16:   &                PEQTHRESH, PERTHRESH 16:   &                PEQTHRESH, PERTHRESH
 17: USE UTILS,ONLY : GETUNIT 17: USE UTILS,ONLY : GETUNIT
 18: IMPLICIT NONE 18: IMPLICIT NONE
 19: DOUBLE PRECISION HIGHESTTS, LOWESTTARG, ETHRESH, BARRIER(NMIN), GMIN, SLOPE, INTERCEPT, CORRELATION, EDISC(NMIN) 19: DOUBLE PRECISION HIGHESTTS, LOWESTTARG, ETHRESH, BARRIER(NMIN), GMIN, SLOPE, INTERCEPT, CORRELATION, EDISC(NMIN)
 20: DOUBLE PRECISION SUMX, SUMY, SUMX2, SUMY2, SUMXY, XTEMP(NMIN), PEQ(NMIN), PFNORM, DUMMY, LOCPFMIN(NMIN) 20: DOUBLE PRECISION SUMX, SUMY, SUMX2, SUMY2, SUMXY, XTEMP(NMIN), PEQ(NMIN), PFNORM, PFMEAN, DUMMY
 21: INTEGER J1, BASIN(NMIN), NBASIN, MINTARG, NGMIN, NCONNECTED, MINID(NMIN) 21: INTEGER J1, BASIN(NMIN), NBASIN, MINTARG, NGMIN, NCONNECTED, MINID(NMIN)
 22: INTEGER BASINS2(NMIN), LUNIT 22: INTEGER BASINS2(NMIN), LUNIT
 23: LOGICAL CHANGED, BASINT(NMIN) 23: LOGICAL CHANGED, BASINT(NMIN)
 24:  24: 
 25: CALL GETBARRIER2(BARRIER,NGMIN,GMIN) 25: CALL GETBARRIER2(BARRIER,NGMIN,GMIN)
 26: PRINT '(A,I6,G20.10)','persistence> global minimum position and energy: ',NGMIN,GMIN 26: PRINT '(A,I6,G20.10)','persistence> global minimum position and energy: ',NGMIN,GMIN
 27: SUMX=0.0D0  27: SUMX=0.0D0 
 28: SUMX2=0.0D0  28: SUMX2=0.0D0 
 29: SUMXY=0.0D0  29: SUMXY=0.0D0 
 30: SUMY=0.0D0  30: SUMY=0.0D0 
 63: DO J1=1,NCONNECTED 63: DO J1=1,NCONNECTED
 64:    PRINT '(I8,3F20.10,I8,F20.10)',J1,EMIN(MINID(J1)),EMIN(MINID(J1))+BARRIER(MINID(J1)),BARRIER(MINID(J1)),MINID(J1), & 64:    PRINT '(I8,3F20.10,I8,F20.10)',J1,EMIN(MINID(J1)),EMIN(MINID(J1))+BARRIER(MINID(J1)),BARRIER(MINID(J1)),MINID(J1), &
 65:   &                        XTEMP(J1) 65:   &                        XTEMP(J1)
 66:    WRITE(LUNIT,'(I8,3F20.10,I8,F20.10)') J1,EMIN(MINID(J1)),EMIN(MINID(J1))+BARRIER(MINID(J1)),BARRIER(MINID(J1)),MINID(J1), & 66:    WRITE(LUNIT,'(I8,3F20.10,I8,F20.10)') J1,EMIN(MINID(J1)),EMIN(MINID(J1))+BARRIER(MINID(J1)),BARRIER(MINID(J1)),MINID(J1), &
 67:   &                        XTEMP(J1) 67:   &                        XTEMP(J1)
 68: ENDDO 68: ENDDO
 69: CLOSE(LUNIT) 69: CLOSE(LUNIT)
 70:  70: 
 71: PFNORM=0.0D0 71: PFNORM=0.0D0
 72: DO J1=1,NMIN 72: DO J1=1,NMIN
 73:    LOCPFMIN(J1) = -EMIN(J1)/TEMPERATURE - FVIBMIN(J1)/2.0D0 - LOG(1.0D0*HORDERMIN(J1)) 73:    PFMIN(J1) = -EMIN(J1)/TEMPERATURE - FVIBMIN(J1)/2.0D0 - LOG(1.0D0*HORDERMIN(J1))
 74:    PFNORM=PFNORM+EXP(LOCPFMIN(J1)-LOCPFMIN(1)) 74:    PFNORM=PFNORM+EXP(PFMIN(J1)-PFMIN(1))
 75: ENDDO 75: ENDDO
 76: PFNORM=LOG(PFNORM)+LOCPFMIN(1) 76: PFNORM=LOG(PFNORM)+PFMIN(1)
 77: DUMMY=0.0D0 77: DUMMY=0.0D0
 78: DO J1=1,NMIN 78: DO J1=1,NMIN
 79:    PEQ(J1)=EXP(LOCPFMIN(J1)-PFNORM) 79:    PEQ(J1)=EXP(PFMIN(J1)-PFNORM)
 80:    DUMMY=DUMMY+PEQ(J1) 80:    DUMMY=DUMMY+PEQ(J1)
 81: ENDDO 81: ENDDO
 82: PRINT '(A,G20.10)','sum of equilibrium occupation probabilities=',DUMMY 82: PRINT '(A,G20.10)','sum of equilibrium occupation probabilities=',DUMMY
 83:  83: 
 84: PRINT '(3(A,G20.10))','data for connected minima with Peq > ',PEQTHRESH,' energy+barrier < ',TSTHRESH,& 84: PRINT '(3(A,G20.10))','data for connected minima with Peq > ',PEQTHRESH,' energy+barrier < ',TSTHRESH,&
 85:                      & ' and deviation > ',PERTHRESH 85:                      & ' and deviation > ',PERTHRESH
 86: LUNIT=GETUNIT() 86: LUNIT=GETUNIT()
 87: OPEN(LUNIT,FILE='min.persist.select',STATUS='UNKNOWN') 87: OPEN(LUNIT,FILE='min.persist.select',STATUS='UNKNOWN')
 88: DO J1=1,NCONNECTED 88: DO J1=1,NCONNECTED
 89:    IF ((PEQ(MINID(J1)).GT.PEQTHRESH).AND.(EMIN(MINID(J1))+BARRIER(MINID(J1)).LT.TSTHRESH).AND.(XTEMP(J1).GT.PERTHRESH)) THEN 89:    IF ((PEQ(MINID(J1)).GT.PEQTHRESH).AND.(EMIN(MINID(J1))+BARRIER(MINID(J1)).LT.TSTHRESH).AND.(XTEMP(J1).GT.PERTHRESH)) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0