hdiff output

r31188/genrigid.f90 2016-09-23 11:30:21.046864277 +0100 r31187/genrigid.f90 2016-09-23 11:30:21.558871101 +0100
753:         DR2(:) = MATMUL(DRMI2,SITESRIGIDBODY(J2,:,J1))753:         DR2(:) = MATMUL(DRMI2,SITESRIGIDBODY(J2,:,J1))
754:         DR3(:) = MATMUL(DRMI3,SITESRIGIDBODY(J2,:,J1))754:         DR3(:) = MATMUL(DRMI3,SITESRIGIDBODY(J2,:,J1))
755:         GR(3*NRIGIDBODY+3*J1-2) = GR(3*NRIGIDBODY+3*J1-2) + DOT_PRODUCT(G(3*J9-2:3*J9),DR1(:))755:         GR(3*NRIGIDBODY+3*J1-2) = GR(3*NRIGIDBODY+3*J1-2) + DOT_PRODUCT(G(3*J9-2:3*J9),DR1(:))
756:         GR(3*NRIGIDBODY+3*J1-1) = GR(3*NRIGIDBODY+3*J1-1) + DOT_PRODUCT(G(3*J9-2:3*J9),DR2(:))756:         GR(3*NRIGIDBODY+3*J1-1) = GR(3*NRIGIDBODY+3*J1-1) + DOT_PRODUCT(G(3*J9-2:3*J9),DR2(:))
757:         GR(3*NRIGIDBODY+3*J1)   = GR(3*NRIGIDBODY+3*J1)   + DOT_PRODUCT(G(3*J9-2:3*J9),DR3(:))757:         GR(3*NRIGIDBODY+3*J1)   = GR(3*NRIGIDBODY+3*J1)   + DOT_PRODUCT(G(3*J9-2:3*J9),DR3(:))
758:      ENDDO758:      ENDDO
759:   ENDDO759:   ENDDO
760: 760: 
761: ! hk286 - testing 6/6/12761: ! hk286 - testing 6/6/12
762:   IF (FREEZERIGIDBODYT .EQV. .TRUE.) THEN762:   IF (FREEZERIGIDBODYT .EQV. .TRUE.) THEN
763:   write(*,*) "Danger: in frozen block" 
764:      GR(3*NRIGIDBODY-2:3*NRIGIDBODY) = 0.0D0763:      GR(3*NRIGIDBODY-2:3*NRIGIDBODY) = 0.0D0
765:      GR(6*NRIGIDBODY-2:6*NRIGIDBODY) = 0.0D0764:      GR(6*NRIGIDBODY-2:6*NRIGIDBODY) = 0.0D0
766:   ENDIF765:   ENDIF
767: 766: 
768: ! hk286 > single atoms767: ! hk286 > single atoms
769: ! vr274 > and lattice768: ! vr274 > and lattice
770:   IF (DEGFREEDOMS > 6 * NRIGIDBODY - NLATTICECOORDS) THEN769:   IF (DEGFREEDOMS > 6 * NRIGIDBODY - NLATTICECOORDS) THEN
771:   write(*,*) "Danger: in single atom block" 
772:      DO J1 = 1, (DEGFREEDOMS - 6*NRIGIDBODY - NLATTICECOORDS)/3770:      DO J1 = 1, (DEGFREEDOMS - 6*NRIGIDBODY - NLATTICECOORDS)/3
773:         J9 = RIGIDSINGLES(J1)771:         J9 = RIGIDSINGLES(J1)
774:         GR(6*NRIGIDBODY + 3*J1-2:6*NRIGIDBODY + 3*J1) = G(3*J9-2:3*J9)772:         GR(6*NRIGIDBODY + 3*J1-2:6*NRIGIDBODY + 3*J1) = G(3*J9-2:3*J9)
775:      ENDDO773:      ENDDO
776:   ENDIF774:   ENDIF
777: 775: 
778:   IF(HAS_LATTICE_COORDS) THEN776:   IF(HAS_LATTICE_COORDS) THEN
779:   write(*,*) "Danger: in lattice coord block" 
780:       CALL GET_LATTICE_MATRIX(XR(DEGFREEDOMS-5:DEGFREEDOMS),MLATTICE)777:       CALL GET_LATTICE_MATRIX(XR(DEGFREEDOMS-5:DEGFREEDOMS),MLATTICE)
781: 778: 
782:       ! vr274> for lattice, go to reduced coordinates779:       ! vr274> for lattice, go to reduced coordinates
783:       DO J1 = 1, NRIGIDBODY780:       DO J1 = 1, NRIGIDBODY
784:           GR(3*J1-2:3*J1) =  matmul(transpose(mlattice), GR(3*J1-2:3*J1))781:           GR(3*J1-2:3*J1) =  matmul(transpose(mlattice), GR(3*J1-2:3*J1))
785:       ENDDO782:       ENDDO
786:       ! vr274> and single atoms783:       ! vr274> and single atoms
787:       IF (DEGFREEDOMS > 6 * NRIGIDBODY + NLATTICECOORDS) THEN784:       IF (DEGFREEDOMS > 6 * NRIGIDBODY + NLATTICECOORDS) THEN
788:           DO J1 = 1, (DEGFREEDOMS - 6*NRIGIDBODY - NLATTICECOORDS)/3785:           DO J1 = 1, (DEGFREEDOMS - 6*NRIGIDBODY - NLATTICECOORDS)/3
789:               J2 = 6*NRIGIDBODY + 3*J1786:               J2 = 6*NRIGIDBODY + 3*J1


