C SJACO8 SOURCE CHAT 05/11/16 21:15:30 5235 ***************************************************** * * Normalisation de phi * **************************************************** C SUBROUTINE SJACO8(IPHI,IPM1) C IMPLICIT INTEGER(I-N) SEGMENT ,MATRIX REAL*8 A(N,N) ENDSEGMENT C POINTEUR IPM1.MATRIX,IPHI.MATRIX, 1 IPHIT.MATRIX,IPRE1.MATRIX,IPRE2.MATRIX C INTEGER N,I,J REAL*8 F,G C CALL SJACO(IPHI,IPHIT) CALL SJACO9 (IPHIT,IPM1,IPRE1) CALL SJACO9 (IPRE1,IPHI,IPRE2) C * SEGACT ,IPHI*MOD,IPRE2 N=IPHI.A(/1) C DO 100 J=1,N G=SQRT(ABS(IPRE2.A(J,J))) DO 200 I=1,N F= IPHI.A(I,J) IPHI.A(I,J)=F/G 200 CONTINUE 100 CONTINUE * SEGDES ,IPHI,IPM1 SEGSUP ,IPHIT,IPRE1,IPRE2 RETURN END