hdiff output

r29905/keywords.f 2016-02-05 15:30:06.862948795 +0000 r29904/keywords.f 2016-02-05 15:30:07.058951464 +0000
133:          DOUBLE PRECISION LJREPBB, LJATTBB, LJREPLL, LJATTLL, LJREPNN, LJATTNN,133:          DOUBLE PRECISION LJREPBB, LJATTBB, LJREPLL, LJATTLL, LJREPNN, LJATTNN,
134:      &   HABLN, HBBLN, HCBLN, HDBLN, EABLN, EBBLN, ECBLN, EDBLN, TABLN, TBBLN, TCBLN, TDBLN134:      &   HABLN, HBBLN, HCBLN, HDBLN, EABLN, EBBLN, ECBLN, EDBLN, TABLN, TBBLN, TCBLN, TDBLN
135:          INTEGER IGBNAB     ! sf344135:          INTEGER IGBNAB     ! sf344
136:          ! LOCAL AMH VARIABLES136:          ! LOCAL AMH VARIABLES
137:          INTEGER NRES_AMH, I_RES, GLY_COUNT137:          INTEGER NRES_AMH, I_RES, GLY_COUNT
138:          ! CHARACTER(LEN=5) TARFL138:          ! CHARACTER(LEN=5) TARFL
139:          ! DOUBLE PRECISION X, Y, Z139:          ! DOUBLE PRECISION X, Y, Z
140:          INTEGER :: GROUPCENTRE140:          INTEGER :: GROUPCENTRE
141:          DOUBLE PRECISION :: GROUPRADIUS,DISTGROUPX2,DISTGROUPY2,DISTGROUPZ2,DISTGROUPCENTRE141:          DOUBLE PRECISION :: GROUPRADIUS,DISTGROUPX2,DISTGROUPY2,DISTGROUPZ2,DISTGROUPCENTRE
142:          CHARACTER (LEN=2) :: FREEZEGROUPTYPE142:          CHARACTER (LEN=2) :: FREEZEGROUPTYPE
143:          LOGICAL :: FREEZEGROUPT, TURNOFFCHECKCHIRALITY, MLPDONE, MLPNORM143:          LOGICAL :: FREEZEGROUPT, TURNOFFCHECKCHIRALITY, MLPDONE
144:          LOGICAL :: RES_IN_LIST144:          LOGICAL :: RES_IN_LIST
145:          DOUBLE PRECISION LPI145:          DOUBLE PRECISION LPI
146:          INTEGER DATA_UNIT146:          INTEGER DATA_UNIT
147:          CHARACTER(LEN=13) :: AAOPTION147:          CHARACTER(LEN=13) :: AAOPTION
148:          CHARACTER(LEN=20) :: AMBERSTR148:          CHARACTER(LEN=20) :: AMBERSTR
149: 149: 
150:          LPI=3.14159265358979323846264338327950288419716939937510D0150:          LPI=3.14159265358979323846264338327950288419716939937510D0
151:          AAA=0151:          AAA=0
152:          AAB=0152:          AAB=0
153:          ABB=0153:          ABB=0
891:          CUDAT=.FALSE.891:          CUDAT=.FALSE.
892:          CUDAPOT=' '892:          CUDAPOT=' '
893:          CUDATIMET=.FALSE.893:          CUDATIMET=.FALSE.
894: 894: 
895: !895: !
896: ! Neural network potential896: ! Neural network potential
897: !897: !
898:          MLP3T=.FALSE.898:          MLP3T=.FALSE.
899:          MLPB3T=.FALSE.899:          MLPB3T=.FALSE.
900:          MLPDONE=.FALSE.900:          MLPDONE=.FALSE.
901:          MLPNORM=.FALSE. 
902:          MLPLAMBDA=0.0D0901:          MLPLAMBDA=0.0D0
903: 902: 
904:          CLSTRINGT=.FALSE.903:          CLSTRINGT=.FALSE.
905:          CLSTRINGTST=.FALSE.904:          CLSTRINGTST=.FALSE.
906:          IF (FILTH2.EQ.0) THEN905:          IF (FILTH2.EQ.0) THEN
907:             OPEN (5,FILE='odata',STATUS='OLD')906:             OPEN (5,FILE='odata',STATUS='OLD')
908:          ELSE907:          ELSE
909:             WRITE(OTEMP,*) FILTH2908:             WRITE(OTEMP,*) FILTH2
910:             WRITE(OSTRING,'(A)') 'odata.' // TRIM(ADJUSTL(OTEMP))909:             WRITE(OSTRING,'(A)') 'odata.' // TRIM(ADJUSTL(OTEMP))
911:             OPEN (5,FILE=OSTRING,STATUS='OLD')910:             OPEN (5,FILE=OSTRING,STATUS='OLD')
3731:             NMLP=MLPHIDDEN*(MLPIN+MLPOUT)+13730:             NMLP=MLPHIDDEN*(MLPIN+MLPOUT)+1
3732:          ELSE3731:          ELSE
3733:             WRITE(*,'(A,4I8,G20.10)') 'MLP3 potential with Nin, Nhidden, Nout, Ndata, lambda=',3732:             WRITE(*,'(A,4I8,G20.10)') 'MLP3 potential with Nin, Nhidden, Nout, Ndata, lambda=',
3734:      &                                   MLPIN,MLPHIDDEN,MLPOUT,MLPDATA,MLPLAMBDA3733:      &                                   MLPIN,MLPHIDDEN,MLPOUT,MLPDATA,MLPLAMBDA
3735:             NMLP=MLPHIDDEN*(MLPIN+MLPOUT)3734:             NMLP=MLPHIDDEN*(MLPIN+MLPOUT)
3736:          ENDIF3735:          ENDIF
3737:          IF (NMLP.NE.NATOMS) THEN3736:          IF (NMLP.NE.NATOMS) THEN
3738:             PRINT '(A,2I8)', 'keywords> ERROR *** NATOMS,NMLP=',NATOMS,NMLP3737:             PRINT '(A,2I8)', 'keywords> ERROR *** NATOMS,NMLP=',NATOMS,NMLP
3739:             STOP3738:             STOP
3740:          ENDIF3739:          ENDIF
 3740:          IF (.NOT.MLPDONE) THEN
 3741:             LUNIT=GETUNIT()
 3742:             OPEN(LUNIT,FILE='MLPdata',STATUS='OLD')
 3743:             ALLOCATE(MLPDAT(MLPDATA,MLPIN),MLPOUTCOME(MLPDATA))
 3744:             DO J1=1,MLPDATA
 3745:                READ(LUNIT,*) MLPDAT(J1,1:MLPIN),MLPOUTCOME(J1)
 3746:                MLPOUTCOME(J1)=MLPOUTCOME(J1)+1 ! to shift the range to 1 to 4
 3747: !           WRITE(*,'(9G15.5,I8)') MLPDAT(J1,1:MLPIN),MLPOUTCOME(J1)
 3748:             ENDDO
 3749:             CLOSE(LUNIT)
 3750:          ENDIF
 3751: ! 
 3752: ! MLPNORM directs OPTIM to rescale the input data columns by dividing each one by the
 3753: ! average of the mean magnitude
 3754: ! 
 3755:       ELSE IF (WORD.EQ.'MLPNORM') THEN
