Esope source files

   1 : C PRODT2    SOURCE    CHAT      05/01/13    02:30:19     5004           
   2 :       SUBROUTINE PRODT2(A,B,C,N)
   3 : C
   4 : C S.P. EFFECTUANT A=TRANSPOSEE(C)*B*C ,OU A(N,N),B(N),C(N,N)
   5 : C     EN FAIT B EST UNE MATRICE DIAGONALE REPRESENTEE 
   6 : C     PAR LE VECTEUR DES VALEURS PROPRES
   7 :       IMPLICIT INTEGER(I-N)
   8 :       REAL*8 A(N,N),B(N),C(N,N)
   9 :       INTEGER N,I,J,K
  10 :       DO J=1,N
  11 :          DO K=1,N
  12 :             A(J,K)=0.D0
  13 :             DO I=1,N
  14 :                A(J,K)=A(J,K)+B(I)*C(J,I)*C(K,I)
  15 :             END DO
  16 :          END DO
  17 :       END DO
  18 :       RETURN
  19 :       END
  20 :