[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / c / sinint64.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:06  mclareni
* Mathlib gen
*
*
#include "gen/pilot.h"
#if !defined(CERNLIB_DOUBLE)
      FUNCTION RSININ(X)
#endif
#if defined(CERNLIB_DOUBLE)
      FUNCTION DSININ(X)
#endif
#include "gen/imp64.inc"
      CHARACTER NAME*(*)
      CHARACTER*80 ERRTXT
#if !defined(CERNLIB_DOUBLE)
      PARAMETER (NAME = 'RCOSIN')
#endif
#if defined(CERNLIB_DOUBLE)
      PARAMETER (NAME = 'RCOSIN/DCOSIN')
#endif
      DIMENSION S(0:15),C(0:15),P(0:28),Q(0:24)

      PARAMETER (Z1 = 1, R8 = Z1/8, R32 = Z1/32)

      PARAMETER (PI = 3.14159 26535 89793 24D0)
      PARAMETER (CE = 0.57721 56649 01532 86D0)
      PARAMETER (PIH = PI/2)

      DATA S( 0) /+1.95222 09759 53071 08D0/
      DATA S( 1) /-0.68840 42321 25715 44D0/
      DATA S( 2) /+0.45518 55132 25584 84D0/
      DATA S( 3) /-0.18045 71236 83877 85D0/
      DATA S( 4) /+0.04104 22133 75859 24D0/
      DATA S( 5) /-0.00595 86169 55588 85D0/
      DATA S( 6) /+0.00060 01427 41414 43D0/
      DATA S( 7) /-0.00004 44708 32910 75D0/
      DATA S( 8) /+0.00000 25300 78230 75D0/
      DATA S( 9) /-0.00000 01141 30759 30D0/
      DATA S(10) /+0.00000 00041 85783 94D0/
      DATA S(11) /-0.00000 00001 27347 06D0/
      DATA S(12) /+0.00000 00000 03267 36D0/
      DATA S(13) /-0.00000 00000 00071 68D0/
      DATA S(14) /+0.00000 00000 00001 36D0/
      DATA S(15) /-0.00000 00000 00000 02D0/

      DATA C( 0) /+1.94054 91464 83554 93D0/
      DATA C( 1) /+0.94134 09132 86521 34D0/
      DATA C( 2) /-0.57984 50342 92992 76D0/
      DATA C( 3) /+0.30915 72011 15927 13D0/
      DATA C( 4) /-0.09161 01792 20771 34D0/
      DATA C( 5) /+0.01644 37407 51546 25D0/
      DATA C( 6) /-0.00197 13091 95216 41D0/
      DATA C( 7) /+0.00016 92538 85083 50D0/
      DATA C( 8) /-0.00001 09393 29573 11D0/
      DATA C( 9) /+0.00000 05522 38574 84D0/
      DATA C(10) /-0.00000 00223 99493 31D0/
      DATA C(11) /+0.00000 00007 46533 25D0/
      DATA C(12) /-0.00000 00000 20818 33D0/
      DATA C(13) /+0.00000 00000 00493 12D0/
      DATA C(14) /-0.00000 00000 00010 05D0/
      DATA C(15) /+0.00000 00000 00000 18D0/

      DATA P( 0) /+0.96074 78397 52035 96D0/
      DATA P( 1) /-0.03711 38962 12398 06D0/
      DATA P( 2) /+0.00194 14398 88991 90D0/
      DATA P( 3) /-0.00017 16598 84251 47D0/
      DATA P( 4) /+0.00002 11263 77532 31D0/
      DATA P( 5) /-0.00000 32716 32567 12D0/
      DATA P( 6) /+0.00000 06006 92116 15D0/
      DATA P( 7) /-0.00000 01258 67944 03D0/
      DATA P( 8) /+0.00000 00293 25634 58D0/
      DATA P( 9) /-0.00000 00074 56959 21D0/
      DATA P(10) /+0.00000 00020 41054 78D0/
      DATA P(11) /-0.00000 00005 95022 30D0/
      DATA P(12) /+0.00000 00001 83229 67D0/
      DATA P(13) /-0.00000 00000 59205 06D0/
      DATA P(14) /+0.00000 00000 19965 17D0/
      DATA P(15) /-0.00000 00000 06995 11D0/
      DATA P(16) /+0.00000 00000 02536 86D0/
      DATA P(17) /-0.00000 00000 00949 29D0/
      DATA P(18) /+0.00000 00000 00365 52D0/
      DATA P(19) /-0.00000 00000 00144 49D0/
      DATA P(20) /+0.00000 00000 00058 51D0/
      DATA P(21) /-0.00000 00000 00024 23D0/
      DATA P(22) /+0.00000 00000 00010 25D0/
      DATA P(23) /-0.00000 00000 00004 42D0/
      DATA P(24) /+0.00000 00000 00001 94D0/
      DATA P(25) /-0.00000 00000 00000 87D0/
      DATA P(26) /+0.00000 00000 00000 39D0/
      DATA P(27) /-0.00000 00000 00000 18D0/
      DATA P(28) /+0.00000 00000 00000 08D0/

      DATA Q( 0) /+0.98604 06569 62382 60D0/
      DATA Q( 1) /-0.01347 17382 08295 21D0/
      DATA Q( 2) /+0.00045 32928 41165 23D0/
      DATA Q( 3) /-0.00003 06728 86516 55D0/
      DATA Q( 4) /+0.00000 31319 91976 01D0/
      DATA Q( 5) /-0.00000 04211 01964 96D0/
      DATA Q( 6) /+0.00000 00690 72448 30D0/
      DATA Q( 7) /-0.00000 00131 83212 90D0/
      DATA Q( 8) /+0.00000 00028 36974 33D0/
      DATA Q( 9) /-0.00000 00006 73292 34D0/
      DATA Q(10) /+0.00000 00001 73396 87D0/
      DATA Q(11) /-0.00000 00000 47869 39D0/
      DATA Q(12) /+0.00000 00000 14032 35D0/
      DATA Q(13) /-0.00000 00000 04334 96D0/
      DATA Q(14) /+0.00000 00000 01402 73D0/
      DATA Q(15) /-0.00000 00000 00473 06D0/
      DATA Q(16) /+0.00000 00000 00165 58D0/
      DATA Q(17) /-0.00000 00000 00059 94D0/
      DATA Q(18) /+0.00000 00000 00022 37D0/
      DATA Q(19) /-0.00000 00000 00008 59D0/
      DATA Q(20) /+0.00000 00000 00003 38D0/
      DATA Q(21) /-0.00000 00000 00001 36D0/
      DATA Q(22) /+0.00000 00000 00000 56D0/
      DATA Q(23) /-0.00000 00000 00000 24D0/
      DATA Q(24) /+0.00000 00000 00000 10D0/

      IF(ABS(X) .LE. 8) THEN
       Y=R8*X
       H=2*Y**2-1
       ALFA=H+H
       B1=0
       B2=0
       DO 1 I = 15,0,-1
       B0=S(I)+ALFA*B1-B2
       B2=B1
    1  B1=B0
       H=Y*(B0-B2)
      ELSE
       R=1/X
       H=128*R**2-1
       ALFA=H+H
       B1=0
       B2=0
       DO 2 I = 28,0,-1
       B0=P(I)+ALFA*B1-B2
       B2=B1
    2  B1=B0
       PP=B0-H*B2
       B1=0
       B2=0
       DO 3 I = 24,0,-1
       B0=Q(I)+ALFA*B1-B2
       B2=B1
    3  B1=B0
       H=SIGN(PIH,X)-R*(R*PP*SIN(X)+(B0-H*B2)*COS(X))
      END IF
      GO TO 9

#if !defined(CERNLIB_DOUBLE)
      ENTRY RCOSIN(X)
#endif
#if defined(CERNLIB_DOUBLE)
      ENTRY DCOSIN(X)
#endif

      IF(X .EQ. 0) THEN
       H=0
       CALL MTLPRT(NAME,'C336.1','ARGUMENT X = 0')
      ELSEIF(ABS(X) .LE. 8) THEN
       H=R32*X**2-1
       ALFA=H+H
       B1=0
       B2=0
       DO 4 I = 15,0,-1
       B0=C(I)+ALFA*B1-B2
       B2=B1
    4  B1=B0
       H=CE+LOG(ABS(X))-B0+H*B2
      ELSE
       R=1/X
       H=128*R**2-1
       ALFA=H+H
       B1=0
       B2=0
       DO 5 I = 28,0,-1
       B0=P(I)+ALFA*B1-B2
       B2=B1
    5  B1=B0
       PP=B0-H*B2
       B1=0
       B2=0
       DO 6 I = 24,0,-1
       B0=Q(I)+ALFA*B1-B2
       B2=B1
    6  B1=B0
       H=R*((B0-H*B2)*SIN(X)-R*PP*COS(X))
      END IF
#if !defined(CERNLIB_DOUBLE)
    9 RSININ=H
#endif
#if defined(CERNLIB_DOUBLE)
    9 DSININ=H
#endif
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/04/01 15:02:06 mclareni
	6	* Mathlib gen
	7	*
	8	*
	9	#include "gen/pilot.h"
	10	#if !defined(CERNLIB_DOUBLE)
	11	FUNCTION RSININ(X)
	12	#endif
	13	#if defined(CERNLIB_DOUBLE)
	14	FUNCTION DSININ(X)
	15	#endif
	16	#include "gen/imp64.inc"
	17	CHARACTER NAME()
	18	CHARACTER*80 ERRTXT
	19	#if !defined(CERNLIB_DOUBLE)
	20	PARAMETER (NAME = 'RCOSIN')
	21	#endif
	22	#if defined(CERNLIB_DOUBLE)
	23	PARAMETER (NAME = 'RCOSIN/DCOSIN')
	24	#endif
	25	DIMENSION S(0:15),C(0:15),P(0:28),Q(0:24)
	26
	27	PARAMETER (Z1 = 1, R8 = Z1/8, R32 = Z1/32)
	28
	29	PARAMETER (PI = 3.14159 26535 89793 24D0)
	30	PARAMETER (CE = 0.57721 56649 01532 86D0)
	31	PARAMETER (PIH = PI/2)
	32
	33	DATA S( 0) /+1.95222 09759 53071 08D0/
	34	DATA S( 1) /-0.68840 42321 25715 44D0/
	35	DATA S( 2) /+0.45518 55132 25584 84D0/
	36	DATA S( 3) /-0.18045 71236 83877 85D0/
	37	DATA S( 4) /+0.04104 22133 75859 24D0/
	38	DATA S( 5) /-0.00595 86169 55588 85D0/
	39	DATA S( 6) /+0.00060 01427 41414 43D0/
	40	DATA S( 7) /-0.00004 44708 32910 75D0/
	41	DATA S( 8) /+0.00000 25300 78230 75D0/
	42	DATA S( 9) /-0.00000 01141 30759 30D0/
	43	DATA S(10) /+0.00000 00041 85783 94D0/
	44	DATA S(11) /-0.00000 00001 27347 06D0/
	45	DATA S(12) /+0.00000 00000 03267 36D0/
	46	DATA S(13) /-0.00000 00000 00071 68D0/
	47	DATA S(14) /+0.00000 00000 00001 36D0/
	48	DATA S(15) /-0.00000 00000 00000 02D0/
	49
	50	DATA C( 0) /+1.94054 91464 83554 93D0/
	51	DATA C( 1) /+0.94134 09132 86521 34D0/
	52	DATA C( 2) /-0.57984 50342 92992 76D0/
	53	DATA C( 3) /+0.30915 72011 15927 13D0/
	54	DATA C( 4) /-0.09161 01792 20771 34D0/
	55	DATA C( 5) /+0.01644 37407 51546 25D0/
	56	DATA C( 6) /-0.00197 13091 95216 41D0/
	57	DATA C( 7) /+0.00016 92538 85083 50D0/
	58	DATA C( 8) /-0.00001 09393 29573 11D0/
	59	DATA C( 9) /+0.00000 05522 38574 84D0/
	60	DATA C(10) /-0.00000 00223 99493 31D0/
	61	DATA C(11) /+0.00000 00007 46533 25D0/
	62	DATA C(12) /-0.00000 00000 20818 33D0/
	63	DATA C(13) /+0.00000 00000 00493 12D0/
	64	DATA C(14) /-0.00000 00000 00010 05D0/
65	DATA C(15) /+0.00000 00000 00000 18D0/
66
67	DATA P( 0) /+0.96074 78397 52035 96D0/
68	DATA P( 1) /-0.03711 38962 12398 06D0/
69	DATA P( 2) /+0.00194 14398 88991 90D0/
70	DATA P( 3) /-0.00017 16598 84251 47D0/
71	DATA P( 4) /+0.00002 11263 77532 31D0/
72	DATA P( 5) /-0.00000 32716 32567 12D0/
73	DATA P( 6) /+0.00000 06006 92116 15D0/
74	DATA P( 7) /-0.00000 01258 67944 03D0/
75	DATA P( 8) /+0.00000 00293 25634 58D0/
76	DATA P( 9) /-0.00000 00074 56959 21D0/
77	DATA P(10) /+0.00000 00020 41054 78D0/
78	DATA P(11) /-0.00000 00005 95022 30D0/
79	DATA P(12) /+0.00000 00001 83229 67D0/
80	DATA P(13) /-0.00000 00000 59205 06D0/
81	DATA P(14) /+0.00000 00000 19965 17D0/
82	DATA P(15) /-0.00000 00000 06995 11D0/
83	DATA P(16) /+0.00000 00000 02536 86D0/
84	DATA P(17) /-0.00000 00000 00949 29D0/
85	DATA P(18) /+0.00000 00000 00365 52D0/
86	DATA P(19) /-0.00000 00000 00144 49D0/
87	DATA P(20) /+0.00000 00000 00058 51D0/
88	DATA P(21) /-0.00000 00000 00024 23D0/
89	DATA P(22) /+0.00000 00000 00010 25D0/
90	DATA P(23) /-0.00000 00000 00004 42D0/
91	DATA P(24) /+0.00000 00000 00001 94D0/
92	DATA P(25) /-0.00000 00000 00000 87D0/
93	DATA P(26) /+0.00000 00000 00000 39D0/
94	DATA P(27) /-0.00000 00000 00000 18D0/
95	DATA P(28) /+0.00000 00000 00000 08D0/
96
97	DATA Q( 0) /+0.98604 06569 62382 60D0/
98	DATA Q( 1) /-0.01347 17382 08295 21D0/
99	DATA Q( 2) /+0.00045 32928 41165 23D0/
100	DATA Q( 3) /-0.00003 06728 86516 55D0/
101	DATA Q( 4) /+0.00000 31319 91976 01D0/
102	DATA Q( 5) /-0.00000 04211 01964 96D0/
103	DATA Q( 6) /+0.00000 00690 72448 30D0/
104	DATA Q( 7) /-0.00000 00131 83212 90D0/
105	DATA Q( 8) /+0.00000 00028 36974 33D0/
106	DATA Q( 9) /-0.00000 00006 73292 34D0/
107	DATA Q(10) /+0.00000 00001 73396 87D0/
108	DATA Q(11) /-0.00000 00000 47869 39D0/
109	DATA Q(12) /+0.00000 00000 14032 35D0/
110	DATA Q(13) /-0.00000 00000 04334 96D0/
111	DATA Q(14) /+0.00000 00000 01402 73D0/
112	DATA Q(15) /-0.00000 00000 00473 06D0/
113	DATA Q(16) /+0.00000 00000 00165 58D0/
114	DATA Q(17) /-0.00000 00000 00059 94D0/
115	DATA Q(18) /+0.00000 00000 00022 37D0/
116	DATA Q(19) /-0.00000 00000 00008 59D0/
117	DATA Q(20) /+0.00000 00000 00003 38D0/
118	DATA Q(21) /-0.00000 00000 00001 36D0/
119	DATA Q(22) /+0.00000 00000 00000 56D0/
120	DATA Q(23) /-0.00000 00000 00000 24D0/
121	DATA Q(24) /+0.00000 00000 00000 10D0/
122
123	IF(ABS(X) .LE. 8) THEN
124	Y=R8*X
125	H=2Y*2-1
126	ALFA=H+H
127	B1=0
128	B2=0
129	DO 1 I = 15,0,-1
130	B0=S(I)+ALFA*B1-B2
131	B2=B1
132	1 B1=B0
133	H=Y*(B0-B2)
134	ELSE
135	R=1/X
136	H=128R*2-1
137	ALFA=H+H
138	B1=0
139	B2=0
140	DO 2 I = 28,0,-1
141	B0=P(I)+ALFA*B1-B2
142	B2=B1
143	2 B1=B0
144	PP=B0-H*B2
145	B1=0
146	B2=0
147	DO 3 I = 24,0,-1
148	B0=Q(I)+ALFA*B1-B2
149	B2=B1
150	3 B1=B0
151	H=SIGN(PIH,X)-R(RPPSIN(X)+(B0-HB2)*COS(X))
152	END IF
153	GO TO 9
154
155	#if !defined(CERNLIB_DOUBLE)
156	ENTRY RCOSIN(X)
157	#endif
158	#if defined(CERNLIB_DOUBLE)
159	ENTRY DCOSIN(X)
160	#endif
161
162	IF(X .EQ. 0) THEN
163	H=0
164	CALL MTLPRT(NAME,'C336.1','ARGUMENT X = 0')
165	ELSEIF(ABS(X) .LE. 8) THEN
166	H=R32X*2-1
167	ALFA=H+H
168	B1=0
169	B2=0
170	DO 4 I = 15,0,-1
171	B0=C(I)+ALFA*B1-B2
172	B2=B1
173	4 B1=B0
174	H=CE+LOG(ABS(X))-B0+H*B2
175	ELSE
176	R=1/X
177	H=128R*2-1
178	ALFA=H+H
179	B1=0
180	B2=0
181	DO 5 I = 28,0,-1
182	B0=P(I)+ALFA*B1-B2
183	B2=B1
184	5 B1=B0
185	PP=B0-H*B2
186	B1=0
187	B2=0
188	DO 6 I = 24,0,-1
189	B0=Q(I)+ALFA*B1-B2
190	B2=B1
191	6 B1=B0
192	H=R((B0-HB2)SIN(X)-RPP*COS(X))
193	END IF
194	#if !defined(CERNLIB_DOUBLE)
195	9 RSININ=H
196	#endif
197	#if defined(CERNLIB_DOUBLE)
198	9 DSININ=H
199	#endif
200	RETURN
201	END