kqsortc.f

Go to the documentation of this file.

      SUBROUTINE kqsortc(A, ORD, N)
!
!==============SORTS THE ARRAY A(I),I=1,2,...,N BY PUTTING THE
!   ASCENDING ORDER VECTOR IN ORD.  THAT IS ASCENDING ORDERED A
!   IS A(ORD(I)),I=1,2,...,N; DESCENDING ORDER A IS A(ORD(N-I+1)),
!   I=1,2,...,N .  THIS SORT RUNS IN TIME PROPORTIONAL TO N LOG N .
!
!
!     ACM QUICKSORT - ALGORITHM #402 - IMPLEMENTED IN FORTRAN BY
!                                 WILLIAM H. VERITY
!                                 COMPUTATION CENTER
!                                 PENNSYLVANIA STATE UNIVERSITY
!                                 UNIVERSITY PARK, PA.  16802
!
      IMPLICIT INTEGER (A-Z)
!
      dimension ord(n),poplst(2,20)
      INTEGER X,XX,Z,ZZ,Y
      CHARACTER *(*) A(n)
!
      ndeep=0
      u1=n
      l1=1
      DO 1  i=1,n
    1 ord(i)=i
    2 IF (u1.GT.l1) GO TO 3
      RETURN
!
    3 l=l1
      u=u1
!
! PART
!
    4 p=l
      q=u
      x=ord(p)
      z=ord(q)
      IF (a(x).LE.a(z)) GO TO 5
      y=x
      x=z
      z=y
      yp=ord(p)
      ord(p)=ord(q)
      ord(q)=yp
    5 IF (u-l.LE.1) GO TO 15
      xx=x
      ix=p
      zz=z
      iz=q
!
! LEFT
!
    6 p=p+1
      IF (p.GE.q) GO TO 7
      x=ord(p)
      IF (a(x).GE.a(xx)) GO TO 8
      GO TO 6
    7 p=q-1
      GO TO 13
!
! RIGHT
!
    8 q=q-1
      IF (q.LE.p) GO TO 9
      z=ord(q)
      IF (a(z).LE.a(zz)) GO TO 10
      GO TO 8
    9 q=p
      p=p-1
      z=x
      x=ord(p)
!
! DIST
!
   10 IF (a(x).LE.a(z)) GO TO 11
      y=x
      x=z
      z=y
      ip=ord(p)
      ord(p)=ord(q)
      ord(q)=ip
   11 IF (a(x).LE.a(xx)) GO TO 12
      xx=x
      ix=p
   12 IF (a(z).GE.a(zz)) GO TO 6
      zz=z
      iz=q
      GO TO 6
!
! OUT
!
   13 CONTINUE
      IF (.NOT.(p.NE.ix.AND.a(x).NE.a(xx))) GO TO 14
      ip=ord(p)
      ord(p)=ord(ix)
      ord(ix)=ip
   14 CONTINUE
      IF (.NOT.(q.NE.iz.AND.a(z).NE.a(zz))) GO TO 15
      iq=ord(q)
      ord(q)=ord(iz)
      ord(iz)=iq
   15 CONTINUE
      IF (u-q.LE.p-l) GO TO 16
      l1=l
      u1=p-1
      l=q+1
      GO TO 17
   16 u1=u
      l1=q+1
      u=p-1
   17 CONTINUE
      IF (u1.LE.l1) GO TO 18
!
! START RECURSIVE CALL
!
      ndeep=ndeep+1
      poplst(1,ndeep)=u
      poplst(2,ndeep)=l
      GO TO 3
   18 IF (u.GT.l) GO TO 4
!
! POP BACK UP IN THE RECURSION LIST
!
      IF (ndeep.EQ.0) GO TO 2
      u=poplst(1,ndeep)
      l=poplst(2,ndeep)
      ndeep=ndeep-1
      GO TO 18
!
! END SORT
! END QSORT
!
      END