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 :
© Cast3M 2003 - All rights reserved.
Disclaimer