Sources Esope | Cast3M

derip2
C DERIP2    SOURCE    CHAT      05/01/12    22:42:57     5004
      SUBROUTINE DERIP2(P, IVAR, DP)
C
C     Procedure de multiplication de 2 polynomes de 2 variables
C     Entree : P(9) : polynome a deriver : 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              IVAR : numero de la variable par rapport a laquelle il
C                       faut deriver
C     Sortie : DP(9) : polynome resultat de la derivation
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION P(9), DP(9), INXD1(9), INXD2(9), C1(9), C2(9)
      DATA INXD1/ 1, 1, 1, 2, 1, 3, 5, 4, 7/
      DATA C1/0., 0., 1., 1., 0., 2., 1., 2., 2./
      DATA INXD2/ 1, 1, 1, 3, 2, 1, 4, 6, 8/
      DATA C2/0., 1., 0., 1., 2., 0., 2., 1., 2./
C
C     Initialisation du polynome resultat
C
      DO 70 I = 1, 9
  70     DP(I) = 0.D0
C
C     Derivation de P
C
      IF (IVAR .EQ. 1) THEN
         DO 80 I = 2, 9
  80        DP(INXD1(I)) = DP(INXD1(I)) + C1(I)*P(I)
      ELSE
         IF (IVAR .EQ. 2) THEN
            DO 90 I = 2, 9
  90           DP(INXD2(I)) = DP(INXD2(I)) + C2(I)*P(I)
 
         ENDIF
      ENDIF
      RETURN
      END