1	/* Function tanf vectorized with SSE4.
2	Copyright (C) 2021-2024 Free Software Foundation, Inc.
3	This file is part of the GNU C Library.
4
5	The GNU C Library is free software; you can redistribute it and/or
6	modify it under the terms of the GNU Lesser General Public
7	License as published by the Free Software Foundation; either
8	version 2.1 of the License, or (at your option) any later version.
9
10	The GNU C Library is distributed in the hope that it will be useful,
11	but WITHOUT ANY WARRANTY; without even the implied warranty of
12	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13	Lesser General Public License for more details.
14
15	You should have received a copy of the GNU Lesser General Public
16	License along with the GNU C Library; if not, see
17	https://www.gnu.org/licenses/. */
18
19	/*
20	* ALGORITHM DESCRIPTION:
21	*
22	* 1) Range reduction to [-Pi/4; +Pi/4] interval
23	* a) Grab sign from source argument and save it.
24	* b) Remove sign using AND 0x7fffffff operation
25	* c) Getting octant Y by 2/Pi multiplication
26	* d) Add "Right Shifter" (0x4B000000) value
27	* e) Treat obtained value as integer for destination sign setting.
28	* Shift first bit of this value to the last (sign) position (S << 31)
29	* f) Change destination sign if source sign is negative
30	* using XOR operation.
31	* g) Subtract "Right Shifter" (0x4B000000) value
32	* h) Subtract Y*(PI/2) from X argument, where PI/2 divided to 4 parts:
33	* X = X - Y*PI1 - Y*PI2 - Y*PI3 - Y*PI4;
34	* 2) Rational polynomial approximation ( at [-Pi/4; +Pi/4] interval)
35	* a) Calculate X^2 = X * X
36	* b) Calculate 2 polynomials:
37	* P = X * (P0 + X^2 * P1);
38	* Q = Q0 + X^2 * (Q1 + x^2 * Q2);
39	* c) Swap P and Q if first bit of obtained value after
40	* Right Shifting is set to 1. Using And, Andnot & Or operations.
41	* d) Divide R = P / Q;
42	* 3) Destination sign setting
43	* a) Set shifted destination sign using XOR operation:
44	* R = XOR( R, S );
45	*
46	*/
47
48	/* Offsets for data table __svml_stan_data_internal
49	*/
50	#define _sInvPI_uisa 0
51	#define _sPI1_uisa 16
52	#define _sPI2_uisa 32
53	#define _sPI3_uisa 48
54	#define _sPI2_ha_uisa 64
55	#define _sPI3_ha_uisa 80
56	#define Th_tbl_uisa 96
57	#define Tl_tbl_uisa 224
58	#define _sPC3_uisa 352
59	#define _sPC5_uisa 368
60	#define _sRangeReductionVal_uisa 384
61	#define _sInvPi 400
62	#define _sSignMask 416
63	#define _sAbsMask 432
64	#define _sRangeVal 448
65	#define _sRShifter 464
66	#define _sOne 480
67	#define _sRangeReductionVal 496
68	#define _sPI1 512
69	#define _sPI2 528
70	#define _sPI3 544
71	#define _sPI4 560
72	#define _sPI1_FMA 576
73	#define _sPI2_FMA 592
74	#define _sPI3_FMA 608
75	#define _sP0 624
76	#define _sP1 640
77	#define _sQ0 656
78	#define _sQ1 672
79	#define _sQ2 688
80	#define _sTwo 704
81	#define _sCoeffs 720
82
83	#include <sysdep.h>
84
85	.section .text.sse4, "ax", @progbits
86	ENTRY(_ZGVbN4v_tanf_sse4)
87	subq $232, %rsp
88	cfi_def_cfa_offset(240)
89	movaps %xmm0, %xmm13
90	movups _sAbsMask+__svml_stan_data_internal(%rip), %xmm12
91
92	/*
93	* Legacy Code
94	* Here HW FMA can be unavailable
95	*/
96	xorl %eax, %eax
97	movaps %xmm12, %xmm4
98	pxor %xmm10, %xmm10
99	movups _sInvPi+__svml_stan_data_internal(%rip), %xmm2
100	andps %xmm13, %xmm4
101	mulps %xmm4, %xmm2
102
103	/* Range reduction */
104	movaps %xmm4, %xmm1
105
106	/*
107	*
108	* Main path (_LA_ and _EP_)
109	*
110	* Octant calculation
111	*/
112	movups _sRShifter+__svml_stan_data_internal(%rip), %xmm3
113
114	/* Large values check */
115	movaps %xmm4, %xmm11
116	movups _sPI1+__svml_stan_data_internal(%rip), %xmm5
117	andnps %xmm13, %xmm12
118	movups _sPI2+__svml_stan_data_internal(%rip), %xmm6
119	addps %xmm3, %xmm2
120	cmpnleps _sRangeReductionVal+__svml_stan_data_internal(%rip), %xmm11
121	movaps %xmm2, %xmm8
122	movups _sPI3+__svml_stan_data_internal(%rip), %xmm7
123	subps %xmm3, %xmm8
124	movmskps %xmm11, %edx
125	movups _sPI4+__svml_stan_data_internal(%rip), %xmm9
126	mulps %xmm8, %xmm5
127	mulps %xmm8, %xmm6
128	mulps %xmm8, %xmm7
129	subps %xmm5, %xmm1
130	mulps %xmm8, %xmm9
131	subps %xmm6, %xmm1
132	movups _sQ2+__svml_stan_data_internal(%rip), %xmm15
133
134	/* Inversion mask and sign calculation */
135	movaps %xmm2, %xmm5
136
137	/* Rational approximation */
138	movups _sP1+__svml_stan_data_internal(%rip), %xmm14
139	pslld $30, %xmm2
140	cmpneqps %xmm10, %xmm2
141	subps %xmm7, %xmm1
142
143	/* Exchanged numerator and denominator if necessary */
144	movaps %xmm2, %xmm0
145	movaps %xmm2, %xmm10
146	pslld $31, %xmm5
147	subps %xmm9, %xmm1
148	movaps %xmm1, %xmm3
149	pxor %xmm12, %xmm5
150	mulps %xmm1, %xmm3
151	mulps %xmm3, %xmm15
152	mulps %xmm3, %xmm14
153	addps _sQ1+__svml_stan_data_internal(%rip), %xmm15
154	addps _sP0+__svml_stan_data_internal(%rip), %xmm14
155	mulps %xmm15, %xmm3
156	mulps %xmm14, %xmm1
157	addps _sQ0+__svml_stan_data_internal(%rip), %xmm3
158	andnps %xmm1, %xmm0
159	andps %xmm3, %xmm10
160	andps %xmm2, %xmm1
161	andnps %xmm3, %xmm2
162	orps %xmm10, %xmm0
163	orps %xmm2, %xmm1
164
165	/* Division */
166	divps %xmm1, %xmm0
167
168	/* Sign setting */
169	pxor %xmm5, %xmm0
170
171	/*
172	*
173	* End of main path (_LA_ and _EP_)
174	*/
175
176	testl %edx, %edx
177
178	/* Go to auxiliary branch */
179	jne L(AUX_BRANCH)
180	# LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm4 xmm11 xmm12 xmm13
181
182	/* Return from auxiliary branch
183	* for out of main path inputs
184	*/
185
186	L(AUX_BRANCH_RETURN):
187	testl %eax, %eax
188
189	/* Go to special inputs processing branch */
190	jne L(SPECIAL_VALUES_BRANCH)
191	# LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm13
192
193	/* Restore registers
194	* and exit the function
195	*/
196
197	L(EXIT):
198	addq $232, %rsp
199	cfi_def_cfa_offset(8)
200	ret
201	cfi_def_cfa_offset(240)
202
203	/* Branch to process
204	* special inputs
205	*/
206
207	L(SPECIAL_VALUES_BRANCH):
208	movups %xmm13, 32(%rsp)
209	movups %xmm0, 48(%rsp)
210	# LOE rbx rbp r12 r13 r14 r15 eax xmm0
211
212	xorl %edx, %edx
213	movq %r12, 16(%rsp)
214	cfi_offset(12, -224)
215	movl %edx, %r12d
216	movq %r13, 8(%rsp)
217	cfi_offset(13, -232)
218	movl %eax, %r13d
219	movq %r14, (%rsp)
220	cfi_offset(14, -240)
221	# LOE rbx rbp r15 r12d r13d
222
223	/* Range mask
224	* bits check
225	*/
226
227	L(RANGEMASK_CHECK):
228	btl %r12d, %r13d
229
230	/* Call scalar math function */
231	jc L(SCALAR_MATH_CALL)
232	# LOE rbx rbp r15 r12d r13d
233
234	/* Special inputs
235	* processing loop
236	*/
237
238	L(SPECIAL_VALUES_LOOP):
239	incl %r12d
240	cmpl $4, %r12d
241
242	/* Check bits in range mask */
243	jl L(RANGEMASK_CHECK)
244	# LOE rbx rbp r15 r12d r13d
245
246	movq 16(%rsp), %r12
247	cfi_restore(12)
248	movq 8(%rsp), %r13
249	cfi_restore(13)
250	movq (%rsp), %r14
251	cfi_restore(14)
252	movups 48(%rsp), %xmm0
253
254	/* Go to exit */
255	jmp L(EXIT)
256	cfi_offset(12, -224)
257	cfi_offset(13, -232)
258	cfi_offset(14, -240)
259	# LOE rbx rbp r12 r13 r14 r15 xmm0
260
261	/* Scalar math function call
262	* to process special input
263	*/
264
265	L(SCALAR_MATH_CALL):
266	movl %r12d, %r14d
267	movss 32(%rsp, %r14, 4), %xmm0
268	call tanf@PLT
269	# LOE rbx rbp r14 r15 r12d r13d xmm0
270
271	movss %xmm0, 48(%rsp, %r14, 4)
272
273	/* Process special inputs in loop */
274	jmp L(SPECIAL_VALUES_LOOP)
275	cfi_restore(12)
276	cfi_restore(13)
277	cfi_restore(14)
278	# LOE rbx rbp r15 r12d r13d
279
280	/* Auxiliary branch
281	* for out of main path inputs
282	*/
283
284	L(AUX_BRANCH):
285	movl $2139095040, %eax
286
287	/*
288	* Get the (2^a / 2pi) mod 1 values from the table.
289	* Because doesn't have I-type gather, we need a trivial cast
290	*/
291	lea __svml_stan_reduction_data_internal(%rip), %r8
292	movups %xmm13, 64(%rsp)
293
294	/*
295	* Also get the significand as an integer
296	* NB: adding in the integer bit is wrong for denorms!
297	* To make this work for denorms we should do something slightly different
298	*/
299	movl $8388607, %r9d
300	movups %xmm12, 80(%rsp)
301	movl $8388608, %r10d
302	movups %xmm11, 96(%rsp)
303
304	/*
305	* Break the P_xxx and m into 16-bit chunks ready for
306	* the long multiplication via 16x16->32 multiplications
307	*/
308	movl $65535, %r11d
309	movd %eax, %xmm3
310	pshufd $0, %xmm3, %xmm2
311	andps %xmm2, %xmm13
312	cmpeqps %xmm2, %xmm13
313	pand %xmm4, %xmm2
314	psrld $23, %xmm2
315	movdqa %xmm2, %xmm12
316	pslld $1, %xmm12
317	paddd %xmm2, %xmm12
318	pslld $2, %xmm12
319	pshufd $1, %xmm12, %xmm10
320	pshufd $2, %xmm12, %xmm11
321	pshufd $3, %xmm12, %xmm14
322	movd %xmm12, %edx
323	movd %xmm10, %ecx
324	movd %xmm11, %esi
325	movd %r9d, %xmm11
326	movd %xmm14, %edi
327	movd 4(%rdx, %r8), %xmm6
328	movd 4(%rcx, %r8), %xmm7
329	movd 4(%rsi, %r8), %xmm3
330	movl $872415232, %r9d
331	movd 4(%rdi, %r8), %xmm5
332	punpckldq %xmm7, %xmm6
333	punpckldq %xmm5, %xmm3
334	movd 8(%rdi, %r8), %xmm10
335	movmskps %xmm13, %eax
336	punpcklqdq %xmm3, %xmm6
337	movd 8(%rdx, %r8), %xmm3
338	movd 8(%rcx, %r8), %xmm2
339	movd 8(%rsi, %r8), %xmm13
340	punpckldq %xmm2, %xmm3
341	punpckldq %xmm10, %xmm13
342	punpcklqdq %xmm13, %xmm3
343	pshufd $0, %xmm11, %xmm13
344	movdqa %xmm3, %xmm2
345	movups %xmm4, 48(%rsp)
346	pand %xmm4, %xmm13
347	movd %r10d, %xmm4
348	psrld $16, %xmm2
349	movd (%rdx, %r8), %xmm9
350
351	/*
352	* We want to incorporate the original sign now too.
353	* Do it here for convenience in getting the right N value,
354	* though we could wait right to the end if we were prepared
355	* to modify the sign of N later too.
356	* So get the appropriate sign mask now (or sooner).
357	*/
358	movl $-2147483648, %edx
359	movd (%rcx, %r8), %xmm8
360
361	/*
362	* Create floating-point high part, implicitly adding integer bit 1
363	* Incorporate overall sign at this stage too.
364	*/
365	movl $1065353216, %ecx
366	movd (%rsi, %r8), %xmm15
367
368	/*
369	* Now round at the 2^-8 bit position for reduction mod pi/2^7
370	* instead of the original 2pi (but still with the same 2pi scaling).
371	* Use a shifter of 2^15 + 2^14.
372	* The N we get is our final version; it has an offset of
373	* 2^8 because of the implicit integer bit, and anyway for negative
374	* starting value it's a 2s complement thing. But we need to mask
375	* off the exponent part anyway so it's fine.
376	*/
377	movl $1195376640, %esi
378	movd (%rdi, %r8), %xmm1
379	movl $511, %r10d
380	movups %xmm0, 112(%rsp)
381	movd %r11d, %xmm0
382	pshufd $0, %xmm4, %xmm12
383	movdqa %xmm2, %xmm4
384	punpckldq %xmm8, %xmm9
385	paddd %xmm12, %xmm13
386	punpckldq %xmm1, %xmm15
387	movdqa %xmm13, %xmm12
388	pshufd $0, %xmm0, %xmm8
389	movdqa %xmm6, %xmm0
390	punpcklqdq %xmm15, %xmm9
391	pand %xmm8, %xmm13
392	movdqa %xmm9, %xmm14
393	pand %xmm8, %xmm9
394	movdqa %xmm13, %xmm10
395	psrld $16, %xmm14
396	movdqu %xmm14, 128(%rsp)
397
398	/* Now do the big multiplication and carry propagation */
399	movdqa %xmm9, %xmm14
400	psrlq $32, %xmm10
401	psrlq $32, %xmm14
402	movdqa %xmm13, %xmm15
403	movdqa %xmm10, %xmm7
404	pmuludq %xmm9, %xmm15
405	psrld $16, %xmm0
406	pmuludq %xmm14, %xmm7
407	movdqu %xmm9, 144(%rsp)
408	psllq $32, %xmm7
409	movdqu .FLT_16(%rip), %xmm9
410	pand %xmm8, %xmm6
411	pand %xmm9, %xmm15
412	psrld $16, %xmm12
413	movdqa %xmm0, %xmm1
414	por %xmm7, %xmm15
415	movdqa %xmm13, %xmm7
416	pand %xmm8, %xmm3
417	movdqu %xmm0, 160(%rsp)
418	movdqa %xmm12, %xmm11
419	movdqu %xmm15, 208(%rsp)
420	psrlq $32, %xmm1
421	pmuludq %xmm0, %xmm7
422	movdqa %xmm6, %xmm5
423	movdqa %xmm10, %xmm15
424	movdqa %xmm12, %xmm0
425	movdqu %xmm14, 176(%rsp)
426	psrlq $32, %xmm11
427	movdqu %xmm1, 192(%rsp)
428	psrlq $32, %xmm5
429	pmuludq %xmm1, %xmm15
430	movdqa %xmm13, %xmm1
431	pmuludq %xmm3, %xmm0
432	pmuludq %xmm6, %xmm1
433	pmuludq %xmm12, %xmm6
434	movdqa %xmm10, %xmm14
435	psrlq $32, %xmm3
436	pmuludq %xmm5, %xmm14
437	pand %xmm9, %xmm1
438	pmuludq %xmm11, %xmm3
439	pmuludq %xmm11, %xmm5
440	psllq $32, %xmm14
441	pand %xmm9, %xmm0
442	psllq $32, %xmm3
443	psrlq $32, %xmm4
444	por %xmm14, %xmm1
445	por %xmm3, %xmm0
446	movdqa %xmm12, %xmm14
447	movdqa %xmm11, %xmm3
448	pmuludq %xmm2, %xmm14
449	pand %xmm9, %xmm7
450	pmuludq %xmm4, %xmm3
451	pmuludq %xmm13, %xmm2
452	pmuludq %xmm10, %xmm4
453	pand %xmm9, %xmm2
454	psllq $32, %xmm4
455	psllq $32, %xmm15
456	pand %xmm9, %xmm14
457	psllq $32, %xmm3
458	por %xmm4, %xmm2
459	por %xmm15, %xmm7
460	por %xmm3, %xmm14
461	psrld $16, %xmm2
462	pand %xmm9, %xmm6
463	psllq $32, %xmm5
464	movdqa %xmm1, %xmm15
465	paddd %xmm2, %xmm14
466	movdqa %xmm7, %xmm2
467	por %xmm5, %xmm6
468	psrld $16, %xmm1
469	pand %xmm8, %xmm2
470	paddd %xmm1, %xmm6
471	movdqu 160(%rsp), %xmm1
472	paddd %xmm6, %xmm2
473	movdqu 192(%rsp), %xmm6
474	psrld $16, %xmm7
475	pmuludq %xmm12, %xmm1
476	pand %xmm8, %xmm15
477	pmuludq %xmm11, %xmm6
478	pmuludq 144(%rsp), %xmm12
479	pmuludq 176(%rsp), %xmm11
480	pand %xmm9, %xmm1
481	psllq $32, %xmm6
482	por %xmm6, %xmm1
483	psrld $16, %xmm0
484	paddd %xmm7, %xmm1
485	paddd %xmm14, %xmm15
486	movdqu 128(%rsp), %xmm7
487	paddd %xmm15, %xmm0
488	pmuludq %xmm7, %xmm13
489	psrlq $32, %xmm7
490	pmuludq %xmm7, %xmm10
491	movdqa %xmm0, %xmm14
492	pand %xmm9, %xmm13
493	movdqu 208(%rsp), %xmm5
494	psrld $16, %xmm14
495	paddd %xmm2, %xmm14
496	movdqa %xmm5, %xmm15
497	movdqa %xmm14, %xmm3
498	pand %xmm8, %xmm15
499	psrld $16, %xmm3
500	paddd %xmm1, %xmm15
501	psllq $32, %xmm10
502	pand %xmm9, %xmm12
503	psllq $32, %xmm11
504	paddd %xmm15, %xmm3
505	por %xmm10, %xmm13
506	por %xmm11, %xmm12
507	psrld $16, %xmm5
508	movdqa %xmm3, %xmm4
509	pand %xmm8, %xmm13
510	paddd %xmm5, %xmm12
511	psrld $16, %xmm4
512	paddd %xmm12, %xmm13
513	paddd %xmm13, %xmm4
514	pand %xmm8, %xmm3
515	pslld $16, %xmm4
516	movd %edx, %xmm9
517	movups 48(%rsp), %xmm15
518	paddd %xmm3, %xmm4
519	pshufd $0, %xmm9, %xmm7
520
521	/* Assemble reduced argument from the pieces */
522	pand %xmm8, %xmm0
523	movd %ecx, %xmm8
524	pand %xmm15, %xmm7
525	pshufd $0, %xmm8, %xmm1
526	movdqa %xmm4, %xmm5
527	psrld $9, %xmm5
528	pxor %xmm7, %xmm1
529	por %xmm1, %xmm5
530	movd %esi, %xmm6
531	pshufd $0, %xmm6, %xmm3
532	movdqa %xmm5, %xmm6
533	movl $262143, %r8d
534
535	/*
536	* Create floating-point low and medium parts, respectively
537	* lo_17, ... lo_0, 0, ..., 0
538	* hi_8, ... hi_0, lo_31, ..., lo_18
539	* then subtract off the implicitly added integer bits,
540	* 2^-46 and 2^-23, respectively.
541	* Put the original sign into all of them at this stage.
542	*/
543	movl $679477248, %edi
544	movd %r10d, %xmm13
545	pslld $16, %xmm14
546	pshufd $0, %xmm13, %xmm1
547	paddd %xmm0, %xmm14
548	movd %r9d, %xmm11
549	pand %xmm4, %xmm1
550	movd %r8d, %xmm9
551	movd %edi, %xmm10
552	pshufd $0, %xmm9, %xmm8
553	pslld $14, %xmm1
554	pshufd $0, %xmm10, %xmm0
555	pand %xmm14, %xmm8
556	pshufd $0, %xmm11, %xmm12
557	psrld $18, %xmm14
558	pxor %xmm7, %xmm0
559	pxor %xmm12, %xmm7
560	por %xmm14, %xmm1
561	pslld $5, %xmm8
562	por %xmm7, %xmm1
563
564	/*
565	* Now multiply those numbers all by 2 pi, reasonably accurately.
566	* The top part uses 2pi = s2pi_lead + s2pi_trail, where
567	* s2pi_lead has 12 significant bits.
568	*/
569	movl $1086918619, %r11d
570
571	/* Split RHi into 12-bit leading and trailing parts. */
572	movl $-4096, %esi
573	por %xmm0, %xmm8
574	movl $1086918656, %edx
575	movl $-1214941318, %ecx
576
577	/*
578	* If the magnitude of the input is <= 2^-20, then
579	* just pass through the input, since no reduction will be needed and
580	* the main path will only work accurately if the reduced argument is
581	* about >= 2^-40 (which it is for all large pi multiples)
582	*/
583	movl $2147483647, %edi
584	addps %xmm3, %xmm6
585	subps %xmm7, %xmm1
586	subps %xmm0, %xmm8
587	movaps %xmm6, %xmm2
588	movd %r11d, %xmm14
589	movd %esi, %xmm4
590	movd %edx, %xmm7
591	movl $897581056, %r8d
592	subps %xmm3, %xmm2
593
594	/* Grab our final N value as an integer, appropriately masked mod 2^8 */
595	movl $255, %r9d
596	subps %xmm2, %xmm5
597
598	/* Now add them up into 2 reasonably aligned pieces */
599	movaps %xmm5, %xmm3
600
601	/*
602	* The output is _VRES_R (high) + _VRES_E (low), and the integer part is _VRES_IND
603	* Set sRp2 = _VRES_R^2 and then resume the original code.
604	* Argument reduction is now finished: x = n * pi/128 + r
605	* where n = iIndex and r = sR (high) + sE (low).
606	* But we have n modulo 256, needed for sin/cos with period 2pi
607	* but we want it modulo 128 since tan has period pi.
608	*/
609	movl $127, %r10d
610	pshufd $0, %xmm14, %xmm2
611	addps %xmm1, %xmm3
612	pshufd $0, %xmm4, %xmm14
613	movd %r8d, %xmm4
614	pshufd $0, %xmm4, %xmm9
615	subps %xmm3, %xmm5
616	movdqa %xmm9, %xmm11
617	addps %xmm5, %xmm1
618	movd %ecx, %xmm5
619	addps %xmm1, %xmm8
620	pshufd $0, %xmm7, %xmm1
621	movdqa %xmm14, %xmm7
622	andps %xmm3, %xmm7
623
624	/*
625	* Do the multiplication as exact top part and "naive" low part.
626	* This still maintains a similar level of offset and doesn't drop
627	* the accuracy much below what we already have.
628	*/
629	movdqa %xmm1, %xmm10
630	pshufd $0, %xmm5, %xmm5
631	subps %xmm7, %xmm3
632	mulps %xmm7, %xmm10
633	mulps %xmm5, %xmm7
634	mulps %xmm3, %xmm1
635	mulps %xmm8, %xmm2
636	mulps %xmm3, %xmm5
637	addps %xmm7, %xmm1
638	addps %xmm5, %xmm2
639	movd %edi, %xmm8
640	addps %xmm2, %xmm1
641
642	/*
643	* Do another stage of compensated summation to get full offset
644	* between the pieces sRedHi + sRedLo.
645	* Depending on the later algorithm, we might avoid this stage.
646	*/
647	movaps %xmm1, %xmm0
648
649	/* Load constants (not all needed at once) */
650	lea _sCoeffs+36+__svml_stan_data_internal(%rip), %rdi
651	pshufd $0, %xmm8, %xmm8
652	addps %xmm10, %xmm0
653	andps %xmm15, %xmm8
654	subps %xmm0, %xmm10
655	cmpltps %xmm8, %xmm11
656	cmpleps %xmm9, %xmm8
657	addps %xmm10, %xmm1
658	andps %xmm15, %xmm8
659	movd %r9d, %xmm15
660	andps %xmm11, %xmm0
661	andps %xmm1, %xmm11
662	pshufd $0, %xmm15, %xmm1
663	movd %r10d, %xmm15
664	pshufd $0, %xmm15, %xmm7
665	pand %xmm1, %xmm6
666	pand %xmm7, %xmm6
667	orps %xmm0, %xmm8
668	movaps %xmm6, %xmm4
669
670	/*
671	* Simply combine the two parts of the reduced argument
672	* since we can afford a few ulps in this case.
673	*/
674	addps %xmm11, %xmm8
675	pslld $2, %xmm4
676	paddd %xmm6, %xmm4
677	pslld $3, %xmm4
678	pshufd $1, %xmm4, %xmm6
679	pshufd $2, %xmm4, %xmm5
680	pshufd $3, %xmm4, %xmm3
681	movd %xmm4, %r11d
682	movd %xmm6, %edx
683	movd %xmm5, %ecx
684	movd %xmm3, %esi
685	movd -32(%r11, %rdi), %xmm15
686	movd -32(%rdx, %rdi), %xmm12
687	movd -32(%rcx, %rdi), %xmm7
688	movd -32(%rsi, %rdi), %xmm13
689	punpckldq %xmm12, %xmm15
690	punpckldq %xmm13, %xmm7
691	movd -28(%rsi, %rdi), %xmm5
692	punpcklqdq %xmm7, %xmm15
693	movd -28(%r11, %rdi), %xmm7
694	movd -28(%rdx, %rdi), %xmm6
695	movd -28(%rcx, %rdi), %xmm4
696	movd -36(%rcx, %rdi), %xmm9
697	movd -36(%r11, %rdi), %xmm1
698	movd -36(%rdx, %rdi), %xmm2
699	movd -24(%rdx, %rdi), %xmm3
700	movd -36(%rsi, %rdi), %xmm10
701	punpckldq %xmm6, %xmm7
702	punpckldq %xmm5, %xmm4
703	movd -24(%r11, %rdi), %xmm6
704	punpckldq %xmm2, %xmm1
705	punpckldq %xmm10, %xmm9
706	punpcklqdq %xmm4, %xmm7
707	movd -16(%r11, %rdi), %xmm4
708	punpckldq %xmm3, %xmm6
709	movd -24(%rcx, %rdi), %xmm10
710	movd -16(%rcx, %rdi), %xmm3
711	movd -24(%rsi, %rdi), %xmm2
712	movd -16(%rsi, %rdi), %xmm13
713	movd -16(%rdx, %rdi), %xmm12
714	punpcklqdq %xmm9, %xmm1
715	movd -20(%rdx, %rdi), %xmm9
716	punpckldq %xmm2, %xmm10
717	movd -20(%r11, %rdi), %xmm5
718	movd -20(%rcx, %rdi), %xmm11
719	movd -20(%rsi, %rdi), %xmm0
720	punpckldq %xmm12, %xmm4
721	punpckldq %xmm13, %xmm3
722	punpcklqdq %xmm10, %xmm6
723	movd -12(%rsi, %rdi), %xmm10
724	punpckldq %xmm9, %xmm5
725	punpckldq %xmm0, %xmm11
726	punpcklqdq %xmm3, %xmm4
727	movd -12(%r11, %rdi), %xmm3
728	movd -12(%rdx, %rdi), %xmm2
729	movd -12(%rcx, %rdi), %xmm9
730	punpcklqdq %xmm11, %xmm5
731	punpckldq %xmm2, %xmm3
732	punpckldq %xmm10, %xmm9
733	movd -8(%rcx, %rdi), %xmm10
734	movd -8(%r11, %rdi), %xmm2
735	movd -8(%rdx, %rdi), %xmm0
736	movd -8(%rsi, %rdi), %xmm11
737	punpckldq %xmm0, %xmm2
738	punpckldq %xmm11, %xmm10
739	movd -4(%rsi, %rdi), %xmm13
740	punpcklqdq %xmm9, %xmm3
741	punpcklqdq %xmm10, %xmm2
742	movd -4(%r11, %rdi), %xmm10
743	movd -4(%rdx, %rdi), %xmm12
744	movd -4(%rcx, %rdi), %xmm9
745	punpckldq %xmm12, %xmm10
746	punpckldq %xmm13, %xmm9
747	punpcklqdq %xmm9, %xmm10
748
749	/*
750	* Compute 2-part reciprocal component
751	* Construct a separate reduced argument modulo pi near pi/2 multiples.
752	* i.e. (pi/2 - x) mod pi, simply by subtracting the reduced argument
753	* from an accurate B_hi + B_lo = (128 - n) pi/128. Force the upper part
754	* of this reduced argument to half-length to simplify accurate
755	* reciprocation later on.
756	*/
757	movdqa %xmm1, %xmm9
758	movd (%r11, %rdi), %xmm13
759	subps %xmm8, %xmm9
760	movd (%rdx, %rdi), %xmm0
761	subps %xmm9, %xmm1
762	punpckldq %xmm0, %xmm13
763	movdqa %xmm14, %xmm0
764	andps %xmm9, %xmm0
765	subps %xmm8, %xmm1
766	subps %xmm0, %xmm9
767	movd (%rcx, %rdi), %xmm12
768	addps %xmm9, %xmm15
769
770	/*
771	* Now compute an approximate reciprocal to mix into the computation
772	* To avoid any danger of nonportability, force it to 12 bits,
773	* though I suspect it always is anyway on current platforms.
774	*/
775	rcpps %xmm0, %xmm9
776	addps %xmm15, %xmm1
777	andps %xmm14, %xmm9
778	mulps %xmm9, %xmm0
779
780	/*
781	* Get a better approximation to 1/sR_hi (not far short of an ulp)
782	* using a third-order polynomial approximation
783	*/
784	movaps %xmm9, %xmm14
785	movd (%rsi, %rdi), %xmm11
786
787	/*
788	* Now compute the error sEr where sRecip_hi = (1/R_hi) * (1 - sEr)
789	* so that we can compensate for it.
790	*/
791	movups _sOne+__svml_stan_data_internal(%rip), %xmm15
792	punpckldq %xmm11, %xmm12
793	movaps %xmm15, %xmm11
794	punpcklqdq %xmm12, %xmm13
795	subps %xmm0, %xmm11
796	mulps %xmm11, %xmm14
797	movups %xmm11, (%rsp)
798	addps %xmm9, %xmm14
799	mulps %xmm11, %xmm11
800	movups %xmm13, 32(%rsp)
801	movups %xmm11, 16(%rsp)
802	movups 112(%rsp), %xmm0
803	movups 96(%rsp), %xmm11
804	movups 80(%rsp), %xmm12
805	movups 64(%rsp), %xmm13
806	# LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm1 xmm2 xmm3 xmm4 xmm5 xmm6 xmm7 xmm8 xmm9 xmm10 xmm11 xmm12 xmm13 xmm14 xmm15
807
808	/*
809	* Compensated sum of dominant component(s)
810	* Compute C0_hi + C1_hi * Z + Recip_hi + Recip_lo = H4 (hi) + H9 (lo)
811	* H1 = C1_hi * Z (exact since C1_hi is 1 bit)
812	*/
813	mulps %xmm8, %xmm4
814	addps 16(%rsp), %xmm15
815
816	/* Finally, multiplex both parts so they are only used in cotangent path */
817	mulps %xmm7, %xmm9
818
819	/*
820	* Higher polynomial terms
821	* Stage 1 (with unlimited parallelism)
822	* P3 = C1_lo + C2 * Z
823	*/
824	mulps %xmm8, %xmm2
825	mulps %xmm15, %xmm14
826	addps %xmm2, %xmm3
827
828	/*
829	* Multiply by sRecip_ok to make sR_lo relative to sR_hi
830	* Since sR_lo is shifted off by about 12 bits, this is accurate enough.
831	*/
832	mulps %xmm14, %xmm1
833
834	/*
835	* Now create a low reciprocal using
836	* (Recip_hi + Er * Recip_ok) * (1 + sR_lo^2 - sR_lo)
837	* =~= Recip_hi + Recip_ok * (Er + sR_lo^2 - sR_lo)
838	*/
839	movaps %xmm1, %xmm15
840	mulps %xmm1, %xmm1
841	subps (%rsp), %xmm15
842
843	/* P4 = C3 + C4 * Z */
844	movups 32(%rsp), %xmm2
845	subps %xmm15, %xmm1
846	mulps %xmm8, %xmm2
847	mulps %xmm1, %xmm14
848	addps %xmm2, %xmm10
849	mulps %xmm14, %xmm7
850
851	/* H2 = high(C0_hi + C1_hi * Z) */
852	movdqa %xmm6, %xmm14
853	addps %xmm4, %xmm14
854
855	/* H4 = high(H2 + Recip_hi) */
856	movaps %xmm14, %xmm1
857
858	/* intermediate in compensated sum */
859	subps %xmm14, %xmm6
860	addps %xmm9, %xmm1
861
862	/* H5 = low(C0_hi + C1_hi * Z) */
863	addps %xmm6, %xmm4
864
865	/* intermediate in compensated sum */
866	subps %xmm1, %xmm9
867
868	/* H7 = low(C0_hi + C1_hi * Z) + Recip_lo */
869	addps %xmm4, %xmm7
870
871	/* H8 = low(H2 + Recip_hi) */
872	addps %xmm9, %xmm14
873
874	/* Z2 = Z^2 */
875	movaps %xmm8, %xmm4
876
877	/* Now H4 + H9 should be that part */
878	addps %xmm14, %xmm7
879	mulps %xmm8, %xmm4
880
881	/* P9 = trail(dominant part) + C0_lo */
882	addps %xmm7, %xmm5
883
884	/*
885	* Stage 2 (with unlimited parallelism)
886	* P6 = C1_lo + C2 * Z + C3 * Z^2 + C4 * Z^3
887	*/
888	mulps %xmm4, %xmm10
889	addps %xmm10, %xmm3
890
891	/* Final accumulation of low part */
892	mulps %xmm3, %xmm8
893
894	/* Merge results from main and large paths: */
895	movaps %xmm11, %xmm3
896	andnps %xmm0, %xmm3
897	addps %xmm8, %xmm5
898	movaps %xmm3, %xmm0
899
900	/* And now the very final summation */
901	addps %xmm5, %xmm1
902
903	/*
904	* The end of implementation (LA with huge args reduction)
905	* End of large arguments path (_HA_, _LA_ and _EP_)
906	*/
907
908	pxor %xmm12, %xmm1
909	andps %xmm11, %xmm1
910	orps %xmm1, %xmm0
911
912	/* Return to main vector processing path */
913	jmp L(AUX_BRANCH_RETURN)
914	# LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm13
915	END(_ZGVbN4v_tanf_sse4)
916
917	.section .rodata, "a"
918	.align 16
919
920	#ifdef __svml_stan_data_internal_typedef
921	typedef unsigned int VUINT32;
922	typedef struct {
923	__declspec(align(16)) VUINT32 _sInvPI_uisa[4][1];
924	__declspec(align(16)) VUINT32 _sPI1_uisa[4][1];
925	__declspec(align(16)) VUINT32 _sPI2_uisa[4][1];
926	__declspec(align(16)) VUINT32 _sPI3_uisa[4][1];
927	__declspec(align(16)) VUINT32 _sPI2_ha_uisa[4][1];
928	__declspec(align(16)) VUINT32 _sPI3_ha_uisa[4][1];
929	__declspec(align(16)) VUINT32 Th_tbl_uisa[32][1];
930	__declspec(align(16)) VUINT32 Tl_tbl_uisa[32][1];
931	__declspec(align(16)) VUINT32 _sPC3_uisa[4][1];
932	__declspec(align(16)) VUINT32 _sPC5_uisa[4][1];
933	__declspec(align(16)) VUINT32 _sRangeReductionVal_uisa[4][1];
934	__declspec(align(16)) VUINT32 _sInvPi[4][1];
935	__declspec(align(16)) VUINT32 _sSignMask[4][1];
936	__declspec(align(16)) VUINT32 _sAbsMask[4][1];
937	__declspec(align(16)) VUINT32 _sRangeVal[4][1];
938	__declspec(align(16)) VUINT32 _sRShifter[4][1];
939	__declspec(align(16)) VUINT32 _sOne[4][1];
940	__declspec(align(16)) VUINT32 _sRangeReductionVal[4][1];
941	__declspec(align(16)) VUINT32 _sPI1[4][1];
942	__declspec(align(16)) VUINT32 _sPI2[4][1];
943	__declspec(align(16)) VUINT32 _sPI3[4][1];
944	__declspec(align(16)) VUINT32 _sPI4[4][1];
945	__declspec(align(16)) VUINT32 _sPI1_FMA[4][1];
946	__declspec(align(16)) VUINT32 _sPI2_FMA[4][1];
947	__declspec(align(16)) VUINT32 _sPI3_FMA[4][1];
948	__declspec(align(16)) VUINT32 _sP0[4][1];
949	__declspec(align(16)) VUINT32 _sP1[4][1];
950	__declspec(align(16)) VUINT32 _sQ0[4][1];
951	__declspec(align(16)) VUINT32 _sQ1[4][1];
952	__declspec(align(16)) VUINT32 _sQ2[4][1];
953	__declspec(align(16)) VUINT32 _sTwo[4][1];
954	__declspec(align(16)) VUINT32 _sCoeffs[128][10][1];
955	} __svml_stan_data_internal;
956	#endif
957	__svml_stan_data_internal:
958	/* UISA */
959	.long 0x4122f983, 0x4122f983, 0x4122f983, 0x4122f983 /* _sInvPI_uisa */
960	.align 16
961	.long 0x3dc90fda, 0x3dc90fda, 0x3dc90fda, 0x3dc90fda /* _sPI1_uisa */
962	.align 16
963	.long 0x31a22168, 0x31a22168, 0x31a22168, 0x31a22168 /* _sPI2_uisa */
964	.align 16
965	.long 0x25c234c5, 0x25c234c5, 0x25c234c5, 0x25c234c5 /* _sPI3_uisa */
966	.align 16
967	.long 0x31a22000, 0x31a22000, 0x31a22000, 0x31a22000 /* _sPI2_ha_uisa */
968	.align 16
969	.long 0x2a34611a, 0x2a34611a, 0x2a34611a, 0x2a34611a /* _sPI3_ha_uisa */
970	/* Th_tbl_uisa for i from 0 to 31 do printsingle(tan(i*Pi/32)); */
971	.align 16
972	.long 0x80000000, 0x3dc9b5dc, 0x3e4bafaf, 0x3e9b5042
973	.long 0x3ed413cd, 0x3f08d5b9, 0x3f2b0dc1, 0x3f521801
974	.long 0x3f800000, 0x3f9bf7ec, 0x3fbf90c7, 0x3fef789e
975	.long 0x401a827a, 0x4052facf, 0x40a0dff7, 0x41227363
976	.long 0xff7fffff, 0xc1227363, 0xc0a0dff7, 0xc052facf
977	.long 0xc01a827a, 0xbfef789e, 0xbfbf90c7, 0xbf9bf7ec
978	.long 0xbf800000, 0xbf521801, 0xbf2b0dc1, 0xbf08d5b9
979	.long 0xbed413cd, 0xbe9b5042, 0xbe4bafaf, 0xbdc9b5dc
980	/* Tl_tbl_uisa for i from 0 to 31 do printsingle(tan(i*Pi/32)-round(tan(i*Pi/32), SG, RN)); */
981	.align 16
982	.long 0x80000000, 0x3145b2da, 0x2f2a62b0, 0xb22a39c2
983	.long 0xb1c0621a, 0xb25ef963, 0x32ab7f99, 0x32ae4285
984	.long 0x00000000, 0x33587608, 0x32169d18, 0xb30c3ec0
985	.long 0xb3cc0622, 0x3390600e, 0x331091dc, 0xb454a046
986	.long 0xf3800000, 0x3454a046, 0xb31091dc, 0xb390600e
987	.long 0x33cc0622, 0x330c3ec0, 0xb2169d18, 0xb3587608
988	.long 0x00000000, 0xb2ae4285, 0xb2ab7f99, 0x325ef963
989	.long 0x31c0621a, 0x322a39c2, 0xaf2a62b0, 0xb145b2da
990	.align 16
991	.long 0x3eaaaaa6, 0x3eaaaaa6, 0x3eaaaaa6, 0x3eaaaaa6 /* _sPC3_uisa */
992	.align 16
993	.long 0x3e08b888, 0x3e08b888, 0x3e08b888, 0x3e08b888 /* _sPC5_uisa */
994	.align 16
995	.long 0x46010000, 0x46010000, 0x46010000, 0x46010000 /* _sRangeReductionVal_uisa */
996	.align 16
997	.long 0x3F22F983, 0x3F22F983, 0x3F22F983, 0x3F22F983 /* _sInvPi */
998	.align 16
999	.long 0x80000000, 0x80000000, 0x80000000, 0x80000000 /* _sSignMask */
1000	.align 16
1001	.long 0x7FFFFFFF, 0x7FFFFFFF, 0x7FFFFFFF, 0x7FFFFFFF /* _sAbsMask */
1002	.align 16
1003	.long 0x7f800000, 0x7f800000, 0x7f800000, 0x7f800000 /* _sRangeVal */
1004	.align 16
1005	.long 0x4B400000, 0x4B400000, 0x4B400000, 0x4B400000 /* _sRShifter */
1006	.align 16
1007	.long 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000 /* _sOne */
1008	.align 16
1009	.long 0x46010000, 0x46010000, 0x46010000, 0x46010000 /* _sRangeVal */
1010	.align 16
1011	.long 0x3FC90000, 0x3FC90000, 0x3FC90000, 0x3FC90000 /* _sPI1 */
1012	.align 16
1013	.long 0x39FDA000, 0x39FDA000, 0x39FDA000, 0x39FDA000 /* _sPI2 */
1014	.align 16
1015	.long 0x33A22000, 0x33A22000, 0x33A22000, 0x33A22000 /* _sPI3 */
1016	.align 16
1017	.long 0x2C34611A, 0x2C34611A, 0x2C34611A, 0x2C34611A /* _sPI4 */
1018	// PI1, PI2, and PI3 when FMA is available
1019	.align 16
1020	.long 0x3FC90FDB, 0x3FC90FDB, 0x3FC90FDB, 0x3FC90FDB /* _sPI1_FMA */
1021	.align 16
1022	.long 0xB33BBD2E, 0xB33BBD2E, 0xB33BBD2E, 0xB33BBD2E /* _sPI2_FMA */
1023	.align 16
1024	.long 0xA6F72CED, 0xA6F72CED, 0xA6F72CED, 0xA6F72CED /* _sPI3_FMA */
1025	.align 16
1026	.long 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC /* _sP0 */
1027	.align 16
1028	.long 0xBDC433B4, 0xBDC433B4, 0xBDC433B4, 0xBDC433B4 /* _sP1 */
1029	.align 16
1030	.long 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC /* _sQ0 */
1031	.align 16
1032	.long 0xBEDBB7AB, 0xBEDBB7AB, 0xBEDBB7AB, 0xBEDBB7AB /* _sQ1 */
1033	.align 16
1034	.long 0x3C1F336B, 0x3C1F336B, 0x3C1F336B, 0x3C1F336B /* _sQ2 */
1035	.align 16
1036	.long 0x40000000, 0x40000000, 0x40000000, 0x40000000 /* _sTwo */
1037	// _sCoeffs Breakpoint B = 0 * pi/128, function tan(B + x)
1038	.align 16
1039	.long 0x3FC90FDB // B' = pi/2 - B (high single)
1040	.long 0xB33BBD2E // B' = pi/2 - B (low single)
1041	.long 0x00000000 // tau (1 for cot path)
1042	.long 0x00000000 // c0 (high single)
1043	.long 0x00000000 // c0 (low single)
1044	.long 0x3F800000 // c1 (high 1 bit)
1045	.long 0x00000000 // c1 (low single)
1046	.long 0x00000000 // c2
1047	.long 0x3EAAACDD // c3
1048	.long 0x00000000 // c4
1049	.long 0x3FC5EB9B // B' = pi/2 - B (high single)
1050	.long 0x32DE638C // B' = pi/2 - B (low single)
1051	.long 0x00000000 // tau (1 for cot path)
1052	.long 0x3CC91A31 // c0 (high single)
1053	.long 0x2F8E8D1A // c0 (low single)
1054	.long 0x3F800000 // c1 (high 1 bit)
1055	.long 0x3A1DFA00 // c1 (low single)
1056	.long 0x3CC9392D // c2
1057	.long 0x3EAB1889 // c3
1058	.long 0x3C885D3B // c4
1059	.long 0x3FC2C75C // B' = pi/2 - B (high single)
1060	.long 0xB2CBBE8A // B' = pi/2 - B (low single)
1061	.long 0x00000000 // tau (1 for cot path)
1062	.long 0x3D49393C // c0 (high single)
1063	.long 0x30A39F5B // c0 (low single)
1064	.long 0x3F800000 // c1 (high 1 bit)
1065	.long 0x3B1E2B00 // c1 (low single)
1066	.long 0x3D49B5D4 // c2
1067	.long 0x3EAC4F10 // c3
1068	.long 0x3CFD9425 // c4
1069	.long 0x3FBFA31C // B' = pi/2 - B (high single)
1070	.long 0x33450FB0 // B' = pi/2 - B (low single)
1071	.long 0x00000000 // tau (1 for cot path)
1072	.long 0x3D9711CE // c0 (high single)
1073	.long 0x314FEB28 // c0 (low single)
1074	.long 0x3F800000 // c1 (high 1 bit)
1075	.long 0x3BB24C00 // c1 (low single)
1076	.long 0x3D97E43A // c2
1077	.long 0x3EAE6A89 // c3
1078	.long 0x3D4D07E0 // c4
1079	.long 0x3FBC7EDD // B' = pi/2 - B (high single)
1080	.long 0xB1800ADD // B' = pi/2 - B (low single)
1081	.long 0x00000000 // tau (1 for cot path)
1082	.long 0x3DC9B5DC // c0 (high single)
1083	.long 0x3145AD86 // c0 (low single)
1084	.long 0x3F800000 // c1 (high 1 bit)
1085	.long 0x3C1EEF20 // c1 (low single)
1086	.long 0x3DCBAAEA // c2
1087	.long 0x3EB14E5E // c3
1088	.long 0x3D858BB2 // c4
1089	.long 0x3FB95A9E // B' = pi/2 - B (high single)
1090	.long 0xB3651267 // B' = pi/2 - B (low single)
1091	.long 0x00000000 // tau (1 for cot path)
1092	.long 0x3DFC98C2 // c0 (high single)
1093	.long 0xB0AE525C // c0 (low single)
1094	.long 0x3F800000 // c1 (high 1 bit)
1095	.long 0x3C793D20 // c1 (low single)
1096	.long 0x3E003845 // c2
1097	.long 0x3EB5271F // c3
1098	.long 0x3DAC669E // c4
1099	.long 0x3FB6365E // B' = pi/2 - B (high single)
1100	.long 0x328BB91C // B' = pi/2 - B (low single)
1101	.long 0x00000000 // tau (1 for cot path)
1102	.long 0x3E17E564 // c0 (high single)
1103	.long 0xB1C5A2E4 // c0 (low single)
1104	.long 0x3F800000 // c1 (high 1 bit)
1105	.long 0x3CB440D0 // c1 (low single)
1106	.long 0x3E1B3D00 // c2
1107	.long 0x3EB9F664 // c3
1108	.long 0x3DD647C0 // c4
1109	.long 0x3FB3121F // B' = pi/2 - B (high single)
1110	.long 0xB30F347D // B' = pi/2 - B (low single)
1111	.long 0x00000000 // tau (1 for cot path)
1112	.long 0x3E31AE4D // c0 (high single)
1113	.long 0xB1F32251 // c0 (low single)
1114	.long 0x3F800000 // c1 (high 1 bit)
1115	.long 0x3CF6A500 // c1 (low single)
1116	.long 0x3E3707DA // c2
1117	.long 0x3EBFA489 // c3
1118	.long 0x3DFBD9C7 // c4
1119	.long 0x3FAFEDDF // B' = pi/2 - B (high single)
1120	.long 0x331BBA77 // B' = pi/2 - B (low single)
1121	.long 0x00000000 // tau (1 for cot path)
1122	.long 0x3E4BAFAF // c0 (high single)
1123	.long 0x2F2A29E0 // c0 (low single)
1124	.long 0x3F800000 // c1 (high 1 bit)
1125	.long 0x3D221018 // c1 (low single)
1126	.long 0x3E53BED0 // c2
1127	.long 0x3EC67E26 // c3
1128	.long 0x3E1568E2 // c4
1129	.long 0x3FACC9A0 // B' = pi/2 - B (high single)
1130	.long 0xB2655A50 // B' = pi/2 - B (low single)
1131	.long 0x00000000 // tau (1 for cot path)
1132	.long 0x3E65F267 // c0 (high single)
1133	.long 0x31B4B1DF // c0 (low single)
1134	.long 0x3F800000 // c1 (high 1 bit)
1135	.long 0x3D4E8B90 // c1 (low single)
1136	.long 0x3E718ACA // c2
1137	.long 0x3ECE7164 // c3
1138	.long 0x3E2DC161 // c4
1139	.long 0x3FA9A560 // B' = pi/2 - B (high single)
1140	.long 0x33719861 // B' = pi/2 - B (low single)
1141	.long 0x00000000 // tau (1 for cot path)
1142	.long 0x3E803FD4 // c0 (high single)
1143	.long 0xB2279E66 // c0 (low single)
1144	.long 0x3F800000 // c1 (high 1 bit)
1145	.long 0x3D807FC8 // c1 (low single)
1146	.long 0x3E884BD4 // c2
1147	.long 0x3ED7812D // c3
1148	.long 0x3E4636EB // c4
1149	.long 0x3FA68121 // B' = pi/2 - B (high single)
1150	.long 0x31E43AAC // B' = pi/2 - B (low single)
1151	.long 0x00000000 // tau (1 for cot path)
1152	.long 0x3E8DB082 // c0 (high single)
1153	.long 0xB132A234 // c0 (low single)
1154	.long 0x3F800000 // c1 (high 1 bit)
1155	.long 0x3D9CD7D0 // c1 (low single)
1156	.long 0x3E988A60 // c2
1157	.long 0x3EE203E3 // c3
1158	.long 0x3E63582C // c4
1159	.long 0x3FA35CE2 // B' = pi/2 - B (high single)
1160	.long 0xB33889B6 // B' = pi/2 - B (low single)
1161	.long 0x00000000 // tau (1 for cot path)
1162	.long 0x3E9B5042 // c0 (high single)
1163	.long 0xB22A3AEE // c0 (low single)
1164	.long 0x3F800000 // c1 (high 1 bit)
1165	.long 0x3DBC7490 // c1 (low single)
1166	.long 0x3EA99AF5 // c2
1167	.long 0x3EEDE107 // c3
1168	.long 0x3E80E9AA // c4
1169	.long 0x3FA038A2 // B' = pi/2 - B (high single)
1170	.long 0x32E4CA7E // B' = pi/2 - B (low single)
1171	.long 0x00000000 // tau (1 for cot path)
1172	.long 0x3EA92457 // c0 (high single)
1173	.long 0x30B80830 // c0 (low single)
1174	.long 0x3F800000 // c1 (high 1 bit)
1175	.long 0x3DDF8200 // c1 (low single)
1176	.long 0x3EBB99E9 // c2
1177	.long 0x3EFB4AA8 // c3
1178	.long 0x3E9182BE // c4
1179	.long 0x3F9D1463 // B' = pi/2 - B (high single)
1180	.long 0xB2C55799 // B' = pi/2 - B (low single)
1181	.long 0x00000000 // tau (1 for cot path)
1182	.long 0x3EB73250 // c0 (high single)
1183	.long 0xB2028823 // c0 (low single)
1184	.long 0x3F800000 // c1 (high 1 bit)
1185	.long 0x3E0318F8 // c1 (low single)
1186	.long 0x3ECEA678 // c2
1187	.long 0x3F053C67 // c3
1188	.long 0x3EA41E53 // c4
1189	.long 0x3F99F023 // B' = pi/2 - B (high single)
1190	.long 0x33484328 // B' = pi/2 - B (low single)
1191	.long 0x00000000 // tau (1 for cot path)
1192	.long 0x3EC5800D // c0 (high single)
1193	.long 0xB214C3C1 // c0 (low single)
1194	.long 0x3F800000 // c1 (high 1 bit)
1195	.long 0x3E185E54 // c1 (low single)
1196	.long 0x3EE2E342 // c2
1197	.long 0x3F0DCA73 // c3
1198	.long 0x3EB8CC21 // c4
1199	.long 0x3F96CBE4 // B' = pi/2 - B (high single)
1200	.long 0xB14CDE2E // B' = pi/2 - B (low single)
1201	.long 0x00000000 // tau (1 for cot path)
1202	.long 0x3ED413CD // c0 (high single)
1203	.long 0xB1C06152 // c0 (low single)
1204	.long 0x3F800000 // c1 (high 1 bit)
1205	.long 0x3E2FB0CC // c1 (low single)
1206	.long 0x3EF876CB // c2
1207	.long 0x3F177807 // c3
1208	.long 0x3ED08437 // c4
1209	.long 0x3F93A7A5 // B' = pi/2 - B (high single)
1210	.long 0xB361DEEE // B' = pi/2 - B (low single)
1211	.long 0x00000000 // tau (1 for cot path)
1212	.long 0x3EE2F439 // c0 (high single)
1213	.long 0xB1F4399E // c0 (low single)
1214	.long 0x3F800000 // c1 (high 1 bit)
1215	.long 0x3E49341C // c1 (low single)
1216	.long 0x3F07C61A // c2
1217	.long 0x3F22560F // c3
1218	.long 0x3EEAA81E // c4
1219	.long 0x3F908365 // B' = pi/2 - B (high single)
1220	.long 0x3292200D // B' = pi/2 - B (low single)
1221	.long 0x00000000 // tau (1 for cot path)
1222	.long 0x3EF22870 // c0 (high single)
1223	.long 0x325271F4 // c0 (low single)
1224	.long 0x3F800000 // c1 (high 1 bit)
1225	.long 0x3E65107A // c1 (low single)
1226	.long 0x3F1429F0 // c2
1227	.long 0x3F2E8AFC // c3
1228	.long 0x3F040498 // c4
1229	.long 0x3F8D5F26 // B' = pi/2 - B (high single)
1230	.long 0xB30C0105 // B' = pi/2 - B (low single)
1231	.long 0x00000000 // tau (1 for cot path)
1232	.long 0x3F00DC0D // c0 (high single)
1233	.long 0xB214AF72 // c0 (low single)
1234	.long 0x3F800000 // c1 (high 1 bit)
1235	.long 0x3E81B994 // c1 (low single)
1236	.long 0x3F218233 // c2
1237	.long 0x3F3C4531 // c3
1238	.long 0x3F149688 // c4
1239	.long 0x3F8A3AE6 // B' = pi/2 - B (high single)
1240	.long 0x331EEDF0 // B' = pi/2 - B (low single)
1241	.long 0x00000000 // tau (1 for cot path)
1242	.long 0x3F08D5B9 // c0 (high single)
1243	.long 0xB25EF98E // c0 (low single)
1244	.long 0x3F800000 // c1 (high 1 bit)
1245	.long 0x3E92478D // c1 (low single)
1246	.long 0x3F2FEDC9 // c2
1247	.long 0x3F4BCD58 // c3
1248	.long 0x3F27AE9E // c4
1249	.long 0x3F8716A7 // B' = pi/2 - B (high single)
1250	.long 0xB2588C6D // B' = pi/2 - B (low single)
1251	.long 0x00000000 // tau (1 for cot path)
1252	.long 0x3F1105AF // c0 (high single)
1253	.long 0x32F045B0 // c0 (low single)
1254	.long 0x3F800000 // c1 (high 1 bit)
1255	.long 0x3EA44EE2 // c1 (low single)
1256	.long 0x3F3F8FDB // c2
1257	.long 0x3F5D3FD0 // c3
1258	.long 0x3F3D0A23 // c4
1259	.long 0x3F83F267 // B' = pi/2 - B (high single)
1260	.long 0x3374CBD9 // B' = pi/2 - B (low single)
1261	.long 0x00000000 // tau (1 for cot path)
1262	.long 0x3F1970C4 // c0 (high single)
1263	.long 0x32904848 // c0 (low single)
1264	.long 0x3F800000 // c1 (high 1 bit)
1265	.long 0x3EB7EFF8 // c1 (low single)
1266	.long 0x3F50907C // c2
1267	.long 0x3F710FEA // c3
1268	.long 0x3F561FED // c4
1269	.long 0x3F80CE28 // B' = pi/2 - B (high single)
1270	.long 0x31FDD672 // B' = pi/2 - B (low single)
1271	.long 0x00000000 // tau (1 for cot path)
1272	.long 0x3F221C37 // c0 (high single)
1273	.long 0xB20C61DC // c0 (low single)
1274	.long 0x3F800000 // c1 (high 1 bit)
1275	.long 0x3ECD4F71 // c1 (low single)
1276	.long 0x3F631DAA // c2
1277	.long 0x3F83B471 // c3
1278	.long 0x3F7281EA // c4
1279	.long 0x3F7B53D1 // B' = pi/2 - B (high single)
1280	.long 0x32955386 // B' = pi/2 - B (low single)
1281	.long 0x00000000 // tau (1 for cot path)
1282	.long 0x3F2B0DC1 // c0 (high single)
1283	.long 0x32AB7EBA // c0 (low single)
1284	.long 0x3F800000 // c1 (high 1 bit)
1285	.long 0x3EE496C2 // c1 (low single)
1286	.long 0x3F776C40 // c2
1287	.long 0x3F9065C1 // c3
1288	.long 0x3F89AFB6 // c4
1289	.long 0x3F750B52 // B' = pi/2 - B (high single)
1290	.long 0x32EB316F // B' = pi/2 - B (low single)
1291	.long 0x00000000 // tau (1 for cot path)
1292	.long 0x3F344BA9 // c0 (high single)
1293	.long 0xB2B8B0EA // c0 (low single)
1294	.long 0x3F800000 // c1 (high 1 bit)
1295	.long 0x3EFDF4F7 // c1 (low single)
1296	.long 0x3F86DCA8 // c2
1297	.long 0x3F9ED53B // c3
1298	.long 0x3F9CBEDE // c4
1299	.long 0x3F6EC2D4 // B' = pi/2 - B (high single)
1300	.long 0xB2BEF0A7 // B' = pi/2 - B (low single)
1301	.long 0x00000000 // tau (1 for cot path)
1302	.long 0x3F3DDCCF // c0 (high single)
1303	.long 0x32D29606 // c0 (low single)
1304	.long 0x40000000 // c1 (high 1 bit)
1305	.long 0xBEE6606F // c1 (low single)
1306	.long 0x3F9325D6 // c2
1307	.long 0x3FAF4E69 // c3
1308	.long 0x3FB3080C // c4
1309	.long 0x3F687A55 // B' = pi/2 - B (high single)
1310	.long 0xB252257B // B' = pi/2 - B (low single)
1311	.long 0x00000000 // tau (1 for cot path)
1312	.long 0x3F47C8CC // c0 (high single)
1313	.long 0xB200F51A // c0 (low single)
1314	.long 0x40000000 // c1 (high 1 bit)
1315	.long 0xBEC82C6C // c1 (low single)
1316	.long 0x3FA0BAE9 // c2
1317	.long 0x3FC2252F // c3
1318	.long 0x3FCD24C7 // c4
1319	.long 0x3F6231D6 // B' = pi/2 - B (high single)
1320	.long 0xB119A6A2 // B' = pi/2 - B (low single)
1321	.long 0x00000000 // tau (1 for cot path)
1322	.long 0x3F521801 // c0 (high single)
1323	.long 0x32AE4178 // c0 (low single)
1324	.long 0x40000000 // c1 (high 1 bit)
1325	.long 0xBEA72938 // c1 (low single)
1326	.long 0x3FAFCC22 // c2
1327	.long 0x3FD7BD4A // c3
1328	.long 0x3FEBB01B // c4
1329	.long 0x3F5BE957 // B' = pi/2 - B (high single)
1330	.long 0x3205522A // B' = pi/2 - B (low single)
1331	.long 0x00000000 // tau (1 for cot path)
1332	.long 0x3F5CD3BE // c0 (high single)
1333	.long 0x31460308 // c0 (low single)
1334	.long 0x40000000 // c1 (high 1 bit)
1335	.long 0xBE8306C5 // c1 (low single)
1336	.long 0x3FC09232 // c2
1337	.long 0x3FF09632 // c3
1338	.long 0x4007DB00 // c4
1339	.long 0x3F55A0D8 // B' = pi/2 - B (high single)
1340	.long 0x329886FF // B' = pi/2 - B (low single)
1341	.long 0x00000000 // tau (1 for cot path)
1342	.long 0x3F68065E // c0 (high single)
1343	.long 0x32670D1A // c0 (low single)
1344	.long 0x40000000 // c1 (high 1 bit)
1345	.long 0xBE36D1D6 // c1 (low single)
1346	.long 0x3FD35007 // c2
1347	.long 0x4006A861 // c3
1348	.long 0x401D4BDA // c4
1349	.long 0x3F4F5859 // B' = pi/2 - B (high single)
1350	.long 0x32EE64E8 // B' = pi/2 - B (low single)
1351	.long 0x00000000 // tau (1 for cot path)
1352	.long 0x3F73BB75 // c0 (high single)
1353	.long 0x32FC908D // c0 (low single)
1354	.long 0x40000000 // c1 (high 1 bit)
1355	.long 0xBDBF94B0 // c1 (low single)
1356	.long 0x3FE8550F // c2
1357	.long 0x40174F67 // c3
1358	.long 0x4036C608 // c4
1359	.long 0x3F490FDB // B' = pi/2 - B (high single)
1360	.long 0xB2BBBD2E // B' = pi/2 - B (low single)
1361	.long 0x3F800000 // tau (1 for cot path)
1362	.long 0xBE8BE60E // c0 (high single)
1363	.long 0x320D8D84 // c0 (low single)
1364	.long 0x3F000000 // c1 (high 1 bit)
1365	.long 0xBDF817B1 // c1 (low single)
1366	.long 0xBD8345EB // c2
1367	.long 0x3D1DFDAC // c3
1368	.long 0xBC52CF6F // c4
1369	.long 0x3F42C75C // B' = pi/2 - B (high single)
1370	.long 0xB24BBE8A // B' = pi/2 - B (low single)
1371	.long 0x3F800000 // tau (1 for cot path)
1372	.long 0xBE87283F // c0 (high single)
1373	.long 0xB268B966 // c0 (low single)
1374	.long 0x3F000000 // c1 (high 1 bit)
1375	.long 0xBDFE6529 // c1 (low single)
1376	.long 0xBD7B1953 // c2
1377	.long 0x3D18E109 // c3
1378	.long 0xBC4570B0 // c4
1379	.long 0x3F3C7EDD // B' = pi/2 - B (high single)
1380	.long 0xB1000ADD // B' = pi/2 - B (low single)
1381	.long 0x3F800000 // tau (1 for cot path)
1382	.long 0xBE827420 // c0 (high single)
1383	.long 0x320B8B4D // c0 (low single)
1384	.long 0x3E800000 // c1 (high 1 bit)
1385	.long 0x3DFB9428 // c1 (low single)
1386	.long 0xBD7002B4 // c2
1387	.long 0x3D142A6C // c3
1388	.long 0xBC3A47FF // c4
1389	.long 0x3F36365E // B' = pi/2 - B (high single)
1390	.long 0x320BB91C // B' = pi/2 - B (low single)
1391	.long 0x3F800000 // tau (1 for cot path)
1392	.long 0xBE7B9282 // c0 (high single)
1393	.long 0xB13383D2 // c0 (low single)
1394	.long 0x3E800000 // c1 (high 1 bit)
1395	.long 0x3DF5D211 // c1 (low single)
1396	.long 0xBD6542B3 // c2
1397	.long 0x3D0FE5E5 // c3
1398	.long 0xBC31FB14 // c4
1399	.long 0x3F2FEDDF // B' = pi/2 - B (high single)
1400	.long 0x329BBA77 // B' = pi/2 - B (low single)
1401	.long 0x3F800000 // tau (1 for cot path)
1402	.long 0xBE724E73 // c0 (high single)
1403	.long 0x3120C3E2 // c0 (low single)
1404	.long 0x3E800000 // c1 (high 1 bit)
1405	.long 0x3DF05283 // c1 (low single)
1406	.long 0xBD5AD45E // c2
1407	.long 0x3D0BAFBF // c3
1408	.long 0xBC27B8BB // c4
1409	.long 0x3F29A560 // B' = pi/2 - B (high single)
1410	.long 0x32F19861 // B' = pi/2 - B (low single)
1411	.long 0x3F800000 // tau (1 for cot path)
1412	.long 0xBE691B44 // c0 (high single)
1413	.long 0x31F18936 // c0 (low single)
1414	.long 0x3E800000 // c1 (high 1 bit)
1415	.long 0x3DEB138B // c1 (low single)
1416	.long 0xBD50B2F7 // c2
1417	.long 0x3D07BE3A // c3
1418	.long 0xBC1E46A7 // c4
1419	.long 0x3F235CE2 // B' = pi/2 - B (high single)
1420	.long 0xB2B889B6 // B' = pi/2 - B (low single)
1421	.long 0x3F800000 // tau (1 for cot path)
1422	.long 0xBE5FF82C // c0 (high single)
1423	.long 0xB170723A // c0 (low single)
1424	.long 0x3E800000 // c1 (high 1 bit)
1425	.long 0x3DE61354 // c1 (low single)
1426	.long 0xBD46DA06 // c2
1427	.long 0x3D0401F8 // c3
1428	.long 0xBC14E013 // c4
1429	.long 0x3F1D1463 // B' = pi/2 - B (high single)
1430	.long 0xB2455799 // B' = pi/2 - B (low single)
1431	.long 0x3F800000 // tau (1 for cot path)
1432	.long 0xBE56E46B // c0 (high single)
1433	.long 0x31E3F001 // c0 (low single)
1434	.long 0x3E800000 // c1 (high 1 bit)
1435	.long 0x3DE15025 // c1 (low single)
1436	.long 0xBD3D4550 // c2
1437	.long 0x3D00462D // c3
1438	.long 0xBC092C98 // c4
1439	.long 0x3F16CBE4 // B' = pi/2 - B (high single)
1440	.long 0xB0CCDE2E // B' = pi/2 - B (low single)
1441	.long 0x3F800000 // tau (1 for cot path)
1442	.long 0xBE4DDF41 // c0 (high single)
1443	.long 0xB1AEA094 // c0 (low single)
1444	.long 0x3E800000 // c1 (high 1 bit)
1445	.long 0x3DDCC85C // c1 (low single)
1446	.long 0xBD33F0BE // c2
1447	.long 0x3CFA23B0 // c3
1448	.long 0xBC01FCF7 // c4
1449	.long 0x3F108365 // B' = pi/2 - B (high single)
1450	.long 0x3212200D // B' = pi/2 - B (low single)
1451	.long 0x3F800000 // tau (1 for cot path)
1452	.long 0xBE44E7F8 // c0 (high single)
1453	.long 0xB1CAA3CB // c0 (low single)
1454	.long 0x3E800000 // c1 (high 1 bit)
1455	.long 0x3DD87A74 // c1 (low single)
1456	.long 0xBD2AD885 // c2
1457	.long 0x3CF3C785 // c3
1458	.long 0xBBF1E348 // c4
1459	.long 0x3F0A3AE6 // B' = pi/2 - B (high single)
1460	.long 0x329EEDF0 // B' = pi/2 - B (low single)
1461	.long 0x3F800000 // tau (1 for cot path)
1462	.long 0xBE3BFDDC // c0 (high single)
1463	.long 0xB132521A // c0 (low single)
1464	.long 0x3E800000 // c1 (high 1 bit)
1465	.long 0x3DD464FC // c1 (low single)
1466	.long 0xBD21F8F1 // c2
1467	.long 0x3CEE3076 // c3
1468	.long 0xBBE6D263 // c4
1469	.long 0x3F03F267 // B' = pi/2 - B (high single)
1470	.long 0x32F4CBD9 // B' = pi/2 - B (low single)
1471	.long 0x3F800000 // tau (1 for cot path)
1472	.long 0xBE33203E // c0 (high single)
1473	.long 0x31FEF5BE // c0 (low single)
1474	.long 0x3E800000 // c1 (high 1 bit)
1475	.long 0x3DD0869C // c1 (low single)
1476	.long 0xBD194E8C // c2
1477	.long 0x3CE8DCA9 // c3
1478	.long 0xBBDADA55 // c4
1479	.long 0x3EFB53D1 // B' = pi/2 - B (high single)
1480	.long 0x32155386 // B' = pi/2 - B (low single)
1481	.long 0x3F800000 // tau (1 for cot path)
1482	.long 0xBE2A4E71 // c0 (high single)
1483	.long 0xB19CFCEC // c0 (low single)
1484	.long 0x3E800000 // c1 (high 1 bit)
1485	.long 0x3DCCDE11 // c1 (low single)
1486	.long 0xBD10D605 // c2
1487	.long 0x3CE382A7 // c3
1488	.long 0xBBC8BD97 // c4
1489	.long 0x3EEEC2D4 // B' = pi/2 - B (high single)
1490	.long 0xB23EF0A7 // B' = pi/2 - B (low single)
1491	.long 0x3F800000 // tau (1 for cot path)
1492	.long 0xBE2187D0 // c0 (high single)
1493	.long 0xB1B7C7F7 // c0 (low single)
1494	.long 0x3E800000 // c1 (high 1 bit)
1495	.long 0x3DC96A2B // c1 (low single)
1496	.long 0xBD088C22 // c2
1497	.long 0x3CDE950E // c3
1498	.long 0xBBB89AD1 // c4
1499	.long 0x3EE231D6 // B' = pi/2 - B (high single)
1500	.long 0xB099A6A2 // B' = pi/2 - B (low single)
1501	.long 0x3F800000 // tau (1 for cot path)
1502	.long 0xBE18CBB7 // c0 (high single)
1503	.long 0xAFE28430 // c0 (low single)
1504	.long 0x3E800000 // c1 (high 1 bit)
1505	.long 0x3DC629CE // c1 (low single)
1506	.long 0xBD006DCD // c2
1507	.long 0x3CDA5A2C // c3
1508	.long 0xBBB0B3D2 // c4
1509	.long 0x3ED5A0D8 // B' = pi/2 - B (high single)
1510	.long 0x321886FF // B' = pi/2 - B (low single)
1511	.long 0x3F800000 // tau (1 for cot path)
1512	.long 0xBE101985 // c0 (high single)
1513	.long 0xB02FB2B8 // c0 (low single)
1514	.long 0x3E800000 // c1 (high 1 bit)
1515	.long 0x3DC31BF3 // c1 (low single)
1516	.long 0xBCF0F04D // c2
1517	.long 0x3CD60BC7 // c3
1518	.long 0xBBA138BA // c4
1519	.long 0x3EC90FDB // B' = pi/2 - B (high single)
1520	.long 0xB23BBD2E // B' = pi/2 - B (low single)
1521	.long 0x3F800000 // tau (1 for cot path)
1522	.long 0xBE07709D // c0 (high single)
1523	.long 0xB18A2A83 // c0 (low single)
1524	.long 0x3E800000 // c1 (high 1 bit)
1525	.long 0x3DC03FA2 // c1 (low single)
1526	.long 0xBCE15096 // c2
1527	.long 0x3CD26472 // c3
1528	.long 0xBB9A1270 // c4
1529	.long 0x3EBC7EDD // B' = pi/2 - B (high single)
1530	.long 0xB0800ADD // B' = pi/2 - B (low single)
1531	.long 0x3F800000 // tau (1 for cot path)
1532	.long 0xBDFDA0CB // c0 (high single)
1533	.long 0x2F14FCA0 // c0 (low single)
1534	.long 0x3E800000 // c1 (high 1 bit)
1535	.long 0x3DBD93F7 // c1 (low single)
1536	.long 0xBCD1F71B // c2
1537	.long 0x3CCEDD2B // c3
1538	.long 0xBB905946 // c4
1539	.long 0x3EAFEDDF // B' = pi/2 - B (high single)
1540	.long 0x321BBA77 // B' = pi/2 - B (low single)
1541	.long 0x3F800000 // tau (1 for cot path)
1542	.long 0xBDEC708C // c0 (high single)
1543	.long 0xB14895C4 // c0 (low single)
1544	.long 0x3E800000 // c1 (high 1 bit)
1545	.long 0x3DBB181E // c1 (low single)
1546	.long 0xBCC2DEA6 // c2
1547	.long 0x3CCB5027 // c3
1548	.long 0xBB7F3969 // c4
1549	.long 0x3EA35CE2 // B' = pi/2 - B (high single)
1550	.long 0xB23889B6 // B' = pi/2 - B (low single)
1551	.long 0x3F800000 // tau (1 for cot path)
1552	.long 0xBDDB4F55 // c0 (high single)
1553	.long 0x30F6437E // c0 (low single)
1554	.long 0x3E800000 // c1 (high 1 bit)
1555	.long 0x3DB8CB52 // c1 (low single)
1556	.long 0xBCB40210 // c2
1557	.long 0x3CC82D45 // c3
1558	.long 0xBB643075 // c4
1559	.long 0x3E96CBE4 // B' = pi/2 - B (high single)
1560	.long 0xB04CDE2E // B' = pi/2 - B (low single)
1561	.long 0x3F800000 // tau (1 for cot path)
1562	.long 0xBDCA3BFF // c0 (high single)
1563	.long 0x311C95EA // c0 (low single)
1564	.long 0x3E800000 // c1 (high 1 bit)
1565	.long 0x3DB6ACDE // c1 (low single)
1566	.long 0xBCA55C5B // c2
1567	.long 0x3CC5BC04 // c3
1568	.long 0xBB63A969 // c4
1569	.long 0x3E8A3AE6 // B' = pi/2 - B (high single)
1570	.long 0x321EEDF0 // B' = pi/2 - B (low single)
1571	.long 0x3F800000 // tau (1 for cot path)
1572	.long 0xBDB93569 // c0 (high single)
1573	.long 0xAFB9ED00 // c0 (low single)
1574	.long 0x3E800000 // c1 (high 1 bit)
1575	.long 0x3DB4BC1F // c1 (low single)
1576	.long 0xBC96E905 // c2
1577	.long 0x3CC2E6F5 // c3
1578	.long 0xBB3E10A6 // c4
1579	.long 0x3E7B53D1 // B' = pi/2 - B (high single)
1580	.long 0x31955386 // B' = pi/2 - B (low single)
1581	.long 0x3F800000 // tau (1 for cot path)
1582	.long 0xBDA83A77 // c0 (high single)
1583	.long 0x316D967A // c0 (low single)
1584	.long 0x3E800000 // c1 (high 1 bit)
1585	.long 0x3DB2F87C // c1 (low single)
1586	.long 0xBC88A31F // c2
1587	.long 0x3CC0E763 // c3
1588	.long 0xBB3F1666 // c4
1589	.long 0x3E6231D6 // B' = pi/2 - B (high single)
1590	.long 0xB019A6A2 // B' = pi/2 - B (low single)
1591	.long 0x3F800000 // tau (1 for cot path)
1592	.long 0xBD974A0D // c0 (high single)
1593	.long 0xB14F365B // c0 (low single)
1594	.long 0x3E800000 // c1 (high 1 bit)
1595	.long 0x3DB1616F // c1 (low single)
1596	.long 0xBC750CD8 // c2
1597	.long 0x3CBEB595 // c3
1598	.long 0xBB22B883 // c4
1599	.long 0x3E490FDB // B' = pi/2 - B (high single)
1600	.long 0xB1BBBD2E // B' = pi/2 - B (low single)
1601	.long 0x3F800000 // tau (1 for cot path)
1602	.long 0xBD866317 // c0 (high single)
1603	.long 0xAFF02140 // c0 (low single)
1604	.long 0x3E800000 // c1 (high 1 bit)
1605	.long 0x3DAFF67D // c1 (low single)
1606	.long 0xBC591CD0 // c2
1607	.long 0x3CBCBEAD // c3
1608	.long 0xBB04BBEC // c4
1609	.long 0x3E2FEDDF // B' = pi/2 - B (high single)
1610	.long 0x319BBA77 // B' = pi/2 - B (low single)
1611	.long 0x3F800000 // tau (1 for cot path)
1612	.long 0xBD6B08FF // c0 (high single)
1613	.long 0xB0EED236 // c0 (low single)
1614	.long 0x3E800000 // c1 (high 1 bit)
1615	.long 0x3DAEB739 // c1 (low single)
1616	.long 0xBC3D6D51 // c2
1617	.long 0x3CBB485D // c3
1618	.long 0xBAFFF5BA // c4
1619	.long 0x3E16CBE4 // B' = pi/2 - B (high single)
1620	.long 0xAFCCDE2E // B' = pi/2 - B (low single)
1621	.long 0x3F800000 // tau (1 for cot path)
1622	.long 0xBD495A6C // c0 (high single)
1623	.long 0xB0A427BD // c0 (low single)
1624	.long 0x3E800000 // c1 (high 1 bit)
1625	.long 0x3DADA345 // c1 (low single)
1626	.long 0xBC21F648 // c2
1627	.long 0x3CB9D1B4 // c3
1628	.long 0xBACB5567 // c4
1629	.long 0x3DFB53D1 // B' = pi/2 - B (high single)
1630	.long 0x31155386 // B' = pi/2 - B (low single)
1631	.long 0x3F800000 // tau (1 for cot path)
1632	.long 0xBD27B856 // c0 (high single)
1633	.long 0xB0F7EE91 // c0 (low single)
1634	.long 0x3E800000 // c1 (high 1 bit)
1635	.long 0x3DACBA4E // c1 (low single)
1636	.long 0xBC06AEE3 // c2
1637	.long 0x3CB8E5DC // c3
1638	.long 0xBAEC00EE // c4
1639	.long 0x3DC90FDB // B' = pi/2 - B (high single)
1640	.long 0xB13BBD2E // B' = pi/2 - B (low single)
1641	.long 0x3F800000 // tau (1 for cot path)
1642	.long 0xBD0620A3 // c0 (high single)
1643	.long 0xB0ECAB40 // c0 (low single)
1644	.long 0x3E800000 // c1 (high 1 bit)
1645	.long 0x3DABFC11 // c1 (low single)
1646	.long 0xBBD7200F // c2
1647	.long 0x3CB79475 // c3
1648	.long 0xBA2B0ADC // c4
1649	.long 0x3D96CBE4 // B' = pi/2 - B (high single)
1650	.long 0xAF4CDE2E // B' = pi/2 - B (low single)
1651	.long 0x3F800000 // tau (1 for cot path)
1652	.long 0xBCC92278 // c0 (high single)
1653	.long 0x302F2E68 // c0 (low single)
1654	.long 0x3E800000 // c1 (high 1 bit)
1655	.long 0x3DAB6854 // c1 (low single)
1656	.long 0xBBA1214F // c2
1657	.long 0x3CB6C1E9 // c3
1658	.long 0x3843C2F3 // c4
1659	.long 0x3D490FDB // B' = pi/2 - B (high single)
1660	.long 0xB0BBBD2E // B' = pi/2 - B (low single)
1661	.long 0x3F800000 // tau (1 for cot path)
1662	.long 0xBC861015 // c0 (high single)
1663	.long 0xAFD68E2E // c0 (low single)
1664	.long 0x3E800000 // c1 (high 1 bit)
1665	.long 0x3DAAFEEB // c1 (low single)
1666	.long 0xBB569F3F // c2
1667	.long 0x3CB6A84E // c3
1668	.long 0xBAC64194 // c4
1669	.long 0x3CC90FDB // B' = pi/2 - B (high single)
1670	.long 0xB03BBD2E // B' = pi/2 - B (low single)
1671	.long 0x3F800000 // tau (1 for cot path)
1672	.long 0xBC060BF3 // c0 (high single)
1673	.long 0x2FE251AE // c0 (low single)
1674	.long 0x3E800000 // c1 (high 1 bit)
1675	.long 0x3DAABFB9 // c1 (low single)
1676	.long 0xBAD67C60 // c2
1677	.long 0x3CB64CA5 // c3
1678	.long 0xBACDE881 // c4
1679	.long 0x00000000 // B' = pi/2 - B (high single)
1680	.long 0x00000000 // B' = pi/2 - B (low single)
1681	.long 0x3F800000 // tau (1 for cot path)
1682	.long 0x00000000 // c0 (high single)
1683	.long 0x00000000 // c0 (low single)
1684	.long 0x3E800000 // c1 (high 1 bit)
1685	.long 0x3DAAAAAB // c1 (low single)
1686	.long 0x00000000 // c2
1687	.long 0x3CB5E28B // c3
1688	.long 0x00000000 // c4
1689	.long 0xBCC90FDB // B' = pi/2 - B (high single)
1690	.long 0x303BBD2E // B' = pi/2 - B (low single)
1691	.long 0x3F800000 // tau (1 for cot path)
1692	.long 0x3C060BF3 // c0 (high single)
1693	.long 0xAFE251AE // c0 (low single)
1694	.long 0x3E800000 // c1 (high 1 bit)
1695	.long 0x3DAABFB9 // c1 (low single)
1696	.long 0x3AD67C60 // c2
1697	.long 0x3CB64CA5 // c3
1698	.long 0x3ACDE881 // c4
1699	.long 0xBD490FDB // B' = pi/2 - B (high single)
1700	.long 0x30BBBD2E // B' = pi/2 - B (low single)
1701	.long 0x3F800000 // tau (1 for cot path)
1702	.long 0x3C861015 // c0 (high single)
1703	.long 0x2FD68E2E // c0 (low single)
1704	.long 0x3E800000 // c1 (high 1 bit)
1705	.long 0x3DAAFEEB // c1 (low single)
1706	.long 0x3B569F3F // c2
1707	.long 0x3CB6A84E // c3
1708	.long 0x3AC64194 // c4
1709	.long 0xBD96CBE4 // B' = pi/2 - B (high single)
1710	.long 0x2F4CDE2E // B' = pi/2 - B (low single)
1711	.long 0x3F800000 // tau (1 for cot path)
1712	.long 0x3CC92278 // c0 (high single)
1713	.long 0xB02F2E68 // c0 (low single)
1714	.long 0x3E800000 // c1 (high 1 bit)
1715	.long 0x3DAB6854 // c1 (low single)
1716	.long 0x3BA1214F // c2
1717	.long 0x3CB6C1E9 // c3
1718	.long 0xB843C2F2 // c4
1719	.long 0xBDC90FDB // B' = pi/2 - B (high single)
1720	.long 0x313BBD2E // B' = pi/2 - B (low single)
1721	.long 0x3F800000 // tau (1 for cot path)
1722	.long 0x3D0620A3 // c0 (high single)
1723	.long 0x30ECAB40 // c0 (low single)
1724	.long 0x3E800000 // c1 (high 1 bit)
1725	.long 0x3DABFC11 // c1 (low single)
1726	.long 0x3BD7200F // c2
1727	.long 0x3CB79475 // c3
1728	.long 0x3A2B0ADC // c4
1729	.long 0xBDFB53D1 // B' = pi/2 - B (high single)
1730	.long 0xB1155386 // B' = pi/2 - B (low single)
1731	.long 0x3F800000 // tau (1 for cot path)
1732	.long 0x3D27B856 // c0 (high single)
1733	.long 0x30F7EE91 // c0 (low single)
1734	.long 0x3E800000 // c1 (high 1 bit)
1735	.long 0x3DACBA4E // c1 (low single)
1736	.long 0x3C06AEE3 // c2
1737	.long 0x3CB8E5DC // c3
1738	.long 0x3AEC00EE // c4
1739	.long 0xBE16CBE4 // B' = pi/2 - B (high single)
1740	.long 0x2FCCDE2E // B' = pi/2 - B (low single)
1741	.long 0x3F800000 // tau (1 for cot path)
1742	.long 0x3D495A6C // c0 (high single)
1743	.long 0x30A427BD // c0 (low single)
1744	.long 0x3E800000 // c1 (high 1 bit)
1745	.long 0x3DADA345 // c1 (low single)
1746	.long 0x3C21F648 // c2
1747	.long 0x3CB9D1B4 // c3
1748	.long 0x3ACB5567 // c4
1749	.long 0xBE2FEDDF // B' = pi/2 - B (high single)
1750	.long 0xB19BBA77 // B' = pi/2 - B (low single)
1751	.long 0x3F800000 // tau (1 for cot path)
1752	.long 0x3D6B08FF // c0 (high single)
1753	.long 0x30EED236 // c0 (low single)
1754	.long 0x3E800000 // c1 (high 1 bit)
1755	.long 0x3DAEB739 // c1 (low single)
1756	.long 0x3C3D6D51 // c2
1757	.long 0x3CBB485D // c3
1758	.long 0x3AFFF5BA // c4
1759	.long 0xBE490FDB // B' = pi/2 - B (high single)
1760	.long 0x31BBBD2E // B' = pi/2 - B (low single)
1761	.long 0x3F800000 // tau (1 for cot path)
1762	.long 0x3D866317 // c0 (high single)
1763	.long 0x2FF02140 // c0 (low single)
1764	.long 0x3E800000 // c1 (high 1 bit)
1765	.long 0x3DAFF67D // c1 (low single)
1766	.long 0x3C591CD0 // c2
1767	.long 0x3CBCBEAD // c3
1768	.long 0x3B04BBEC // c4
1769	.long 0xBE6231D6 // B' = pi/2 - B (high single)
1770	.long 0x3019A6A2 // B' = pi/2 - B (low single)
1771	.long 0x3F800000 // tau (1 for cot path)
1772	.long 0x3D974A0D // c0 (high single)
1773	.long 0x314F365B // c0 (low single)
1774	.long 0x3E800000 // c1 (high 1 bit)
1775	.long 0x3DB1616F // c1 (low single)
1776	.long 0x3C750CD8 // c2
1777	.long 0x3CBEB595 // c3
1778	.long 0x3B22B883 // c4
1779	.long 0xBE7B53D1 // B' = pi/2 - B (high single)
1780	.long 0xB1955386 // B' = pi/2 - B (low single)
1781	.long 0x3F800000 // tau (1 for cot path)
1782	.long 0x3DA83A77 // c0 (high single)
1783	.long 0xB16D967A // c0 (low single)
1784	.long 0x3E800000 // c1 (high 1 bit)
1785	.long 0x3DB2F87C // c1 (low single)
1786	.long 0x3C88A31F // c2
1787	.long 0x3CC0E763 // c3
1788	.long 0x3B3F1666 // c4
1789	.long 0xBE8A3AE6 // B' = pi/2 - B (high single)
1790	.long 0xB21EEDF0 // B' = pi/2 - B (low single)
1791	.long 0x3F800000 // tau (1 for cot path)
1792	.long 0x3DB93569 // c0 (high single)
1793	.long 0x2FB9ED00 // c0 (low single)
1794	.long 0x3E800000 // c1 (high 1 bit)
1795	.long 0x3DB4BC1F // c1 (low single)
1796	.long 0x3C96E905 // c2
1797	.long 0x3CC2E6F5 // c3
1798	.long 0x3B3E10A6 // c4
1799	.long 0xBE96CBE4 // B' = pi/2 - B (high single)
1800	.long 0x304CDE2E // B' = pi/2 - B (low single)
1801	.long 0x3F800000 // tau (1 for cot path)
1802	.long 0x3DCA3BFF // c0 (high single)
1803	.long 0xB11C95EA // c0 (low single)
1804	.long 0x3E800000 // c1 (high 1 bit)
1805	.long 0x3DB6ACDE // c1 (low single)
1806	.long 0x3CA55C5B // c2
1807	.long 0x3CC5BC04 // c3
1808	.long 0x3B63A969 // c4
1809	.long 0xBEA35CE2 // B' = pi/2 - B (high single)
1810	.long 0x323889B6 // B' = pi/2 - B (low single)
1811	.long 0x3F800000 // tau (1 for cot path)
1812	.long 0x3DDB4F55 // c0 (high single)
1813	.long 0xB0F6437E // c0 (low single)
1814	.long 0x3E800000 // c1 (high 1 bit)
1815	.long 0x3DB8CB52 // c1 (low single)
1816	.long 0x3CB40210 // c2
1817	.long 0x3CC82D45 // c3
1818	.long 0x3B643075 // c4
1819	.long 0xBEAFEDDF // B' = pi/2 - B (high single)
1820	.long 0xB21BBA77 // B' = pi/2 - B (low single)
1821	.long 0x3F800000 // tau (1 for cot path)
1822	.long 0x3DEC708C // c0 (high single)
1823	.long 0x314895C4 // c0 (low single)
1824	.long 0x3E800000 // c1 (high 1 bit)
1825	.long 0x3DBB181E // c1 (low single)
1826	.long 0x3CC2DEA6 // c2
1827	.long 0x3CCB5027 // c3
1828	.long 0x3B7F3969 // c4
1829	.long 0xBEBC7EDD // B' = pi/2 - B (high single)
1830	.long 0x30800ADD // B' = pi/2 - B (low single)
1831	.long 0x3F800000 // tau (1 for cot path)
1832	.long 0x3DFDA0CB // c0 (high single)
1833	.long 0xAF14FCA0 // c0 (low single)
1834	.long 0x3E800000 // c1 (high 1 bit)
1835	.long 0x3DBD93F7 // c1 (low single)
1836	.long 0x3CD1F71B // c2
1837	.long 0x3CCEDD2B // c3
1838	.long 0x3B905946 // c4
1839	.long 0xBEC90FDB // B' = pi/2 - B (high single)
1840	.long 0x323BBD2E // B' = pi/2 - B (low single)
1841	.long 0x3F800000 // tau (1 for cot path)
1842	.long 0x3E07709D // c0 (high single)
1843	.long 0x318A2A83 // c0 (low single)
1844	.long 0x3E800000 // c1 (high 1 bit)
1845	.long 0x3DC03FA2 // c1 (low single)
1846	.long 0x3CE15096 // c2
1847	.long 0x3CD26472 // c3
1848	.long 0x3B9A1270 // c4
1849	.long 0xBED5A0D8 // B' = pi/2 - B (high single)
1850	.long 0xB21886FF // B' = pi/2 - B (low single)
1851	.long 0x3F800000 // tau (1 for cot path)
1852	.long 0x3E101985 // c0 (high single)
1853	.long 0x302FB2B8 // c0 (low single)
1854	.long 0x3E800000 // c1 (high 1 bit)
1855	.long 0x3DC31BF3 // c1 (low single)
1856	.long 0x3CF0F04D // c2
1857	.long 0x3CD60BC7 // c3
1858	.long 0x3BA138BA // c4
1859	.long 0xBEE231D6 // B' = pi/2 - B (high single)
1860	.long 0x3099A6A2 // B' = pi/2 - B (low single)
1861	.long 0x3F800000 // tau (1 for cot path)
1862	.long 0x3E18CBB7 // c0 (high single)
1863	.long 0x2FE28430 // c0 (low single)
1864	.long 0x3E800000 // c1 (high 1 bit)
1865	.long 0x3DC629CE // c1 (low single)
1866	.long 0x3D006DCD // c2
1867	.long 0x3CDA5A2C // c3
1868	.long 0x3BB0B3D2 // c4
1869	.long 0xBEEEC2D4 // B' = pi/2 - B (high single)
1870	.long 0x323EF0A7 // B' = pi/2 - B (low single)
1871	.long 0x3F800000 // tau (1 for cot path)
1872	.long 0x3E2187D0 // c0 (high single)
1873	.long 0x31B7C7F7 // c0 (low single)
1874	.long 0x3E800000 // c1 (high 1 bit)
1875	.long 0x3DC96A2B // c1 (low single)
1876	.long 0x3D088C22 // c2
1877	.long 0x3CDE950E // c3
1878	.long 0x3BB89AD1 // c4
1879	.long 0xBEFB53D1 // B' = pi/2 - B (high single)
1880	.long 0xB2155386 // B' = pi/2 - B (low single)
1881	.long 0x3F800000 // tau (1 for cot path)
1882	.long 0x3E2A4E71 // c0 (high single)
1883	.long 0x319CFCEC // c0 (low single)
1884	.long 0x3E800000 // c1 (high 1 bit)
1885	.long 0x3DCCDE11 // c1 (low single)
1886	.long 0x3D10D605 // c2
1887	.long 0x3CE382A7 // c3
1888	.long 0x3BC8BD97 // c4
1889	.long 0xBF03F267 // B' = pi/2 - B (high single)
1890	.long 0xB2F4CBD9 // B' = pi/2 - B (low single)
1891	.long 0x3F800000 // tau (1 for cot path)
1892	.long 0x3E33203E // c0 (high single)
1893	.long 0xB1FEF5BE // c0 (low single)
1894	.long 0x3E800000 // c1 (high 1 bit)
1895	.long 0x3DD0869C // c1 (low single)
1896	.long 0x3D194E8C // c2
1897	.long 0x3CE8DCA9 // c3
1898	.long 0x3BDADA55 // c4
1899	.long 0xBF0A3AE6 // B' = pi/2 - B (high single)
1900	.long 0xB29EEDF0 // B' = pi/2 - B (low single)
1901	.long 0x3F800000 // tau (1 for cot path)
1902	.long 0x3E3BFDDC // c0 (high single)
1903	.long 0x3132521A // c0 (low single)
1904	.long 0x3E800000 // c1 (high 1 bit)
1905	.long 0x3DD464FC // c1 (low single)
1906	.long 0x3D21F8F1 // c2
1907	.long 0x3CEE3076 // c3
1908	.long 0x3BE6D263 // c4
1909	.long 0xBF108365 // B' = pi/2 - B (high single)
1910	.long 0xB212200D // B' = pi/2 - B (low single)
1911	.long 0x3F800000 // tau (1 for cot path)
1912	.long 0x3E44E7F8 // c0 (high single)
1913	.long 0x31CAA3CB // c0 (low single)
1914	.long 0x3E800000 // c1 (high 1 bit)
1915	.long 0x3DD87A74 // c1 (low single)
1916	.long 0x3D2AD885 // c2
1917	.long 0x3CF3C785 // c3
1918	.long 0x3BF1E348 // c4
1919	.long 0xBF16CBE4 // B' = pi/2 - B (high single)
1920	.long 0x30CCDE2E // B' = pi/2 - B (low single)
1921	.long 0x3F800000 // tau (1 for cot path)
1922	.long 0x3E4DDF41 // c0 (high single)
1923	.long 0x31AEA094 // c0 (low single)
1924	.long 0x3E800000 // c1 (high 1 bit)
1925	.long 0x3DDCC85C // c1 (low single)
1926	.long 0x3D33F0BE // c2
1927	.long 0x3CFA23B0 // c3
1928	.long 0x3C01FCF7 // c4
1929	.long 0xBF1D1463 // B' = pi/2 - B (high single)
1930	.long 0x32455799 // B' = pi/2 - B (low single)
1931	.long 0x3F800000 // tau (1 for cot path)
1932	.long 0x3E56E46B // c0 (high single)
1933	.long 0xB1E3F001 // c0 (low single)
1934	.long 0x3E800000 // c1 (high 1 bit)
1935	.long 0x3DE15025 // c1 (low single)
1936	.long 0x3D3D4550 // c2
1937	.long 0x3D00462D // c3
1938	.long 0x3C092C98 // c4
1939	.long 0xBF235CE2 // B' = pi/2 - B (high single)
1940	.long 0x32B889B6 // B' = pi/2 - B (low single)
1941	.long 0x3F800000 // tau (1 for cot path)
1942	.long 0x3E5FF82C // c0 (high single)
1943	.long 0x3170723A // c0 (low single)
1944	.long 0x3E800000 // c1 (high 1 bit)
1945	.long 0x3DE61354 // c1 (low single)
1946	.long 0x3D46DA06 // c2
1947	.long 0x3D0401F8 // c3
1948	.long 0x3C14E013 // c4
1949	.long 0xBF29A560 // B' = pi/2 - B (high single)
1950	.long 0xB2F19861 // B' = pi/2 - B (low single)
1951	.long 0x3F800000 // tau (1 for cot path)
1952	.long 0x3E691B44 // c0 (high single)
1953	.long 0xB1F18936 // c0 (low single)
1954	.long 0x3E800000 // c1 (high 1 bit)
1955	.long 0x3DEB138B // c1 (low single)
1956	.long 0x3D50B2F7 // c2
1957	.long 0x3D07BE3A // c3
1958	.long 0x3C1E46A7 // c4
1959	.long 0xBF2FEDDF // B' = pi/2 - B (high single)
1960	.long 0xB29BBA77 // B' = pi/2 - B (low single)
1961	.long 0x3F800000 // tau (1 for cot path)
1962	.long 0x3E724E73 // c0 (high single)
1963	.long 0xB120C3E2 // c0 (low single)
1964	.long 0x3E800000 // c1 (high 1 bit)
1965	.long 0x3DF05283 // c1 (low single)
1966	.long 0x3D5AD45E // c2
1967	.long 0x3D0BAFBF // c3
1968	.long 0x3C27B8BB // c4
1969	.long 0xBF36365E // B' = pi/2 - B (high single)
1970	.long 0xB20BB91C // B' = pi/2 - B (low single)
1971	.long 0x3F800000 // tau (1 for cot path)
1972	.long 0x3E7B9282 // c0 (high single)
1973	.long 0x313383D2 // c0 (low single)
1974	.long 0x3E800000 // c1 (high 1 bit)
1975	.long 0x3DF5D211 // c1 (low single)
1976	.long 0x3D6542B3 // c2
1977	.long 0x3D0FE5E5 // c3
1978	.long 0x3C31FB14 // c4
1979	.long 0xBF3C7EDD // B' = pi/2 - B (high single)
1980	.long 0x31000ADD // B' = pi/2 - B (low single)
1981	.long 0x3F800000 // tau (1 for cot path)
1982	.long 0x3E827420 // c0 (high single)
1983	.long 0xB20B8B4D // c0 (low single)
1984	.long 0x3E800000 // c1 (high 1 bit)
1985	.long 0x3DFB9428 // c1 (low single)
1986	.long 0x3D7002B4 // c2
1987	.long 0x3D142A6C // c3
1988	.long 0x3C3A47FF // c4
1989	.long 0xBF42C75C // B' = pi/2 - B (high single)
1990	.long 0x324BBE8A // B' = pi/2 - B (low single)
1991	.long 0x3F800000 // tau (1 for cot path)
1992	.long 0x3E87283F // c0 (high single)
1993	.long 0x3268B966 // c0 (low single)
1994	.long 0x3F000000 // c1 (high 1 bit)
1995	.long 0xBDFE6529 // c1 (low single)
1996	.long 0x3D7B1953 // c2
1997	.long 0x3D18E109 // c3
1998	.long 0x3C4570B0 // c4
1999	.long 0xBF490FDB // B' = pi/2 - B (high single)
2000	.long 0x32BBBD2E // B' = pi/2 - B (low single)
2001	.long 0x00000000 // tau (1 for cot path)
2002	.long 0xBF800000 // c0 (high single)
2003	.long 0x2B410000 // c0 (low single)
2004	.long 0x40000000 // c1 (high 1 bit)
2005	.long 0xB3000000 // c1 (low single)
2006	.long 0xC0000000 // c2
2007	.long 0x402AB7C8 // c3
2008	.long 0xC05561DB // c4
2009	.long 0xBF4F5859 // B' = pi/2 - B (high single)
2010	.long 0xB2EE64E8 // B' = pi/2 - B (low single)
2011	.long 0x00000000 // tau (1 for cot path)
2012	.long 0xBF73BB75 // c0 (high single)
2013	.long 0xB2FC908D // c0 (low single)
2014	.long 0x40000000 // c1 (high 1 bit)
2015	.long 0xBDBF94B0 // c1 (low single)
2016	.long 0xBFE8550F // c2
2017	.long 0x40174F67 // c3
2018	.long 0xC036C608 // c4
2019	.long 0xBF55A0D8 // B' = pi/2 - B (high single)
2020	.long 0xB29886FF // B' = pi/2 - B (low single)
2021	.long 0x00000000 // tau (1 for cot path)
2022	.long 0xBF68065E // c0 (high single)
2023	.long 0xB2670D1A // c0 (low single)
2024	.long 0x40000000 // c1 (high 1 bit)
2025	.long 0xBE36D1D6 // c1 (low single)
2026	.long 0xBFD35007 // c2
2027	.long 0x4006A861 // c3
2028	.long 0xC01D4BDA // c4
2029	.long 0xBF5BE957 // B' = pi/2 - B (high single)
2030	.long 0xB205522A // B' = pi/2 - B (low single)
2031	.long 0x00000000 // tau (1 for cot path)
2032	.long 0xBF5CD3BE // c0 (high single)
2033	.long 0xB1460308 // c0 (low single)
2034	.long 0x40000000 // c1 (high 1 bit)
2035	.long 0xBE8306C5 // c1 (low single)
2036	.long 0xBFC09232 // c2
2037	.long 0x3FF09632 // c3
2038	.long 0xC007DB00 // c4
2039	.long 0xBF6231D6 // B' = pi/2 - B (high single)
2040	.long 0x3119A6A2 // B' = pi/2 - B (low single)
2041	.long 0x00000000 // tau (1 for cot path)
2042	.long 0xBF521801 // c0 (high single)
2043	.long 0xB2AE4178 // c0 (low single)
2044	.long 0x40000000 // c1 (high 1 bit)
2045	.long 0xBEA72938 // c1 (low single)
2046	.long 0xBFAFCC22 // c2
2047	.long 0x3FD7BD4A // c3
2048	.long 0xBFEBB01B // c4
2049	.long 0xBF687A55 // B' = pi/2 - B (high single)
2050	.long 0x3252257B // B' = pi/2 - B (low single)
2051	.long 0x00000000 // tau (1 for cot path)
2052	.long 0xBF47C8CC // c0 (high single)
2053	.long 0x3200F51A // c0 (low single)
2054	.long 0x40000000 // c1 (high 1 bit)
2055	.long 0xBEC82C6C // c1 (low single)
2056	.long 0xBFA0BAE9 // c2
2057	.long 0x3FC2252F // c3
2058	.long 0xBFCD24C7 // c4
2059	.long 0xBF6EC2D4 // B' = pi/2 - B (high single)
2060	.long 0x32BEF0A7 // B' = pi/2 - B (low single)
2061	.long 0x00000000 // tau (1 for cot path)
2062	.long 0xBF3DDCCF // c0 (high single)
2063	.long 0xB2D29606 // c0 (low single)
2064	.long 0x40000000 // c1 (high 1 bit)
2065	.long 0xBEE6606F // c1 (low single)
2066	.long 0xBF9325D6 // c2
2067	.long 0x3FAF4E69 // c3
2068	.long 0xBFB3080C // c4
2069	.long 0xBF750B52 // B' = pi/2 - B (high single)
2070	.long 0xB2EB316F // B' = pi/2 - B (low single)
2071	.long 0x00000000 // tau (1 for cot path)
2072	.long 0xBF344BA9 // c0 (high single)
2073	.long 0x32B8B0EA // c0 (low single)
2074	.long 0x3F800000 // c1 (high 1 bit)
2075	.long 0x3EFDF4F7 // c1 (low single)
2076	.long 0xBF86DCA8 // c2
2077	.long 0x3F9ED53B // c3
2078	.long 0xBF9CBEDE // c4
2079	.long 0xBF7B53D1 // B' = pi/2 - B (high single)
2080	.long 0xB2955386 // B' = pi/2 - B (low single)
2081	.long 0x00000000 // tau (1 for cot path)
2082	.long 0xBF2B0DC1 // c0 (high single)
2083	.long 0xB2AB7EBA // c0 (low single)
2084	.long 0x3F800000 // c1 (high 1 bit)
2085	.long 0x3EE496C2 // c1 (low single)
2086	.long 0xBF776C40 // c2
2087	.long 0x3F9065C1 // c3
2088	.long 0xBF89AFB6 // c4
2089	.long 0xBF80CE28 // B' = pi/2 - B (high single)
2090	.long 0xB1FDD672 // B' = pi/2 - B (low single)
2091	.long 0x00000000 // tau (1 for cot path)
2092	.long 0xBF221C37 // c0 (high single)
2093	.long 0x320C61DC // c0 (low single)
2094	.long 0x3F800000 // c1 (high 1 bit)
2095	.long 0x3ECD4F71 // c1 (low single)
2096	.long 0xBF631DAA // c2
2097	.long 0x3F83B471 // c3
2098	.long 0xBF7281EA // c4
2099	.long 0xBF83F267 // B' = pi/2 - B (high single)
2100	.long 0xB374CBD9 // B' = pi/2 - B (low single)
2101	.long 0x00000000 // tau (1 for cot path)
2102	.long 0xBF1970C4 // c0 (high single)
2103	.long 0xB2904848 // c0 (low single)
2104	.long 0x3F800000 // c1 (high 1 bit)
2105	.long 0x3EB7EFF8 // c1 (low single)
2106	.long 0xBF50907C // c2
2107	.long 0x3F710FEA // c3
2108	.long 0xBF561FED // c4
2109	.long 0xBF8716A7 // B' = pi/2 - B (high single)
2110	.long 0x32588C6D // B' = pi/2 - B (low single)
2111	.long 0x00000000 // tau (1 for cot path)
2112	.long 0xBF1105AF // c0 (high single)
2113	.long 0xB2F045B0 // c0 (low single)
2114	.long 0x3F800000 // c1 (high 1 bit)
2115	.long 0x3EA44EE2 // c1 (low single)
2116	.long 0xBF3F8FDB // c2
2117	.long 0x3F5D3FD0 // c3
2118	.long 0xBF3D0A23 // c4
2119	.long 0xBF8A3AE6 // B' = pi/2 - B (high single)
2120	.long 0xB31EEDF0 // B' = pi/2 - B (low single)
2121	.long 0x00000000 // tau (1 for cot path)
2122	.long 0xBF08D5B9 // c0 (high single)
2123	.long 0x325EF98E // c0 (low single)
2124	.long 0x3F800000 // c1 (high 1 bit)
2125	.long 0x3E92478D // c1 (low single)
2126	.long 0xBF2FEDC9 // c2
2127	.long 0x3F4BCD58 // c3
2128	.long 0xBF27AE9E // c4
2129	.long 0xBF8D5F26 // B' = pi/2 - B (high single)
2130	.long 0x330C0105 // B' = pi/2 - B (low single)
2131	.long 0x00000000 // tau (1 for cot path)
2132	.long 0xBF00DC0D // c0 (high single)
2133	.long 0x3214AF72 // c0 (low single)
2134	.long 0x3F800000 // c1 (high 1 bit)
2135	.long 0x3E81B994 // c1 (low single)
2136	.long 0xBF218233 // c2
2137	.long 0x3F3C4531 // c3
2138	.long 0xBF149688 // c4
2139	.long 0xBF908365 // B' = pi/2 - B (high single)
2140	.long 0xB292200D // B' = pi/2 - B (low single)
2141	.long 0x00000000 // tau (1 for cot path)
2142	.long 0xBEF22870 // c0 (high single)
2143	.long 0xB25271F4 // c0 (low single)
2144	.long 0x3F800000 // c1 (high 1 bit)
2145	.long 0x3E65107A // c1 (low single)
2146	.long 0xBF1429F0 // c2
2147	.long 0x3F2E8AFC // c3
2148	.long 0xBF040498 // c4
2149	.long 0xBF93A7A5 // B' = pi/2 - B (high single)
2150	.long 0x3361DEEE // B' = pi/2 - B (low single)
2151	.long 0x00000000 // tau (1 for cot path)
2152	.long 0xBEE2F439 // c0 (high single)
2153	.long 0x31F4399E // c0 (low single)
2154	.long 0x3F800000 // c1 (high 1 bit)
2155	.long 0x3E49341C // c1 (low single)
2156	.long 0xBF07C61A // c2
2157	.long 0x3F22560F // c3
2158	.long 0xBEEAA81E // c4
2159	.long 0xBF96CBE4 // B' = pi/2 - B (high single)
2160	.long 0x314CDE2E // B' = pi/2 - B (low single)
2161	.long 0x00000000 // tau (1 for cot path)
2162	.long 0xBED413CD // c0 (high single)
2163	.long 0x31C06152 // c0 (low single)
2164	.long 0x3F800000 // c1 (high 1 bit)
2165	.long 0x3E2FB0CC // c1 (low single)
2166	.long 0xBEF876CB // c2
2167	.long 0x3F177807 // c3
2168	.long 0xBED08437 // c4
2169	.long 0xBF99F023 // B' = pi/2 - B (high single)
2170	.long 0xB3484328 // B' = pi/2 - B (low single)
2171	.long 0x00000000 // tau (1 for cot path)
2172	.long 0xBEC5800D // c0 (high single)
2173	.long 0x3214C3C1 // c0 (low single)
2174	.long 0x3F800000 // c1 (high 1 bit)
2175	.long 0x3E185E54 // c1 (low single)
2176	.long 0xBEE2E342 // c2
2177	.long 0x3F0DCA73 // c3
2178	.long 0xBEB8CC21 // c4
2179	.long 0xBF9D1463 // B' = pi/2 - B (high single)
2180	.long 0x32C55799 // B' = pi/2 - B (low single)
2181	.long 0x00000000 // tau (1 for cot path)
2182	.long 0xBEB73250 // c0 (high single)
2183	.long 0x32028823 // c0 (low single)
2184	.long 0x3F800000 // c1 (high 1 bit)
2185	.long 0x3E0318F8 // c1 (low single)
2186	.long 0xBECEA678 // c2
2187	.long 0x3F053C67 // c3
2188	.long 0xBEA41E53 // c4
2189	.long 0xBFA038A2 // B' = pi/2 - B (high single)
2190	.long 0xB2E4CA7E // B' = pi/2 - B (low single)
2191	.long 0x00000000 // tau (1 for cot path)
2192	.long 0xBEA92457 // c0 (high single)
2193	.long 0xB0B80830 // c0 (low single)
2194	.long 0x3F800000 // c1 (high 1 bit)
2195	.long 0x3DDF8200 // c1 (low single)
2196	.long 0xBEBB99E9 // c2
2197	.long 0x3EFB4AA8 // c3
2198	.long 0xBE9182BE // c4
2199	.long 0xBFA35CE2 // B' = pi/2 - B (high single)
2200	.long 0x333889B6 // B' = pi/2 - B (low single)
2201	.long 0x00000000 // tau (1 for cot path)
2202	.long 0xBE9B5042 // c0 (high single)
2203	.long 0x322A3AEE // c0 (low single)
2204	.long 0x3F800000 // c1 (high 1 bit)
2205	.long 0x3DBC7490 // c1 (low single)
2206	.long 0xBEA99AF5 // c2
2207	.long 0x3EEDE107 // c3
2208	.long 0xBE80E9AA // c4
2209	.long 0xBFA68121 // B' = pi/2 - B (high single)
2210	.long 0xB1E43AAC // B' = pi/2 - B (low single)
2211	.long 0x00000000 // tau (1 for cot path)
2212	.long 0xBE8DB082 // c0 (high single)
2213	.long 0x3132A234 // c0 (low single)
2214	.long 0x3F800000 // c1 (high 1 bit)
2215	.long 0x3D9CD7D0 // c1 (low single)
2216	.long 0xBE988A60 // c2
2217	.long 0x3EE203E3 // c3
2218	.long 0xBE63582C // c4
2219	.long 0xBFA9A560 // B' = pi/2 - B (high single)
2220	.long 0xB3719861 // B' = pi/2 - B (low single)
2221	.long 0x00000000 // tau (1 for cot path)
2222	.long 0xBE803FD4 // c0 (high single)
2223	.long 0x32279E66 // c0 (low single)
2224	.long 0x3F800000 // c1 (high 1 bit)
2225	.long 0x3D807FC8 // c1 (low single)
2226	.long 0xBE884BD4 // c2
2227	.long 0x3ED7812D // c3
2228	.long 0xBE4636EB // c4
2229	.long 0xBFACC9A0 // B' = pi/2 - B (high single)
2230	.long 0x32655A50 // B' = pi/2 - B (low single)
2231	.long 0x00000000 // tau (1 for cot path)
2232	.long 0xBE65F267 // c0 (high single)
2233	.long 0xB1B4B1DF // c0 (low single)
2234	.long 0x3F800000 // c1 (high 1 bit)
2235	.long 0x3D4E8B90 // c1 (low single)
2236	.long 0xBE718ACA // c2
2237	.long 0x3ECE7164 // c3
2238	.long 0xBE2DC161 // c4
2239	.long 0xBFAFEDDF // B' = pi/2 - B (high single)
2240	.long 0xB31BBA77 // B' = pi/2 - B (low single)
2241	.long 0x00000000 // tau (1 for cot path)
2242	.long 0xBE4BAFAF // c0 (high single)
2243	.long 0xAF2A29E0 // c0 (low single)
2244	.long 0x3F800000 // c1 (high 1 bit)
2245	.long 0x3D221018 // c1 (low single)
2246	.long 0xBE53BED0 // c2
2247	.long 0x3EC67E26 // c3
2248	.long 0xBE1568E2 // c4
2249	.long 0xBFB3121F // B' = pi/2 - B (high single)
2250	.long 0x330F347D // B' = pi/2 - B (low single)
2251	.long 0x00000000 // tau (1 for cot path)
2252	.long 0xBE31AE4D // c0 (high single)
2253	.long 0x31F32251 // c0 (low single)
2254	.long 0x3F800000 // c1 (high 1 bit)
2255	.long 0x3CF6A500 // c1 (low single)
2256	.long 0xBE3707DA // c2
2257	.long 0x3EBFA489 // c3
2258	.long 0xBDFBD9C7 // c4
2259	.long 0xBFB6365E // B' = pi/2 - B (high single)
2260	.long 0xB28BB91C // B' = pi/2 - B (low single)
2261	.long 0x00000000 // tau (1 for cot path)
2262	.long 0xBE17E564 // c0 (high single)
2263	.long 0x31C5A2E4 // c0 (low single)
2264	.long 0x3F800000 // c1 (high 1 bit)
2265	.long 0x3CB440D0 // c1 (low single)
2266	.long 0xBE1B3D00 // c2
2267	.long 0x3EB9F664 // c3
2268	.long 0xBDD647C0 // c4
2269	.long 0xBFB95A9E // B' = pi/2 - B (high single)
2270	.long 0x33651267 // B' = pi/2 - B (low single)
2271	.long 0x00000000 // tau (1 for cot path)
2272	.long 0xBDFC98C2 // c0 (high single)
2273	.long 0x30AE525C // c0 (low single)
2274	.long 0x3F800000 // c1 (high 1 bit)
2275	.long 0x3C793D20 // c1 (low single)
2276	.long 0xBE003845 // c2
2277	.long 0x3EB5271F // c3
2278	.long 0xBDAC669E // c4
2279	.long 0xBFBC7EDD // B' = pi/2 - B (high single)
2280	.long 0x31800ADD // B' = pi/2 - B (low single)
2281	.long 0x00000000 // tau (1 for cot path)
2282	.long 0xBDC9B5DC // c0 (high single)
2283	.long 0xB145AD86 // c0 (low single)
2284	.long 0x3F800000 // c1 (high 1 bit)
2285	.long 0x3C1EEF20 // c1 (low single)
2286	.long 0xBDCBAAEA // c2
2287	.long 0x3EB14E5E // c3
2288	.long 0xBD858BB2 // c4
2289	.long 0xBFBFA31C // B' = pi/2 - B (high single)
2290	.long 0xB3450FB0 // B' = pi/2 - B (low single)
2291	.long 0x00000000 // tau (1 for cot path)
2292	.long 0xBD9711CE // c0 (high single)
2293	.long 0xB14FEB28 // c0 (low single)
2294	.long 0x3F800000 // c1 (high 1 bit)
2295	.long 0x3BB24C00 // c1 (low single)
2296	.long 0xBD97E43A // c2
2297	.long 0x3EAE6A89 // c3
2298	.long 0xBD4D07E0 // c4
2299	.long 0xBFC2C75C // B' = pi/2 - B (high single)
2300	.long 0x32CBBE8A // B' = pi/2 - B (low single)
2301	.long 0x00000000 // tau (1 for cot path)
2302	.long 0xBD49393C // c0 (high single)
2303	.long 0xB0A39F5B // c0 (low single)
2304	.long 0x3F800000 // c1 (high 1 bit)
2305	.long 0x3B1E2B00 // c1 (low single)
2306	.long 0xBD49B5D4 // c2
2307	.long 0x3EAC4F10 // c3
2308	.long 0xBCFD9425 // c4
2309	.long 0xBFC5EB9B // B' = pi/2 - B (high single)
2310	.long 0xB2DE638C // B' = pi/2 - B (low single)
2311	.long 0x00000000 // tau (1 for cot path)
2312	.long 0xBCC91A31 // c0 (high single)
2313	.long 0xAF8E8D1A // c0 (low single)
2314	.long 0x3F800000 // c1 (high 1 bit)
2315	.long 0x3A1DFA00 // c1 (low single)
2316	.long 0xBCC9392D // c2
2317	.long 0x3EAB1889 // c3
2318	.long 0xBC885D3B // c4
2319	.align 16
2320	.type __svml_stan_data_internal, @object
2321	.size __svml_stan_data_internal, .-__svml_stan_data_internal
2322	.space 16, 0x00
2323	.align 16
2324
2325	#ifdef __svml_stan_reduction_data_internal_typedef
2326	typedef unsigned int VUINT32;
2327	typedef struct {
2328	__declspec(align(16)) VUINT32 _sPtable[256][3][1];
2329	} __svml_stan_reduction_data_internal;
2330	#endif
2331	__svml_stan_reduction_data_internal:
2332	/* P_hi P_med P_lo */
2333	.long 0x00000000, 0x00000000, 0x00000000 /* 0 */
2334	.long 0x00000000, 0x00000000, 0x00000000 /* 1 */
2335	.long 0x00000000, 0x00000000, 0x00000000 /* 2 */
2336	.long 0x00000000, 0x00000000, 0x00000000 /* 3 */
2337	.long 0x00000000, 0x00000000, 0x00000000 /* 4 */
2338	.long 0x00000000, 0x00000000, 0x00000000 /* 5 */
2339	.long 0x00000000, 0x00000000, 0x00000000 /* 6 */
2340	.long 0x00000000, 0x00000000, 0x00000000 /* 7 */
2341	.long 0x00000000, 0x00000000, 0x00000000 /* 8 */
2342	.long 0x00000000, 0x00000000, 0x00000000 /* 9 */
2343	.long 0x00000000, 0x00000000, 0x00000000 /* 10 */
2344	.long 0x00000000, 0x00000000, 0x00000000 /* 11 */
2345	.long 0x00000000, 0x00000000, 0x00000000 /* 12 */
2346	.long 0x00000000, 0x00000000, 0x00000000 /* 13 */
2347	.long 0x00000000, 0x00000000, 0x00000000 /* 14 */
2348	.long 0x00000000, 0x00000000, 0x00000000 /* 15 */
2349	.long 0x00000000, 0x00000000, 0x00000000 /* 16 */
2350	.long 0x00000000, 0x00000000, 0x00000000 /* 17 */
2351	.long 0x00000000, 0x00000000, 0x00000000 /* 18 */
2352	.long 0x00000000, 0x00000000, 0x00000000 /* 19 */
2353	.long 0x00000000, 0x00000000, 0x00000000 /* 20 */
2354	.long 0x00000000, 0x00000000, 0x00000000 /* 21 */
2355	.long 0x00000000, 0x00000000, 0x00000000 /* 22 */
2356	.long 0x00000000, 0x00000000, 0x00000000 /* 23 */
2357	.long 0x00000000, 0x00000000, 0x00000000 /* 24 */
2358	.long 0x00000000, 0x00000000, 0x00000000 /* 25 */
2359	.long 0x00000000, 0x00000000, 0x00000000 /* 26 */
2360	.long 0x00000000, 0x00000000, 0x00000000 /* 27 */
2361	.long 0x00000000, 0x00000000, 0x00000000 /* 28 */
2362	.long 0x00000000, 0x00000000, 0x00000000 /* 29 */
2363	.long 0x00000000, 0x00000000, 0x00000000 /* 30 */
2364	.long 0x00000000, 0x00000000, 0x00000000 /* 31 */
2365	.long 0x00000000, 0x00000000, 0x00000000 /* 32 */
2366	.long 0x00000000, 0x00000000, 0x00000000 /* 33 */
2367	.long 0x00000000, 0x00000000, 0x00000000 /* 34 */
2368	.long 0x00000000, 0x00000000, 0x00000000 /* 35 */
2369	.long 0x00000000, 0x00000000, 0x00000000 /* 36 */
2370	.long 0x00000000, 0x00000000, 0x00000000 /* 37 */
2371	.long 0x00000000, 0x00000000, 0x00000000 /* 38 */
2372	.long 0x00000000, 0x00000000, 0x00000000 /* 39 */
2373	.long 0x00000000, 0x00000000, 0x00000000 /* 40 */
2374	.long 0x00000000, 0x00000000, 0x00000000 /* 41 */
2375	.long 0x00000000, 0x00000000, 0x00000000 /* 42 */
2376	.long 0x00000000, 0x00000000, 0x00000000 /* 43 */
2377	.long 0x00000000, 0x00000000, 0x00000000 /* 44 */
2378	.long 0x00000000, 0x00000000, 0x00000000 /* 45 */
2379	.long 0x00000000, 0x00000000, 0x00000000 /* 46 */
2380	.long 0x00000000, 0x00000000, 0x00000000 /* 47 */
2381	.long 0x00000000, 0x00000000, 0x00000000 /* 48 */
2382	.long 0x00000000, 0x00000000, 0x00000000 /* 49 */
2383	.long 0x00000000, 0x00000000, 0x00000000 /* 50 */
2384	.long 0x00000000, 0x00000000, 0x00000000 /* 51 */
2385	.long 0x00000000, 0x00000000, 0x00000000 /* 52 */
2386	.long 0x00000000, 0x00000000, 0x00000000 /* 53 */
2387	.long 0x00000000, 0x00000000, 0x00000000 /* 54 */
2388	.long 0x00000000, 0x00000000, 0x00000000 /* 55 */
2389	.long 0x00000000, 0x00000000, 0x00000000 /* 56 */
2390	.long 0x00000000, 0x00000000, 0x00000001 /* 57 */
2391	.long 0x00000000, 0x00000000, 0x00000002 /* 58 */
2392	.long 0x00000000, 0x00000000, 0x00000005 /* 59 */
2393	.long 0x00000000, 0x00000000, 0x0000000A /* 60 */
2394	.long 0x00000000, 0x00000000, 0x00000014 /* 61 */
2395	.long 0x00000000, 0x00000000, 0x00000028 /* 62 */
2396	.long 0x00000000, 0x00000000, 0x00000051 /* 63 */
2397	.long 0x00000000, 0x00000000, 0x000000A2 /* 64 */
2398	.long 0x00000000, 0x00000000, 0x00000145 /* 65 */
2399	.long 0x00000000, 0x00000000, 0x0000028B /* 66 */
2400	.long 0x00000000, 0x00000000, 0x00000517 /* 67 */
2401	.long 0x00000000, 0x00000000, 0x00000A2F /* 68 */
2402	.long 0x00000000, 0x00000000, 0x0000145F /* 69 */
2403	.long 0x00000000, 0x00000000, 0x000028BE /* 70 */
2404	.long 0x00000000, 0x00000000, 0x0000517C /* 71 */
2405	.long 0x00000000, 0x00000000, 0x0000A2F9 /* 72 */
2406	.long 0x00000000, 0x00000000, 0x000145F3 /* 73 */
2407	.long 0x00000000, 0x00000000, 0x00028BE6 /* 74 */
2408	.long 0x00000000, 0x00000000, 0x000517CC /* 75 */
2409	.long 0x00000000, 0x00000000, 0x000A2F98 /* 76 */
2410	.long 0x00000000, 0x00000000, 0x00145F30 /* 77 */
2411	.long 0x00000000, 0x00000000, 0x0028BE60 /* 78 */
2412	.long 0x00000000, 0x00000000, 0x00517CC1 /* 79 */
2413	.long 0x00000000, 0x00000000, 0x00A2F983 /* 80 */
2414	.long 0x00000000, 0x00000000, 0x0145F306 /* 81 */
2415	.long 0x00000000, 0x00000000, 0x028BE60D /* 82 */
2416	.long 0x00000000, 0x00000000, 0x0517CC1B /* 83 */
2417	.long 0x00000000, 0x00000000, 0x0A2F9836 /* 84 */
2418	.long 0x00000000, 0x00000000, 0x145F306D /* 85 */
2419	.long 0x00000000, 0x00000000, 0x28BE60DB /* 86 */
2420	.long 0x00000000, 0x00000000, 0x517CC1B7 /* 87 */
2421	.long 0x00000000, 0x00000000, 0xA2F9836E /* 88 */
2422	.long 0x00000000, 0x00000001, 0x45F306DC /* 89 */
2423	.long 0x00000000, 0x00000002, 0x8BE60DB9 /* 90 */
2424	.long 0x00000000, 0x00000005, 0x17CC1B72 /* 91 */
2425	.long 0x00000000, 0x0000000A, 0x2F9836E4 /* 92 */
2426	.long 0x00000000, 0x00000014, 0x5F306DC9 /* 93 */
2427	.long 0x00000000, 0x00000028, 0xBE60DB93 /* 94 */
2428	.long 0x00000000, 0x00000051, 0x7CC1B727 /* 95 */
2429	.long 0x00000000, 0x000000A2, 0xF9836E4E /* 96 */
2430	.long 0x00000000, 0x00000145, 0xF306DC9C /* 97 */
2431	.long 0x00000000, 0x0000028B, 0xE60DB939 /* 98 */
2432	.long 0x00000000, 0x00000517, 0xCC1B7272 /* 99 */
2433	.long 0x00000000, 0x00000A2F, 0x9836E4E4 /* 100 */
2434	.long 0x00000000, 0x0000145F, 0x306DC9C8 /* 101 */
2435	.long 0x00000000, 0x000028BE, 0x60DB9391 /* 102 */
2436	.long 0x00000000, 0x0000517C, 0xC1B72722 /* 103 */
2437	.long 0x00000000, 0x0000A2F9, 0x836E4E44 /* 104 */
2438	.long 0x00000000, 0x000145F3, 0x06DC9C88 /* 105 */
2439	.long 0x00000000, 0x00028BE6, 0x0DB93910 /* 106 */
2440	.long 0x00000000, 0x000517CC, 0x1B727220 /* 107 */
2441	.long 0x00000000, 0x000A2F98, 0x36E4E441 /* 108 */
2442	.long 0x00000000, 0x00145F30, 0x6DC9C882 /* 109 */
2443	.long 0x00000000, 0x0028BE60, 0xDB939105 /* 110 */
2444	.long 0x00000000, 0x00517CC1, 0xB727220A /* 111 */
2445	.long 0x00000000, 0x00A2F983, 0x6E4E4415 /* 112 */
2446	.long 0x00000000, 0x0145F306, 0xDC9C882A /* 113 */
2447	.long 0x00000000, 0x028BE60D, 0xB9391054 /* 114 */
2448	.long 0x00000000, 0x0517CC1B, 0x727220A9 /* 115 */
2449	.long 0x00000000, 0x0A2F9836, 0xE4E44152 /* 116 */
2450	.long 0x00000000, 0x145F306D, 0xC9C882A5 /* 117 */
2451	.long 0x00000000, 0x28BE60DB, 0x9391054A /* 118 */
2452	.long 0x00000000, 0x517CC1B7, 0x27220A94 /* 119 */
2453	.long 0x00000000, 0xA2F9836E, 0x4E441529 /* 120 */
2454	.long 0x00000001, 0x45F306DC, 0x9C882A53 /* 121 */
2455	.long 0x00000002, 0x8BE60DB9, 0x391054A7 /* 122 */
2456	.long 0x00000005, 0x17CC1B72, 0x7220A94F /* 123 */
2457	.long 0x0000000A, 0x2F9836E4, 0xE441529F /* 124 */
2458	.long 0x00000014, 0x5F306DC9, 0xC882A53F /* 125 */
2459	.long 0x00000028, 0xBE60DB93, 0x91054A7F /* 126 */
2460	.long 0x00000051, 0x7CC1B727, 0x220A94FE /* 127 */
2461	.long 0x000000A2, 0xF9836E4E, 0x441529FC /* 128 */
2462	.long 0x00000145, 0xF306DC9C, 0x882A53F8 /* 129 */
2463	.long 0x0000028B, 0xE60DB939, 0x1054A7F0 /* 130 */
2464	.long 0x00000517, 0xCC1B7272, 0x20A94FE1 /* 131 */
2465	.long 0x00000A2F, 0x9836E4E4, 0x41529FC2 /* 132 */
2466	.long 0x0000145F, 0x306DC9C8, 0x82A53F84 /* 133 */
2467	.long 0x000028BE, 0x60DB9391, 0x054A7F09 /* 134 */
2468	.long 0x0000517C, 0xC1B72722, 0x0A94FE13 /* 135 */
2469	.long 0x0000A2F9, 0x836E4E44, 0x1529FC27 /* 136 */
2470	.long 0x000145F3, 0x06DC9C88, 0x2A53F84E /* 137 */
2471	.long 0x00028BE6, 0x0DB93910, 0x54A7F09D /* 138 */
2472	.long 0x000517CC, 0x1B727220, 0xA94FE13A /* 139 */
2473	.long 0x000A2F98, 0x36E4E441, 0x529FC275 /* 140 */
2474	.long 0x00145F30, 0x6DC9C882, 0xA53F84EA /* 141 */
2475	.long 0x0028BE60, 0xDB939105, 0x4A7F09D5 /* 142 */
2476	.long 0x00517CC1, 0xB727220A, 0x94FE13AB /* 143 */
2477	.long 0x00A2F983, 0x6E4E4415, 0x29FC2757 /* 144 */
2478	.long 0x0145F306, 0xDC9C882A, 0x53F84EAF /* 145 */
2479	.long 0x028BE60D, 0xB9391054, 0xA7F09D5F /* 146 */
2480	.long 0x0517CC1B, 0x727220A9, 0x4FE13ABE /* 147 */
2481	.long 0x0A2F9836, 0xE4E44152, 0x9FC2757D /* 148 */
2482	.long 0x145F306D, 0xC9C882A5, 0x3F84EAFA /* 149 */
2483	.long 0x28BE60DB, 0x9391054A, 0x7F09D5F4 /* 150 */
2484	.long 0x517CC1B7, 0x27220A94, 0xFE13ABE8 /* 151 */
2485	.long 0xA2F9836E, 0x4E441529, 0xFC2757D1 /* 152 */
2486	.long 0x45F306DC, 0x9C882A53, 0xF84EAFA3 /* 153 */
2487	.long 0x8BE60DB9, 0x391054A7, 0xF09D5F47 /* 154 */
2488	.long 0x17CC1B72, 0x7220A94F, 0xE13ABE8F /* 155 */
2489	.long 0x2F9836E4, 0xE441529F, 0xC2757D1F /* 156 */
2490	.long 0x5F306DC9, 0xC882A53F, 0x84EAFA3E /* 157 */
2491	.long 0xBE60DB93, 0x91054A7F, 0x09D5F47D /* 158 */
2492	.long 0x7CC1B727, 0x220A94FE, 0x13ABE8FA /* 159 */
2493	.long 0xF9836E4E, 0x441529FC, 0x2757D1F5 /* 160 */
2494	.long 0xF306DC9C, 0x882A53F8, 0x4EAFA3EA /* 161 */
2495	.long 0xE60DB939, 0x1054A7F0, 0x9D5F47D4 /* 162 */
2496	.long 0xCC1B7272, 0x20A94FE1, 0x3ABE8FA9 /* 163 */
2497	.long 0x9836E4E4, 0x41529FC2, 0x757D1F53 /* 164 */
2498	.long 0x306DC9C8, 0x82A53F84, 0xEAFA3EA6 /* 165 */
2499	.long 0x60DB9391, 0x054A7F09, 0xD5F47D4D /* 166 */
2500	.long 0xC1B72722, 0x0A94FE13, 0xABE8FA9A /* 167 */
2501	.long 0x836E4E44, 0x1529FC27, 0x57D1F534 /* 168 */
2502	.long 0x06DC9C88, 0x2A53F84E, 0xAFA3EA69 /* 169 */
2503	.long 0x0DB93910, 0x54A7F09D, 0x5F47D4D3 /* 170 */
2504	.long 0x1B727220, 0xA94FE13A, 0xBE8FA9A6 /* 171 */
2505	.long 0x36E4E441, 0x529FC275, 0x7D1F534D /* 172 */
2506	.long 0x6DC9C882, 0xA53F84EA, 0xFA3EA69B /* 173 */
2507	.long 0xDB939105, 0x4A7F09D5, 0xF47D4D37 /* 174 */
2508	.long 0xB727220A, 0x94FE13AB, 0xE8FA9A6E /* 175 */
2509	.long 0x6E4E4415, 0x29FC2757, 0xD1F534DD /* 176 */
2510	.long 0xDC9C882A, 0x53F84EAF, 0xA3EA69BB /* 177 */
2511	.long 0xB9391054, 0xA7F09D5F, 0x47D4D377 /* 178 */
2512	.long 0x727220A9, 0x4FE13ABE, 0x8FA9A6EE /* 179 */
2513	.long 0xE4E44152, 0x9FC2757D, 0x1F534DDC /* 180 */
2514	.long 0xC9C882A5, 0x3F84EAFA, 0x3EA69BB8 /* 181 */
2515	.long 0x9391054A, 0x7F09D5F4, 0x7D4D3770 /* 182 */
2516	.long 0x27220A94, 0xFE13ABE8, 0xFA9A6EE0 /* 183 */
2517	.long 0x4E441529, 0xFC2757D1, 0xF534DDC0 /* 184 */
2518	.long 0x9C882A53, 0xF84EAFA3, 0xEA69BB81 /* 185 */
2519	.long 0x391054A7, 0xF09D5F47, 0xD4D37703 /* 186 */
2520	.long 0x7220A94F, 0xE13ABE8F, 0xA9A6EE06 /* 187 */
2521	.long 0xE441529F, 0xC2757D1F, 0x534DDC0D /* 188 */
2522	.long 0xC882A53F, 0x84EAFA3E, 0xA69BB81B /* 189 */
2523	.long 0x91054A7F, 0x09D5F47D, 0x4D377036 /* 190 */
2524	.long 0x220A94FE, 0x13ABE8FA, 0x9A6EE06D /* 191 */
2525	.long 0x441529FC, 0x2757D1F5, 0x34DDC0DB /* 192 */
2526	.long 0x882A53F8, 0x4EAFA3EA, 0x69BB81B6 /* 193 */
2527	.long 0x1054A7F0, 0x9D5F47D4, 0xD377036D /* 194 */
2528	.long 0x20A94FE1, 0x3ABE8FA9, 0xA6EE06DB /* 195 */
2529	.long 0x41529FC2, 0x757D1F53, 0x4DDC0DB6 /* 196 */
2530	.long 0x82A53F84, 0xEAFA3EA6, 0x9BB81B6C /* 197 */
2531	.long 0x054A7F09, 0xD5F47D4D, 0x377036D8 /* 198 */
2532	.long 0x0A94FE13, 0xABE8FA9A, 0x6EE06DB1 /* 199 */
2533	.long 0x1529FC27, 0x57D1F534, 0xDDC0DB62 /* 200 */
2534	.long 0x2A53F84E, 0xAFA3EA69, 0xBB81B6C5 /* 201 */
2535	.long 0x54A7F09D, 0x5F47D4D3, 0x77036D8A /* 202 */
2536	.long 0xA94FE13A, 0xBE8FA9A6, 0xEE06DB14 /* 203 */
2537	.long 0x529FC275, 0x7D1F534D, 0xDC0DB629 /* 204 */
2538	.long 0xA53F84EA, 0xFA3EA69B, 0xB81B6C52 /* 205 */
2539	.long 0x4A7F09D5, 0xF47D4D37, 0x7036D8A5 /* 206 */
2540	.long 0x94FE13AB, 0xE8FA9A6E, 0xE06DB14A /* 207 */
2541	.long 0x29FC2757, 0xD1F534DD, 0xC0DB6295 /* 208 */
2542	.long 0x53F84EAF, 0xA3EA69BB, 0x81B6C52B /* 209 */
2543	.long 0xA7F09D5F, 0x47D4D377, 0x036D8A56 /* 210 */
2544	.long 0x4FE13ABE, 0x8FA9A6EE, 0x06DB14AC /* 211 */
2545	.long 0x9FC2757D, 0x1F534DDC, 0x0DB62959 /* 212 */
2546	.long 0x3F84EAFA, 0x3EA69BB8, 0x1B6C52B3 /* 213 */
2547	.long 0x7F09D5F4, 0x7D4D3770, 0x36D8A566 /* 214 */
2548	.long 0xFE13ABE8, 0xFA9A6EE0, 0x6DB14ACC /* 215 */
2549	.long 0xFC2757D1, 0xF534DDC0, 0xDB629599 /* 216 */
2550	.long 0xF84EAFA3, 0xEA69BB81, 0xB6C52B32 /* 217 */
2551	.long 0xF09D5F47, 0xD4D37703, 0x6D8A5664 /* 218 */
2552	.long 0xE13ABE8F, 0xA9A6EE06, 0xDB14ACC9 /* 219 */
2553	.long 0xC2757D1F, 0x534DDC0D, 0xB6295993 /* 220 */
2554	.long 0x84EAFA3E, 0xA69BB81B, 0x6C52B327 /* 221 */
2555	.long 0x09D5F47D, 0x4D377036, 0xD8A5664F /* 222 */
2556	.long 0x13ABE8FA, 0x9A6EE06D, 0xB14ACC9E /* 223 */
2557	.long 0x2757D1F5, 0x34DDC0DB, 0x6295993C /* 224 */
2558	.long 0x4EAFA3EA, 0x69BB81B6, 0xC52B3278 /* 225 */
2559	.long 0x9D5F47D4, 0xD377036D, 0x8A5664F1 /* 226 */
2560	.long 0x3ABE8FA9, 0xA6EE06DB, 0x14ACC9E2 /* 227 */
2561	.long 0x757D1F53, 0x4DDC0DB6, 0x295993C4 /* 228 */
2562	.long 0xEAFA3EA6, 0x9BB81B6C, 0x52B32788 /* 229 */
2563	.long 0xD5F47D4D, 0x377036D8, 0xA5664F10 /* 230 */
2564	.long 0xABE8FA9A, 0x6EE06DB1, 0x4ACC9E21 /* 231 */
2565	.long 0x57D1F534, 0xDDC0DB62, 0x95993C43 /* 232 */
2566	.long 0xAFA3EA69, 0xBB81B6C5, 0x2B327887 /* 233 */
2567	.long 0x5F47D4D3, 0x77036D8A, 0x5664F10E /* 234 */
2568	.long 0xBE8FA9A6, 0xEE06DB14, 0xACC9E21C /* 235 */
2569	.long 0x7D1F534D, 0xDC0DB629, 0x5993C439 /* 236 */
2570	.long 0xFA3EA69B, 0xB81B6C52, 0xB3278872 /* 237 */
2571	.long 0xF47D4D37, 0x7036D8A5, 0x664F10E4 /* 238 */
2572	.long 0xE8FA9A6E, 0xE06DB14A, 0xCC9E21C8 /* 239 */
2573	.long 0xD1F534DD, 0xC0DB6295, 0x993C4390 /* 240 */
2574	.long 0xA3EA69BB, 0x81B6C52B, 0x32788720 /* 241 */
2575	.long 0x47D4D377, 0x036D8A56, 0x64F10E41 /* 242 */
2576	.long 0x8FA9A6EE, 0x06DB14AC, 0xC9E21C82 /* 243 */
2577	.long 0x1F534DDC, 0x0DB62959, 0x93C43904 /* 244 */
2578	.long 0x3EA69BB8, 0x1B6C52B3, 0x27887208 /* 245 */
2579	.long 0x7D4D3770, 0x36D8A566, 0x4F10E410 /* 246 */
2580	.long 0xFA9A6EE0, 0x6DB14ACC, 0x9E21C820 /* 247 */
2581	.long 0xF534DDC0, 0xDB629599, 0x3C439041 /* 248 */
2582	.long 0xEA69BB81, 0xB6C52B32, 0x78872083 /* 249 */
2583	.long 0xD4D37703, 0x6D8A5664, 0xF10E4107 /* 250 */
2584	.long 0xA9A6EE06, 0xDB14ACC9, 0xE21C820F /* 251 */
2585	.long 0x534DDC0D, 0xB6295993, 0xC439041F /* 252 */
2586	.long 0xA69BB81B, 0x6C52B327, 0x8872083F /* 253 */
2587	.long 0x4D377036, 0xD8A5664F, 0x10E4107F /* 254 */
2588	.long 0x9A6EE06D, 0xB14ACC9E, 0x21C820FF /* 255 */
2589	.align 16
2590	.type __svml_stan_reduction_data_internal, @object
2591	.size __svml_stan_reduction_data_internal, .-__svml_stan_reduction_data_internal
2592	.align 16
2593
2594	.FLT_16:
2595	.long 0xffffffff, 0x00000000, 0xffffffff, 0x00000000
2596	.type .FLT_16, @object
2597	.size .FLT_16, 16
2598