Télécharger la plaquette Cast3M


fr en

Télécharger vmult4.eso

Retour à la liste

Numérotation des lignes : 

vmult4

  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