Télécharger vmult2.eso

Retour à la liste

Numérotation des lignes :

vmult2
  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.  
  42.  

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