| 21886bb6 |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame