Télécharger la plaquette Cast3M


fr en

Télécharger vmult1.eso

Retour à la liste

Numérotation des lignes : 

vmult1

  1. C VMULT1    SOURCE    CHAT      05/01/13    04:08:24     5004
  2.       SUBROUTINE VMULT1(A,AA,B)
  3. C=======================================================================
  4. C     CALCULE LE PRODUIT  AA = TRPOSE ( B ) * A * B
  5. C    ENTREES
  6. C     A(9)   = MATRICE STOCKEE PAR COLONNE  ( SYMMMETRIQUE
  7. C              A11=A(1) A12=A(4)  A21=A(2) A22=A(5) A33=A(9) NON NULS
  8. C     B(3,3) = MATRICE DE PASSAGE           ( NON SYMMETRIQUE
  9. C    SORTIES
  10. C     AA(9)  =   MATRICE  STOCKEE PAR COLONNE ( SYMMMETRIQUE
  11. C=======================================================================
  12.       IMPLICIT INTEGER(I-N)
  13.       IMPLICIT REAL*8(A-H,O-Z)
  14.       DIMENSION A(*),AA(*),B(3,*)
  15. C
  16.       B11=B(1,1)
  17.       B12=B(1,2)
  18.       B13=B(1,3)
  19.       B21=B(2,1)
  20.       B22=B(2,2)
  21.       B23=B(2,3)
  22.       B31=B(3,1)
  23.       B32=B(3,2)
  24.       B33=B(3,3)
  25. C
  26.       C1=B11*A(1)+B21*A(4)
  27.       C2=B11*A(4)+B21*A(5)
  28.       C3=B12*A(1)+B22*A(4)
  29.       C4=B12*A(4)+B22*A(5)
  30.       C5=B13*A(1)+B23*A(4)
  31.       C6=B13*A(4)+B23*A(5)
  32.       C7=                  B31*A(9)
  33.       C8=                  B32*A(9)
  34.       C9=                  B33*A(9)
  35. C
  36.       AA(1)=B11*C1+B21*C2+B31*C7
  37.       AA(4)=B12*C1+B22*C2+B32*C7
  38.       AA(5)=B12*C3+B22*C4+B32*C8
  39.       AA(7)=B13*C1+B23*C2+B33*C7
  40.       AA(8)=B12*C5+B22*C6+B32*C9
  41.       AA(9)=B13*C5+B23*C6+B33*C9
  42.       AA(2)=AA(4)
  43.       AA(3)=AA(7)
  44.       AA(6)=AA(8)
  45.       RETURN
  46.       END
  47.  
  48.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales