hdiff output

r30531/stealthy.f90 2016-05-30 13:30:04.330755647 +0100 r30530/stealthy.f90 2016-05-30 13:30:04.522758209 +0100
 22:        TIDU=0 22:        TIDU=0
 23:     END IF 23:     END IF
 24:  24: 
 25:     IF (STEST) THEN 25:     IF (STEST) THEN
 26:        HESS=0 26:        HESS=0
 27:     END IF 27:     END IF
 28:  28: 
 29: !   get the initial box 29: !   get the initial box
 30:  30: 
 31: !   calculate the wave vector 31: !   calculate the wave vector
 32:     KX(1)=2*PI/BULK_BOXVEC(1) 32:     KX(1)=360/BULK_BOXVEC(1)
 33:     KX(2)=0 33:     KX(2)=0
 34:     KX(3)=0 34:     KX(3)=0
 35:     KY(1)=0 35:     KY(1)=0
 36:     KY(2)=2*PI/BULK_BOXVEC(2) 36:     KY(2)=360/BULK_BOXVEC(2)
 37:     KY(3)=0 37:     KY(3)=0
 38:     KZ(1)=0 38:     KZ(1)=0
 39:     KZ(2)=0 39:     KZ(2)=0
 40:     KZ(3)=2*PI/BULK_BOXVEC(3) 40:     KZ(3)=360/BULK_BOXVEC(3)
 41:     VF=BULK_BOXVEC(1)*BULK_BOXVEC(2)*BULK_BOXVEC(3) 41:     VF=BULK_BOXVEC(1)*BULK_BOXVEC(2)*BULK_BOXVEC(3)
 42: !   WRITE(*,*) KI(1), KI(2), KI(3) 42: !   WRITE(*,*) KI(1), KI(2), KI(3)
 43:  43: 
 44: !   periodic boundary condition 44: !   periodic boundary condition
 45:     IF (BULKT) THEN 45:     IF (BULKT) THEN
 46:        DO J1=1, NATOMS 46:        DO J1=1, NATOMS
 47:           J2=3*(J1-1) 47:           J2=3*(J1-1)
 48:           XS(J2+1)=XS(J2+1) - BULK_BOXVEC(1)*DNINT(XS(J2+1)/BULK_BOXVEC(1)) 48:           XS(J2+1)=XS(J2+1) - BULK_BOXVEC(1)*DNINT(XS(J2+1)/BULK_BOXVEC(1))
 49:           XS(J2+2)=XS(J2+2) - BULK_BOXVEC(2)*DNINT(XS(J2+2)/BULK_BOXVEC(2)) 49:           XS(J2+2)=XS(J2+2) - BULK_BOXVEC(2)*DNINT(XS(J2+2)/BULK_BOXVEC(2))
 50:           XS(J2+3)=XS(J2+3) - BULK_BOXVEC(3)*DNINT(XS(J2+3)/BULK_BOXVEC(3)) 50:           XS(J2+3)=XS(J2+3) - BULK_BOXVEC(3)*DNINT(XS(J2+3)/BULK_BOXVEC(3))
 66:              C=0 66:              C=0
 67:              S=0 67:              S=0
 68:              KN(1)=KX(1)*M1 + KY(1)*M2 + KZ(1)*M3 68:              KN(1)=KX(1)*M1 + KY(1)*M2 + KZ(1)*M3
 69:              KN(2)=KX(2)*M1 + KY(2)*M2 + KZ(2)*M3 69:              KN(2)=KX(2)*M1 + KY(2)*M2 + KZ(2)*M3
 70:              KN(3)=KX(3)*M1 + KY(3)*M2 + KZ(3)*M3 70:              KN(3)=KX(3)*M1 + KY(3)*M2 + KZ(3)*M3
 71:  71: 
 72: !            calculation of Re[n(k)] and Im[n(k)] 72: !            calculation of Re[n(k)] and Im[n(k)]
 73:              DO I1=1, NATOMS 73:              DO I1=1, NATOMS
 74:                 I2=3*(I1-1) 74:                 I2=3*(I1-1)
 75:                 KR=KN(1)*XS(I2+1) + KN(2)*XS(I2+2) + KN(3)*XS(I2+3) 75:                 KR=KN(1)*XS(I2+1) + KN(2)*XS(I2+2) + KN(3)*XS(I2+3)
 76:                 CKR(I1)=COS(KR) 76:                 CKR(I1)=COSD(KR)
 77:                 SKR(I1)=SIN(KR) 77:                 SKR(I1)=SIND(KR)
 78:                 !if (M1==0.AND.M2==-3) write(*,'(T1,4(I,1X),F,1X,F)') M1, M2, M3, I1, CKR(I1), SKR(I1) 78:                 !if (M1==0.AND.M2==-3) write(*,'(T1,4(I,1X),F,1X,F)') M1, M2, M3, I1, CKR(I1), SKR(I1)
 79:  79: 
 80:                 C=C+CKR(I1) 80:                 C=C+CKR(I1)
 81:                 S=S+SKR(I1) 81:                 S=S+SKR(I1)
 82: !                if (M1==0.AND.M2==-3) write(*,'(T1,F,1X,F)') C, S 82: !                if (M1==0.AND.M2==-3) write(*,'(T1,F,1X,F)') C, S
 83:              END DO 83:              END DO
 84:              !write(*,*) M1, M2, M3 84:              !write(*,*) M1, M2, M3
 85:              !write(*,*) CKR(1:6) 85:              !write(*,*) CKR(1:6)
 86:              !write(*,*) SKR(1:6) 86:              !write(*,*) SKR(1:6)
 87:              !if (M1==0.AND.M2==-8) then 87:              !if (M1==0.AND.M2==-8) then
