Esope source files

   1 : C SODIRI    SOURCE    CHAT      05/01/13    03:21:20     5004           
   2 :         SUBROUTINE SODIRI(D,X,B,MM,LL,KERRE)
   3 : *
   4 : *
   5 : *      resolution d'un du systeme A.X=B
   6 : *      par la methode du pivot de Gauss
   7 : *
   8 :       IMPLICIT INTEGER(I-N)
   9 :        IMPLICIT REAL*8(A-H,O-Z)
  10 : -INC CCOPTIO
  11 : *
  12 :         REAL*8   X(*), B(*)
  13 :       DIMENSION D(1)                                                            
  14 : *            DO 12 I=1,LL
  15 : *            DO 13 J=I,LL
  16 : *              D(I+J)=A(I,J)
  17 : *13          CONTINUE
  18 : *12          CONTINUE
  19 : 
  20 :       KERRE=0
  21 :       PREC=1.e-12
  22 :       DO 200 N=1,MM
  23 :       NN=LL*(N-1)+N
  24 :       DIAG=D(NN)
  25 :       IF(ABS(DIAG).LE.PREC) THEN
  26 :       KERRE=49
  27 :       RETURN
  28 :       ENDIF
  29 :       DIAG=1./DIAG
  30 :       DO 100 J=1,MM 
  31 :       NJ=LL*(N-1)+J
  32 :   100 D(NJ)=-D(NJ)*DIAG
  33 :       DO 150 I=1,MM
  34 :       IN=LL*(I-1)+N
  35 :       IF(I.EQ.N) GO TO 150
  36 :       DO 140 J=1,MM
  37 :       IF(J.EQ.N) GO TO 140
  38 :       NJ=LL*(N-1)+J
  39 :       IJ=LL*(I-1)+J
  40 :       D(IJ)=D(IJ)+D(IN)*D(NJ) 
  41 :   140 CONTINUE 
  42 :   150 D(IN)=D(IN)*DIAG
  43 :   200 D(NN)=DIAG
  44 :         DO 20 I=1,MM
  45 :           X(I)=0
  46 :            DO 21 J=1,MM
  47 :             IND=LL*(I-1)+J
  48 :             X(I)=X(I)+D(IND)*B(J)   
  49 : 21         CONTINUE
  50 :         
  51 : 20      CONTINUE
  52 :       RETURN                                                                    
  53 :       END                                                                       
  54 :  
  55 :