Télécharger promah.eso

Retour à la liste

Numérotation des lignes :

promah
  1. C PROMAH SOURCE CHAT 05/01/13 02:33:10 5004
  2. c Sous-Programme du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE PROMAH (H,X, Y)
  5. C --------------------------------------------------------------------------
  6. C Calcul des coefficients de la matrice Y = M.X, en fonction de ceux de la
  7. C matrice 6x6 M et de la matrice colonne X a 6 lignes pour M de la forme :
  8. C M11 = H(2)+H(3) , M22 = H(3)+H(1) , M33 = H(1)+H(2)
  9. C M12 = M21 = -H(3) , M13 = M31 = -H(2) , M23 = M32 = -H(1)
  10. C M44 = H(4) , M55 = H(5) , M66 = H(6)
  11. C les autres coefficients etant nuls.
  12. C --------------------------------------------------------------------------
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8 (A-H, O-Z)
  15. DIMENSION H(1:6),X(1:6),Y(1:6)
  16. Y(1) = (H(2)+H(3))*X(1)-H(3)*X(2)-H(2)*X(3)
  17. Y(2) = (H(3)+H(1))*X(2)-H(1)*X(3)-H(3)*X(1)
  18. Y(3) = (H(1)+H(2))*X(3)-H(2)*X(1)-H(1)*X(2)
  19. DO L = 4,6
  20. Y(L) = H(L)*X(L)
  21. END DO
  22. RETURN
  23. END
  24.  
  25.  
  26.  

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