3741:          LUNIT=GETUNIT()3756:          LUNIT=GETUNIT()
3742:          OPEN(LUNIT,FILE='MLPdata',STATUS='OLD')3757:          OPEN(LUNIT,FILE='MLPdata',STATUS='OLD')
 3758:          IF (ALLOCATED(MLPDAT)) DEALLOCATE(MLPDAT)
 3759:          IF (ALLOCATED(MLPOUTCOME)) DEALLOCATE(MLPOUTCOME)
3743:          ALLOCATE(MLPDAT(MLPDATA,MLPIN),MLPOUTCOME(MLPDATA),MLPMEAN(MLPIN))3760:          ALLOCATE(MLPDAT(MLPDATA,MLPIN),MLPOUTCOME(MLPDATA),MLPMEAN(MLPIN))
3744:          MLPMEAN(1:MLPIN)=0.0D03761:          MLPMEAN(1:MLPIN)=0.0D0
3745:          DO J1=1,MLPDATA3762:          DO J1=1,MLPDATA
3746:             READ(LUNIT,*) MLPDAT(J1,1:MLPIN),MLPOUTCOME(J1)3763:             READ(LUNIT,*) MLPDAT(J1,1:MLPIN),MLPOUTCOME(J1)
3747:             MLPOUTCOME(J1)=MLPOUTCOME(J1)+1 ! to shift the range to 1 to 43764:             MLPOUTCOME(J1)=MLPOUTCOME(J1)+1 ! to shift the range to start from 1 instead of zero
3748: !           WRITE(*,'(9G15.5,I8)') MLPDAT(J1,1:MLPIN),MLPOUTCOME(J1) 
3749:             DO J2=1,MLPIN3765:             DO J2=1,MLPIN
3750:                MLPMEAN(J2)=MLPMEAN(J2)+ABS(MLPDAT(J1,J2))3766:                MLPMEAN(J2)=MLPMEAN(J2)+ABS(MLPDAT(J1,J2))
3751:             ENDDO3767:             ENDDO
3752:          ENDDO3768:          ENDDO
3753:          CLOSE(LUNIT)3769:          CLOSE(LUNIT)
3754:          IF (MLPNORM) THEN3770:          MLPMEAN(1:MLPIN)=MLPMEAN(1:MLPIN)/MLPDATA
3755:             MLPMEAN(1:MLPIN)=MLPMEAN(1:MLPIN)/MLPDATA3771:          PRINT '(A)','keyword> Rescaling inputs by mean absolute values:'
3756:             WRITE(MYUNIT,'(A)') 'keyword> Rescaling inputs by mean absolute values:'3772:          PRINT '(6G20.10)',MLPMEAN(1:MLPIN)
3757:             WRITE(MYUNIT,'(6G20.10)') MLPMEAN(1:MLPIN)3773:          DO J1=1,MLPIN
3758:             DO J1=1,MLPIN3774:             MLPDAT(1:MLPDATA,J1)=MLPDAT(1:MLPDATA,J1)/MLPMEAN(J1)
3759:                MLPDAT(1:MLPDATA,J1)=MLPDAT(1:MLPDATA,J1)/MLPMEAN(J1)3775:          ENDDO
3760:             ENDDO 
3761:          ENDIF 
3762:          DEALLOCATE(MLPMEAN) 
3763:          MLPDONE=.TRUE.3776:          MLPDONE=.TRUE.
3764: ! 3777: !        PRINT '(A)','keyword> rescaled MLPdata:'
3765: ! MLPNORM directs OPTIM to rescale the input data columns by dividing each one by the3778: !        DO J1=1,MLPDATA
3766: ! average of the mean magnitude3779: !           WRITE(*,'(50G15.5)') MLPDAT(J1,1:MLPIN)
3767: ! 3780: !        ENDDO
3768:       ELSE IF (WORD.EQ.'MLPNORM') THEN3781:          DEALLOCATE(MLPMEAN)
3769:          MLPNORM=.TRUE.3782: 
3770:          IF (MLPDONE) THEN 
3771:             LUNIT=GETUNIT() 
3772:             OPEN(LUNIT,FILE='MLPdata',STATUS='OLD') 
3773:             ALLOCATE(MLPMEAN(MLPIN)) 
3774:             MLPMEAN(1:MLPIN)=0.0D0 
3775:             DO J1=1,MLPDATA 
3776:                READ(LUNIT,*) MLPDAT(J1,1:MLPIN),MLPOUTCOME(J1) 
3777:                MLPOUTCOME(J1)=MLPOUTCOME(J1)+1 ! to shift the range to start from 1 instead of zero 
3778:                DO J2=1,MLPIN 
3779:                   MLPMEAN(J2)=MLPMEAN(J2)+ABS(MLPDAT(J1,J2)) 
3780:                ENDDO 
3781:             ENDDO 
3782:             CLOSE(LUNIT) 
3783:             MLPMEAN(1:MLPIN)=MLPMEAN(1:MLPIN)/MLPDATA 
3784:             PRINT '(A)','keyword> Rescaling inputs by mean absolute values:' 
3785:             PRINT '(6G20.10)',MLPMEAN(1:MLPIN) 
3786:             DO J1=1,MLPIN 
3787:                MLPDAT(1:MLPDATA,J1)=MLPDAT(1:MLPDATA,J1)/MLPMEAN(J1) 
3788:             ENDDO 
3789:             DEALLOCATE(MLPMEAN) 
3790:          ENDIF 
3791: ! 3783: ! 
3792: ! MODE n  specifies the eigenvector to follow                  - default n=03784: ! MODE n  specifies the eigenvector to follow                  - default n=0
3793: ! 3785: ! 
3794: 3786: 
3795:          ELSE IF (WORD.EQ.'MODE') THEN3787:          ELSE IF (WORD.EQ.'MODE') THEN
3796:             CALL READI(IVEC)3788:             CALL READI(IVEC)
3797:             IF (NITEMS.GT.2) THEN3789:             IF (NITEMS.GT.2) THEN
3798:                CALL READI(IVEC2)3790:                CALL READI(IVEC2)
3799:             ELSE3791:             ELSE
3800:                ! IVEC2=IVEC3792:                ! IVEC2=IVEC


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0