Dummy subroutines to avoid files with no code in
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / c / gamma128.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1999/05/18 15:55:34  fca
6 * AliRoot sources
7 *
8 * Revision 1.1.1.1  1996/04/01 15:01:54  mclareni
9 * Mathlib gen
10 *
11 *
12 #include "gen/pilot.h"
13 #if defined(CERNLIB_QUAD)
14 #if !defined(CERNLIB_DOUBLE)
15       FUNCTION DGAMMA(X)
16 #endif
17 #if defined(CERNLIB_DOUBLE)
18       FUNCTION QGAMMA(X)
19 #endif
20
21 #include "gen/imp128.inc"
22       CHARACTER NAME*(*)
23       CHARACTER*80 ERRTXT
24 #if !defined(CERNLIB_DOUBLE)
25       PARAMETER (NAME = 'DGAMMA')
26       DIMENSION C(0:26)
27
28       DATA C( 0) / 3.65738 77250 83382 43849 88068 52039D0/
29       DATA C( 1) / 1.95754 34566 61268 26928 33742 32877D0/
30       DATA C( 2) / 0.33829 71138 26160 38915 58510 72725D0/
31       DATA C( 3) / 0.04208 95127 65575 49198 51083 97028D0/
32       DATA C( 4) / 0.00428 76504 82129 08770 04289 08116D0/
33       DATA C( 5) / 0.00036 52121 69294 61767 02198 21538D0/
34       DATA C( 6) / 0.00002 74006 42226 42200 27170 66599D0/
35       DATA C( 7) / 0.00000 18124 02333 65124 44603 04653D0/
36       DATA C( 8) / 0.00000 01096 57758 65997 06993 06003D0/
37       DATA C( 9) / 0.00000 00059 87184 04552 00046 94981D0/
38       DATA C(10) / 0.00000 00003 07690 80535 24777 70970D0/
39       DATA C(11) / 0.00000 00000 14317 93029 61915 76369D0/
40       DATA C(12) / 0.00000 00000 00651 08773 34803 70077D0/
41       DATA C(13) / 0.00000 00000 00025 95849 89822 27998D0/
42       DATA C(14) / 0.00000 00000 00001 10789 38922 59522D0/
43       DATA C(15) / 0.00000 00000 00000 03547 43620 16354D0/
44       DATA C(16) / 0.00000 00000 00000 00168 86075 02335D0/
45       DATA C(17) / 0.00000 00000 00000 00002 73543 57562D0/
46       DATA C(18) / 0.00000 00000 00000 00000 30297 74054D0/
47       DATA C(19) /-0.00000 00000 00000 00000 00571 22030D0/
48       DATA C(20) / 0.00000 00000 00000 00000 00090 76859D0/
49       DATA C(21) /-0.00000 00000 00000 00000 00005 04522D0/
50       DATA C(22) / 0.00000 00000 00000 00000 00000 40496D0/
51       DATA C(23) /-0.00000 00000 00000 00000 00000 02792D0/
52       DATA C(24) / 0.00000 00000 00000 00000 00000 00204D0/
53       DATA C(25) /-0.00000 00000 00000 00000 00000 00015D0/
54       DATA C(26) / 0.00000 00000 00000 00000 00000 00001D0/
55
56 #endif
57 #if defined(CERNLIB_DOUBLE)
58       PARAMETER (NAME = 'QGAMMA')
59       DIMENSION C(0:27)
60
61       DATA C( 0) / 3.65738 77250 83382 43849 88068 52039 34Q0/
62       DATA C( 1) / 1.95754 34566 61268 26928 33742 32876 84Q0/
63       DATA C( 2) / 0.33829 71138 26160 38915 58510 72725 40Q0/
64       DATA C( 3) / 0.04208 95127 65575 49198 51083 97027 90Q0/
65       DATA C( 4) / 0.00428 76504 82129 08770 04289 08116 23Q0/
66       DATA C( 5) / 0.00036 52121 69294 61767 02198 21538 33Q0/
67       DATA C( 6) / 0.00002 74006 42226 42200 27170 66598 98Q0/
68       DATA C( 7) / 0.00000 18124 02333 65124 44603 04653 06Q0/
69       DATA C( 8) / 0.00000 01096 57758 65997 06993 06002 52Q0/
70       DATA C( 9) / 0.00000 00059 87184 04552 00046 94980 50Q0/
71       DATA C(10) / 0.00000 00003 07690 80535 24777 70969 71Q0/
72       DATA C(11) / 0.00000 00000 14317 93029 61915 76369 12Q0/
73       DATA C(12) / 0.00000 00000 00651 08773 34803 70076 65Q0/
74       DATA C(13) / 0.00000 00000 00025 95849 89822 27997 94Q0/
75       DATA C(14) / 0.00000 00000 00001 10789 38922 59521 65Q0/
76       DATA C(15) / 0.00000 00000 00000 03547 43620 16353 80Q0/
77       DATA C(16) / 0.00000 00000 00000 00168 86075 02335 19Q0/
78       DATA C(17) / 0.00000 00000 00000 00002 73543 57561 54Q0/
79       DATA C(18) / 0.00000 00000 00000 00000 30297 74054 44Q0/
80       DATA C(19) /-0.00000 00000 00000 00000 00571 22030 32Q0/
81       DATA C(20) / 0.00000 00000 00000 00000 00090 76858 53Q0/
82       DATA C(21) /-0.00000 00000 00000 00000 00005 04522 44Q0/
83       DATA C(22) / 0.00000 00000 00000 00000 00000 40496 46Q0/
84       DATA C(23) /-0.00000 00000 00000 00000 00000 02791 84Q0/
85       DATA C(24) / 0.00000 00000 00000 00000 00000 00203 70Q0/
86       DATA C(25) /-0.00000 00000 00000 00000 00000 00014 54Q0/
87       DATA C(26) / 0.00000 00000 00000 00000 00000 00001 06Q0/
88       DATA C(27) /-0.00000 00000 00000 00000 00000 00000 07Q0/
89 #endif
90       U=X
91       IF(U .LE. 0) THEN
92        H=0
93        WRITE(ERRTXT,101) U
94        CALL MTLPRT(NAME,'C302.1',ERRTXT)
95        GO TO 9
96       ENDIF
97       F=1
98       IF(U .LT. 3) THEN
99        DO 1 I = 1,INT(4-U)
100        F=F/U
101     1  U=U+1
102       ELSE
103        DO 2 I = 1,INT(U-3)
104        U=U-1
105     2  F=F*U
106       END IF
107       H=U+U-7
108       ALFA=H+H
109       B1=0
110       B2=0
111 #if !defined(CERNLIB_DOUBLE)
112       DO 3 I = 26,0,-1
113 #endif
114 #if defined(CERNLIB_DOUBLE)
115       DO 3 I = 27,0,-1
116 #endif
117       B0=C(I)+ALFA*B1-B2
118       B2=B1
119     3 B1=B0
120       H=F*(B0-H*B2)
121 #if !defined(CERNLIB_DOUBLE)
122     9 DGAMMA=H
123 #endif
124 #if defined(CERNLIB_DOUBLE)
125     9 QGAMMA=H
126 #endif
127       RETURN
128   101 FORMAT('ARGUMENT IS NEGATIVE = ',1P,E15.1)
129       END
130 #else
131       SUBROUTINE GAMMA128_DUMMY
132       END
133 #endif