Sources Esope | Cast3M

prodt3
C PRODT3    SOURCE    CB215821  17/06/06    21:15:10     9449           
      SUBROUTINE PRODT3(A,B,C,I,J)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
C     EFFECTUANT A=TRANSPOSEE(C)*B*C ,OU A(I,I),B(J,J),C(3,I)
 
      DIMENSION A(I,1),B(J,1),C(3,1)
      DO 1 L=1,I
        DO 3 K=1,I
          A(K,L)=0.D0
   3    CONTINUE
        DO 1 M=1,J
          S=0.D0
          DO 2 N=1,J
            S=S+B(M,N)*C(N,L)
   2      CONTINUE
          DO 1 K=1,I
            A(K,L)=A(K,L)+C(M,K)*S
   1  CONTINUE
      RETURN
      END