hdiff output

r29772/commons.f90 2016-01-19 18:30:13.128938275 +0000 r29771/commons.f90 2016-01-19 18:30:14.316954065 +0000
 98:      &        BARRIERSHORT, FREEZE, RATESHORT, DUMMYRUNT, REWEIGHTT, REGROUPFREET, RFMULTIT, REGROUPFREEABT, READMINT, & 98:      &        BARRIERSHORT, FREEZE, RATESHORT, DUMMYRUNT, REWEIGHTT, REGROUPFREET, RFMULTIT, REGROUPFREEABT, READMINT, &
 99:      &        DUMPGROUPST, FREEPAIRT, KSHORTESTPATHST, KSHORT_FULL_PRINTT, DIJINITFLYT, BHINTERPT, ICINTERPT, & 99:      &        DUMPGROUPST, FREEPAIRT, KSHORTESTPATHST, KSHORT_FULL_PRINTT, DIJINITFLYT, BHINTERPT, ICINTERPT, &
100:      &        DUMMYTST, DOCKT, DSTAGE(6), USEPAIRST, LOWESTFRQT, BISECTT, NGTDISCONNECTALL, ANGLEAXIS2, TFOLDT, &100:      &        DUMMYTST, DOCKT, DSTAGE(6), USEPAIRST, LOWESTFRQT, BISECTT, NGTDISCONNECTALL, ANGLEAXIS2, TFOLDT, &
101:      &        SLURMT, INDEXCOSTFUNCTION, CVT, DOST, IMFRQT, CLOSEFILEST, PULLT, FRICTIONT, ATOMMATCHFULL, &101:      &        SLURMT, INDEXCOSTFUNCTION, CVT, DOST, IMFRQT, CLOSEFILEST, PULLT, FRICTIONT, ATOMMATCHFULL, &
102:      &        INTCONSTRAINTT, CHECKCONINT, INTLJT, INTERPCOSTFUNCTION, REMOVEUNCONNECTEDT, ATOMMATCHDIST, &102:      &        INTCONSTRAINTT, CHECKCONINT, INTLJT, INTERPCOSTFUNCTION, REMOVEUNCONNECTEDT, ATOMMATCHDIST, &
103:      &        DBPT, DBPTDT, DMBLPYT, EFIELDT, MSSTOCKT, NTIPT, PAHAT, PAPT, PATCHYDT, STOCKAAT, RBAAT, RBSYMT, TRAPT, SILANET, &103:      &        DBPT, DBPTDT, DMBLPYT, EFIELDT, MSSTOCKT, NTIPT, PAHAT, PAPT, PATCHYDT, STOCKAAT, RBAAT, RBSYMT, TRAPT, SILANET, &
104:      &        OHCELLT, INTFREEZET, LPERMDIST, PBST, RANDOMMETRICT, SSHT, ALLTST, USERPOTT, CHECKMINT, &104:      &        OHCELLT, INTFREEZET, LPERMDIST, PBST, RANDOMMETRICT, SSHT, ALLTST, USERPOTT, CHECKMINT, &
105:      &        CHECKTST, CHECKSPT, FROMLOWESTT, ADDMINXYZT, MACHINE, RATESCYCLET, NOINVERSION, NEWCONNECTIONST, NIMET, NIHEAM7T, &105:      &        CHECKTST, CHECKSPT, FROMLOWESTT, ADDMINXYZT, MACHINE, RATESCYCLET, NOINVERSION, NEWCONNECTIONST, NIMET, NIHEAM7T, &
106:      &        NIH2LEPST, DISTANCET, RATETARGETT, TARGETHIT, ALLOWABT, MICROTHERMT, RFKMCT, REGROUPKMCT, ONEREGROUPT, PHI4MODT, &106:      &        NIH2LEPST, DISTANCET, RATETARGETT, TARGETHIT, ALLOWABT, MICROTHERMT, RFKMCT, REGROUPKMCT, ONEREGROUPT, PHI4MODT, &
107:      &        PERSISTT, REGROUPPERSISTT, NOLABELST, SHANNONT, MAKEPAIRS, SKIPPAIRST, PERSISTAPPROXT, ALLCOMPONENTST, &107:      &        PERSISTT, REGROUPPERSISTT, NOLABELST, SHANNONT, MAKEPAIRS, SKIPPAIRST, PERSISTAPPROXT, ALLCOMPONENTST, &
108:      &        SHANNONRT, SHANNONZT, CUDAT, DIJPRUNET, PRINTSUMMARYT108:      &        SHANNONRT, SHANNONZT, CUDAT, DIJPRUNET
109: 109: 
110:       LOGICAL, ALLOCATABLE :: SHIFTABLE(:)110:       LOGICAL, ALLOCATABLE :: SHIFTABLE(:)
111:       CHARACTER(LEN=80) COORDSLIGANDSTR, COORDSCOMPLEXSTR, COORDSPROTEINSTR111:       CHARACTER(LEN=80) COORDSLIGANDSTR, COORDSCOMPLEXSTR, COORDSPROTEINSTR
112:       CHARACTER(LEN=80) EXEC,EXECGMIN112:       CHARACTER(LEN=80) EXEC,EXECGMIN
113:       CHARACTER(LEN=80) PATHNAME, MINNAME, ADDMINXYZNAME, ALLCOMPS113:       CHARACTER(LEN=80) PATHNAME, MINNAME, ADDMINXYZNAME, ALLCOMPS
114:       CHARACTER(LEN=150) COPYFILES114:       CHARACTER(LEN=150) COPYFILES
115:       CHARACTER(LEN=80) USEPAIRSFILE115:       CHARACTER(LEN=80) USEPAIRSFILE
116:       CHARACTER(LEN=80) MAKEPAIRSFILE116:       CHARACTER(LEN=80) MAKEPAIRSFILE
117:       CHARACTER(LEN=2) DIRECTION117:       CHARACTER(LEN=2) DIRECTION
118:       CHARACTER(LEN=5) UNCONNECTEDS118:       CHARACTER(LEN=5) UNCONNECTEDS


