PDF/spdf/afgetfv.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:27:31  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:28:59  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13        FUNCTION AFGETFV(X,FVL)
  14 C  NOUVEAU PROGRAMME D'INTERPOLATION UTILISANT UNE ROUTINE DE MATH. RECIPES
  15        DIMENSION FVL(32)
  16 #include "pdf/expdp.inc"
  17      +                XI(32),WI(32),XX(33)
  18        COMMON/W5051I9/XI,WI,XX,NTERMS
  19        DIMENSION A(4),B(4)
  20        N=4
  21        EPS=1.E-7
  22        XAM=XX(1)-EPS
  23        XAP=XX(1)+EPS
  24 C      IF(X.LT.XAM) PRINT*,' X = ',X
  25        IF(X.GT.XAM.AND.X.LT.XAP) GO TO 50
  26        GO TO 80
  27    50  Y=FVL(1)
  28        GO TO 77
  29    80  IF(X.LT.XX(2)) GO TO 51
  30        IF(X.GT.XX(30)) GO TO 61
  31        DO 1 I=3,30
  32        IF(X.GT.XX(I)) GO TO 1
  33        A(1)=XX(I-2)
  34        A(2)=XX(I-1)
  35        A(3)=XX(I)
  36        A(4)=XX(I+1)
  37        B(1)=FVL(I-2)
  38        B(2)=FVL(I-1)
  39        B(3)=FVL(I)
  40        B(4)=FVL(I+1)
  41        GO TO 70
  42    1   CONTINUE
  43   61   A(1)=XX(29)
  44        A(2)=XX(30)
  45        A(3)=XX(31)
  46        A(4)=XX(32)
  47        B(1)=FVL(29)
  48        B(2)=FVL(30)
  49        B(3)=FVL(31)
  50        B(4)=FVL(32)
  51        GO  TO 70
  52   51   A(1)=XX(1)
  53        A(2)=XX(2)
  54        A(3)=XX(3)
  55        A(4)=XX(4)
  56        B(1)=FVL(1)
  57        B(2)=FVL(2)
  58        B(3)=FVL(3)
  59        B(4)=FVL(4)
  60 C 70   IF(X.GT..2.AND.X.LT..8) THEN
  61 C            CALL AFPOLIN(A,B,N,X,Y,DY)
  62 C      ELSE
  63 C            CALL AFRATIN(A,B,N,X,Y,DY)
  64 C      ENDIF
  65   70   CONTINUE
  66              CALL AFPOLIN(A,B,N,X,Y,DY)
  67   77   AFGETFV=Y
  68        RETURN
  69        END