hdiff output

r32324/percsphere.f90 2017-04-24 11:30:12.512708767 +0100 r32323/percsphere.f90 2017-04-24 11:30:12.732711621 +0100
155: ! Set up ROTVEC as the angle-axis rotation vector to take JR2 to JR1155: ! Set up ROTVEC as the angle-axis rotation vector to take JR2 to JR1
156:         ROTVEC(1) = P(3*JR1-1) * P(3*JR2)   - P(3*JR2-1) * P(3*JR1)156:         ROTVEC(1) = P(3*JR1-1) * P(3*JR2)   - P(3*JR2-1) * P(3*JR1)
157:         ROTVEC(2) = P(3*JR1)   * P(3*JR2-2) - P(3*JR2)   * P(3*JR1-2)157:         ROTVEC(2) = P(3*JR1)   * P(3*JR2-2) - P(3*JR2)   * P(3*JR1-2)
158:         ROTVEC(3) = P(3*JR1-2) * P(3*JR2-1) - P(3*JR2-2) * P(3*JR1-1)158:         ROTVEC(3) = P(3*JR1-2) * P(3*JR2-1) - P(3*JR2-2) * P(3*JR1-1)
159:         DUMMY = DOT_PRODUCT(ROTVEC,ROTVEC)159:         DUMMY = DOT_PRODUCT(ROTVEC,ROTVEC)
160:         ROTVEC = THETA * ROTVEC / DSQRT(DUMMY)160:         ROTVEC = THETA * ROTVEC / DSQRT(DUMMY)
161: 161: 
162: ! Test rotation of JR2 onto JR1162: ! Test rotation of JR2 onto JR1
163:         CALL RMDRVT(ROTVEC, RMI, DRMI, DRMI, DRMI, .FALSE.)163:         CALL RMDRVT(ROTVEC, RMI, DRMI, DRMI, DRMI, .FALSE.)
164:         TEMP = MATMUL(RMI,P(3*JR2-2:3*JR2))164:         TEMP = MATMUL(RMI,P(3*JR2-2:3*JR2))
165: !        DUMMY=(TEMP(3*JR2-2)-P(3*JR1-2))**2+(TEMP(3*JR2-1)-P(3*JR1-1))**2+(TEMP(3*JR2)-P(3*JR1))**2165:         DUMMY=(TEMP(3*JR2-2)-P(3*JR1-2))**2+(TEMP(3*JR2-1)-P(3*JR1-1))**2+(TEMP(3*JR2)-P(3*JR1))**2
166:         DUMMY = (TEMP(1) - P(3*JR1-2))**2 + (TEMP(2) - P(3*JR1-2))**2 + (TEMP(3) - P(3*JR1-2))**2 
167:         IF (DUMMY > PERCCUT) THEN166:         IF (DUMMY > PERCCUT) THEN
168: ! Reverse the angle if we've rotated the wrong way167: ! Reverse the angle if we've rotated the wrong way
169:            ROTVEC = -1.0D0 * ROTVEC168:            ROTVEC = -1.0D0 * ROTVEC
170:         ENDIF169:         ENDIF
171: 170: 
172: ! Reduce the rotation a bit, so we don't put bodies on top of each other171: ! Reduce the rotation a bit, so we don't put bodies on top of each other
173:         ROTVEC = ROTVEC * 0.975D0172:         ROTVEC = ROTVEC * 0.975D0
174:         CALL RMDRVT(ROTVEC, RMI, DRMI, DRMI, DRMI, .FALSE.)173:         CALL RMDRVT(ROTVEC, RMI, DRMI, DRMI, DRMI, .FALSE.)
175: 174: 
176: ! Rotate all the bodies connected to the root of cluster J3175: ! Rotate all the bodies connected to the root of cluster J3


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0