## hdiff output

r33478/connect_unconnected.f90 2017-11-14 18:30:27.151949005 +0000 r33477/connect_unconnected.f90 2017-11-14 18:30:27.407952402 +0000
277:    EREF=EMIN(REFMIN)        277:    EREF=EMIN(REFMIN)
278:    ALLOCATE(EUNDIFF(NUNCONN))278:    ALLOCATE(EUNDIFF(NUNCONN))
279:    DO J1=1,NUNCONN,1279:    DO J1=1,NUNCONN,1
280:       EUNDIFF(J1)=ABS(EUNCON(J1)-EREF)  !use absolute difference280:       EUNDIFF(J1)=ABS(EUNCON(J1)-EREF)  !use absolute difference
281:    ENDDO281:    ENDDO
282:    !CALL SORT(NUNCONN,NUNCONN,EUNDIFF,UNCMIN)  !sort UNCMIN by the energy differences282:    !CALL SORT(NUNCONN,NUNCONN,EUNDIFF,UNCMIN)  !sort UNCMIN by the energy differences
283:    CALL DSORT(EUNDIFF,UNCMIN,NUNCONN,2)283:    CALL DSORT(EUNDIFF,UNCMIN,NUNCONN,2)
284:    !now find the pairs to connect284:    !now find the pairs to connect
285:    J7=0285:    J7=0
286:    DO J4=1,PAIRSTODO286:    DO J4=1,PAIRSTODO
287: 55    IF (J7.LT.NUNCONN) THEN 287: 55    J7=J7+1
288:          J7=J7+1
289:       ELSE
290:          PRINT '(A)', 'connectclosestE> exhausted the list of unconnected minima for this reference'
291:          STOP
292:       ENDIF
293:       DO J3=1,NPAIRDONE,1288:       DO J3=1,NPAIRDONE,1
294:          IF ((PAIR1(J3).EQ.UNCMIN(J7)).AND.(PAIR2(J3).EQ.REFMIN)) THEN289:          IF ((PAIR1(J3).EQ.UNCMIN(J7)).AND.(PAIR2(J3).EQ.REFMIN)) THEN
295:             GOTO 55290:             GOTO 55
296:          ELSE IF ((PAIR1(J3).EQ.REFMIN).AND.(PAIR2(J3).EQ.UNCMIN(J7))) THEN291:          ELSE IF ((PAIR1(J3).EQ.REFMIN).AND.(PAIR2(J3).EQ.UNCMIN(J7))) THEN
297:             GOTO 55292:             GOTO 55
298:          ENDIF293:          ENDIF
299:       ENDDO294:       ENDDO
300:       DMIN1(J4)=UNCMIN(J7)295:       DMIN1(J4)=UNCMIN(J7)
301:       DMIN2(J4)=REFMIN296:       DMIN2(J4)=REFMIN
302:       PRINT '(A,I8,A,F15.8)','connectclosestE> Minimum 1:',DMIN1(J4),' Energy:',EMIN(DMIN1(J4))297:       PRINT '(A,I8,A,F15.8)','connectclosestE> Minimum 1:',DMIN1(J4),' Energy:',EMIN(DMIN1(J4))
313:       CALL MINPERMDIST(SETX,REFX,NATOMS,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT,TWOD,DISTANCE, &308:       CALL MINPERMDIST(SETX,REFX,NATOMS,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT,TWOD,DISTANCE, &
314: &                         DIST2,RIGIDBODY,RMAT,.FALSE.)309: &                         DIST2,RIGIDBODY,RMAT,.FALSE.)
315:       UNCONDIST(J5) = DISTANCE                           310:       UNCONDIST(J5) = DISTANCE
316:    ENDDO311:    ENDDO
317:   !CALL SORT(NUNCONN,NUNCONN,UNCONDIST,UNCMIN)312:   !CALL SORT(NUNCONN,NUNCONN,UNCONDIST,UNCMIN)
318:    CALL DSORT(UNCONDIST,UNCMIN,NUNCONN,2)313:    CALL DSORT(UNCONDIST,UNCMIN,NUNCONN,2)
319:    !now find the pairs to connect314:    !now find the pairs to connect
320:    J7=0315:    J7=0
321:    DO J4=1,PAIRSTODO316:    DO J4=1,PAIRSTODO
322: 56    IF (J7.LT.NUNCONN) THEN317: 56    J7=J7+1
323:          J7=J7+1
324:       ELSE
325:          PRINT '(A)', 'connectclosestE> exhausted the list of unconnected minima for this reference'
326:          STOP
327:       ENDIF
328:       DO J3=1,NPAIRDONE,1318:       DO J3=1,NPAIRDONE,1
329:          IF ((PAIR1(J3).EQ.UNCMIN(J7)).AND.(PAIR2(J3).EQ.REFMIN)) THEN319:          IF ((PAIR1(J3).EQ.UNCMIN(J7)).AND.(PAIR2(J3).EQ.REFMIN)) THEN
330:             GOTO 56320:             GOTO 56
331:          ELSE IF ((PAIR1(J3).EQ.REFMIN).AND.(PAIR2(J3).EQ.UNCMIN(J7))) THEN321:          ELSE IF ((PAIR1(J3).EQ.REFMIN).AND.(PAIR2(J3).EQ.UNCMIN(J7))) THEN
332:             GOTO 56322:             GOTO 56
333:          ENDIF323:          ENDIF
334:       ENDDO324:       ENDDO
335:       DMIN1(J4)=UNCMIN(J7)325:       DMIN1(J4)=UNCMIN(J7)
336:       DMIN2(J4)=REFMIN326:       DMIN2(J4)=REFMIN
337:       PRINT '(A,I8,A,F15.8)','connectclosestD> Minimum 1:',DMIN1(J4),' Energy:',EMIN(DMIN1(J4))327:       PRINT '(A,I8,A,F15.8)','connectclosestD> Minimum 1:',DMIN1(J4),' Energy:',EMIN(DMIN1(J4))

legend