r2sort.f

Go to the documentation of this file.

      SUBROUTINE r2sort(IJ,RIJ2)
*
*
*       Sorting of square chain distances.
*       ---------------------------------
*
      include 'commonc.h'
      REAL*8 rij2(nmx,nmx)
      INTEGER ij(nmx)
*
*
      rijmin=1.e30
      fmax = 0.0
      DO i=1,n-1
      DO j=i+1,n
      rij2(i,j) = (x(3*i-2)-x(3*j-2))**2+(x(3*i-1)-x(3*j-1))**2+
     &            (x(3*i)-x(3*j))**2
      rij2(j,i)=rij2(i,j)
*     IF(RIJ2(I,J).LT.RIJMIN)THEN
*     RIJMIN=RIJ2(I,J)
*       Save dominant pair instead of smallest distance (08/09).
      IF ((m(i)+m(j))/rij2(i,j).GT.fmax)THEN
      fmax=(m(i)+m(j))/rij2(i,j)
      ij(1)=i
      ij(2)=j
      END IF
      END DO
      END DO
*
      l=2
      DO i=1,n
      IF(i.NE.ij(1).AND.i.NE.ij(2))THEN
      l=l+1
      ij(l)=i
      END IF
      END DO
*
      RETURN
      END