hdiff output

r32260/lopermdist.f90 2017-04-05 16:30:12.999026310 +0100 r32259/lopermdist.f90 2017-04-05 16:30:13.283029992 +0100
 39: IMPLICIT NONE 39: IMPLICIT NONE
 40:  40: 
 41: INTEGER, PARAMETER :: MAXIMUMTRIES=10 41: INTEGER, PARAMETER :: MAXIMUMTRIES=10
 42: INTEGER NATOMS, NPERM, PATOMS, NRB, OPNUM,  NORBIT1, NORBIT2, NCHOOSE2, NCHOOSE1, NTRIES, NORBITB1, NORBITB2 42: INTEGER NATOMS, NPERM, PATOMS, NRB, OPNUM,  NORBIT1, NORBIT2, NCHOOSE2, NCHOOSE1, NTRIES, NORBITB1, NORBITB2
 43: INTEGER J3, J4, NDUMMY, LPERM(NATOMS), J1, J2, NOTHER, LPERMBEST(NATOMS), NCHOOSEB1, NCHOOSEB2, & 43: INTEGER J3, J4, NDUMMY, LPERM(NATOMS), J1, J2, NOTHER, LPERMBEST(NATOMS), NCHOOSEB1, NCHOOSEB2, &
 44:         LPERMBESTATOM(NATOMS) 44:         LPERMBESTATOM(NATOMS)
 45: DOUBLE PRECISION DIST2, COORDSA(3*NATOMS), COORDSB(3*NATOMS), DISTANCE, DUMMYA(3*NATOMS), & 45: DOUBLE PRECISION DIST2, COORDSA(3*NATOMS), COORDSB(3*NATOMS), DISTANCE, DUMMYA(3*NATOMS), &
 46:   &              BESTA(3*NATOMS), DUMMYB(3*NATOMS), DUMMY(3*NATOMS), DIST, DSUM 46:   &              BESTA(3*NATOMS), DUMMYB(3*NATOMS), DUMMY(3*NATOMS), DIST, DSUM
 47: DOUBLE PRECISION BOXLX,BOXLY,BOXLZ,WORSTRAD,RMAT(3,3),ENERGY, VNEW(3*NATOMS), DX, DY, DZ, RMS, DBEST, XBEST(3*NATOMS) 47: DOUBLE PRECISION BOXLX,BOXLY,BOXLZ,WORSTRAD,RMAT(3,3),ENERGY, VNEW(3*NATOMS), DX, DY, DZ, RMS, DBEST, XBEST(3*NATOMS)
 48: DOUBLE PRECISION CMXA, CMXB, CMXC, QBEST(4), SITESA(3*NTSITES), SITESB(3*NTSITES) 48: DOUBLE PRECISION CMXA, CMXB, CMXC, QBEST(4), SITESA(3*NTSITES), SITESB(3*NTSITES)
 49: DOUBLE PRECISION ROTA(3,3), ROTINVA(3,3), ROTB(3,3), ROTINVB(3,3), RMATBEST(3,3), TMAT(3,3) 49: DOUBLE PRECISION ROTA(3,3), ROTINVA(3,3), ROTB(3,3), ROTINVB(3,3), RMATBEST(3,3), TMAT(3,3), LPC2, LPC22
 50: DOUBLE PRECISION PVEC(3), RTEMP1(3,3), RTEMP2(3,3) 50: DOUBLE PRECISION PVEC(3), RTEMP1(3,3), RTEMP2(3,3)
 51: LOGICAL DEBUG, TWOD, RIGID, BULKT, PITEST, AOK, BOK, ADDED, PERMUTABLE(NATOMS) 51: LOGICAL DEBUG, TWOD, RIGID, BULKT, PITEST, AOK, BOK, ADDED, PERMUTABLE(NATOMS)
 52: DOUBLE PRECISION PDUMMYA(3*NATOMS), PDUMMYB(3*NATOMS), LDISTANCE, DUMMYC(3*NATOMS), XDUMMY, DUMMYD(3*NATOMS), & 52: DOUBLE PRECISION PDUMMYA(3*NATOMS), PDUMMYB(3*NATOMS), LDISTANCE, DUMMYC(3*NATOMS), XDUMMY, DUMMYD(3*NATOMS), &
 53:    &             LDBEST(NPERMGROUP), LDBESTATOM 53:    &             LDBEST(NPERMGROUP), LDBESTATOM
 54: DOUBLE PRECISION SPDUMMYA(3*NATOMS), SPDUMMYB(3*NATOMS), AINIT, BINIT 54: DOUBLE PRECISION SPDUMMYA(3*NATOMS), SPDUMMYB(3*NATOMS), AINIT, BINIT
 55: INTEGER NEWPERM(NATOMS), ALLPERM(NATOMS), SAVEPERM(NATOMS) 55: INTEGER NEWPERM(NATOMS), ALLPERM(NATOMS), SAVEPERM(NATOMS)
 56: DOUBLE PRECISION TIME0, TIME1 56: DOUBLE PRECISION TIME0, TIME1
 57: DOUBLE PRECISION, ALLOCATABLE :: TEMPA(:), TEMPB(:) 57: DOUBLE PRECISION, ALLOCATABLE :: TEMPA(:), TEMPB(:)
 58: CHARACTER(LEN=5) ZSYMSAVE 58: CHARACTER(LEN=5) ZSYMSAVE
 59: COMMON /SYS/ ZSYMSAVE 59: COMMON /SYS/ ZSYMSAVE
 68:       PRINT '(A)',' lopermdist> WARNING *** RMS for structure A is outside tolerance' 68:       PRINT '(A)',' lopermdist> WARNING *** RMS for structure A is outside tolerance'
 69:    ENDIF 69:    ENDIF
 70:    IF (CHRMMT) CALL UPDATENBONDS(COORDSB) 70:    IF (CHRMMT) CALL UPDATENBONDS(COORDSB)
 71:    CALL POTENTIAL(COORDSB,BINIT,VNEW,.TRUE.,.FALSE.,RMS,.FALSE.,.FALSE.) 71:    CALL POTENTIAL(COORDSB,BINIT,VNEW,.TRUE.,.FALSE.,RMS,.FALSE.,.FALSE.)
 72:    PRINT '(2(A,G25.15))',' initial energy for structure B=             ',BINIT,' RMS=',RMS 72:    PRINT '(2(A,G25.15))',' initial energy for structure B=             ',BINIT,' RMS=',RMS
 73:    IF (RMS-MAX(GMAX,CONVR).GT.1.0D-6) THEN 73:    IF (RMS-MAX(GMAX,CONVR).GT.1.0D-6) THEN
 74:       PRINT '(A)',' lopermdist> WARNING *** RMS for structure B is outside tolerance' 74:       PRINT '(A)',' lopermdist> WARNING *** RMS for structure B is outside tolerance'
 75:    ENDIF 75:    ENDIF
 76: ENDIF 76: ENDIF
 77:  77: 
  78: LPC2=LOCALPERMCUT**2
  79: LPC22=LOCALPERMCUT2**2
 78: DBEST=1.0D100 80: DBEST=1.0D100
 79: PERMUTABLE(1:NATOMS)=.FALSE. 81: PERMUTABLE(1:NATOMS)=.FALSE.
 80: NDUMMY=1 82: NDUMMY=1
 81: DO J1=1,NPERMGROUP 83: DO J1=1,NPERMGROUP
 82:    DO J2=1,NPERMSIZE(J1) 84:    DO J2=1,NPERMSIZE(J1)
 83:       PERMUTABLE(PERMGROUP(NDUMMY+J2-1))=.TRUE. 85:       PERMUTABLE(PERMGROUP(NDUMMY+J2-1))=.TRUE.
 84:       INGROUP(PERMGROUP(NDUMMY+J2-1))=J1 86:       INGROUP(PERMGROUP(NDUMMY+J2-1))=J1
 85:    ENDDO 87:    ENDDO
 86:    NDUMMY=NDUMMY+NPERMSIZE(J1) 88:    NDUMMY=NDUMMY+NPERMSIZE(J1)
 87: ENDDO 89: ENDDO
