Sources Esope | Cast3M

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