Esope source files

   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 :