Esope source files

   1 : C SJACO9    SOURCE    CHAT      05/11/16    21:15:33     5235           
   2 : *****************************************************
   3 : *
   4 : * Multiplie la matrice A par B 
   5 : *
   6 : ****************************************************     
   7 : C
   8 :           SUBROUTINE SJACO9(IPA,IPB,IPC)
   9 : C
  10 :       IMPLICIT INTEGER(I-N)
  11 :           SEGMENT ,MATRIX
  12 :            REAL*8 A(N,N)
  13 :           ENDSEGMENT
  14 : C
  15 :           POINTEUR IPA.MATRIX,IPB.MATRIX,IPC.MATRIX
  16 : C         
  17 :           INTEGER N,I,J,K
  18 :           REAL*8 S
  19 : C
  20 : *          SEGACT ,IPA,IPB
  21 : C
  22 :           N=IPA.A(/1)
  23 : C
  24 :           SEGINI ,IPC
  25 : C
  26 :           DO 100 I=1,N
  27 : C
  28 :              DO 200 J=1,N
  29 :                S=0.D0
  30 : C
  31 :                DO 300 K=1,N
  32 : C
  33 :                  S=S+(IPA.A(I,K)*IPB.A(K,J))
  34 : C
  35 : 300            CONTINUE
  36 : C
  37 :              IPC.A(I,J)=S
  38 : C
  39 : 200          CONTINUE
  40 : C
  41 : 100       CONTINUE
  42 : C          
  43 : *          SEGDES ,IPA,IPB,IPC
  44 :           RETURN
  45 :           END
  46 :  
  47 :