199:    LDBESTATOM=1.0D100201:    LDBESTATOM=1.0D100
200:    NOTHER=0202:    NOTHER=0
201:    DO J2=1,NATOMS203:    DO J2=1,NATOMS
202:       IF (TRIED(J2).EQ.1) THEN204:       IF (TRIED(J2).EQ.1) THEN
203:          NOTHER=NOTHER+1205:          NOTHER=NOTHER+1
204:          DLIST(NOTHER)=J2206:          DLIST(NOTHER)=J2
205:       ENDIF207:       ENDIF
206:    ENDDO208:    ENDDO
207:    ADDED=.FALSE.209:    ADDED=.FALSE.
208:    outer2: DO J2=1,NATOMS210:    outer2: DO J2=1,NATOMS
209:       IF (DMEAN(J2).GT.LOCALPERMCUT2) THEN211:       IF (DMEAN(J2).GT.LPC22) THEN
210: !        PRINT '(A)',' lopermdist> No more atoms within cutoff'212: !        PRINT '(A)',' lopermdist> No more atoms within cutoff'
211:          GOTO 91213:          GOTO 91
212:       ENDIF214:       ENDIF
213:       IF (TRIED(SORTLIST(J2)).EQ.0) THEN215:       IF (TRIED(SORTLIST(J2)).EQ.0) THEN
214:          ADDED=.TRUE.216:          ADDED=.TRUE.
215:          NOTHER=NOTHER+1217:          NOTHER=NOTHER+1
216:          IF (NOTHER+PATOMS.GT.NATOMS) THEN218:          IF (NOTHER+PATOMS.GT.NATOMS) THEN
217:             PRINT '(A,I6)', &219:             PRINT '(A,I6)', &
218:   & ' lopermdist> ERROR *** number of neighbours plus number of permutable atoms exceeds total for group ',J1220:   & ' lopermdist> ERROR *** number of neighbours plus number of permutable atoms exceeds total for group ',J1
219:             STOP221:             STOP


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0