hdiff output

r31925/getallmin.f 2017-02-16 10:30:17.423543997 +0000 r31924/getallmin.f 2017-02-16 10:30:17.871549975 +0000
 22: C 22: C
 23: C  This subroutine analyses a min.data.info file in the format generated  23: C  This subroutine analyses a min.data.info file in the format generated 
 24: C  the OPTIM keyword BHINTERP or BISECTT. 24: C  the OPTIM keyword BHINTERP or BISECTT.
 25: C 25: C
 26:       SUBROUTINE GETALLMIN(MINS,MINF) 26:       SUBROUTINE GETALLMIN(MINS,MINF)
 27:       USE PORFUNCS 27:       USE PORFUNCS
 28:       USE COMMONS 28:       USE COMMONS
 29:       USE UTILS,ONLY : GETUNIT 29:       USE UTILS,ONLY : GETUNIT
 30:       IMPLICIT NONE 30:       IMPLICIT NONE
 31:  31: 
 32:       INTEGER J1, J2, J3, MIN1, MIN2, J4, LUNIT 32:       INTEGER J1, J2, ISTAT, J3, MIN1, MIN2, J4, LUNIT
 33:       DOUBLE PRECISION LOCALPOINTS(3*NATOMS), NEWEMIN, DISTANCE, RMAT(3,3), LOCALPOINTS2(3*NATOMS) 33:       DOUBLE PRECISION LOCALPOINTS(3*NATOMS), NEWEMIN, DISTANCE, RMAT(3,3), LOCALPOINTS2(3*NATOMS)
 34:       DOUBLE PRECISION DIST2  34:       DOUBLE PRECISION DIST2 
 35:       DOUBLE PRECISION NEWFVIBMIN, NEWPOINTSMIN(3*NATOMS), NEWIXMIN,  NEWIYMIN, NEWIZMIN, TSEGUESS, NEWMINCURVE, NEWMINFRQ2, 35:       DOUBLE PRECISION NEWFVIBMIN, NEWPOINTSMIN(3*NATOMS), NEWIXMIN,  NEWIYMIN, NEWIZMIN, TSEGUESS, NEWMINCURVE, NEWMINFRQ2,
 36:      &                 TSFVIBGUESS, ETSDUMMY, DUMMY, FRICTIONFAC 36:      &                 TSFVIBGUESS, ETSDUMMY, DUMMY, FRICTIONFAC
 37:       INTEGER NEWHORDERMIN, NEWMIN, MINS, MINF, NMINSAVE 37:       INTEGER NEWHORDERMIN, NEWMIN, MINS, MINF, NMINSAVE
 38:       LOGICAL MINISOLD, MATCHED, REWRITE 38:       LOGICAL MINISOLD, MATCHED, REWRITE
 39:  39: 
 40:       LUNIT=GETUNIT() 40:       LUNIT=GETUNIT()
 41:       IF (MACHINE) THEN 41:       IF (MACHINE) THEN
 42:          OPEN(LUNIT,FILE='min.data.info',STATUS='OLD',form='unformatted') 42:          OPEN(LUNIT,FILE='min.data.info',STATUS='OLD',form='unformatted')
