Esope source files

   1 : C PRODT     SOURCE    CHAT      05/01/13    02:30:24     5004           
   2 :       SUBROUTINE PRODT(A,B,C,I,J)                                               
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C S.P. EFFECTUANT A=TRANSPOSEE(C)*B*C ,OU A(I,I),B(J,J),C(J,I)---(A COMP        
   6 :       DIMENSIONA(I,1),B(J,1),C(J,1)                                             
   7 :       DO 1 L=1,I                                                                
   8 :       DO 3 K=1,I                                                                
   9 :       A(K,L)=0.D0                                                               
  10 :    3  CONTINUE                                                                  
  11 :       DO 1 M=1,J                                                                
  12 :       S=0.D0                                                                    
  13 :       DO 2 N=1,J                                                                
  14 :       S=S+B(M,N)*C(N,L)                                                         
  15 :    2  CONTINUE                                                                  
  16 :       DO 1 K=1,I                                                                
  17 :       A(K,L)=A(K,L)+C(M,K)*S                                                    
  18 :    1  CONTINUE                                                                  
  19 :       RETURN                                                                    
  20 :       END                                                                       
  21 :