hdiff output

r29521/mie_field.f90 2016-03-16 18:33:15.486874977 +0000 r29520/mie_field.f90 2016-03-16 18:33:15.710877281 +0000
 59:   IF (MIEF_CUTT) THEN 59:   IF (MIEF_CUTT) THEN
 60:      ALLOCATE(MIEF_U_RCUT(NSPECIES(0))) 60:      ALLOCATE(MIEF_U_RCUT(NSPECIES(0)))
 61:      ALLOCATE(MIEF_DUDR_RCUT(NSPECIES(0))) 61:      ALLOCATE(MIEF_DUDR_RCUT(NSPECIES(0)))
 62:      PREF = DBLE(MIEF_N)/DBLE(MIEF_N-MIEF_M)*& 62:      PREF = DBLE(MIEF_N)/DBLE(MIEF_N-MIEF_M)*&
 63:           (DBLE(MIEF_N)/DBLE(MIEF_M))**& 63:           (DBLE(MIEF_N)/DBLE(MIEF_M))**&
 64:           (DBLE(MIEF_M)/DBLE(MIEF_N-MIEF_M)) 64:           (DBLE(MIEF_M)/DBLE(MIEF_N-MIEF_M))
 65:      DO I=1,NSPECIES(0) 65:      DO I=1,NSPECIES(0)
 66:         DUMMY = MIEF_SIG(I)/MIEF_RCUT 66:         DUMMY = MIEF_SIG(I)/MIEF_RCUT
 67:         REP_TERM=DUMMY**MIEF_N 67:         REP_TERM=DUMMY**MIEF_N
 68:         ATT_TERM=DUMMY**MIEF_M 68:         ATT_TERM=DUMMY**MIEF_M
 69:         ! ds656> Had a fucking mistake here for ages!!!!!  69:         MIEF_U_RCUT=PREF*MIEF_EPS(I)*(REP_TERM-ATT_TERM)
 70:         !MIEF_U_RCUT=PREF*MIEF_EPS(I)*(REP_TERM-ATT_TERM) 70:         MIEF_DUDR_RCUT=PREF*MIEF_EPS(I)*( -DBLE(MIEF_N)*REP_TERM + &
 71:         !MIEF_DUDR_RCUT=PREF*MIEF_EPS(I)*( -DBLE(MIEF_N)*REP_TERM + & 
 72:         !     DBLE(MIEF_M)*ATT_TERM ) / MIEF_RCUT 
 73:         ! ds656> Spot the difference... why was the compiler quiet?!? 
 74:         MIEF_U_RCUT(I)=PREF*MIEF_EPS(I)*(REP_TERM-ATT_TERM) 
 75:         MIEF_DUDR_RCUT(I)=PREF*MIEF_EPS(I)*( -DBLE(MIEF_N)*REP_TERM + & 
 76:              DBLE(MIEF_M)*ATT_TERM ) / MIEF_RCUT 71:              DBLE(MIEF_M)*ATT_TERM ) / MIEF_RCUT
 77:      ENDDO 72:      ENDDO
 78:   ENDIF 73:   ENDIF
 79:   ! 74:   !
 80:   RETURN 75:   RETURN
 81:   ! 76:   !
 82: END SUBROUTINE MIEF_INI 77: END SUBROUTINE MIEF_INI
 83: ! 78: !
 84: SUBROUTINE MIEF(X,GRAD,EREAL,GRADT, STRESST) 79: SUBROUTINE MIEF(X,GRAD,EREAL,GRADT, STRESST)
 85:   ! 80:   !
 96:   LOGICAL, INTENT(IN) :: GRADT, STRESST 91:   LOGICAL, INTENT(IN) :: GRADT, STRESST
 97:   ! 92:   !
 98:   INTEGER :: I, J, K, L, I3, ITYPE 93:   INTEGER :: I, J, K, L, I3, ITYPE
 99:   DOUBLE PRECISION :: PREF, DX(3),DIST2,IDIST,SIG,EPS,& 94:   DOUBLE PRECISION :: PREF, DX(3),DIST2,IDIST,SIG,EPS,&
100:        DUMMY, ATT_TERM, REP_TERM, DIST 95:        DUMMY, ATT_TERM, REP_TERM, DIST
101:   ! 96:   !
102:   PREF = DBLE(MIEF_N)/DBLE(MIEF_N-MIEF_M)*& 97:   PREF = DBLE(MIEF_N)/DBLE(MIEF_N-MIEF_M)*&
103:        (DBLE(MIEF_N)/DBLE(MIEF_M))**& 98:        (DBLE(MIEF_N)/DBLE(MIEF_M))**&
104:        (DBLE(MIEF_M)/DBLE(MIEF_N-MIEF_M)) 99:        (DBLE(MIEF_M)/DBLE(MIEF_N-MIEF_M))
105:   !100:   !
106:   DO I=1,NATOMS ! Loop over ALL atoms101:   DO I=1,NATOMS ! Loop over atoms
107:      !102:      !
108:      I3=3*(I-1)103:      I3=3*(I-1)
109:      ITYPE=INVATOMLISTS(I,1)104:      ITYPE=INVATOMLISTS(I,1)
110:      EPS=MIEF_EPS(ITYPE)*PREF105:      EPS=MIEF_EPS(ITYPE)*PREF
111:      SIG=MIEF_SIG(ITYPE)106:      SIG=MIEF_SIG(ITYPE)
112:      !107:      !
113:      DO J=1,MIEF_NSITES ! Loope over Mie sites108:      DO J=1,MIEF_NSITES ! Loope over Mie sites
114:         !109:         !
115:         DIST2=0.0110:         DIST2=0.0
116:         DO K=1,3111:         DO K=1,3


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0