1 : C VMULT2 SOURCE CHAT 05/01/13 04:08:27 5004 2 : SUBROUTINE VMULT2(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 A31,A32 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 : C1=B11*A(3)+B21*A(6) 27 : C2=B12*A(3)+B22*A(6) 28 : C3=B13*A(3)+B23*A(6) 29 : C 30 : AA(1)=B31*C1 31 : AA(2)=B32*C1 32 : AA(3)=B33*C1 33 : AA(4)=B31*C2 34 : AA(5)=B32*C2 35 : AA(6)=B33*C2 36 : AA(7)=B31*C3 37 : AA(8)=B32*C3 38 : AA(9)=B33*C3 39 : RETURN 40 : END 41 :
© Cast3M 2003 - All rights reserved.
Disclaimer