246: 130      CONTINUE246: 130      CONTINUE
247:       ENDDO247:       ENDDO
248: 110   CLOSE(LUNIT)248: 110   CLOSE(LUNIT)
249:       IF (DIJINITT) THEN249:       IF (DIJINITT) THEN
250:          IF (NMIN.GT.NMINSAVE) THEN ! write new pairdist and pairlist files250:          IF (NMIN.GT.NMINSAVE) THEN ! write new pairdist and pairlist files
251:             LUNIT=GETUNIT()251:             LUNIT=GETUNIT()
252:             OPEN(UNIT=LUNIT,FILE='pairdist',STATUS='UNKNOWN',POSITION='APPEND')252:             OPEN(UNIT=LUNIT,FILE='pairdist',STATUS='UNKNOWN',POSITION='APPEND')
253:             DO J3=NMINSAVE+1,NMIN253:             DO J3=NMINSAVE+1,NMIN
254:                WRITE(LUNIT,'(10G20.10)') (PAIRDIST(J3,J4),J4=1,PAIRDISTMAX)254:                WRITE(LUNIT,'(10G20.10)') (PAIRDIST(J3,J4),J4=1,PAIRDISTMAX)
255:             ENDDO255:             ENDDO
256:             CALL FLUSH(LUNIT)256:             CALL FLUSH(LUNIT,ISTAT)
257:             CLOSE(LUNIT)257:             CLOSE(LUNIT)
258:             OPEN(UNIT=LUNIT,FILE='pairlist',STATUS='UNKNOWN',POSITION='APPEND')258:             OPEN(UNIT=LUNIT,FILE='pairlist',STATUS='UNKNOWN',POSITION='APPEND')
259:             DO J3=NMINSAVE+1,NMIN259:             DO J3=NMINSAVE+1,NMIN
260:                WRITE(LUNIT,'(10I10)') (PAIRLIST(J3,J4),J4=1,PAIRDISTMAX)260:                WRITE(LUNIT,'(10I10)') (PAIRLIST(J3,J4),J4=1,PAIRDISTMAX)
261:             ENDDO261:             ENDDO
262:             CALL FLUSH(LUNIT)262:             CALL FLUSH(LUNIT,ISTAT)
263:             CLOSE(LUNIT)263:             CLOSE(LUNIT)
264:          ENDIF264:          ENDIF
265:       ENDIF265:       ENDIF
266: !266: !
267: !  If we found new minima between the original ones then adjust the barrier for MINS and MINF267: !  If we found new minima between the original ones then adjust the barrier for MINS and MINF
268: !  according to the sum of distances.268: !  according to the sum of distances.
269: !269: !
270:       IF ((NMIN-NMINSAVE.GT.0).AND.(BISECTT)) THEN270:       IF ((NMIN-NMINSAVE.GT.0).AND.(BISECTT)) THEN
271:          MATCHED=.FALSE.271:          MATCHED=.FALSE.
272:          DO J3=1,NTS272:          DO J3=1,NTS


r31925/getallpaths.f 2017-02-16 10:30:17.647546987 +0000 r31924/getallpaths.f 2017-02-16 10:30:18.091552910 +0000
 22: C 22: C
 23: C  This subroutine analyses a path.info in the new min-sad-min min-sad-min format 23: C  This subroutine analyses a path.info in the new min-sad-min min-sad-min format
 24: C  as generated with OPTIM keyword DUMPALLPATHS. 24: C  as generated with OPTIM keyword DUMPALLPATHS.
 25: C 25: C
 26:       SUBROUTINE GETALLPATHS 26:       SUBROUTINE GETALLPATHS
 27:       USE PORFUNCS 27:       USE PORFUNCS
 28:       USE COMMONS 28:       USE COMMONS
 29:       USE UTILS,ONLY : GETUNIT 29:       USE UTILS,ONLY : GETUNIT
 30:       IMPLICIT NONE 30:       IMPLICIT NONE
 31:  31: 
 32:       INTEGER J1, J2, NMINOLD, TSNUMBER, J3, NCOUNT, NMINSAVE, NTSSAVE, J4, J5, LUNIT, PLUSMIN 32:       INTEGER J1, J2, ISTAT, NMINOLD, TSNUMBER, J3, NCOUNT, NMINSAVE, NTSSAVE, J4, J5, LUNIT, PLUSMIN
 33:       DOUBLE PRECISION LOCALPOINTS(NOPT), ENERGY, NEWEMIN, NEWETS, DISTANCE, RMAT(3,3), 33:       DOUBLE PRECISION LOCALPOINTS(NOPT), ENERGY, NEWEMIN, NEWETS, DISTANCE, RMAT(3,3),
 34:      1                 LPOINTSTS(NOPT), LPLUS(NOPT), LMINUS(NOPT), LOCALPOINTS2(NOPT) 34:      1                 LPOINTSTS(NOPT), LPLUS(NOPT), LMINUS(NOPT), LOCALPOINTS2(NOPT)
 35:       DOUBLE PRECISION DUMMY, DIST2, ELAPSED, TNEW 35:       DOUBLE PRECISION DUMMY, DIST2, ELAPSED, TNEW
 36:       DOUBLE PRECISION NEWFVIBMIN, NEWFVIBTS, NEWNEGEIG, NEWPOINTSMIN(NOPT), NEWPOINTSMINPLUS(NOPT), EPLUS, 36:       DOUBLE PRECISION NEWFVIBMIN, NEWFVIBTS, NEWNEGEIG, NEWPOINTSMIN(NOPT), NEWPOINTSMINPLUS(NOPT), EPLUS,
 37:      1                 NEWPOINTSTS(NOPT), NEWIXMIN,  NEWIYMIN, NEWIZMIN, IXPLUS, IYPLUS, IZPLUS, 37:      1                 NEWPOINTSTS(NOPT), NEWIXMIN,  NEWIYMIN, NEWIZMIN, IXPLUS, IYPLUS, IZPLUS,
 38:      2                 NEWIXTS,  NEWIYTS, NEWIZTS, IXMINUS, IYMINUS, IZMINUS, FRICTIONFAC, TEMPD(PAIRDISTMAX) 38:      2                 NEWIXTS,  NEWIYTS, NEWIZTS, IXMINUS, IYMINUS, IZMINUS, FRICTIONFAC, TEMPD(PAIRDISTMAX)
 39:       INTEGER NEWHORDERMIN, NEWHORDERTS, NEWMIN, NEWTS, NTRIPLES, TEMPL(PAIRDISTMAX), NFRQS, NVARS 39:       INTEGER NEWHORDERMIN, NEWHORDERTS, NEWMIN, NEWTS, NTRIPLES, TEMPL(PAIRDISTMAX), NFRQS, NVARS
 40:       LOGICAL TSISOLD, FAILED, MINPOLD, MINMOLD, BADTRIPLE, TESTOP, TESTNAME 40:       LOGICAL TSISOLD, FAILED, MINPOLD, MINMOLD, BADTRIPLE, TESTOP, TESTNAME
 41:       CHARACTER(LEN=1) DUMMYSTRING 41:       CHARACTER(LEN=1) DUMMYSTRING
 42:  42: 
