[u/mrichter/AliRoot.git] / PDF / spdf / aurpi1.F

*
* $Id$
*
* $Log$
* Revision 1.4  2000/09/18 10:02:35  hristov
* Makefile added to PDF8
*
* Revision 1.1.1.2  1996/10/30 08:27:36  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:03  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
C
      SUBROUTINE AURPI1(I,NDRV,X,S,ANS)
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION
#endif
#if defined(CERNLIB_SINGLE)
      REAL
#endif
     +       CALCPIO(8,20,25)
      COMMON/W5051I4/CALCPIO
      REAL   F1(25),F2(25)
      DATA DELTA/.10/
      ANS=0.
      IF(X.GT.0.9985) RETURN
      IF(I.EQ.3.AND.X.GT.0.95) RETURN
      IF(I.EQ.8.AND.X.GT.0.95) RETURN
      IS=S/DELTA+1
      IS1=IS+1
      DO 1 L=1,25
      KL=L+NDRV*25
      F1(L)=CALCPIO(I,IS,KL)
      F2(L)=CALCPIO(I,IS1,KL)
    1 CONTINUE
      A1=AUGETFV(X,F1)
      A2=AUGETFV(X,F2)
      S1=(IS-1)*DELTA
      S2=S1+DELTA
      ANS=A1*(S-S2)/(S1-S2)+A2*(S-S1)/(S2-S1)
      RETURN
      END
Commit	Line	Data
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:36 cernlib
	9	* Version 7.04
	10	*
	11	* Revision 1.1.1.1 1996/04/12 15:29:03 plothow
	12	* Version 7.01
	13	*
	14	*
	15	#include "pdf/pilot.h"
	16	C
	17	SUBROUTINE AURPI1(I,NDRV,X,S,ANS)
7ef50f50	18	#if defined(CERNLIB_DOUBLE)
	19	DOUBLE PRECISION
	20	#endif
	21	#if defined(CERNLIB_SINGLE)
	22	REAL
	23	#endif
21886bb6	24	+ CALCPIO(8,20,25)
	25	COMMON/W5051I4/CALCPIO
	26	REAL F1(25),F2(25)
	27	DATA DELTA/.10/
	28	ANS=0.
	29	IF(X.GT.0.9985) RETURN
	30	IF(I.EQ.3.AND.X.GT.0.95) RETURN
	31	IF(I.EQ.8.AND.X.GT.0.95) RETURN
	32	IS=S/DELTA+1
	33	IS1=IS+1
	34	DO 1 L=1,25
	35	KL=L+NDRV*25
	36	F1(L)=CALCPIO(I,IS,KL)
	37	F2(L)=CALCPIO(I,IS1,KL)
	38	1 CONTINUE
	39	A1=AUGETFV(X,F1)
	40	A2=AUGETFV(X,F2)
	41	S1=(IS-1)*DELTA
	42	S2=S1+DELTA
	43	ANS=A1(S-S2)/(S1-S2)+A2(S-S1)/(S2-S1)
	44	RETURN
	45	END