Sources Esope | Cast3M

prodt
C PRODT     SOURCE    CHAT      05/01/13    02:30:24     5004
      SUBROUTINE PRODT(A,B,C,I,J)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C S.P. EFFECTUANT A=TRANSPOSEE(C)*B*C ,OU A(I,I),B(J,J),C(J,I)---(A COMP
      DIMENSIONA(I,1),B(J,1),C(J,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