Esope source files

   1 : C DERIP2    SOURCE    CHAT      05/01/12    22:42:57     5004           
   2 :       SUBROUTINE DERIP2(P, IVAR, DP)
   3 : C
   4 : C     Procedure de multiplication de 2 polynomes de 2 variables
   5 : C     Entree : P(9) : polynome a deriver : les coefficients sont
   6 : C                       ranges comme suit :
   7 : C               indice : 1 2  3    4    5    6      7       8        9
   8 : C               terme :  1 T2 T1 T1*T2 T2^2 T1^2 T1*T2^2 T1^2*T2 T1^2*T2^2
   9 : C              IVAR : numero de la variable par rapport a laquelle il 
  10 : C                       faut deriver
  11 : C     Sortie : DP(9) : polynome resultat de la derivation
  12 : C
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8(A-H,O-Z)
  15 :       DIMENSION P(9), DP(9), INXD1(9), INXD2(9), C1(9), C2(9)
  16 :       DATA INXD1/ 1, 1, 1, 2, 1, 3, 5, 4, 7/
  17 :       DATA C1/0., 0., 1., 1., 0., 2., 1., 2., 2./
  18 :       DATA INXD2/ 1, 1, 1, 3, 2, 1, 4, 6, 8/
  19 :       DATA C2/0., 1., 0., 1., 2., 0., 2., 1., 2./
  20 : C
  21 : C     Initialisation du polynome resultat
  22 : C
  23 :       DO 70 I = 1, 9
  24 :   70     DP(I) = 0.D0
  25 : C
  26 : C     Derivation de P
  27 : C
  28 :       IF (IVAR .EQ. 1) THEN
  29 :          DO 80 I = 2, 9
  30 :   80        DP(INXD1(I)) = DP(INXD1(I)) + C1(I)*P(I)
  31 :       ELSE
  32 :          IF (IVAR .EQ. 2) THEN
  33 :             DO 90 I = 2, 9
  34 :   90           DP(INXD2(I)) = DP(INXD2(I)) + C2(I)*P(I)
  35 : 
  36 :          ENDIF
  37 :       ENDIF
  38 :       RETURN
  39 :       END
  40 :  
  41 :