147:       HESS(I2+3, I2+3)=HESS(I2+3, I2+3) - 2*K(3)*K(3)/V * CF147:       HESS(I2+3, I2+3)=HESS(I2+3, I2+3) - 2*K(3)*K(3)/V * CF
148: 148: 
149:       DO J1=I1+1, NATOMS149:       DO J1=I1+1, NATOMS
150:          J2=3*(J1-1)150:          J2=3*(J1-1)
151: 151: 
152:          RIJ(1)=X(I2+1) - X(J2+1)152:          RIJ(1)=X(I2+1) - X(J2+1)
153:          RIJ(2)=X(I2+2) - X(J2+2)153:          RIJ(2)=X(I2+2) - X(J2+2)
154:          RIJ(3)=X(I2+3) - X(J2+3)154:          RIJ(3)=X(I2+3) - X(J2+3)
155: 155: 
156:          KRIJ=K(1)*RIJ(1) + K(2)*RIJ(2) + K(3)*RIJ(3)156:          KRIJ=K(1)*RIJ(1) + K(2)*RIJ(2) + K(3)*RIJ(3)
157:          CSKR= COS(KRIJ)*SCA157:          CSKR= COSD(KRIJ)*SCA
158: 158: 
159:          HESS(I2+1, J2+1)=HESS(I2+1, J2+1) + 2*K(1)*K(1)/V * CSKR159:          HESS(I2+1, J2+1)=HESS(I2+1, J2+1) + 2*K(1)*K(1)/V * CSKR
160:          HESS(I2+1, J2+2)=HESS(I2+1, J2+2) + 2*K(1)*K(2)/V * CSKR160:          HESS(I2+1, J2+2)=HESS(I2+1, J2+2) + 2*K(1)*K(2)/V * CSKR
161:          HESS(I2+1, J2+3)=HESS(I2+1, J2+3) + 2*K(1)*K(3)/V * CSKR161:          HESS(I2+1, J2+3)=HESS(I2+1, J2+3) + 2*K(1)*K(3)/V * CSKR
162: 162: 
163:          HESS(I2+2, J2+1)=HESS(I2+2, J2+1) + 2*K(2)*K(1)/V * CSKR163:          HESS(I2+2, J2+1)=HESS(I2+2, J2+1) + 2*K(2)*K(1)/V * CSKR
164:          HESS(I2+2, J2+2)=HESS(I2+2, J2+2) + 2*K(2)*K(2)/V * CSKR164:          HESS(I2+2, J2+2)=HESS(I2+2, J2+2) + 2*K(2)*K(2)/V * CSKR
165:          HESS(I2+2, J2+3)=HESS(I2+2, J2+3) + 2*K(2)*K(3)/V * CSKR165:          HESS(I2+2, J2+3)=HESS(I2+2, J2+3) + 2*K(2)*K(3)/V * CSKR
166: 166: 
167:          HESS(I2+3, J2+1)=HESS(I2+3, J2+1) + 2*K(3)*K(1)/V * CSKR167:          HESS(I2+3, J2+1)=HESS(I2+3, J2+1) + 2*K(3)*K(1)/V * CSKR


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0