Sources Esope | Cast3M

prodt2
C PRODT2    SOURCE    CHAT      05/01/13    02:30:19     5004
      SUBROUTINE PRODT2(A,B,C,N)
C
C S.P. EFFECTUANT A=TRANSPOSEE(C)*B*C ,OU A(N,N),B(N),C(N,N)
C     EN FAIT B EST UNE MATRICE DIAGONALE REPRESENTEE
C     PAR LE VECTEUR DES VALEURS PROPRES
      IMPLICIT INTEGER(I-N)
      REAL*8 A(N,N),B(N),C(N,N)
      INTEGER N,I,J,K
      DO J=1,N
         DO K=1,N
            A(J,K)=0.D0
            DO I=1,N
               A(J,K)=A(J,K)+B(I)*C(J,I)*C(K,I)
            END DO
         END DO
      END DO
      RETURN
      END