Esope source files

   1 : C SJACO6    SOURCE    CHAT      05/11/16    21:15:23     5235           
   2 : *****************************************************
   3 : *
   4 : * affecte 1 à T si la convergence est réalisee
   5 : *
   6 : *****************************************************
   7 : C
   8 :           SUBROUTINE SJACO6(IPK1,IPK2,IPM1,IPM2,T,N,K)
   9 : C
  10 :       IMPLICIT INTEGER(I-N)
  11 :           SEGMENT ,MATRIX
  12 :            REAL*8 A(N,N)
  13 :           ENDSEGMENT
  14 : C
  15 :           POINTEUR IPK1.MATRIX,IPK2.MATRIX,IPM1.MATRIX,IPM2.MATRIX
  16 : C         
  17 :           INTEGER N,T,I
  18 :           REAL*8 F,S,TOL
  19 : C
  20 : *          SEGACT ,IPK1,IPK2,IPM1,IPM2
  21 : C
  22 :           T=1
  23 : *          S=10.D-12
  24 :           S=1.D-12
  25 : C
  26 :           DO 100 I=1,N
  27 : C
  28 :             TOL=S*( IPK1.A(I,I)/IPM1.A(I,I) )
  29 :             F=ABS( IPK1.A(I,I)/IPM1.A(I,I)
  30 :      1           - IPK2.A(I,I)/IPM2.A(I,I) ) 
  31 : C
  32 :             IF (F.GT.TOL) THEN      
  33 :               T=0
  34 :               GOTO 10
  35 :             ENDIF
  36 : C
  37 : 100       CONTINUE
  38 : C
  39 : 10      continue
  40 : *         SEGDES ,IPK1,IPK2,IPM1,IPM2
  41 : C          
  42 :           RETURN
  43 :           END 
  44 :  
  45 :