iparti.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
integer function	iparti (i, j, tab, n2)

Function/Subroutine Documentation

iparti()

integer function iparti	(	integer	i,
		integer	j,
			tab,
		integer	n2 )

Definition at line 31 of file iparti.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER I, J, N2
C     REAL
        my_real
     .     tab(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER MILIEU, U, V
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
        INTEGER ICOMPC
        iparti=i
        IF(j>i) THEN
          milieu=(i+j)/2
          u=i
          v=j+1
C
          IF(icompc(j,milieu,tab)<0) THEN
            CALL echang(j,milieu,tab)
            CALL echang(n2+j,n2+milieu,tab)
          ENDIF
          IF(icompc(j,i,tab)<0) THEN
            CALL echang(j,i,tab)
            CALL echang(n2+j,n2+i,tab)
          ENDIF
          IF(icompc(i,milieu,tab)<0) THEN
            CALL echang(i,milieu,tab)
            CALL echang(n2+i,n2+milieu,tab)
          ENDIF
C
 1000     IF(u<v) THEN
C
 1100       CONTINUE
            u=u+1
            IF(icompc(u,i,tab)<0) GO TO 1100
C
 1200       CONTINUE
            v=v-1
            IF(icompc(v,i,tab)>0) GO TO 1200
C
            CALL echang(u,v,tab)
            CALL echang(n2+u,n2+v,tab)
            GO TO 1000
          ENDIF
          IF(u>v) CALL echang(u,v,tab)
          IF(u>v) CALL echang(n2+u,n2+v,tab)
          u=min0(u,v)
          CALL echang(i,u,tab)
          CALL echang(n2+i,n2+u,tab)
          iparti=u
        ENDIF
        RETURN