800:                  ENDIF800:                  ENDIF
801: !801: !
802: ! Since entries have changed we'd better rewrite these files rather802: ! Since entries have changed we'd better rewrite these files rather
803: ! than just append to them. 803: ! than just append to them. 
804: !804: !
805:                   LUNIT=GETUNIT()805:                   LUNIT=GETUNIT()
806:                   OPEN(UNIT=LUNIT,FILE='pairdist',STATUS='UNKNOWN')806:                   OPEN(UNIT=LUNIT,FILE='pairdist',STATUS='UNKNOWN')
807:                   DO J3=1,NMIN807:                   DO J3=1,NMIN
808:                      WRITE(LUNIT,'(10G20.10)') (PAIRDIST(J3,J4),J4=1,PAIRDISTMAX)808:                      WRITE(LUNIT,'(10G20.10)') (PAIRDIST(J3,J4),J4=1,PAIRDISTMAX)
809:                   ENDDO809:                   ENDDO
810:                   CALL FLUSH(LUNIT)810:                   CALL FLUSH(LUNIT,ISTAT)
811:                   CLOSE(LUNIT)811:                   CLOSE(LUNIT)
812:                   OPEN(UNIT=LUNIT,FILE='pairlist',STATUS='UNKNOWN')812:                   OPEN(UNIT=LUNIT,FILE='pairlist',STATUS='UNKNOWN')
813:                   DO J3=1,NMIN813:                   DO J3=1,NMIN
814:                      WRITE(LUNIT,'(10I10)') (PAIRLIST(J3,J4),J4=1,PAIRDISTMAX)814:                      WRITE(LUNIT,'(10I10)') (PAIRLIST(J3,J4),J4=1,PAIRDISTMAX)
815:                   ENDDO815:                   ENDDO
816:                   CALL FLUSH(LUNIT)816:                   CALL FLUSH(LUNIT,ISTAT)
817:                   CLOSE(LUNIT)817:                   CLOSE(LUNIT)
818:                ENDIF818:                ENDIF
819:             ENDIF819:             ENDIF
820: C820: C
821: C  Update ts pointers.821: C  Update ts pointers.
822: C822: C
823:             POINTERP(NTS)=-1823:             POINTERP(NTS)=-1
824:             POINTERM(NTS)=-1824:             POINTERM(NTS)=-1
825:             IF (TOPPOINTER(PLUS(NTS)).GT.0) POINTERP(NTS)=TOPPOINTER(PLUS(NTS))825:             IF (TOPPOINTER(PLUS(NTS)).GT.0) POINTERP(NTS)=TOPPOINTER(PLUS(NTS))
826:             IF (TOPPOINTER(MINUS(NTS)).GT.0) POINTERM(NTS)=TOPPOINTER(MINUS(NTS))826:             IF (TOPPOINTER(MINUS(NTS)).GT.0) POINTERM(NTS)=TOPPOINTER(MINUS(NTS))


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0