hdiff output

r33402/alignutils.f90 2017-10-17 20:30:16.811307902 +0100 r33401/alignutils.f90 2017-10-17 20:30:17.251312828 +0100
366: IMPLICIT NONE366: IMPLICIT NONE
367: INTEGER, INTENT(IN) :: NCOORDS367: INTEGER, INTENT(IN) :: NCOORDS
368: DOUBLE PRECISION, INTENT(IN) :: COORDSA(3*NCOORDS), COORDSB(3*NCOORDS)368: DOUBLE PRECISION, INTENT(IN) :: COORDSA(3*NCOORDS), COORDSB(3*NCOORDS)
369: 369: 
370: DOUBLE PRECISION, INTENT(OUT) :: DISP(3), DIST370: DOUBLE PRECISION, INTENT(OUT) :: DISP(3), DIST
371: 371: 
372: INTEGER J1372: INTEGER J1
373: DOUBLE PRECISION XM, YM, ZM373: DOUBLE PRECISION XM, YM, ZM
374: 374: 
375: ! Calculate average displacement375: ! Calculate average displacement
376: DISP(1:3) = 0.D0 
377: DO J1=1,NCOORDS376: DO J1=1,NCOORDS
378:     XM = COORDSB(3*J1-2) - COORDSA(3*J1-2)377:     XM = COORDSB(3*J1-2) - COORDSA(3*J1-2)
379:     YM = COORDSB(3*J1-1) - COORDSA(3*J1-1)378:     YM = COORDSB(3*J1-1) - COORDSA(3*J1-1)
380:     DISP(1) = DISP(1) + XM - BULK_BOXVEC(1)*NINT(XM/BULK_BOXVEC(1))379:     DISP(1) = DISP(1) + XM - BULK_BOXVEC(1)*NINT(XM/BULK_BOXVEC(1))
381:     DISP(2) = DISP(2) + YM - BULK_BOXVEC(2)*NINT(YM/BULK_BOXVEC(2))380:     DISP(2) = DISP(2) + YM - BULK_BOXVEC(2)*NINT(YM/BULK_BOXVEC(2))
382: ENDDO381: ENDDO
383: 382: 
384: IF (TWOD) THEN383: IF (TWOD) THEN
385:     DISP(3) = 0.D0384:     DISP(3) = 0.D0
386: ELSE385: ELSE


r33402/fastbulk.f90 2017-10-17 20:30:17.031309876 +0100 r33401/fastbulk.f90 2017-10-17 20:30:17.471315783 +0100
255: INTEGER, INTENT(IN) :: NCOORDS255: INTEGER, INTENT(IN) :: NCOORDS
256: DOUBLE PRECISION, INTENT(IN) :: BOXLX,BOXLY,BOXLZ256: DOUBLE PRECISION, INTENT(IN) :: BOXLX,BOXLY,BOXLZ
257: DOUBLE PRECISION, INTENT(OUT) :: KERNELWIDTH257: DOUBLE PRECISION, INTENT(OUT) :: KERNELWIDTH
258: INTEGER, INTENT(OUT) :: NWAVE,NFSPACE258: INTEGER, INTENT(OUT) :: NWAVE,NFSPACE
259: 259: 
260: DOUBLE PRECISION MAXWAVEK260: DOUBLE PRECISION MAXWAVEK
261: 261: 
262: NATOMS=NCOORDS262: NATOMS=NCOORDS
263: KERNELWIDTH = (BOXLX*BOXLY*BOXLZ/NATOMS)**(1.D0/3.D0) / 3.D0263: KERNELWIDTH = (BOXLX*BOXLY*BOXLZ/NATOMS)**(1.D0/3.D0) / 3.D0
264: MAXWAVEK = 1.5 / KERNELWIDTH264: MAXWAVEK = 1.5 / KERNELWIDTH
265: NWAVE = CEILING(1.3 * NATOMS**(1.D0/3.D0), 4)265: NWAVE = CEILING(2*3.14159265359/MIN(BOXLX,BOXLY,BOXLZ)*MAXWAVEK, 4)
266: 266: 
267: 267: 
268: IF((2*NWAVE+1).LE.200) THEN268: IF((2*NWAVE+1).LE.200) THEN
269:     NFSPACE = FASTLEN(4*NWAVE+3)269:     NFSPACE = FASTLEN(4*NWAVE+3)
270: ELSE270: ELSE
271:     ! PROBABLY NOT THE BEST WAY TO CALCULATE THIS!271:     ! PROBABLY NOT THE BEST WAY TO CALCULATE THIS!
272:     NFSPACE = 2**CEILING(LOG(4.D0*NWAVE+3.D0)/LOG(2.D0),4)272:     NFSPACE = 2**CEILING(LOG(4.D0*NWAVE+3.D0)/LOG(2.D0),4)
273: ENDIF273: ENDIF
274: 274: 
275: END SUBROUTINE CALCDEFAULTS275: END SUBROUTINE CALCDEFAULTS


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0