1 : C VMULT3 SOURCE CHAT 05/01/13 04:08:30 5004 2 : SUBROUTINE VMULT3(A,AA,B) 3 : C======================================================================= 4 : C CALCULE LE PRODUIT AA = TRPOSE( B ) * A * B 5 : C ENTREES 6 : C A(9) = MATRICE NON SYMETRIQUE STOCKEE PAR COLONNE 7 : C A13,A23 SONT NON NULS SEULEMENT 8 : C B(3,3) = MATRICE DE PASSAGE 9 : C SORTIES 10 : C AA(9) = MATRICE NON SYMETRIQUE STOCKEE PAR COLONNE 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 : AA(1)=B11*A(7)*B31+B21*A(8)*B31 27 : AA(2)=B12*A(7)*B31+B22*A(8)*B31 28 : AA(3)=B13*A(7)*B31+B23*A(8)*B31 29 : AA(4)=B11*A(7)*B32+B21*A(8)*B32 30 : AA(5)=B12*A(7)*B32+B22*A(8)*B32 31 : AA(6)=B13*A(7)*B32+B23*A(8)*B32 32 : AA(7)=B11*A(7)*B33+B21*A(8)*B33 33 : AA(8)=B12*A(7)*B33+B22*A(8)*B33 34 : AA(9)=B13*A(7)*B33+B23*A(8)*B33 35 : RETURN 36 : END 37 :
© Cast3M 2003 - All rights reserved.
Disclaimer