hdiff output

r31924/rbperm.f90 2017-02-16 09:30:11.555578418 +0000 r31923/rbperm.f90 2017-02-16 09:30:12.175586785 +0000
 32: !     centre of coordinates of COORDSA will be the same as for COORDSB. 32: !     centre of coordinates of COORDSA will be the same as for COORDSB.
 33: ! 33: !
 34: !     ---------------------------------------------------------------------------------------------- 34: !     ----------------------------------------------------------------------------------------------
 35:   35:  
 36:       SUBROUTINE RBMINPERMDIST(COORDSB,COORDSA,DISTANCE,DIST2,QBEST,RMATBEST,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT,SITESB,SITESA) 36:       SUBROUTINE RBMINPERMDIST(COORDSB,COORDSA,DISTANCE,DIST2,QBEST,RMATBEST,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT,SITESB,SITESA)
 37:  37: 
 38: !     returns DISTANCE as the actual distance, rather than the squared distance 38: !     returns DISTANCE as the actual distance, rather than the squared distance
 39:  39: 
 40:       USE COMMONS, ONLY : NATOMS, NRBSITES  40:       USE COMMONS, ONLY : NATOMS, NRBSITES 
 41:       USE KEY,ONLY : NTSITES, NPERMGROUP, NPERMSIZE, PERMGROUP, NSETS, GEOMDIFFTOL, BESTPERM, EFIELDT, PAPT, PTSTSTT, & 41:       USE KEY,ONLY : NTSITES, NPERMGROUP, NPERMSIZE, PERMGROUP, NSETS, GEOMDIFFTOL, BESTPERM, EFIELDT, PAPT, PTSTSTT, &
 42:      &               RBOPS, NRBGROUP, RBSYMT, DBPTDT, MULTISITEPYT, PYGPERIODICT  42:      &               RBOPS, NRBGROUP, RBSYMT, DBPTDT, MULTISITEPYT 
 43:  43: 
 44:       IMPLICIT NONE 44:       IMPLICIT NONE
 45:  45: 
 46:       INTEGER, PARAMETER :: MAXIMUMTRIES = 100 46:       INTEGER, PARAMETER :: MAXIMUMTRIES = 100
 47:       INTEGER            :: NPERM, PATOMS, NTRIES, NSIZE, J1, J2, J3, NRB !jwrm2> unused, INFO, JMAX, LOCMAX(1) 47:       INTEGER            :: NPERM, PATOMS, NTRIES, NSIZE, J1, J2, J3, NRB !jwrm2> unused, INFO, JMAX, LOCMAX(1)
 48:       INTEGER            :: INVERT, NORBIT1, NORBIT2, PERM(NATOMS), NCHOOSE2, NDUMMY, LPERM(NATOMS), NCHOOSE1 48:       INTEGER            :: INVERT, NORBIT1, NORBIT2, PERM(NATOMS), NCHOOSE2, NDUMMY, LPERM(NATOMS), NCHOOSE1
 49:       INTEGER            :: NEWPERM(NATOMS), ALLPERM(NATOMS) 49:       INTEGER            :: NEWPERM(NATOMS), ALLPERM(NATOMS)
 50:       DOUBLE PRECISION   :: COORDSA(3*NATOMS), COORDSB(3*NATOMS), DISTANCE, DIST2 !jwrm2> unused, DISTWP, TEMPA(9*NATOMS)   50:       DOUBLE PRECISION   :: COORDSA(3*NATOMS), COORDSB(3*NATOMS), DISTANCE, DIST2 !jwrm2> unused, DISTWP, TEMPA(9*NATOMS)  
 51:       DOUBLE PRECISION   :: DUMMYA(3*NATOMS), DUMMYB(3*NATOMS), DUMMY(3*NATOMS) !jwrm2> unused, DUMMYWP(3*NATOMS) 51:       DOUBLE PRECISION   :: DUMMYA(3*NATOMS), DUMMYB(3*NATOMS), DUMMY(3*NATOMS) !jwrm2> unused, DUMMYWP(3*NATOMS)
 52:       DOUBLE PRECISION   :: XA(3*NTSITES),  XB(3*NTSITES), XTMP(3*NATOMS) !jwrm2> unused, XBS(3*NTSITES) 52:       DOUBLE PRECISION   :: XA(3*NTSITES),  XB(3*NTSITES), XTMP(3*NATOMS) !jwrm2> unused, XBS(3*NTSITES)
 68:       COORDSCOMA(:) = COORDSA(1:3*NATOMS/2) 68:       COORDSCOMA(:) = COORDSA(1:3*NATOMS/2)
 69:       COORDSCOMB(:) = COORDSB(1:3*NATOMS/2) 69:       COORDSCOMB(:) = COORDSB(1:3*NATOMS/2)
 70:  70: 
 71:       IF (DEBUG) THEN 71:       IF (DEBUG) THEN
 72:          WRITE(*,'(3F20.10)')COORDSA(1:3*NATOMS) 72:          WRITE(*,'(3F20.10)')COORDSA(1:3*NATOMS)
 73:          WRITE(*,*) 73:          WRITE(*,*)
 74:          WRITE(*,'(3F20.10)')COORDSB(1:3*NATOMS) 74:          WRITE(*,'(3F20.10)')COORDSB(1:3*NATOMS)
 75:       ENDIF 75:       ENDIF
 76:  76: 
 77: !     do not call minpermdistrbCOM for multisite PY  77: !     do not call minpermdistrbCOM for multisite PY 
 78:       IF(MULTISITEPYT.OR.PYGPERIODICT) GOTO 100 78:       IF(MULTISITEPYT) GOTO 100
 79:         NATOMS = NATOMS/2 79:         NATOMS = NATOMS/2
 80:  80: 
 81:         CALL MINPERMDISTRBCOM(COORDSCOMB,COORDSCOMA,DISTANCE,DIST2,QBEST,RMATBEST,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT) 81:         CALL MINPERMDISTRBCOM(COORDSCOMB,COORDSCOMA,DISTANCE,DIST2,QBEST,RMATBEST,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT)
 82:  82: 
 83:         NATOMS = 2*NATOMS 83:         NATOMS = 2*NATOMS
 84:  84: 
 85:         IF (SQRT(DISTANCE) <= GEOMDIFFTOL) THEN 85:         IF (SQRT(DISTANCE) <= GEOMDIFFTOL) THEN
 86:          DISTANCE = SQRT(DISTANCE) 86:          DISTANCE = SQRT(DISTANCE)
 87:          IF (DEBUG) PRINT '(A)',' rbpermdist> minpermdistrbcom suggests identical' 87:          IF (DEBUG) PRINT '(A)',' rbpermdist> minpermdistrbcom suggests identical'
 88:          CMAX = 0.0D0; CMAY = 0.0D0; CMAZ = 0.0D0 88:          CMAX = 0.0D0; CMAY = 0.0D0; CMAZ = 0.0D0


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0