hdiff output

r32791/keywords.f 2017-06-14 12:30:14.735644311 +0100 r32790/keywords.f 2017-06-14 12:30:15.287651563 +0100
904:          DIAGSCALE=1.0D0904:          DIAGSCALE=1.0D0
905:          IF (NITEMS.GT.1) CALL READI(NCONNMIN)905:          IF (NITEMS.GT.1) CALL READI(NCONNMIN)
906:          IF (NITEMS.GT.2) CALL READI(NDIAGEIG)906:          IF (NITEMS.GT.2) CALL READI(NDIAGEIG)
907:          IF (NITEMS.GT.3) CALL READF(DIAGSCALE)907:          IF (NITEMS.GT.3) CALL READF(DIAGSCALE)
908: C908: C
909: C  DIJINIT specifies a Dijkstra analysis to try and construct an initial path.909: C  DIJINIT specifies a Dijkstra analysis to try and construct an initial path.
910: C  We only need files start and finish containing the end point coordinates for910: C  We only need files start and finish containing the end point coordinates for
911: C  DIJINITSTART - this keyword creates new min.A, min.B, points.min, points.ts,911: C  DIJINITSTART - this keyword creates new min.A, min.B, points.min, points.ts,
912: C  min.data and ts.data files, so it is potentially dangerous!912: C  min.data and ts.data files, so it is potentially dangerous!
913: C  DIJINITCONT specifies an initial path run where the above files already exist.913: C  DIJINITCONT specifies an initial path run where the above files already exist.
914: C  In a DIJINITSTART run 914: C  In a DIJINITSTART run C  file min.A is created exist with entries 1 1 and file min.B
915: C  file min.A is created with entries 1 1 and file min.B 
916: C  with entries 1 2. min.data is set up with the two entries for915: C  with entries 1 2. min.data is set up with the two entries for
917: C  these minima. The coordinates are obtained via odata.start and odata.finish and put in records916: C  these minima. The coordinates are obtained via odata.start and odata.finish and put in records
918: C  1 and 2 in the points.min file. 917: C  1 and 2 in the points.min file. 
919: C  PAIRDISTMAX is the number of nearest neighbours to save in PAIRLIST judged by918: C  PAIRDISTMAX is the number of nearest neighbours to save in PAIRLIST judged by
920: C  the chosen distance or interpolation metric.919: C  the chosen distance or interpolation metric.
921: C920: C
922:       ELSE IF ((WORD.EQ.'DIJINITCONT').OR.(WORD.EQ.'DIJPRUNE')) THEN921:       ELSE IF ((WORD.EQ.'DIJINITCONT').OR.(WORD.EQ.'DIJPRUNE')) THEN
923:          IF (WORD.EQ.'DIJPRUNE') DIJPRUNET=.TRUE.922:          IF (WORD.EQ.'DIJPRUNE') DIJPRUNET=.TRUE.
924:          IF (NITEMS.GT.1) THEN923:          IF (NITEMS.GT.1) THEN
925:                CALL READU(WW)924:                CALL READU(WW)


