Esope source files

   1 : C MULTP2    SOURCE    CHAT      05/01/13    01:55:35     5004           
   2 : 
   3 :       SUBROUTINE MULTP2(P, Q, R)
   4 : C
   5 : C     Procedure de multiplication de 2 polynomes de 2 variables
   6 : C     Entree : P(9) : premier polynome : les coefficients sont
   7 : C                       ranges comme suit :
   8 : C               indice : 1 2  3    4    5    6      7       8        9
   9 : C               terme :  1 T2 T1 T1*T2 T2^2 T1^2 T1*T2^2 T1^2*T2 T1^2*T2^2
  10 : C              Q(9) : deuxieme polynome
  11 : C     Sortie : R(9) : polynome resultat de la multiplication
  12 : C
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8(A-H,O-Z)
  15 :       DIMENSION P(9), Q(9), R(9), INXM(4, 4)
  16 :       DATA INXM/ 1,   2,   3,   4,
  17 :      1           2,   5,   4,   7,
  18 :      2           3,   4,   6,   8,
  19 :      3           4,   7,   8,   9/
  20 : C
  21 : C     Initialisation du polynome resultat
  22 : C
  23 :       DO 50 I = 1, 9
  24 :   50     R(I) = 0.D0
  25 : C
  26 : C     Multiplication de P et Q
  27 : C
  28 :       DO 60 I = 1, 4
  29 :          DO 60 J = 1, 4
  30 :   60        R(INXM(I,J)) = R(INXM(I,J)) + P(I)*Q(J)
  31 :       RETURN
  32 :       END
  33 :  
  34 :  
  35 :