r29772/keywords.f 2016-01-19 18:30:13.372941518 +0000 r29771/keywords.f 2016-01-19 18:30:14.540957043 +0000
263:       SHANNONRT=.FALSE.263:       SHANNONRT=.FALSE.
264:       SHANNONZT=.FALSE.264:       SHANNONZT=.FALSE.
265:       NPEQ=100265:       NPEQ=100
266:       MICROTHERMT=.FALSE.266:       MICROTHERMT=.FALSE.
267:       DOST=.FALSE.267:       DOST=.FALSE.
268:       CHECKCONNECTIONST=.FALSE.268:       CHECKCONNECTIONST=.FALSE.
269:       NEWCONNECTIONST=.FALSE.269:       NEWCONNECTIONST=.FALSE.
270:       CONNMINSTART=1270:       CONNMINSTART=1
271:       CLOSEFILEST=.FALSE.271:       CLOSEFILEST=.FALSE.
272:       PULLT=.FALSE.272:       PULLT=.FALSE.
273:       PRINTSUMMARYT=.FALSE. 
274: !     DC430 >273: !     DC430 >
275:       RBAAT  = .FALSE.274:       RBAAT  = .FALSE.
276:       RBSYMT = .FALSE.275:       RBSYMT = .FALSE.
277:       FRICTIONT=.FALSE.276:       FRICTIONT=.FALSE.
278:       GAMMAFRICTION=0.0D0277:       GAMMAFRICTION=0.0D0
279:       REMOVEUNCONNECTEDT=.FALSE.278:       REMOVEUNCONNECTEDT=.FALSE.
280:       UNCONNECTEDS='AB'279:       UNCONNECTEDS='AB'
281:       OHCELLT=.FALSE.280:       OHCELLT=.FALSE.
282:       ATOMMATCHDIST=.FALSE.281:       ATOMMATCHDIST=.FALSE.
283:       ATOMMATCHFULL=.FALSE.282:       ATOMMATCHFULL=.FALSE.
1769:          CALL READF(PLANCK)1768:          CALL READF(PLANCK)
1770:       ELSE IF (WORD.EQ.'PSCALE') THEN1769:       ELSE IF (WORD.EQ.'PSCALE') THEN
1771:          CALL READF(PSCALE)1770:          CALL READF(PSCALE)
1772: !1771: !
1773: ! PHI4MOD for mean field phi^4 model1772: ! PHI4MOD for mean field phi^4 model
1774: !1773: !
1775:       ELSE IF (WORD.EQ.'PHI4MOD') THEN1774:       ELSE IF (WORD.EQ.'PHI4MOD') THEN
1776:          PHI4MODT=.TRUE.1775:          PHI4MODT=.TRUE.
1777:          IF (NITEMS.GT.1) CALL READF(JPARAM)1776:          IF (NITEMS.GT.1) CALL READF(JPARAM)
1778:          PRINT '(A,G20.10)',' keywords> PHI4 mean field model with J=',JPARAM1777:          PRINT '(A,G20.10)',' keywords> PHI4 mean field model with J=',JPARAM
1779: !1778: 
1780: ! Print some summary info about the kinetic transition network to 
1781: ! standard out 
1782: ! 
1783:       ELSE IF (WORD.EQ.'PRINTSUMMARY') THEN 
1784:           PRINTSUMMARYT=.TRUE. 
1785: C1779: C
1786: C  Try connecting minima that have the largest values of equilibrium occupation probability times1780: C  Try connecting minima that have the largest values of equilibrium occupation probability times
1787: C  waiting time to minima on the fastest path.1781: C  waiting time to minima on the fastest path.
1788: C  Don't bother - UNTRAP with free energy sorting should do the job.1782: C  Don't bother - UNTRAP with free energy sorting should do the job.
1789: C1783: C
1790: !     ELSE IF (WORD.EQ.'PTAU') THEN1784: !     ELSE IF (WORD.EQ.'PTAU') THEN
1791: !        PTAUT=.TRUE.1785: !        PTAUT=.TRUE.
1792: !        DIJKSTRAT=.TRUE.1786: !        DIJKSTRAT=.TRUE.
1793:       ELSE IF (WORD.EQ.'PULL') THEN1787:       ELSE IF (WORD.EQ.'PULL') THEN
1794:          PULLT=.TRUE.1788:          PULLT=.TRUE.


