hdiff output

r31192/genrigid.f90 2016-09-23 11:30:12.662752519 +0100 r31191/genrigid.f90 2016-09-23 11:30:12.918755932 +0100
1693:   IR     = 0  ! IR+1 is the index of the first coordinate for the current rigid body1693:   IR     = 0  ! IR+1 is the index of the first coordinate for the current rigid body
1694:   IC     = 01694:   IC     = 0
1695:   ART = ATOMRIGIDCOORDT1695:   ART = ATOMRIGIDCOORDT
1696: 1696: 
1697: ! Transform coordinates1697: ! Transform coordinates
1698:   IF ( ATOMRIGIDCOORDT .EQV. .TRUE. ) THEN1698:   IF ( ATOMRIGIDCOORDT .EQV. .TRUE. ) THEN
1699:      CALL TRANSFORMCTORIGID (X, XRIGIDCOORDS)1699:      CALL TRANSFORMCTORIGID (X, XRIGIDCOORDS)
1700:   ELSE1700:   ELSE
1701:      XRIGIDCOORDS(1:DEGFREEDOMS) = X(1:DEGFREEDOMS)1701:      XRIGIDCOORDS(1:DEGFREEDOMS) = X(1:DEGFREEDOMS)
1702:   ENDIF1702:   ENDIF
1703:   1703:      
1704:   DO J1 = 1, NRIGIDBODY 
1705:  
1706:      J3 = 3*J11704:      J3 = 3*J1
1707:      J5 = OFFSET + J31705:      J5 = OFFSET + J3
1708:      P  = XRIGIDCOORDS(J5-2:J5)1706:      P  = XRIGIDCOORDS(J5-2:J5)
1709:      ! KBLOCK is the moment of inertia matrix for this rigid body.1707:      ! KBLOCK is the moment of inertia matrix for this rigid body.
1710:      KBLOCK(:,:) = 0.D0     1708:      KBLOCK(:,:) = 0.D0     
1711:      ! Get the rotation matrix RMI that corresponds to P (doesn't calculate any derivatives)1709:      ! Get the rotation matrix RMI that corresponds to P (doesn't calculate any derivatives)
1712:      CALL RMDFAS(P, RMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, .FALSE., .FALSE.)1710:      CALL RMDFAS(P, RMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, DRMI, .FALSE., .FALSE.)
1713: 1711: 
1714:      TMASS = 0.0D0 ! The total mass of the rigid body1712:      TMASS = 0.0D0 ! The total mass of the rigid body
1715:      DO J2 = 1, NSITEPERBODY(J1)1713:      DO J2 = 1, NSITEPERBODY(J1)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0