r31188/ncutils.f90 2016-09-23 11:30:20.794860919 +0100 r31187/ncutils.f90 2016-09-23 11:30:21.306867742 +0100
2219:         ENDIF2219:         ENDIF
2220:         PREVIOUSTS(1:NOPT)=TS(MI(DUMMY%I)%DATA%CTS(DUMMY%J))%DATA%X(1:NOPT)2220:         PREVIOUSTS(1:NOPT)=TS(MI(DUMMY%I)%DATA%CTS(DUMMY%J))%DATA%X(1:NOPT)
2221:         DUMMY=>DUMMY%NEXT2221:         DUMMY=>DUMMY%NEXT
2222:         I=I+12222:         I=I+1
2223:      ENDDO2223:      ENDDO
2224:      CALL FLUSH(88)2224:      CALL FLUSH(88)
2225:      CLOSE(88)2225:      CLOSE(88)
2226: 2226: 
2227:      END SUBROUTINE MAKEPATHINFO2227:      END SUBROUTINE MAKEPATHINFO
2228: 2228: 
2229: ! -------------------------------------------------------------------------------------------------------------------- 
2230:  
2231: !2229: !
2232: !  Dump min to path.info2230: !  Dump min to path.info
2233: !2231: !
2234:      SUBROUTINE MAKEMINPATHINFO(MINIMUM)2232:      SUBROUTINE MAKEMINPATHINFO(MINIMUM)
2235:      USE SYMINF 2233:      USE SYMINF 
2236:      USE MODHESS2234:      USE MODHESS
2237:      USE MODCHARMM2235:      USE MODCHARMM
2238:      USE MODUNRES2236:      USE MODUNRES
2239:      USE KEY, ONLY: FILTH, FILTHSTR, UNRST, TWOD, BULKT, MACHINE, NOFRQS, AMBERT, NABT, AMHT, SEQ, TARFL, NRES_AMH_TEMP, SDT, &2237:      USE KEY, ONLY: FILTH, FILTHSTR, UNRST, TWOD, BULKT, MACHINE, NOFRQS, AMBERT, NABT, AMHT, SEQ, TARFL, NRES_AMH_TEMP, SDT, &
2240:           AMBER12T, RBAAT, MACROCYCLET, GTHOMSONT, TTM3T, BOWMANT, HESSDUMPT, INSTANTONSTARTDUMPT,FREEZE,NONFREEZE, METRICTENSOR ! hk2862238:           AMBER12T, RBAAT, MACROCYCLET, GTHOMSONT, TTM3T, BOWMANT, HESSDUMPT, INSTANTONSTARTDUMPT,FREEZE,NONFREEZE, METRICTENSOR ! hk286


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0