r29772/main.F 2016-01-19 18:30:13.628944920 +0000 r29771/main.F 2016-01-19 18:30:14.768960073 +0000
532:             PRINT '(A)','main> Job submission using srun is not implemented for keyword DOCK. '532:             PRINT '(A)','main> Job submission using srun is not implemented for keyword DOCK. '
533:             STOP533:             STOP
534:          END IF534:          END IF
535:          CALL DOCK535:          CALL DOCK
536:       END IF536:       END IF
537:       IF (SIST) THEN537:       IF (SIST) THEN
538:          CALL SETUP_SIS538:          CALL SETUP_SIS
539:       ELSE539:       ELSE
540:          CALL SETUP540:          CALL SETUP
541:       ENDIF541:       ENDIF
542:       IF (PRINTSUMMARYT) CALL PRINTSUMMARY 
543:       IF (CHECKMINT) THEN542:       IF (CHECKMINT) THEN
544:          IF (CHECKSPS.LE.0) CHECKSPS=1543:          IF (CHECKSPS.LE.0) CHECKSPS=1
545:          IF (CHECKSPF.LE.0) CHECKSPF=NMIN544:          IF (CHECKSPF.LE.0) CHECKSPF=NMIN
546:          PRINT '(2(A,I10))','Checking minima from ',CHECKSPS,' to ',CHECKSPF545:          PRINT '(2(A,I10))','Checking minima from ',CHECKSPS,' to ',CHECKSPF
547:          INQUIRE(FILE='odata.checksp',EXIST=YESNOT)546:          INQUIRE(FILE='odata.checksp',EXIST=YESNOT)
548:          IF (.NOT.YESNOT) THEN547:          IF (.NOT.YESNOT) THEN
549:             WRITE(*,'(A)') 'main> ERROR - no odata.checksp file present'548:             WRITE(*,'(A)') 'main> ERROR - no odata.checksp file present'
550:             STOP549:             STOP
551:          ENDIF550:          ENDIF
552:          IF (GETMINFRQST) THEN551:          IF (GETMINFRQST) THEN


