Sources Esope | Cast3M

multp2
C MULTP2    SOURCE    CHAT      05/01/13    01:55:35     5004
 
      SUBROUTINE MULTP2(P, Q, R)
C
C     Procedure de multiplication de 2 polynomes de 2 variables
C     Entree : P(9) : premier polynome : les coefficients sont
C                       ranges comme suit :
C               indice : 1 2  3    4    5    6      7       8        9
C               terme :  1 T2 T1 T1*T2 T2^2 T1^2 T1*T2^2 T1^2*T2 T1^2*T2^2
C              Q(9) : deuxieme polynome
C     Sortie : R(9) : polynome resultat de la multiplication
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION P(9), Q(9), R(9), INXM(4, 4)
      DATA INXM/ 1,   2,   3,   4,
     1           2,   5,   4,   7,
     2           3,   4,   6,   8,
     3           4,   7,   8,   9/
C
C     Initialisation du polynome resultat
C
      DO 50 I = 1, 9
  50     R(I) = 0.D0
C
C     Multiplication de P et Q
C
      DO 60 I = 1, 4
         DO 60 J = 1, 4
  60        R(INXM(I,J)) = R(INXM(I,J)) + P(I)*Q(J)
      RETURN
      END