r32791/KMCcommit.f 2017-06-14 12:30:14.347639167 +0100 r32790/KMCcommit.f 2017-06-14 12:30:15.007647873 +0100
776: C  Check row normalisation.776: C  Check row normalisation.
777: C777: C
778:          DO J1=1,NMIN778:          DO J1=1,NMIN
779:             DUMMY=0.0D0779:             DUMMY=0.0D0
780:             DO J2=1,NCOL(J1)780:             DO J2=1,NCOL(J1)
781:                DUMMY=DUMMY+DMAT0(J2,J1)781:                DUMMY=DUMMY+DMAT0(J2,J1)
782:                IF (DEBUG) WRITE(*,'(A,3I6,3G20.10)') 'J1,J2,NVAL,DMAT0,sum=',J1,J2,NVAL(J2,J1),DMAT0(J2,J1),DUMMY782:                IF (DEBUG) WRITE(*,'(A,3I6,3G20.10)') 'J1,J2,NVAL,DMAT0,sum=',J1,J2,NVAL(J2,J1),DMAT0(J2,J1),DUMMY
783:             ENDDO783:             ENDDO
784:             IF (DEBUG.AND.(NCOL(J1).GT.0)) WRITE(*,'(A,2I6,3G20.10)') 'J1,ncol,sum=',J1,NCOL(J1),DUMMY784:             IF (DEBUG.AND.(NCOL(J1).GT.0)) WRITE(*,'(A,2I6,3G20.10)') 'J1,ncol,sum=',J1,NCOL(J1),DUMMY
785:             IF ((NCOL(J1).GT.0).AND.(ABS(DUMMY-1.0D0).GT.1.0D-10)) THEN785:             IF ((NCOL(J1).GT.0).AND.(ABS(DUMMY-1.0D0).GT.1.0D-10)) THEN
786:             PRINT*,'MAKED> ERROR - J1,NCOL(J1),DUMMY=',J1,NCOL(J1),DUMMY786:             PRINT*,'KMCcommit> ERROR - J1,NCOL(J1),DUMMY=',J1,NCOL(J1),DUMMY
787:                   STOP787:                   STOP
788:             ENDIF788:             ENDIF
789:          ENDDO789:          ENDDO
790:       ELSE790:       ELSE
791: C791: C
792: C  Reorder DMAT0 so that the connected minima for minimum J1792: C  Reorder DMAT0 so that the connected minima for minimum J1
793: C  appear in ascending order. Could speed up PFOLD calculation in 793: C  appear in ascending order. Could speed up PFOLD calculation in 
794: C  new compressed row storage format. Doesn't seem to make much difference for LJ55.794: C  new compressed row storage format. Doesn't seem to make much difference for LJ55.
795: C795: C
796:          DO J1=1,NMIN796:          DO J1=1,NMIN
908: C  Check row normalisation.908: C  Check row normalisation.
909: C909: C
910:       DO J1=1,NMIN910:       DO J1=1,NMIN
911:          DUMMY=0.0D0911:          DUMMY=0.0D0
912:          DO J2=1,NCOL(J1)912:          DO J2=1,NCOL(J1)
913:             DUMMY=DUMMY+DMAT0(NOFF(J1)+J2)913:             DUMMY=DUMMY+DMAT0(NOFF(J1)+J2)
914:             IF (DEBUG) WRITE(*,'(A,3I8,3G20.10)') 'J1,J2,NVAL,DMAT0,sum=',J1,J2,NVAL(NOFF(J1)+J2),DMAT0(NOFF(J1)+J2),DUMMY914:             IF (DEBUG) WRITE(*,'(A,3I8,3G20.10)') 'J1,J2,NVAL,DMAT0,sum=',J1,J2,NVAL(NOFF(J1)+J2),DMAT0(NOFF(J1)+J2),DUMMY
915:          ENDDO915:          ENDDO
916:          IF (DEBUG.AND.(NCOL(J1).GT.0)) WRITE(*,'(A,2I8,3G20.10)') 'J1,ncol,sum=',J1,NCOL(J1),DUMMY916:          IF (DEBUG.AND.(NCOL(J1).GT.0)) WRITE(*,'(A,2I8,3G20.10)') 'J1,ncol,sum=',J1,NCOL(J1),DUMMY
917:          IF ((NCOL(J1).GT.0).AND.(ABS(DUMMY-1.0D0).GT.1.0D-10)) THEN917:          IF ((NCOL(J1).GT.0).AND.(ABS(DUMMY-1.0D0).GT.1.0D-10)) THEN
918:          PRINT*,'MAKEDLIN> ERROR - J1,NCOL(J1),DUMMY=',J1,NCOL(J1),DUMMY918:          PRINT*,'KMCcommit> ERROR - J1,NCOL(J1),DUMMY=',J1,NCOL(J1),DUMMY
919:                STOP919:                STOP
920:          ENDIF920:          ENDIF
921:       ENDDO921:       ENDDO
922: 922: 
923:       RETURN923:       RETURN
924:       END924:       END


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0