r29772/Makefile 2016-01-19 18:30:12.876934925 +0000 r29771/Makefile 2016-01-19 18:30:14.092951088 +0000
  7:         mindouble.o tsdouble.o Pfold.o NGT.o GT.o getallpaths.o cycle2.o dprand.o sdprnd.o \  7:         mindouble.o tsdouble.o Pfold.o NGT.o GT.o getallpaths.o cycle2.o dprand.o sdprnd.o \
  8:         Dijkstra.o getdpair.o Dijinit.o calcorder.o getrpair.o connectd.o getspair.o \  8:         Dijkstra.o getdpair.o Dijinit.o calcorder.o getrpair.o connectd.o getspair.o \
  9:         getpair.o pairdouble.o donedouble.o newmindist.o mergedb.o getupair.o \  9:         getpair.o pairdouble.o donedouble.o newmindist.o mergedb.o getupair.o \
 10:         GT2.o GT2data.o GT2FibonacciHeap.o GT2DLL.o GT2input.o GT2FreeMemory.o \ 10:         GT2.o GT2data.o GT2FibonacciHeap.o GT2DLL.o GT2input.o GT2FreeMemory.o \
 11:         minperm.o minpermdist.o rigidbodymod.o mathsconstants.o quaternionmatch.o \ 11:         minperm.o minpermdist.o rigidbodymod.o mathsconstants.o quaternionmatch.o \
 12:         regroupfree.o regroup.o dsort.o savestate.o orderodata.o regroupfree2.o \ 12:         regroupfree.o regroup.o dsort.o savestate.o orderodata.o regroupfree2.o \
 13:         probacc.o newconn.o getfreepair.o getfreebarrier.o kshortestpaths.o Dijinitfly.o \ 13:         probacc.o newconn.o getfreepair.o getfreebarrier.o kshortestpaths.o Dijinitfly.o \
 14:         getallmin.o myorient.o getusepair.o NGTmem.o NGTrealloc.o NGTrenorm.o \ 14:         getallmin.o myorient.o getusepair.o NGTmem.o NGTrealloc.o NGTrenorm.o \
 15:         NGTremoveid.o NGTremovei.o swapnode.o mymerge.o rbinertia.o rigidb.o diagonalise2.o reweight.o \ 15:         NGTremoveid.o NGTremovei.o swapnode.o mymerge.o rbinertia.o rigidb.o diagonalise2.o reweight.o \
 16:         rbperm.o minpermdistrbcom.o Cv.o DOS.o bulkmindist.o frictionfac.o make_conpot.o setup_sis.o \ 16:         rbperm.o minpermdistrbcom.o Cv.o DOS.o bulkmindist.o frictionfac.o make_conpot.o setup_sis.o \
 17:         rateconst_setup.o checkspodata.o addminxyz.o printsummary.o \ 17:         rateconst_setup.o checkspodata.o addminxyz.o \
 18:         checkTS.o remove_unconnected.o getupairmetric.o \ 18:         checkTS.o remove_unconnected.o getupairmetric.o \
 19:         NGT_crstorage.o NGTremovei_crstorage.o NGTremoveid_crstorage.o NGTrenorm_crstorage.o \ 19:         NGT_crstorage.o NGTremovei_crstorage.o NGTremoveid_crstorage.o NGTrenorm_crstorage.o \
 20:         NGTrealloc_crstorage.o lopermdist.o getmetric.o rotations.o vec3.o newconnodata.o Gthomson.o \ 20:         NGTrealloc_crstorage.o lopermdist.o getmetric.o rotations.o vec3.o newconnodata.o Gthomson.o \
 21:         microtherm.o getbarrier2.o persistence.o regrouppersist.o shannon.o get_exact_persistence.o 21:         microtherm.o getbarrier2.o persistence.o regrouppersist.o shannon.o get_exact_persistence.o
 22:  22: 
 23:  23: 
 24: # note that nag64/5.1-216 fails for large memory, but version 365 works 24: # note that nag64/5.1-216 fails for large memory, but version 365 works
 25: # large memory for pgi and ifort requires the -mcmodel flag! 25: # large memory for pgi and ifort requires the -mcmodel flag!
 26:  26: 
 27: # WARNING - points.min and point.ts created with ifort executables cannot be 27: # WARNING - points.min and point.ts created with ifort executables cannot be
235: rateconst_setup.o: commons.o porfuncs.o235: rateconst_setup.o: commons.o porfuncs.o
236: checkTS.o: commons.o porfuncs.o236: checkTS.o: commons.o porfuncs.o
237: lopermdist.o: commons.o 237: lopermdist.o: commons.o 
238: getmetric.o: commons.o porfuncs.o238: getmetric.o: commons.o porfuncs.o
239: myorient.o: commons.o 239: myorient.o: commons.o 
240: microtherm.o: commons.o utils.o240: microtherm.o: commons.o utils.o
241: getbarrier2.o: commons.o241: getbarrier2.o: commons.o
242: persistence.o: commons.o242: persistence.o: commons.o
243: shannon.o: commons.o porfuncs.o savestate.o utils.o243: shannon.o: commons.o porfuncs.o savestate.o utils.o
244: get_exact_persistence.o: commons.o244: get_exact_persistence.o: commons.o
245: printsummary.o:        commons.o porfuncs.o 


