Télécharger vmult4.eso

Retour à la liste

Numérotation des lignes :

  1. C VMULT4 SOURCE CHAT 05/01/13 04:08:33 5004
  2. SUBROUTINE VMULT4(A,AA,B)
  3. C=======================================================================
  4. C CALCULE LE PRODUIT AA = TRANSPOSE ( B ) * A * B
  5. C ENTREES
  6. C A(3,3) NON SYMETRIQUE
  7. C B(3,3) NON SYMETRIQUE
  8. C SORTIES
  9. C AA(9) NON SYMETRIQUE STOCKEE PAR COLONNE
  10. C=======================================================================
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8(A-H,O-Z)
  13. PARAMETER(XZER=0.D0)
  14. DIMENSION A(3,*),AA(*),B(3,*)
  15. DIMENSION T(3)
  16. C
  17. IJ=1
  18. DO 40 J=1,3
  19. DO 20 I1=1,3
  20. C=XZER
  21. DO 10 J1=1,3
  22. 10 C=C+A(I1,J1)*B(J1,J)
  23. 20 T(I1)=C
  24. DO 40 I=1,3
  25. C=XZER
  26. DO 30 J1=1,3
  27. 30 C=C+T(J1)*B(J1,I)
  28. AA(IJ)=C
  29. IJ=IJ+1
  30. 40 CONTINUE
  31. RETURN
  32. END
  33.  
  34.  

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