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