hdiff output

r33401/alignutils.f90 2017-10-17 20:30:12.375247340 +0100 r33400/alignutils.f90 2017-10-17 20:30:12.835253519 +0100
369: IMPLICIT NONE369: IMPLICIT NONE
370: INTEGER, INTENT(IN) :: NCOORDS370: INTEGER, INTENT(IN) :: NCOORDS
371: DOUBLE PRECISION, INTENT(IN) :: COORDSA(3*NCOORDS), COORDSB(3*NCOORDS)371: DOUBLE PRECISION, INTENT(IN) :: COORDSA(3*NCOORDS), COORDSB(3*NCOORDS)
372: 372: 
373: DOUBLE PRECISION, INTENT(OUT) :: DISP(3), DIST373: DOUBLE PRECISION, INTENT(OUT) :: DISP(3), DIST
374: 374: 
375: INTEGER J1375: INTEGER J1
376: DOUBLE PRECISION XM, YM, ZM376: DOUBLE PRECISION XM, YM, ZM
377: 377: 
378: ! Calculate average displacement378: ! Calculate average displacement
379: DISP(1:3) = 0.D0 
380: DO J1=1,NCOORDS379: DO J1=1,NCOORDS
381:     XM = COORDSB(3*J1-2) - COORDSA(3*J1-2)380:     XM = COORDSB(3*J1-2) - COORDSA(3*J1-2)
382:     YM = COORDSB(3*J1-1) - COORDSA(3*J1-1)381:     YM = COORDSB(3*J1-1) - COORDSA(3*J1-1)
383:     DISP(1) = DISP(1) + XM - BOXLX*NINT(XM/BOXLX)382:     DISP(1) = DISP(1) + XM - BOXLX*NINT(XM/BOXLX)
384:     DISP(2) = DISP(2) + YM - BOXLY*NINT(YM/BOXLY)383:     DISP(2) = DISP(2) + YM - BOXLY*NINT(YM/BOXLY)
385: ENDDO384: ENDDO
386: 385: 
387: IF (TWOD) THEN386: IF (TWOD) THEN
388:     DISP(3) = 0.D0387:     DISP(3) = 0.D0
389: ELSE388: ELSE


r33401/fastbulk.f90 2017-10-17 20:30:12.603250402 +0100 r33400/fastbulk.f90 2017-10-17 20:30:13.055256473 +0100
259: INTEGER, INTENT(IN) :: NCOORDS259: INTEGER, INTENT(IN) :: NCOORDS
260: DOUBLE PRECISION, INTENT(IN) :: BOXLX,BOXLY,BOXLZ260: DOUBLE PRECISION, INTENT(IN) :: BOXLX,BOXLY,BOXLZ
261: DOUBLE PRECISION, INTENT(OUT) :: KERNELWIDTH261: DOUBLE PRECISION, INTENT(OUT) :: KERNELWIDTH
262: INTEGER, INTENT(OUT) :: NWAVE,NFSPACE262: INTEGER, INTENT(OUT) :: NWAVE,NFSPACE
263: 263: 
264: DOUBLE PRECISION MAXWAVEK264: DOUBLE PRECISION MAXWAVEK
265: 265: 
266: NATOMS=NCOORDS266: NATOMS=NCOORDS
267: KERNELWIDTH = (BOXLX*BOXLY*BOXLZ/NATOMS)**(1.D0/3.D0) / 3.D0267: KERNELWIDTH = (BOXLX*BOXLY*BOXLZ/NATOMS)**(1.D0/3.D0) / 3.D0
268: MAXWAVEK = 1.5 / KERNELWIDTH268: MAXWAVEK = 1.5 / KERNELWIDTH
269: NWAVE = CEILING(1.3 * NATOMS**(1.D0/3.D0), 4)269: NWAVE = CEILING(2*3.14159265359/MIN(BOXLX,BOXLY,BOXLZ)*MAXWAVEK, 4)
270: 270: 
271: 271: 
272: IF((2*NWAVE+1).LE.200) THEN272: IF((2*NWAVE+1).LE.200) THEN
273:     NFSPACE = FASTLEN(4*NWAVE+3)273:     NFSPACE = FASTLEN(4*NWAVE+3)
274: ELSE274: ELSE
275:     ! PROBABLY NOT THE BEST WAY TO CALCULATE THIS!275:     ! PROBABLY NOT THE BEST WAY TO CALCULATE THIS!
276:     NFSPACE = 2**CEILING(LOG(4.D0*NWAVE+3.D0)/LOG(2.D0),4)276:     NFSPACE = 2**CEILING(LOG(4.D0*NWAVE+3.D0)/LOG(2.D0),4)
277: ENDIF277: ENDIF
278: 278: 
279: END SUBROUTINE CALCDEFAULTS279: END SUBROUTINE CALCDEFAULTS


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0