hdiff output

r33460/commons.f90 2017-11-10 15:30:16.834118657 +0000 r33459/commons.f90 2017-11-10 15:30:17.446126712 +0000
658:       INTEGER, ALLOCATABLE ::  LJADDNN(:,:)658:       INTEGER, ALLOCATABLE ::  LJADDNN(:,:)
659: 659: 
660:       INTEGER, DIMENSION(:,:), ALLOCATABLE :: BONDS !for QCIAMBER660:       INTEGER, DIMENSION(:,:), ALLOCATABLE :: BONDS !for QCIAMBER
661: 661: 
662: !OPEP interface662: !OPEP interface
663:       LOGICAL :: OPEPT, OPEP_RNAT663:       LOGICAL :: OPEPT, OPEP_RNAT
664: 664: 
665: !AMBER mutational steps665: !AMBER mutational steps
666:       LOGICAL :: AMBERMUTATIONT, AMBERMUTRIGIDT666:       LOGICAL :: AMBERMUTATIONT, AMBERMUTRIGIDT
667:       INTEGER :: MUTUNIT,NMUTATION,MUTATIONFREQ,MUTTESTSTEPS,AMBERMUTFF,AMBERMUTIGB,MUTENERGY,MUTTERMID667:       INTEGER :: MUTUNIT,NMUTATION,MUTATIONFREQ,MUTTESTSTEPS,AMBERMUTFF,AMBERMUTIGB,MUTENERGY,MUTTERMID
668:       DOUBLE PRECISION :: BESTMUTSCORE 
669:       INTEGER :: BESTMUTN, NSEQSTORED 
670:       CHARACTER(LEN=4), ALLOCATABLE :: BESTMUTSEQ(:) , MUTSEQ(:,:) 
671:       LOGICAL, ALLOCATABLE :: SEQREJECTED(:) 
672: 668: 
673: !Orbital variables669: !Orbital variables
674:       LOGICAL :: ORBITALS670:       LOGICAL :: ORBITALS
675:       INTEGER :: NROTS, NORBS, ORBVAREXPONENT671:       INTEGER :: NROTS, NORBS, ORBVAREXPONENT
676:       DOUBLE PRECISION, ALLOCATABLE :: R2INTS(:,:), DIPINTS(:,:,:)672:       DOUBLE PRECISION, ALLOCATABLE :: R2INTS(:,:), DIPINTS(:,:,:)
677: END MODULE COMMONS673: END MODULE COMMONS


r33460/mc.F 2017-11-10 15:30:17.174123121 +0000 r33459/mc.F 2017-11-10 15:30:17.714130254 +0000
537:          ENDIF537:          ENDIF
538: ! hk286538: ! hk286
539:       ENDDO539:       ENDDO
540:       EPSSPHERE=EPSSAVE540:       EPSSPHERE=EPSSAVE
541: 541: 
542: ! kr366 - save score for initial quench 542: ! kr366 - save score for initial quench 
543:       IF (AMBERMUTATIONT) THEN543:       IF (AMBERMUTATIONT) THEN
544:          CALL MUTATION_E(MUTSCORE,COORDS(1:3*NATOMS,1),MUTENERGY,MUTTERMID)544:          CALL MUTATION_E(MUTSCORE,COORDS(1:3*NATOMS,1),MUTENERGY,MUTTERMID)
545:          WRITE(MUTSCOREUNIT, '(A,F20.10)') '  0  ',MUTSCORE545:          WRITE(MUTSCOREUNIT, '(A,F20.10)') '  0  ',MUTSCORE
546:          CALL FLUSH(MUTSCOREUNIT)546:          CALL FLUSH(MUTSCOREUNIT)
547:          BESTMUTSEQ(:) = MUTSEQ(NSEQSTORED, :) 
548:          BESTMUTN = NSEQSTORED 
549:          BESTMUTSCORE = MUTSCORE 
550:       ENDIF547:       ENDIF
551: 548: 
552: #endif549: #endif
553: 550: 
554:       IF (THOMSONT .AND. (.NOT. GTHOMSONT)) THEN551:       IF (THOMSONT .AND. (.NOT. GTHOMSONT)) THEN
555: !552: !
556: ! Scale maximum step size for the Thomson problem according to the mean nearest-neighbour 553: ! Scale maximum step size for the Thomson problem according to the mean nearest-neighbour 
557: ! distance after the first quench.554: ! distance after the first quench.
558: !555: !
559:          DUMMY1=0.0D0556:          DUMMY1=0.0D0
1783:                 WRITE(MYRUNIT,'(I10,F15.5)') NQ(JP),RMSD1780:                 WRITE(MYRUNIT,'(I10,F15.5)') NQ(JP),RMSD
1784:                 CALL FLUSH(MYRUNIT)1781:                 CALL FLUSH(MYRUNIT)
1785:              ENDIF1782:              ENDIF
1786:           ENDIF1783:           ENDIF
1787: 1784: 
1788: 1785: 
1789:           IF (AMBERMUTATIONT) THEN1786:           IF (AMBERMUTATIONT) THEN
1790:              CALL MUTATION_E(MUTSCORE,COORDS(:,1),MUTENERGY,MUTTERMID)1787:              CALL MUTATION_E(MUTSCORE,COORDS(:,1),MUTENERGY,MUTTERMID)
1791:              WRITE(MUTSCOREUNIT, '(I8,F20.10)') J1,MUTSCORE1788:              WRITE(MUTSCOREUNIT, '(I8,F20.10)') J1,MUTSCORE
1792:              CALL FLUSH(MUTSCOREUNIT)1789:              CALL FLUSH(MUTSCOREUNIT)
1793:              IF (MUTSCORE.LT.BESTMUTSCORE) THEN 
1794:                 IF (NSEQSTORED.EQ.BESTMUTN) THEN 
1795:                    BESTMUTSCORE = MUTSCORE 
1796:                 ELSE 
1797:                    BESTMUTSCORE = MUTSCORE 
1798:                    BESTMUTSEQ(:) = MUTSEQ(NSEQSTORED , :) 
1799:                    BESTMUTN = NSEQSTORED 
1800:                 ENDIF 
1801:              ENDIF 
1802:           ENDIF1790:           ENDIF
1803: 1791: 
1804: 1792: 
1805: ! DAESTAT keyword prints two files :1793: ! DAESTAT keyword prints two files :
1806: ! stat.all which contains all quenches and their energies1794: ! stat.all which contains all quenches and their energies
1807: ! stat.acc which contains only accepted quenches and their energies1795: ! stat.acc which contains only accepted quenches and their energies
1808: ! used to analyse how many new minima are being found1796: ! used to analyse how many new minima are being found
1809:  1797:  
1810:             IF (DAESTAT) THEN1798:             IF (DAESTAT) THEN
1811:                PRINT*,'DAESTAT block in mc.f not implemented'1799:                PRINT*,'DAESTAT block in mc.f not implemented'


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0