Esope source files

   1 : C SIMP3     SOURCE    CHAT      05/01/13    03:18:45     5004           
   2 :       SUBROUTINE SIMP3(A,MP,NP,I1,K1,IP,KP)
   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)
  13 : C
  14 :       PIV=1.D0/A(IP+1,KP+1)
  15 :       IF(I1.GE.0.D0)THEN
  16 :         DO 12 II=1,I1+1
  17 :           IF(II-1.NE.IP)THEN
  18 :             A(II,KP+1)=A(II,KP+1)*PIV
  19 :             DO 11 KK=1,K1+1
  20 :               IF(KK-1.NE.KP)THEN
  21 :                 A(II,KK)=A(II,KK)-A(IP+1,KK)*A(II,KP+1)
  22 :               ENDIF
  23 :  11         CONTINUE
  24 :           ENDIF
  25 :  12     CONTINUE
  26 :       ENDIF
  27 : C
  28 :       DO 13 KK=1,K1+1
  29 :         IF(KK-1.NE.KP)A(IP+1,KK)=-A(IP+1,KK)*PIV
  30 :  13   CONTINUE
  31 :       A(IP+1,KP+1)=PIV
  32 : C
  33 :       RETURN
  34 :       END
  35 :