hdiff output

r30556/nnutils.f90 2016-06-03 19:30:16.142716144 +0100 r30555/nnutils.f90 2016-06-03 19:30:16.430719935 +0100
348:                             ! Probably both interpolations are bad and the connect run will fail, 348:                             ! Probably both interpolations are bad and the connect run will fail, 
349:                             ! but we may as well use the better of the two.349:                             ! but we may as well use the better of the two.
350:                             IF(DEBUG) WRITE(*,*) "Energy got worse. Resetting band"350:                             IF(DEBUG) WRITE(*,*) "Energy got worse. Resetting band"
351:                             XYZ = SAVE_BAND351:                             XYZ = SAVE_BAND
352:                         ELSE352:                         ELSE
353:                             IF(DEBUG) WRITE(*,*) "Energy improved after reinterpolation. Keeping new band"353:                             IF(DEBUG) WRITE(*,*) "Energy improved after reinterpolation. Keeping new band"
354:                             IF (EWORST_NEW.LT.COLL_TOL) THEN354:                             IF (EWORST_NEW.LT.COLL_TOL) THEN
355:                                 IF(DEBUG) WRITE(*,*) "Energy has fallen below tolerance. Hopefully all collisions have been averted"355:                                 IF(DEBUG) WRITE(*,*) "Energy has fallen below tolerance. Hopefully all collisions have been averted"
356:                                 EXIT356:                                 EXIT
357:                             ENDIF357:                             ENDIF
358:                             EWORST = EWORST_NEW 
359:                         ENDIF358:                         ENDIF
360:                     ENDDO359:                     ENDDO
361:                 ENDIF  ! End of IF (EWORST .GT. COLL_TOL)360:                 ENDIF  ! End of IF (EWORST .GT. COLL_TOL)
362:             ENDIF  ! End of IF (AVOID_COLLISIONS)361:             ENDIF  ! End of IF (AVOID_COLLISIONS)
363: 362: 
364: 363: 
365: 364: 
366: ! sf344> now we have the interpolated coordinates. For MULTISITEPY, check for overlap365: ! sf344> now we have the interpolated coordinates. For MULTISITEPY, check for overlap
367: !        and if there's any (very likely), move the offending bodies randomly around. 366: !        and if there's any (very likely), move the offending bodies randomly around. 
368:             IF(MULTISITEPYT) THEN    ! sn402 : added this IF statement367:             IF(MULTISITEPYT) THEN    ! sn402 : added this IF statement
936:         ! and last of these sets correspond to the start and endpoints, and the sets inbetween correspond935:         ! and last of these sets correspond to the start and endpoints, and the sets inbetween correspond
937:         ! to the images.936:         ! to the images.
938: 937: 
939:         IMPLICIT NONE938:         IMPLICIT NONE
940: 939: 
941:         LOGICAL, INTENT(IN) :: LONG_WAY940:         LOGICAL, INTENT(IN) :: LONG_WAY
942:         DOUBLE PRECISION,ALLOCATABLE :: QTN(:,:), PTN(:,:)941:         DOUBLE PRECISION,ALLOCATABLE :: QTN(:,:), PTN(:,:)
943:         DOUBLE PRECISION THETA, THETAH, ST, CT, P(3), FCT942:         DOUBLE PRECISION THETA, THETAH, ST, CT, P(3), FCT
944:         ! J1 is the index of the rigid body we are interested in. RBOFFSET is the point within the list943:         ! J1 is the index of the rigid body we are interested in. RBOFFSET is the point within the list
945:         ! of coordinates for this image at which the angle-axis components begin to be listed.944:         ! of coordinates for this image at which the angle-axis components begin to be listed.
946:         ! RBOFFSET is the index of the COORDS array at which the angle-axis coordinates start. The945:         ! RBOFFSET is different depending on which rigid-body representation is being used (GENRIGID or
947:         ! expression for this value is different depending on which rigid-body representation is being used946:         ! RBAA) so must be passed as a parameter.
948:         ! (GENRIGID or RBAA) so it must be passed as a parameter. 
949:         INTEGER J1, K, J, I, RBOFFSET947:         INTEGER J1, K, J, I, RBOFFSET
950: 948: 
951:         ! Lists of parallel and perpendicular quaternions for body J1 in all the images.949:         ! Lists of parallel and perpendicular quaternions for body J1 in all the images.
952:         ALLOCATE(QTN(NIMAGE+2,4))950:         ALLOCATE(QTN(NIMAGE+2,4))
953:         ALLOCATE(PTN(NIMAGE+2,4))951:         ALLOCATE(PTN(NIMAGE+2,4))
954: 952: 
955: !       CONVERT FROM AA -> QUATERNION FOR THE INITIAL AND FINAL FRAMES953: !       CONVERT FROM AA -> QUATERNION FOR THE INITIAL AND FINAL FRAMES
956:         DO K = 1, (NIMAGE+2), (NIMAGE+1)954:         DO K = 1, (NIMAGE+2), (NIMAGE+1)
957:            ! P is the angle-axis vector for rigid body J1 in the initial or final frame955:            ! P is the angle-axis vector for rigid body J1 in the initial or final frame
958:            J = (K-1)*NOPT + RBOFFSET + 3*(J1-1)956:            J = (K-1)*NOPT + RBOFFSET + 3*(J1-1)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0