hdiff output

r33568/congrad.f90 2017-12-07 17:30:13.728796893 +0000 r33567/congrad.f90 2017-12-07 17:30:13.988800352 +0000
218: !218: !
219: ! Squared distance between atoms A and B for theta=0 - distance in image 2219: ! Squared distance between atoms A and B for theta=0 - distance in image 2
220: !220: !
221:       DSQ2=G2(1)**2 + G2(2)**2 + G2(3)**2221:       DSQ2=G2(1)**2 + G2(2)**2 + G2(3)**2
222: !222: !
223: ! Squared distance between atoms A and B for theta=Pi/2 - distance in image 1223: ! Squared distance between atoms A and B for theta=Pi/2 - distance in image 1
224: !224: !
225:       DSQ1=G1(1)**2 + G1(2)**2 + G1(3)**2225:       DSQ1=G1(1)**2 + G1(2)**2 + G1(3)**2
226:       DCUT=NREPCUT(J2)**2226:       DCUT=NREPCUT(J2)**2
227:       IF ((DSQ1.GT.DCUT).AND.(DSQ2.GT.DCUT)) CYCLE ! don't look for an internal minimum if both repulsions outside cutoff227:       IF ((DSQ1.GT.DCUT).AND.(DSQ2.GT.DCUT)) CYCLE ! don't look for an internal minimum if both repulsions outside cutoff
228:       IF (ABS(NREPI(J2)-NREPJ(J2)).LT.QCIINTREPMINSEP) THEN ! don't check for internal minimum in distance - atoms too close for chain crossing.228:       IF (ABS(NREPI(J2)-NREPJ(J2)).GT.QCIINTREPMINSEP) THEN ! don't check for internal minimum in distance - atoms too close for chain crossing.
229:          r1apr2bmr2amr1bsq=0.0D0229:          r1apr2bmr2amr1bsq=0.0D0
230:       ELSE230:       ELSE
231: !        WRITE(*,'(A,I6,A,I6,A,2G20.10,A,G20.10)') 'distances in ',J1-1,' and ',J1,' are ',SQRT(DSQ1),SQRT(DSQ2),' cutoff=',NREPCUT(J2)231: !        WRITE(*,'(A,I6,A,I6,A,2G20.10,A,G20.10)') 'distances in ',J1-1,' and ',J1,' are ',SQRT(DSQ1),SQRT(DSQ2),' cutoff=',NREPCUT(J2)
232: !        WRITE(*,'(A,I6,A,2I6,A,6G15.7)') 'image ',J1-1,' atoms ',NREPI(J2),NREPJ(J2),' coords ',XYZ(NI1+1:NI1+3),XYZ(NJ1+1:NJ1+3)232: !        WRITE(*,'(A,I6,A,2I6,A,6G15.7)') 'image ',J1-1,' atoms ',NREPI(J2),NREPJ(J2),' coords ',XYZ(NI1+1:NI1+3),XYZ(NJ1+1:NJ1+3)
233: !        WRITE(*,'(A,I6,A,2I6,A,6G15.7)') 'image ',J1  ,' atoms ',NREPI(J2),NREPJ(J2),' coords ',XYZ(NI2+1:NI2+3),XYZ(NJ2+1:NJ2+3)233: !        WRITE(*,'(A,I6,A,2I6,A,6G15.7)') 'image ',J1  ,' atoms ',NREPI(J2),NREPJ(J2),' coords ',XYZ(NI2+1:NI2+3),XYZ(NJ2+1:NJ2+3)
234:          r1amr1bdr2amr2b=G1(1)*G2(1)+G1(2)*G2(2)+G1(3)*G2(3)234:          r1amr1bdr2amr2b=G1(1)*G2(1)+G1(2)*G2(2)+G1(3)*G2(3)
235:          r1apr2bmr2amr1bsq=DSQ1+DSQ2-2.0D0*r1amr1bdr2amr2b235:          r1apr2bmr2amr1bsq=DSQ1+DSQ2-2.0D0*r1amr1bdr2amr2b
236:       ENDIF236:       ENDIF
237: !237: !
238: ! Is the denominator of the d^2 internal minimum term close to zero?238: ! Is the denominator of the d^2 internal minimum term close to zero?


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0