r29772/printsummary.f90 2016-01-19 18:30:13.848947845 +0000 r29771/printsummary.f90 2016-01-19 18:30:14.960962625 +0000
  1: !   PATHSAMPLE: A driver for OPTIM to create stationary point databases using discrete path sampling and perform kinetic analysis  1: svn: E195012: Unable to find repository location for 'svn+ssh://svn.ch.private.cam.ac.uk/groups/wales/trunk/PATHSAMPLE/source/printsummary.f90' in revision 29771
  2: !   Copyright (C) 1999-2009 David J. Wales 
  3: !   This file is part of PATHSAMPLE. 
  4: ! 
  5: !   PATHSAMPLE is free software; you can redistribute it and/or modify 
  6: !   it under the terms of the GNU General Public License as published by 
  7: !   the Free Software Foundation; either version 2 of the License, or 
  8: !   (at your option) any later version. 
  9: ! 
 10: !   PATHSAMPLE is distributed in the hope that it will be useful, 
 11: !   but WITHOUT ANY WARRANTY; without even the implied warranty of 
 12: !   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the 
 13: !   GNU General Public License for more details. 
 14: ! 
 15: !   You should have received a copy of the GNU General Public License 
 16: !   along with this program; if not, write to the Free Software 
 17: !   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA 
 18: ! 
 19:  
 20: !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
 21: ! 
 22: ! Printing some summary info about the KTN, to standard out  
 23:  
 24:       SUBROUTINE PRINTSUMMARY 
 25:       USE COMMONS, ONLY : EMIN, ETS, NMIN, NTS, NCONN, PLUS, MINUS, NCONNMIN, TSTHRESH, MAXBARRIER, MINBARRIER, & 
 26:                           KPLUS, KMINUS 
 27:       IMPLICIT NONE 
 28:  
 29:       INTEGER :: NDEAD, MINTARG, NLEFTTS, NCOUNT, NDONE, J1, NLEFTMIN, MAXTSINDEX, MINTSINDEX, CURRTARG 
 30:       INTEGER :: NDISTA(NMIN) 
 31:       DOUBLE PRECISION :: MAXTS, MINTS, GMIN, CURRENERGY 
 32:       DOUBLE PRECISION :: CUT_UNDERFLOW=-300.0D0 ! set here but NOT used in the current coding of checkTS 
 33:       LOGICAL :: CHANGED 
 34:       LOGICAL :: DEADTS(NTS), INCC(NMIN) 
 35:  
 36:       WRITE(*,'(A)')'*************************************************************************************************' 
 37:       WRITE(*,'(A)')'printsummary> Network summary information' 
 38:       WRITE(*,'(A)')'using thresholds:' 
 39:       WRITE(*,'(A,I8)')'NCONNMIN=',NCONNMIN 
 40:       WRITE(*,'(A,G20.10)')'MAXBARRIER=',MAXBARRIER 
 41:       WRITE(*,'(A,G20.10)')'MINBARRIER=',MINBARRIER 
 42:       WRITE(*,'(A,G20.10)')'TSTHRESH or MAXTSENERGY=',TSTHRESH 
 43:       WRITE(*,'(A)')'and NOT including degenerate rearrangements.' 
 44:  
 45:       CALL GETNCONN ! GETNCONN only counts non-degenerate connections 
 46:       NDEAD=0 
 47:       IF (NCONNMIN.GE.0) THEN 
 48:          DO J1=1,NMIN 
 49:             IF (NCONN(J1).LE.NCONNMIN) THEN 
 50:                NDEAD=NDEAD+1 
 51:             ENDIF  
 52:          ENDDO 
 53:       ENDIF 
 54:       WRITE(*,'(A,I8,A,I8)')'There are ',NMIN-NDEAD,' valid minima out of ',NMIN 
 55:  
 56:       GMIN=HUGE(1.0D0) 
 57:       DO J1=1,NMIN 
 58:          IF (EMIN(J1).LT.GMIN .AND. (NCONN(J1).GT.NCONNMIN)) THEN 
 59:             GMIN=EMIN(J1) 
 60:             MINTARG=J1 
 61:          ENDIF 
 62:       ENDDO 
 63:       WRITE(*,'(A,I8,A,G20.10)')'Valid global minimum at index ',MINTARG,' with energy ',GMIN 
 64:  
 65:       NLEFTTS=0 
 66:       DEADTS(1:NTS)=.FALSE. 
 67:       DO J1=1,NTS 
 68:          CALL CHECKTS(ETS(J1),EMIN(PLUS(J1)),EMIN(MINUS(J1)),KPLUS(J1),KMINUS(J1),NCONN(PLUS(J1)),NCONN(MINUS(J1)), & 
 69:                       PLUS(J1),MINUS(J1),.FALSE.,CUT_UNDERFLOW,DEADTS(J1)) 
 70:          IF (.NOT.DEADTS(J1)) NLEFTTS=NLEFTTS+1 
 71:       ENDDO 
 72:       WRITE(*,'(A,I8,A,I8)')'There are ',NLEFTTS,' valid transition states out of ',NTS 
 73:  
 74:       IF ((NMIN-NDEAD.EQ.0).OR.(NLEFTTS.EQ.0)) RETURN ! sanity check 
 75:  
 76:       MINTS=HUGE(1.0D0) 
 77:       MAXTS=-HUGE(1.0D0) 
 78:       DO J1=1,NTS 
 79:          IF (ETS(J1).LE.MINTS .AND. (.NOT.DEADTS(J1))) THEN 
 80:             MINTS=ETS(J1) 
 81:             MINTSINDEX=J1 
 82:          ENDIF 
 83:          IF (ETS(J1).GE.MAXTS .AND. (.NOT.DEADTS(J1))) THEN 
 84:             MAXTS=ETS(J1) 
 85:             MAXTSINDEX=J1 
 86:          ENDIF 
 87:       ENDDO 
 88:       WRITE(*,'(A,I8,A,G20.10)')'Highest valid transition state at index ',MAXTSINDEX,' with energy ',MAXTS 
 89:       WRITE(*,'(A,I8,A,G20.10)')'Lowest valid transition state at index ',MINTSINDEX,' with energy ',MINTS 
 90:  
 91:       WRITE(*,'(A)')'Connected components for valid minima and TSs:' 
 92:       WRITE(*,'(A)')'CC #        # minima    index of lowest min    energy of lowest min' 
 93:       NCOUNT=0 
 94:       NDONE=0 
 95:       CURRTARG=MINTARG ! start with the valid global minimum  
 96:       CURRENERGY=GMIN 
 97:       INCC(1:NMIN)=.FALSE. 
 98:       DO  
 99:          NCOUNT=NCOUNT+1 
