Sources Esope | Cast3M

sodiri
C SODIRI    SOURCE    CHAT      05/01/13    03:21:20     5004
        SUBROUTINE SODIRI(D,X,B,MM,LL,KERRE)
*
*
*      resolution d'un du systeme A.X=B
*      par la methode du pivot de Gauss
*
      IMPLICIT INTEGER(I-N)
       IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
*
        REAL*8   X(*), B(*)
      DIMENSION D(1)
*            DO 12 I=1,LL
*            DO 13 J=I,LL
*              D(I+J)=A(I,J)
*13          CONTINUE
*12          CONTINUE
 
      KERRE=0
      PREC=1.e-12
      DO 200 N=1,MM
      NN=LL*(N-1)+N
      DIAG=D(NN)
      IF(ABS(DIAG).LE.PREC) THEN
      KERRE=49
      RETURN
      ENDIF
      DIAG=1./DIAG
      DO 100 J=1,MM
      NJ=LL*(N-1)+J
  100 D(NJ)=-D(NJ)*DIAG
      DO 150 I=1,MM
      IN=LL*(I-1)+N
      IF(I.EQ.N) GO TO 150
      DO 140 J=1,MM
      IF(J.EQ.N) GO TO 140
      NJ=LL*(N-1)+J
      IJ=LL*(I-1)+J
      D(IJ)=D(IJ)+D(IN)*D(NJ)
  140 CONTINUE
  150 D(IN)=D(IN)*DIAG
  200 D(NN)=DIAG
        DO 20 I=1,MM
          X(I)=0
           DO 21 J=1,MM
            IND=LL*(I-1)+J
            X(I)=X(I)+D(IND)*B(J)
21         CONTINUE
 
20      CONTINUE
      RETURN
      END