Esope source files

   1 : C SJACO7    SOURCE    CHAT      05/11/16    21:15:27     5235           
   2 : *****************************************************
   3 : *
   4 : * affecte 1 à T si la convergence est réalisee
   5 : *
   6 : *****************************************************
   7 : C
   8 :           SUBROUTINE SJACO7(IPK1,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
  16 :           REAL*8 F,S
  17 : C         
  18 :           INTEGER N,T,I,J 
  19 : C
  20 : *          SEGACT ,IPK1
  21 : C
  22 :           T=1
  23 :           S=(1.D-12)
  24 : C
  25 :           DO 100 J=1,N
  26 : C
  27 :              DO 150 I=(J+1),N
  28 :                F=SQRT((IPK1.A(I,J)**2)/(IPK1.A(I,I)*IPK1.A(J,J)))
  29 :                IF (F.GT.S) THEN
  30 :                  T=0
  31 :                  GOTO 10
  32 :                ENDIF
  33 : 150          CONTINUE
  34 : C
  35 : 100       CONTINUE
  36 : C
  37 : 10      continue
  38 : *        SEGDES ,IPK1
  39 :           RETURN
  40 :           END
  41 :  
  42 :