Esope source files

   1 : C SIMP2     SOURCE    CHAT      05/01/13    03:18:41     5004           
   2 :       SUBROUTINE SIMP2(A,M,N,MP,NP,L2,NL2,IP,KP,Q1,EPS)
   3 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC        
   4 : C                                                                     C
   5 : C     METHODE DU SIMPLEX : ESCLAVE DE SIMPLX                          C
   6 : C                                                                     C
   7 : C     P.PEGON 31/8/92                                                 C
   8 : C                                                                     C
   9 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8(A-H,O-Z)
  12 :       DIMENSION A(MP,NP),L2(MP)
  13 : C
  14 :       IP=0
  15 :       IF(NL2.LT.1)RETURN
  16 : C
  17 :       DO 11 I=1,NL2
  18 :         IF(A(L2(I)+1,KP+1).LT.-EPS)GOTO 2
  19 :  11   CONTINUE
  20 :       RETURN
  21 : C
  22 :  2    Q1=-A(L2(I)+1,1)/A(L2(I)+1,KP+1)
  23 :       IP=L2(I)
  24 :       IF(I+1.GT.NL2)RETURN
  25 :       DO 13 I=I+1,NL2
  26 :         II=L2(I)
  27 :         IF(A(II+1,KP+1).LT.-EPS)THEN
  28 :           Q=-A(II+1,1)/A(II+1,KP+1)
  29 :           IF(Q.LT.Q1)THEN
  30 :             IP=II
  31 :             Q1=Q
  32 :           ELSEIF(Q.EQ.Q1)THEN
  33 :             DO 12 K=1,N
  34 :               QP=-A(IP+1,K+1)/A(IP+1,KP+1)
  35 :               Q0=-A(II+1,K+1)/A(II+1,KP+1)
  36 :               IF(Q0.NE.QP)GOTO 6
  37 :  12         CONTINUE
  38 :  6          IF(Q0.LT.QP)IP=II
  39 :           ENDIF
  40 :         ENDIF
  41 :  13   CONTINUE
  42 : C
  43 :       RETURN
  44 :       END
  45 :