Esope source files

   1 : C SJACO8    SOURCE    CHAT      05/11/16    21:15:30     5235           
   2 : *****************************************************
   3 : *
   4 : * Normalisation de phi
   5 : *
   6 : ****************************************************     
   7 : C
   8 :           SUBROUTINE SJACO8(IPHI,IPM1)
   9 : C
  10 :       IMPLICIT INTEGER(I-N)
  11 :           SEGMENT ,MATRIX
  12 :            REAL*8 A(N,N)
  13 :           ENDSEGMENT
  14 : C
  15 :           POINTEUR IPM1.MATRIX,IPHI.MATRIX,
  16 :      1    IPHIT.MATRIX,IPRE1.MATRIX,IPRE2.MATRIX
  17 : C         
  18 :           INTEGER N,I,J
  19 :           REAL*8 F,G
  20 : C  
  20 :           CALL SJACO(IPHI,IPHIT)
  20 :           CALL SJACO9 (IPHIT,IPM1,IPRE1)
  20 :           CALL SJACO9 (IPRE1,IPHI,IPRE2)
  21 : C
  22 : *          SEGACT ,IPHI*MOD,IPRE2
  23 :           N=IPHI.A(/1)
  24 : C
  25 :           DO 100 J=1,N
  26 :           G=SQRT(ABS(IPRE2.A(J,J)))
  27 :              DO 200 I=1,N
  28 :                F= IPHI.A(I,J)
  29 :                IPHI.A(I,J)=F/G
  30 : 200          CONTINUE
  31 : 100       CONTINUE
  32 : *          SEGDES ,IPHI,IPM1
  33 :           SEGSUP ,IPHIT,IPRE1,IPRE2
  34 :           RETURN
  35 :           END
  36 :  
  37 :  
  38 :  
  39 :