100:          NDISTA(1:NMIN)=1000000 
101:          NDISTA(CURRTARG)=0 
102:  
103: 5        CHANGED=.FALSE. 
104:          DO J1=1,NTS 
105:             IF (DEADTS(J1)) CYCLE 
106:             IF (NDISTA(MINUS(J1))+1.LT.NDISTA(PLUS(J1))) THEN 
107:                CHANGED=.TRUE. 
108:                NDISTA(PLUS(J1))=NDISTA(MINUS(J1))+1 
109:             ENDIF 
110:             IF (NDISTA(PLUS(J1))+1.LT.NDISTA(MINUS(J1))) THEN 
111:                CHANGED=.TRUE. 
112:                NDISTA(MINUS(J1))=NDISTA(PLUS(J1))+1 
113:             ENDIF 
114:          ENDDO 
115:          IF (CHANGED) GOTO 5 
116:  
117:          NLEFTMIN=0 
118:          DO J1=1,NMIN 
119:             IF (INCC(J1)) CYCLE 
120:             IF ((NDISTA(J1).LT.1000000).AND.(NCONN(J1).GT.NCONNMIN)) THEN 
121:                NLEFTMIN=NLEFTMIN+1 
122:                NDONE=NDONE+1 
123:                INCC(J1)=.TRUE. 
124:             ENDIF 
125:          ENDDO 
126:  
127:          WRITE(*,'(2(I8,4X),I8,15X,G20.10)')NCOUNT,NLEFTMIN,CURRTARG,CURRENERGY 
128:          IF (NDONE.EQ.NMIN-NDEAD) EXIT 
129:  
130:          ! update CURRTARG and CURRENERGY: find the lowest minimum not yet in a CC 
131:          CURRENERGY=HUGE(1.0D0) 
132:          DO J1=1,NMIN 
133:             IF (INCC(J1)) CYCLE 
134:             IF (EMIN(J1).LT.CURRENERGY .AND. (NCONN(J1).GT.NCONNMIN)) THEN 
135:                CURRENERGY=EMIN(J1) 
136:                CURRTARG=J1 
137:             ENDIF 
138:          ENDDO 
139:       ENDDO 
140:  
141:       WRITE(*,'(A)')'*************************************************************************************************' 
142:  
143:       RETURN 
144:       END SUBROUTINE PRINTSUMMARY 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0