Sources Esope | Cast3M

simp2
C SIMP2     SOURCE    CHAT      05/01/13    03:18:41     5004
      SUBROUTINE SIMP2(A,M,N,MP,NP,L2,NL2,IP,KP,Q1,EPS)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C                                                                     C
C     METHODE DU SIMPLEX : ESCLAVE DE SIMPLX                          C
C                                                                     C
C     P.PEGON 31/8/92                                                 C
C                                                                     C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION A(MP,NP),L2(MP)
C
      IP=0
      IF(NL2.LT.1)RETURN
C
      DO 11 I=1,NL2
        IF(A(L2(I)+1,KP+1).LT.-EPS)GOTO 2
 11   CONTINUE
      RETURN
C
 2    Q1=-A(L2(I)+1,1)/A(L2(I)+1,KP+1)
      IP=L2(I)
      IF(I+1.GT.NL2)RETURN
      DO 13 I=I+1,NL2
        II=L2(I)
        IF(A(II+1,KP+1).LT.-EPS)THEN
          Q=-A(II+1,1)/A(II+1,KP+1)
          IF(Q.LT.Q1)THEN
            IP=II
            Q1=Q
          ELSEIF(Q.EQ.Q1)THEN
            DO 12 K=1,N
              QP=-A(IP+1,K+1)/A(IP+1,KP+1)
              Q0=-A(II+1,K+1)/A(II+1,KP+1)
              IF(Q0.NE.QP)GOTO 6
 12         CONTINUE
 6          IF(Q0.LT.QP)IP=II
          ENDIF
        ENDIF
 13   CONTINUE
C
      RETURN
      END