/* Function tanf vectorized with SSE4. Copyright (C) 2021-2022 Free Software Foundation, Inc. This file is part of the GNU C Library. The GNU C Library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version. The GNU C Library is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details. You should have received a copy of the GNU Lesser General Public License along with the GNU C Library; if not, see https://www.gnu.org/licenses/. */ /* * ALGORITHM DESCRIPTION: * * 1) Range reduction to [-Pi/4; +Pi/4] interval * a) Grab sign from source argument and save it. * b) Remove sign using AND 0x7fffffff operation * c) Getting octant Y by 2/Pi multiplication * d) Add "Right Shifter" (0x4B000000) value * e) Treat obtained value as integer for destination sign setting. * Shift first bit of this value to the last (sign) position (S << 31) * f) Change destination sign if source sign is negative * using XOR operation. * g) Subtract "Right Shifter" (0x4B000000) value * h) Subtract Y*(PI/2) from X argument, where PI/2 divided to 4 parts: * X = X - Y*PI1 - Y*PI2 - Y*PI3 - Y*PI4; * 2) Rational polynomial approximation ( at [-Pi/4; +Pi/4] interval) * a) Calculate X^2 = X * X * b) Calculate 2 polynomials: * P = X * (P0 + X^2 * P1); * Q = Q0 + X^2 * (Q1 + x^2 * Q2); * c) Swap P and Q if first bit of obtained value after * Right Shifting is set to 1. Using And, Andnot & Or operations. * d) Divide R = P / Q; * 3) Destination sign setting * a) Set shifted destination sign using XOR operation: * R = XOR( R, S ); * */ /* Offsets for data table __svml_stan_data_internal */ #define _sInvPI_uisa 0 #define _sPI1_uisa 16 #define _sPI2_uisa 32 #define _sPI3_uisa 48 #define _sPI2_ha_uisa 64 #define _sPI3_ha_uisa 80 #define Th_tbl_uisa 96 #define Tl_tbl_uisa 224 #define _sPC3_uisa 352 #define _sPC5_uisa 368 #define _sRangeReductionVal_uisa 384 #define _sInvPi 400 #define _sSignMask 416 #define _sAbsMask 432 #define _sRangeVal 448 #define _sRShifter 464 #define _sOne 480 #define _sRangeReductionVal 496 #define _sPI1 512 #define _sPI2 528 #define _sPI3 544 #define _sPI4 560 #define _sPI1_FMA 576 #define _sPI2_FMA 592 #define _sPI3_FMA 608 #define _sP0 624 #define _sP1 640 #define _sQ0 656 #define _sQ1 672 #define _sQ2 688 #define _sTwo 704 #define _sCoeffs 720 #include .section .text.sse4, "ax", @progbits ENTRY(_ZGVbN4v_tanf_sse4) subq $232, %rsp cfi_def_cfa_offset(240) movaps %xmm0, %xmm13 movups _sAbsMask+__svml_stan_data_internal(%rip), %xmm12 /* * Legacy Code * Here HW FMA can be unavailable */ xorl %eax, %eax movaps %xmm12, %xmm4 pxor %xmm10, %xmm10 movups _sInvPi+__svml_stan_data_internal(%rip), %xmm2 andps %xmm13, %xmm4 mulps %xmm4, %xmm2 /* Range reduction */ movaps %xmm4, %xmm1 /* * * Main path (_LA_ and _EP_) * * Octant calculation */ movups _sRShifter+__svml_stan_data_internal(%rip), %xmm3 /* Large values check */ movaps %xmm4, %xmm11 movups _sPI1+__svml_stan_data_internal(%rip), %xmm5 andnps %xmm13, %xmm12 movups _sPI2+__svml_stan_data_internal(%rip), %xmm6 addps %xmm3, %xmm2 cmpnleps _sRangeReductionVal+__svml_stan_data_internal(%rip), %xmm11 movaps %xmm2, %xmm8 movups _sPI3+__svml_stan_data_internal(%rip), %xmm7 subps %xmm3, %xmm8 movmskps %xmm11, %edx movups _sPI4+__svml_stan_data_internal(%rip), %xmm9 mulps %xmm8, %xmm5 mulps %xmm8, %xmm6 mulps %xmm8, %xmm7 subps %xmm5, %xmm1 mulps %xmm8, %xmm9 subps %xmm6, %xmm1 movups _sQ2+__svml_stan_data_internal(%rip), %xmm15 /* Inversion mask and sign calculation */ movaps %xmm2, %xmm5 /* Rational approximation */ movups _sP1+__svml_stan_data_internal(%rip), %xmm14 pslld $30, %xmm2 cmpneqps %xmm10, %xmm2 subps %xmm7, %xmm1 /* Exchanged numerator and denominator if necessary */ movaps %xmm2, %xmm0 movaps %xmm2, %xmm10 pslld $31, %xmm5 subps %xmm9, %xmm1 movaps %xmm1, %xmm3 pxor %xmm12, %xmm5 mulps %xmm1, %xmm3 mulps %xmm3, %xmm15 mulps %xmm3, %xmm14 addps _sQ1+__svml_stan_data_internal(%rip), %xmm15 addps _sP0+__svml_stan_data_internal(%rip), %xmm14 mulps %xmm15, %xmm3 mulps %xmm14, %xmm1 addps _sQ0+__svml_stan_data_internal(%rip), %xmm3 andnps %xmm1, %xmm0 andps %xmm3, %xmm10 andps %xmm2, %xmm1 andnps %xmm3, %xmm2 orps %xmm10, %xmm0 orps %xmm2, %xmm1 /* Division */ divps %xmm1, %xmm0 /* Sign setting */ pxor %xmm5, %xmm0 /* * * End of main path (_LA_ and _EP_) */ testl %edx, %edx /* Go to auxilary branch */ jne L(AUX_BRANCH) # LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm4 xmm11 xmm12 xmm13 /* Return from auxilary branch * for out of main path inputs */ L(AUX_BRANCH_RETURN): testl %eax, %eax /* Go to special inputs processing branch */ jne L(SPECIAL_VALUES_BRANCH) # LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm13 /* Restore registers * and exit the function */ L(EXIT): addq $232, %rsp cfi_def_cfa_offset(8) ret cfi_def_cfa_offset(240) /* Branch to process * special inputs */ L(SPECIAL_VALUES_BRANCH): movups %xmm13, 32(%rsp) movups %xmm0, 48(%rsp) # LOE rbx rbp r12 r13 r14 r15 eax xmm0 xorl %edx, %edx movq %r12, 16(%rsp) cfi_offset(12, -224) movl %edx, %r12d movq %r13, 8(%rsp) cfi_offset(13, -232) movl %eax, %r13d movq %r14, (%rsp) cfi_offset(14, -240) # LOE rbx rbp r15 r12d r13d /* Range mask * bits check */ L(RANGEMASK_CHECK): btl %r12d, %r13d /* Call scalar math function */ jc L(SCALAR_MATH_CALL) # LOE rbx rbp r15 r12d r13d /* Special inputs * processing loop */ L(SPECIAL_VALUES_LOOP): incl %r12d cmpl $4, %r12d /* Check bits in range mask */ jl L(RANGEMASK_CHECK) # LOE rbx rbp r15 r12d r13d movq 16(%rsp), %r12 cfi_restore(12) movq 8(%rsp), %r13 cfi_restore(13) movq (%rsp), %r14 cfi_restore(14) movups 48(%rsp), %xmm0 /* Go to exit */ jmp L(EXIT) cfi_offset(12, -224) cfi_offset(13, -232) cfi_offset(14, -240) # LOE rbx rbp r12 r13 r14 r15 xmm0 /* Scalar math fucntion call * to process special input */ L(SCALAR_MATH_CALL): movl %r12d, %r14d movss 32(%rsp, %r14, 4), %xmm0 call tanf@PLT # LOE rbx rbp r14 r15 r12d r13d xmm0 movss %xmm0, 48(%rsp, %r14, 4) /* Process special inputs in loop */ jmp L(SPECIAL_VALUES_LOOP) cfi_restore(12) cfi_restore(13) cfi_restore(14) # LOE rbx rbp r15 r12d r13d /* Auxilary branch * for out of main path inputs */ L(AUX_BRANCH): movl $2139095040, %eax /* * Get the (2^a / 2pi) mod 1 values from the table. * Because doesn't have I-type gather, we need a trivial cast */ lea __svml_stan_reduction_data_internal(%rip), %r8 movups %xmm13, 64(%rsp) /* * Also get the significand as an integer * NB: adding in the integer bit is wrong for denorms! * To make this work for denorms we should do something slightly different */ movl $8388607, %r9d movups %xmm12, 80(%rsp) movl $8388608, %r10d movups %xmm11, 96(%rsp) /* * Break the P_xxx and m into 16-bit chunks ready for * the long multiplication via 16x16->32 multiplications */ movl $65535, %r11d movd %eax, %xmm3 pshufd $0, %xmm3, %xmm2 andps %xmm2, %xmm13 cmpeqps %xmm2, %xmm13 pand %xmm4, %xmm2 psrld $23, %xmm2 movdqa %xmm2, %xmm12 pslld $1, %xmm12 paddd %xmm2, %xmm12 pslld $2, %xmm12 pshufd $1, %xmm12, %xmm10 pshufd $2, %xmm12, %xmm11 pshufd $3, %xmm12, %xmm14 movd %xmm12, %edx movd %xmm10, %ecx movd %xmm11, %esi movd %r9d, %xmm11 movd %xmm14, %edi movd 4(%rdx, %r8), %xmm6 movd 4(%rcx, %r8), %xmm7 movd 4(%rsi, %r8), %xmm3 movl $872415232, %r9d movd 4(%rdi, %r8), %xmm5 punpckldq %xmm7, %xmm6 punpckldq %xmm5, %xmm3 movd 8(%rdi, %r8), %xmm10 movmskps %xmm13, %eax punpcklqdq %xmm3, %xmm6 movd 8(%rdx, %r8), %xmm3 movd 8(%rcx, %r8), %xmm2 movd 8(%rsi, %r8), %xmm13 punpckldq %xmm2, %xmm3 punpckldq %xmm10, %xmm13 punpcklqdq %xmm13, %xmm3 pshufd $0, %xmm11, %xmm13 movdqa %xmm3, %xmm2 movups %xmm4, 48(%rsp) pand %xmm4, %xmm13 movd %r10d, %xmm4 psrld $16, %xmm2 movd (%rdx, %r8), %xmm9 /* * We want to incorporate the original sign now too. * Do it here for convenience in getting the right N value, * though we could wait right to the end if we were prepared * to modify the sign of N later too. * So get the appropriate sign mask now (or sooner). */ movl $-2147483648, %edx movd (%rcx, %r8), %xmm8 /* * Create floating-point high part, implicitly adding integer bit 1 * Incorporate overall sign at this stage too. */ movl $1065353216, %ecx movd (%rsi, %r8), %xmm15 /* * Now round at the 2^-8 bit position for reduction mod pi/2^7 * instead of the original 2pi (but still with the same 2pi scaling). * Use a shifter of 2^15 + 2^14. * The N we get is our final version; it has an offset of * 2^8 because of the implicit integer bit, and anyway for negative * starting value it's a 2s complement thing. But we need to mask * off the exponent part anyway so it's fine. */ movl $1195376640, %esi movd (%rdi, %r8), %xmm1 movl $511, %r10d movups %xmm0, 112(%rsp) movd %r11d, %xmm0 pshufd $0, %xmm4, %xmm12 movdqa %xmm2, %xmm4 punpckldq %xmm8, %xmm9 paddd %xmm12, %xmm13 punpckldq %xmm1, %xmm15 movdqa %xmm13, %xmm12 pshufd $0, %xmm0, %xmm8 movdqa %xmm6, %xmm0 punpcklqdq %xmm15, %xmm9 pand %xmm8, %xmm13 movdqa %xmm9, %xmm14 pand %xmm8, %xmm9 movdqa %xmm13, %xmm10 psrld $16, %xmm14 movdqu %xmm14, 128(%rsp) /* Now do the big multiplication and carry propagation */ movdqa %xmm9, %xmm14 psrlq $32, %xmm10 psrlq $32, %xmm14 movdqa %xmm13, %xmm15 movdqa %xmm10, %xmm7 pmuludq %xmm9, %xmm15 psrld $16, %xmm0 pmuludq %xmm14, %xmm7 movdqu %xmm9, 144(%rsp) psllq $32, %xmm7 movdqu .FLT_16(%rip), %xmm9 pand %xmm8, %xmm6 pand %xmm9, %xmm15 psrld $16, %xmm12 movdqa %xmm0, %xmm1 por %xmm7, %xmm15 movdqa %xmm13, %xmm7 pand %xmm8, %xmm3 movdqu %xmm0, 160(%rsp) movdqa %xmm12, %xmm11 movdqu %xmm15, 208(%rsp) psrlq $32, %xmm1 pmuludq %xmm0, %xmm7 movdqa %xmm6, %xmm5 movdqa %xmm10, %xmm15 movdqa %xmm12, %xmm0 movdqu %xmm14, 176(%rsp) psrlq $32, %xmm11 movdqu %xmm1, 192(%rsp) psrlq $32, %xmm5 pmuludq %xmm1, %xmm15 movdqa %xmm13, %xmm1 pmuludq %xmm3, %xmm0 pmuludq %xmm6, %xmm1 pmuludq %xmm12, %xmm6 movdqa %xmm10, %xmm14 psrlq $32, %xmm3 pmuludq %xmm5, %xmm14 pand %xmm9, %xmm1 pmuludq %xmm11, %xmm3 pmuludq %xmm11, %xmm5 psllq $32, %xmm14 pand %xmm9, %xmm0 psllq $32, %xmm3 psrlq $32, %xmm4 por %xmm14, %xmm1 por %xmm3, %xmm0 movdqa %xmm12, %xmm14 movdqa %xmm11, %xmm3 pmuludq %xmm2, %xmm14 pand %xmm9, %xmm7 pmuludq %xmm4, %xmm3 pmuludq %xmm13, %xmm2 pmuludq %xmm10, %xmm4 pand %xmm9, %xmm2 psllq $32, %xmm4 psllq $32, %xmm15 pand %xmm9, %xmm14 psllq $32, %xmm3 por %xmm4, %xmm2 por %xmm15, %xmm7 por %xmm3, %xmm14 psrld $16, %xmm2 pand %xmm9, %xmm6 psllq $32, %xmm5 movdqa %xmm1, %xmm15 paddd %xmm2, %xmm14 movdqa %xmm7, %xmm2 por %xmm5, %xmm6 psrld $16, %xmm1 pand %xmm8, %xmm2 paddd %xmm1, %xmm6 movdqu 160(%rsp), %xmm1 paddd %xmm6, %xmm2 movdqu 192(%rsp), %xmm6 psrld $16, %xmm7 pmuludq %xmm12, %xmm1 pand %xmm8, %xmm15 pmuludq %xmm11, %xmm6 pmuludq 144(%rsp), %xmm12 pmuludq 176(%rsp), %xmm11 pand %xmm9, %xmm1 psllq $32, %xmm6 por %xmm6, %xmm1 psrld $16, %xmm0 paddd %xmm7, %xmm1 paddd %xmm14, %xmm15 movdqu 128(%rsp), %xmm7 paddd %xmm15, %xmm0 pmuludq %xmm7, %xmm13 psrlq $32, %xmm7 pmuludq %xmm7, %xmm10 movdqa %xmm0, %xmm14 pand %xmm9, %xmm13 movdqu 208(%rsp), %xmm5 psrld $16, %xmm14 paddd %xmm2, %xmm14 movdqa %xmm5, %xmm15 movdqa %xmm14, %xmm3 pand %xmm8, %xmm15 psrld $16, %xmm3 paddd %xmm1, %xmm15 psllq $32, %xmm10 pand %xmm9, %xmm12 psllq $32, %xmm11 paddd %xmm15, %xmm3 por %xmm10, %xmm13 por %xmm11, %xmm12 psrld $16, %xmm5 movdqa %xmm3, %xmm4 pand %xmm8, %xmm13 paddd %xmm5, %xmm12 psrld $16, %xmm4 paddd %xmm12, %xmm13 paddd %xmm13, %xmm4 pand %xmm8, %xmm3 pslld $16, %xmm4 movd %edx, %xmm9 movups 48(%rsp), %xmm15 paddd %xmm3, %xmm4 pshufd $0, %xmm9, %xmm7 /* Assemble reduced argument from the pieces */ pand %xmm8, %xmm0 movd %ecx, %xmm8 pand %xmm15, %xmm7 pshufd $0, %xmm8, %xmm1 movdqa %xmm4, %xmm5 psrld $9, %xmm5 pxor %xmm7, %xmm1 por %xmm1, %xmm5 movd %esi, %xmm6 pshufd $0, %xmm6, %xmm3 movdqa %xmm5, %xmm6 movl $262143, %r8d /* * Create floating-point low and medium parts, respectively * lo_17, ... lo_0, 0, ..., 0 * hi_8, ... hi_0, lo_31, ..., lo_18 * then subtract off the implicitly added integer bits, * 2^-46 and 2^-23, respectively. * Put the original sign into all of them at this stage. */ movl $679477248, %edi movd %r10d, %xmm13 pslld $16, %xmm14 pshufd $0, %xmm13, %xmm1 paddd %xmm0, %xmm14 movd %r9d, %xmm11 pand %xmm4, %xmm1 movd %r8d, %xmm9 movd %edi, %xmm10 pshufd $0, %xmm9, %xmm8 pslld $14, %xmm1 pshufd $0, %xmm10, %xmm0 pand %xmm14, %xmm8 pshufd $0, %xmm11, %xmm12 psrld $18, %xmm14 pxor %xmm7, %xmm0 pxor %xmm12, %xmm7 por %xmm14, %xmm1 pslld $5, %xmm8 por %xmm7, %xmm1 /* * Now multiply those numbers all by 2 pi, reasonably accurately. * The top part uses 2pi = s2pi_lead + s2pi_trail, where * s2pi_lead has 12 significant bits. */ movl $1086918619, %r11d /* Split RHi into 12-bit leading and trailing parts. */ movl $-4096, %esi por %xmm0, %xmm8 movl $1086918656, %edx movl $-1214941318, %ecx /* * If the magnitude of the input is <= 2^-20, then * just pass through the input, since no reduction will be needed and * the main path will only work accurately if the reduced argument is * about >= 2^-40 (which it is for all large pi multiples) */ movl $2147483647, %edi addps %xmm3, %xmm6 subps %xmm7, %xmm1 subps %xmm0, %xmm8 movaps %xmm6, %xmm2 movd %r11d, %xmm14 movd %esi, %xmm4 movd %edx, %xmm7 movl $897581056, %r8d subps %xmm3, %xmm2 /* Grab our final N value as an integer, appropriately masked mod 2^8 */ movl $255, %r9d subps %xmm2, %xmm5 /* Now add them up into 2 reasonably aligned pieces */ movaps %xmm5, %xmm3 /* * The output is _VRES_R (high) + _VRES_E (low), and the integer part is _VRES_IND * Set sRp2 = _VRES_R^2 and then resume the original code. * Argument reduction is now finished: x = n * pi/128 + r * where n = iIndex and r = sR (high) + sE (low). * But we have n modulo 256, needed for sin/cos with period 2pi * but we want it modulo 128 since tan has period pi. */ movl $127, %r10d pshufd $0, %xmm14, %xmm2 addps %xmm1, %xmm3 pshufd $0, %xmm4, %xmm14 movd %r8d, %xmm4 pshufd $0, %xmm4, %xmm9 subps %xmm3, %xmm5 movdqa %xmm9, %xmm11 addps %xmm5, %xmm1 movd %ecx, %xmm5 addps %xmm1, %xmm8 pshufd $0, %xmm7, %xmm1 movdqa %xmm14, %xmm7 andps %xmm3, %xmm7 /* * Do the multiplication as exact top part and "naive" low part. * This still maintains a similar level of offset and doesn't drop * the accuracy much below what we already have. */ movdqa %xmm1, %xmm10 pshufd $0, %xmm5, %xmm5 subps %xmm7, %xmm3 mulps %xmm7, %xmm10 mulps %xmm5, %xmm7 mulps %xmm3, %xmm1 mulps %xmm8, %xmm2 mulps %xmm3, %xmm5 addps %xmm7, %xmm1 addps %xmm5, %xmm2 movd %edi, %xmm8 addps %xmm2, %xmm1 /* * Do another stage of compensated summation to get full offset * between the pieces sRedHi + sRedLo. * Depending on the later algorithm, we might avoid this stage. */ movaps %xmm1, %xmm0 /* Load constants (not all needed at once) */ lea _sCoeffs+36+__svml_stan_data_internal(%rip), %rdi pshufd $0, %xmm8, %xmm8 addps %xmm10, %xmm0 andps %xmm15, %xmm8 subps %xmm0, %xmm10 cmpltps %xmm8, %xmm11 cmpleps %xmm9, %xmm8 addps %xmm10, %xmm1 andps %xmm15, %xmm8 movd %r9d, %xmm15 andps %xmm11, %xmm0 andps %xmm1, %xmm11 pshufd $0, %xmm15, %xmm1 movd %r10d, %xmm15 pshufd $0, %xmm15, %xmm7 pand %xmm1, %xmm6 pand %xmm7, %xmm6 orps %xmm0, %xmm8 movaps %xmm6, %xmm4 /* * Simply combine the two parts of the reduced argument * since we can afford a few ulps in this case. */ addps %xmm11, %xmm8 pslld $2, %xmm4 paddd %xmm6, %xmm4 pslld $3, %xmm4 pshufd $1, %xmm4, %xmm6 pshufd $2, %xmm4, %xmm5 pshufd $3, %xmm4, %xmm3 movd %xmm4, %r11d movd %xmm6, %edx movd %xmm5, %ecx movd %xmm3, %esi movd -32(%r11, %rdi), %xmm15 movd -32(%rdx, %rdi), %xmm12 movd -32(%rcx, %rdi), %xmm7 movd -32(%rsi, %rdi), %xmm13 punpckldq %xmm12, %xmm15 punpckldq %xmm13, %xmm7 movd -28(%rsi, %rdi), %xmm5 punpcklqdq %xmm7, %xmm15 movd -28(%r11, %rdi), %xmm7 movd -28(%rdx, %rdi), %xmm6 movd -28(%rcx, %rdi), %xmm4 movd -36(%rcx, %rdi), %xmm9 movd -36(%r11, %rdi), %xmm1 movd -36(%rdx, %rdi), %xmm2 movd -24(%rdx, %rdi), %xmm3 movd -36(%rsi, %rdi), %xmm10 punpckldq %xmm6, %xmm7 punpckldq %xmm5, %xmm4 movd -24(%r11, %rdi), %xmm6 punpckldq %xmm2, %xmm1 punpckldq %xmm10, %xmm9 punpcklqdq %xmm4, %xmm7 movd -16(%r11, %rdi), %xmm4 punpckldq %xmm3, %xmm6 movd -24(%rcx, %rdi), %xmm10 movd -16(%rcx, %rdi), %xmm3 movd -24(%rsi, %rdi), %xmm2 movd -16(%rsi, %rdi), %xmm13 movd -16(%rdx, %rdi), %xmm12 punpcklqdq %xmm9, %xmm1 movd -20(%rdx, %rdi), %xmm9 punpckldq %xmm2, %xmm10 movd -20(%r11, %rdi), %xmm5 movd -20(%rcx, %rdi), %xmm11 movd -20(%rsi, %rdi), %xmm0 punpckldq %xmm12, %xmm4 punpckldq %xmm13, %xmm3 punpcklqdq %xmm10, %xmm6 movd -12(%rsi, %rdi), %xmm10 punpckldq %xmm9, %xmm5 punpckldq %xmm0, %xmm11 punpcklqdq %xmm3, %xmm4 movd -12(%r11, %rdi), %xmm3 movd -12(%rdx, %rdi), %xmm2 movd -12(%rcx, %rdi), %xmm9 punpcklqdq %xmm11, %xmm5 punpckldq %xmm2, %xmm3 punpckldq %xmm10, %xmm9 movd -8(%rcx, %rdi), %xmm10 movd -8(%r11, %rdi), %xmm2 movd -8(%rdx, %rdi), %xmm0 movd -8(%rsi, %rdi), %xmm11 punpckldq %xmm0, %xmm2 punpckldq %xmm11, %xmm10 movd -4(%rsi, %rdi), %xmm13 punpcklqdq %xmm9, %xmm3 punpcklqdq %xmm10, %xmm2 movd -4(%r11, %rdi), %xmm10 movd -4(%rdx, %rdi), %xmm12 movd -4(%rcx, %rdi), %xmm9 punpckldq %xmm12, %xmm10 punpckldq %xmm13, %xmm9 punpcklqdq %xmm9, %xmm10 /* * Compute 2-part reciprocal component * Construct a separate reduced argument modulo pi near pi/2 multiples. * i.e. (pi/2 - x) mod pi, simply by subtracting the reduced argument * from an accurate B_hi + B_lo = (128 - n) pi/128. Force the upper part * of this reduced argument to half-length to simplify accurate * reciprocation later on. */ movdqa %xmm1, %xmm9 movd (%r11, %rdi), %xmm13 subps %xmm8, %xmm9 movd (%rdx, %rdi), %xmm0 subps %xmm9, %xmm1 punpckldq %xmm0, %xmm13 movdqa %xmm14, %xmm0 andps %xmm9, %xmm0 subps %xmm8, %xmm1 subps %xmm0, %xmm9 movd (%rcx, %rdi), %xmm12 addps %xmm9, %xmm15 /* * Now compute an approximate reciprocal to mix into the computation * To avoid any danger of nonportability, force it to 12 bits, * though I suspect it always is anyway on current platforms. */ rcpps %xmm0, %xmm9 addps %xmm15, %xmm1 andps %xmm14, %xmm9 mulps %xmm9, %xmm0 /* * Get a better approximation to 1/sR_hi (not far short of an ulp) * using a third-order polynomial approximation */ movaps %xmm9, %xmm14 movd (%rsi, %rdi), %xmm11 /* * Now compute the error sEr where sRecip_hi = (1/R_hi) * (1 - sEr) * so that we can compensate for it. */ movups _sOne+__svml_stan_data_internal(%rip), %xmm15 punpckldq %xmm11, %xmm12 movaps %xmm15, %xmm11 punpcklqdq %xmm12, %xmm13 subps %xmm0, %xmm11 mulps %xmm11, %xmm14 movups %xmm11, (%rsp) addps %xmm9, %xmm14 mulps %xmm11, %xmm11 movups %xmm13, 32(%rsp) movups %xmm11, 16(%rsp) movups 112(%rsp), %xmm0 movups 96(%rsp), %xmm11 movups 80(%rsp), %xmm12 movups 64(%rsp), %xmm13 # LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm1 xmm2 xmm3 xmm4 xmm5 xmm6 xmm7 xmm8 xmm9 xmm10 xmm11 xmm12 xmm13 xmm14 xmm15 /* * Compensated sum of dominant component(s) * Compute C0_hi + C1_hi * Z + Recip_hi + Recip_lo = H4 (hi) + H9 (lo) * H1 = C1_hi * Z (exact since C1_hi is 1 bit) */ mulps %xmm8, %xmm4 addps 16(%rsp), %xmm15 /* Finally, multiplex both parts so they are only used in cotangent path */ mulps %xmm7, %xmm9 /* * Higher polynomial terms * Stage 1 (with unlimited parallelism) * P3 = C1_lo + C2 * Z */ mulps %xmm8, %xmm2 mulps %xmm15, %xmm14 addps %xmm2, %xmm3 /* * Multiply by sRecip_ok to make sR_lo relative to sR_hi * Since sR_lo is shifted off by about 12 bits, this is accurate enough. */ mulps %xmm14, %xmm1 /* * Now create a low reciprocal using * (Recip_hi + Er * Recip_ok) * (1 + sR_lo^2 - sR_lo) * =~= Recip_hi + Recip_ok * (Er + sR_lo^2 - sR_lo) */ movaps %xmm1, %xmm15 mulps %xmm1, %xmm1 subps (%rsp), %xmm15 /* P4 = C3 + C4 * Z */ movups 32(%rsp), %xmm2 subps %xmm15, %xmm1 mulps %xmm8, %xmm2 mulps %xmm1, %xmm14 addps %xmm2, %xmm10 mulps %xmm14, %xmm7 /* H2 = high(C0_hi + C1_hi * Z) */ movdqa %xmm6, %xmm14 addps %xmm4, %xmm14 /* H4 = high(H2 + Recip_hi) */ movaps %xmm14, %xmm1 /* intermediate in compensated sum */ subps %xmm14, %xmm6 addps %xmm9, %xmm1 /* H5 = low(C0_hi + C1_hi * Z) */ addps %xmm6, %xmm4 /* intermediate in compensated sum */ subps %xmm1, %xmm9 /* H7 = low(C0_hi + C1_hi * Z) + Recip_lo */ addps %xmm4, %xmm7 /* H8 = low(H2 + Recip_hi) */ addps %xmm9, %xmm14 /* Z2 = Z^2 */ movaps %xmm8, %xmm4 /* Now H4 + H9 should be that part */ addps %xmm14, %xmm7 mulps %xmm8, %xmm4 /* P9 = trail(dominant part) + C0_lo */ addps %xmm7, %xmm5 /* * Stage 2 (with unlimited parallelism) * P6 = C1_lo + C2 * Z + C3 * Z^2 + C4 * Z^3 */ mulps %xmm4, %xmm10 addps %xmm10, %xmm3 /* Final accumulation of low part */ mulps %xmm3, %xmm8 /* Merge results from main and large paths: */ movaps %xmm11, %xmm3 andnps %xmm0, %xmm3 addps %xmm8, %xmm5 movaps %xmm3, %xmm0 /* And now the very final summation */ addps %xmm5, %xmm1 /* * The end of implementation (LA with huge args reduction) * End of large arguments path (_HA_, _LA_ and _EP_) */ pxor %xmm12, %xmm1 andps %xmm11, %xmm1 orps %xmm1, %xmm0 /* Return to main vector processing path */ jmp L(AUX_BRANCH_RETURN) # LOE rbx rbp r12 r13 r14 r15 eax xmm0 xmm13 END(_ZGVbN4v_tanf_sse4) .section .rodata, "a" .align 16 #ifdef __svml_stan_data_internal_typedef typedef unsigned int VUINT32; typedef struct { __declspec(align(16)) VUINT32 _sInvPI_uisa[4][1]; __declspec(align(16)) VUINT32 _sPI1_uisa[4][1]; __declspec(align(16)) VUINT32 _sPI2_uisa[4][1]; __declspec(align(16)) VUINT32 _sPI3_uisa[4][1]; __declspec(align(16)) VUINT32 _sPI2_ha_uisa[4][1]; __declspec(align(16)) VUINT32 _sPI3_ha_uisa[4][1]; __declspec(align(16)) VUINT32 Th_tbl_uisa[32][1]; __declspec(align(16)) VUINT32 Tl_tbl_uisa[32][1]; __declspec(align(16)) VUINT32 _sPC3_uisa[4][1]; __declspec(align(16)) VUINT32 _sPC5_uisa[4][1]; __declspec(align(16)) VUINT32 _sRangeReductionVal_uisa[4][1]; __declspec(align(16)) VUINT32 _sInvPi[4][1]; __declspec(align(16)) VUINT32 _sSignMask[4][1]; __declspec(align(16)) VUINT32 _sAbsMask[4][1]; __declspec(align(16)) VUINT32 _sRangeVal[4][1]; __declspec(align(16)) VUINT32 _sRShifter[4][1]; __declspec(align(16)) VUINT32 _sOne[4][1]; __declspec(align(16)) VUINT32 _sRangeReductionVal[4][1]; __declspec(align(16)) VUINT32 _sPI1[4][1]; __declspec(align(16)) VUINT32 _sPI2[4][1]; __declspec(align(16)) VUINT32 _sPI3[4][1]; __declspec(align(16)) VUINT32 _sPI4[4][1]; __declspec(align(16)) VUINT32 _sPI1_FMA[4][1]; __declspec(align(16)) VUINT32 _sPI2_FMA[4][1]; __declspec(align(16)) VUINT32 _sPI3_FMA[4][1]; __declspec(align(16)) VUINT32 _sP0[4][1]; __declspec(align(16)) VUINT32 _sP1[4][1]; __declspec(align(16)) VUINT32 _sQ0[4][1]; __declspec(align(16)) VUINT32 _sQ1[4][1]; __declspec(align(16)) VUINT32 _sQ2[4][1]; __declspec(align(16)) VUINT32 _sTwo[4][1]; __declspec(align(16)) VUINT32 _sCoeffs[128][10][1]; } __svml_stan_data_internal; #endif __svml_stan_data_internal: /* UISA */ .long 0x4122f983, 0x4122f983, 0x4122f983, 0x4122f983 /* _sInvPI_uisa */ .align 16 .long 0x3dc90fda, 0x3dc90fda, 0x3dc90fda, 0x3dc90fda /* _sPI1_uisa */ .align 16 .long 0x31a22168, 0x31a22168, 0x31a22168, 0x31a22168 /* _sPI2_uisa */ .align 16 .long 0x25c234c5, 0x25c234c5, 0x25c234c5, 0x25c234c5 /* _sPI3_uisa */ .align 16 .long 0x31a22000, 0x31a22000, 0x31a22000, 0x31a22000 /* _sPI2_ha_uisa */ .align 16 .long 0x2a34611a, 0x2a34611a, 0x2a34611a, 0x2a34611a /* _sPI3_ha_uisa */ /* Th_tbl_uisa for i from 0 to 31 do printsingle(tan(i*Pi/32)); */ .align 16 .long 0x80000000, 0x3dc9b5dc, 0x3e4bafaf, 0x3e9b5042 .long 0x3ed413cd, 0x3f08d5b9, 0x3f2b0dc1, 0x3f521801 .long 0x3f800000, 0x3f9bf7ec, 0x3fbf90c7, 0x3fef789e .long 0x401a827a, 0x4052facf, 0x40a0dff7, 0x41227363 .long 0xff7fffff, 0xc1227363, 0xc0a0dff7, 0xc052facf .long 0xc01a827a, 0xbfef789e, 0xbfbf90c7, 0xbf9bf7ec .long 0xbf800000, 0xbf521801, 0xbf2b0dc1, 0xbf08d5b9 .long 0xbed413cd, 0xbe9b5042, 0xbe4bafaf, 0xbdc9b5dc /* Tl_tbl_uisa for i from 0 to 31 do printsingle(tan(i*Pi/32)-round(tan(i*Pi/32), SG, RN)); */ .align 16 .long 0x80000000, 0x3145b2da, 0x2f2a62b0, 0xb22a39c2 .long 0xb1c0621a, 0xb25ef963, 0x32ab7f99, 0x32ae4285 .long 0x00000000, 0x33587608, 0x32169d18, 0xb30c3ec0 .long 0xb3cc0622, 0x3390600e, 0x331091dc, 0xb454a046 .long 0xf3800000, 0x3454a046, 0xb31091dc, 0xb390600e .long 0x33cc0622, 0x330c3ec0, 0xb2169d18, 0xb3587608 .long 0x00000000, 0xb2ae4285, 0xb2ab7f99, 0x325ef963 .long 0x31c0621a, 0x322a39c2, 0xaf2a62b0, 0xb145b2da .align 16 .long 0x3eaaaaa6, 0x3eaaaaa6, 0x3eaaaaa6, 0x3eaaaaa6 /* _sPC3_uisa */ .align 16 .long 0x3e08b888, 0x3e08b888, 0x3e08b888, 0x3e08b888 /* _sPC5_uisa */ .align 16 .long 0x46010000, 0x46010000, 0x46010000, 0x46010000 /* _sRangeReductionVal_uisa */ .align 16 .long 0x3F22F983, 0x3F22F983, 0x3F22F983, 0x3F22F983 /* _sInvPi */ .align 16 .long 0x80000000, 0x80000000, 0x80000000, 0x80000000 /* _sSignMask */ .align 16 .long 0x7FFFFFFF, 0x7FFFFFFF, 0x7FFFFFFF, 0x7FFFFFFF /* _sAbsMask */ .align 16 .long 0x7f800000, 0x7f800000, 0x7f800000, 0x7f800000 /* _sRangeVal */ .align 16 .long 0x4B400000, 0x4B400000, 0x4B400000, 0x4B400000 /* _sRShifter */ .align 16 .long 0x3f800000, 0x3f800000, 0x3f800000, 0x3f800000 /* _sOne */ .align 16 .long 0x46010000, 0x46010000, 0x46010000, 0x46010000 /* _sRangeVal */ .align 16 .long 0x3FC90000, 0x3FC90000, 0x3FC90000, 0x3FC90000 /* _sPI1 */ .align 16 .long 0x39FDA000, 0x39FDA000, 0x39FDA000, 0x39FDA000 /* _sPI2 */ .align 16 .long 0x33A22000, 0x33A22000, 0x33A22000, 0x33A22000 /* _sPI3 */ .align 16 .long 0x2C34611A, 0x2C34611A, 0x2C34611A, 0x2C34611A /* _sPI4 */ // PI1, PI2, and PI3 when FMA is available .align 16 .long 0x3FC90FDB, 0x3FC90FDB, 0x3FC90FDB, 0x3FC90FDB /* _sPI1_FMA */ .align 16 .long 0xB33BBD2E, 0xB33BBD2E, 0xB33BBD2E, 0xB33BBD2E /* _sPI2_FMA */ .align 16 .long 0xA6F72CED, 0xA6F72CED, 0xA6F72CED, 0xA6F72CED /* _sPI3_FMA */ .align 16 .long 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC /* _sP0 */ .align 16 .long 0xBDC433B4, 0xBDC433B4, 0xBDC433B4, 0xBDC433B4 /* _sP1 */ .align 16 .long 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC, 0x3F7FFFFC /* _sQ0 */ .align 16 .long 0xBEDBB7AB, 0xBEDBB7AB, 0xBEDBB7AB, 0xBEDBB7AB /* _sQ1 */ .align 16 .long 0x3C1F336B, 0x3C1F336B, 0x3C1F336B, 0x3C1F336B /* _sQ2 */ .align 16 .long 0x40000000, 0x40000000, 0x40000000, 0x40000000 /* _sTwo */ // _sCoeffs Breakpoint B = 0 * pi/128, function tan(B + x) .align 16 .long 0x3FC90FDB // B' = pi/2 - B (high single) .long 0xB33BBD2E // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x00000000 // c0 (high single) .long 0x00000000 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x00000000 // c1 (low single) .long 0x00000000 // c2 .long 0x3EAAACDD // c3 .long 0x00000000 // c4 .long 0x3FC5EB9B // B' = pi/2 - B (high single) .long 0x32DE638C // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3CC91A31 // c0 (high single) .long 0x2F8E8D1A // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3A1DFA00 // c1 (low single) .long 0x3CC9392D // c2 .long 0x3EAB1889 // c3 .long 0x3C885D3B // c4 .long 0x3FC2C75C // B' = pi/2 - B (high single) .long 0xB2CBBE8A // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3D49393C // c0 (high single) .long 0x30A39F5B // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3B1E2B00 // c1 (low single) .long 0x3D49B5D4 // c2 .long 0x3EAC4F10 // c3 .long 0x3CFD9425 // c4 .long 0x3FBFA31C // B' = pi/2 - B (high single) .long 0x33450FB0 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3D9711CE // c0 (high single) .long 0x314FEB28 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3BB24C00 // c1 (low single) .long 0x3D97E43A // c2 .long 0x3EAE6A89 // c3 .long 0x3D4D07E0 // c4 .long 0x3FBC7EDD // B' = pi/2 - B (high single) .long 0xB1800ADD // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3DC9B5DC // c0 (high single) .long 0x3145AD86 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3C1EEF20 // c1 (low single) .long 0x3DCBAAEA // c2 .long 0x3EB14E5E // c3 .long 0x3D858BB2 // c4 .long 0x3FB95A9E // B' = pi/2 - B (high single) .long 0xB3651267 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3DFC98C2 // c0 (high single) .long 0xB0AE525C // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3C793D20 // c1 (low single) .long 0x3E003845 // c2 .long 0x3EB5271F // c3 .long 0x3DAC669E // c4 .long 0x3FB6365E // B' = pi/2 - B (high single) .long 0x328BB91C // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E17E564 // c0 (high single) .long 0xB1C5A2E4 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3CB440D0 // c1 (low single) .long 0x3E1B3D00 // c2 .long 0x3EB9F664 // c3 .long 0x3DD647C0 // c4 .long 0x3FB3121F // B' = pi/2 - B (high single) .long 0xB30F347D // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E31AE4D // c0 (high single) .long 0xB1F32251 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3CF6A500 // c1 (low single) .long 0x3E3707DA // c2 .long 0x3EBFA489 // c3 .long 0x3DFBD9C7 // c4 .long 0x3FAFEDDF // B' = pi/2 - B (high single) .long 0x331BBA77 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E4BAFAF // c0 (high single) .long 0x2F2A29E0 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D221018 // c1 (low single) .long 0x3E53BED0 // c2 .long 0x3EC67E26 // c3 .long 0x3E1568E2 // c4 .long 0x3FACC9A0 // B' = pi/2 - B (high single) .long 0xB2655A50 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E65F267 // c0 (high single) .long 0x31B4B1DF // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D4E8B90 // c1 (low single) .long 0x3E718ACA // c2 .long 0x3ECE7164 // c3 .long 0x3E2DC161 // c4 .long 0x3FA9A560 // B' = pi/2 - B (high single) .long 0x33719861 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E803FD4 // c0 (high single) .long 0xB2279E66 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D807FC8 // c1 (low single) .long 0x3E884BD4 // c2 .long 0x3ED7812D // c3 .long 0x3E4636EB // c4 .long 0x3FA68121 // B' = pi/2 - B (high single) .long 0x31E43AAC // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E8DB082 // c0 (high single) .long 0xB132A234 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D9CD7D0 // c1 (low single) .long 0x3E988A60 // c2 .long 0x3EE203E3 // c3 .long 0x3E63582C // c4 .long 0x3FA35CE2 // B' = pi/2 - B (high single) .long 0xB33889B6 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3E9B5042 // c0 (high single) .long 0xB22A3AEE // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3DBC7490 // c1 (low single) .long 0x3EA99AF5 // c2 .long 0x3EEDE107 // c3 .long 0x3E80E9AA // c4 .long 0x3FA038A2 // B' = pi/2 - B (high single) .long 0x32E4CA7E // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3EA92457 // c0 (high single) .long 0x30B80830 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3DDF8200 // c1 (low single) .long 0x3EBB99E9 // c2 .long 0x3EFB4AA8 // c3 .long 0x3E9182BE // c4 .long 0x3F9D1463 // B' = pi/2 - B (high single) .long 0xB2C55799 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3EB73250 // c0 (high single) .long 0xB2028823 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E0318F8 // c1 (low single) .long 0x3ECEA678 // c2 .long 0x3F053C67 // c3 .long 0x3EA41E53 // c4 .long 0x3F99F023 // B' = pi/2 - B (high single) .long 0x33484328 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3EC5800D // c0 (high single) .long 0xB214C3C1 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E185E54 // c1 (low single) .long 0x3EE2E342 // c2 .long 0x3F0DCA73 // c3 .long 0x3EB8CC21 // c4 .long 0x3F96CBE4 // B' = pi/2 - B (high single) .long 0xB14CDE2E // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3ED413CD // c0 (high single) .long 0xB1C06152 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E2FB0CC // c1 (low single) .long 0x3EF876CB // c2 .long 0x3F177807 // c3 .long 0x3ED08437 // c4 .long 0x3F93A7A5 // B' = pi/2 - B (high single) .long 0xB361DEEE // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3EE2F439 // c0 (high single) .long 0xB1F4399E // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E49341C // c1 (low single) .long 0x3F07C61A // c2 .long 0x3F22560F // c3 .long 0x3EEAA81E // c4 .long 0x3F908365 // B' = pi/2 - B (high single) .long 0x3292200D // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3EF22870 // c0 (high single) .long 0x325271F4 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E65107A // c1 (low single) .long 0x3F1429F0 // c2 .long 0x3F2E8AFC // c3 .long 0x3F040498 // c4 .long 0x3F8D5F26 // B' = pi/2 - B (high single) .long 0xB30C0105 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F00DC0D // c0 (high single) .long 0xB214AF72 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E81B994 // c1 (low single) .long 0x3F218233 // c2 .long 0x3F3C4531 // c3 .long 0x3F149688 // c4 .long 0x3F8A3AE6 // B' = pi/2 - B (high single) .long 0x331EEDF0 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F08D5B9 // c0 (high single) .long 0xB25EF98E // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E92478D // c1 (low single) .long 0x3F2FEDC9 // c2 .long 0x3F4BCD58 // c3 .long 0x3F27AE9E // c4 .long 0x3F8716A7 // B' = pi/2 - B (high single) .long 0xB2588C6D // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F1105AF // c0 (high single) .long 0x32F045B0 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EA44EE2 // c1 (low single) .long 0x3F3F8FDB // c2 .long 0x3F5D3FD0 // c3 .long 0x3F3D0A23 // c4 .long 0x3F83F267 // B' = pi/2 - B (high single) .long 0x3374CBD9 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F1970C4 // c0 (high single) .long 0x32904848 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EB7EFF8 // c1 (low single) .long 0x3F50907C // c2 .long 0x3F710FEA // c3 .long 0x3F561FED // c4 .long 0x3F80CE28 // B' = pi/2 - B (high single) .long 0x31FDD672 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F221C37 // c0 (high single) .long 0xB20C61DC // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3ECD4F71 // c1 (low single) .long 0x3F631DAA // c2 .long 0x3F83B471 // c3 .long 0x3F7281EA // c4 .long 0x3F7B53D1 // B' = pi/2 - B (high single) .long 0x32955386 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F2B0DC1 // c0 (high single) .long 0x32AB7EBA // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EE496C2 // c1 (low single) .long 0x3F776C40 // c2 .long 0x3F9065C1 // c3 .long 0x3F89AFB6 // c4 .long 0x3F750B52 // B' = pi/2 - B (high single) .long 0x32EB316F // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F344BA9 // c0 (high single) .long 0xB2B8B0EA // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EFDF4F7 // c1 (low single) .long 0x3F86DCA8 // c2 .long 0x3F9ED53B // c3 .long 0x3F9CBEDE // c4 .long 0x3F6EC2D4 // B' = pi/2 - B (high single) .long 0xB2BEF0A7 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F3DDCCF // c0 (high single) .long 0x32D29606 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBEE6606F // c1 (low single) .long 0x3F9325D6 // c2 .long 0x3FAF4E69 // c3 .long 0x3FB3080C // c4 .long 0x3F687A55 // B' = pi/2 - B (high single) .long 0xB252257B // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F47C8CC // c0 (high single) .long 0xB200F51A // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBEC82C6C // c1 (low single) .long 0x3FA0BAE9 // c2 .long 0x3FC2252F // c3 .long 0x3FCD24C7 // c4 .long 0x3F6231D6 // B' = pi/2 - B (high single) .long 0xB119A6A2 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F521801 // c0 (high single) .long 0x32AE4178 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBEA72938 // c1 (low single) .long 0x3FAFCC22 // c2 .long 0x3FD7BD4A // c3 .long 0x3FEBB01B // c4 .long 0x3F5BE957 // B' = pi/2 - B (high single) .long 0x3205522A // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F5CD3BE // c0 (high single) .long 0x31460308 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBE8306C5 // c1 (low single) .long 0x3FC09232 // c2 .long 0x3FF09632 // c3 .long 0x4007DB00 // c4 .long 0x3F55A0D8 // B' = pi/2 - B (high single) .long 0x329886FF // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F68065E // c0 (high single) .long 0x32670D1A // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBE36D1D6 // c1 (low single) .long 0x3FD35007 // c2 .long 0x4006A861 // c3 .long 0x401D4BDA // c4 .long 0x3F4F5859 // B' = pi/2 - B (high single) .long 0x32EE64E8 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0x3F73BB75 // c0 (high single) .long 0x32FC908D // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBDBF94B0 // c1 (low single) .long 0x3FE8550F // c2 .long 0x40174F67 // c3 .long 0x4036C608 // c4 .long 0x3F490FDB // B' = pi/2 - B (high single) .long 0xB2BBBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE8BE60E // c0 (high single) .long 0x320D8D84 // c0 (low single) .long 0x3F000000 // c1 (high 1 bit) .long 0xBDF817B1 // c1 (low single) .long 0xBD8345EB // c2 .long 0x3D1DFDAC // c3 .long 0xBC52CF6F // c4 .long 0x3F42C75C // B' = pi/2 - B (high single) .long 0xB24BBE8A // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE87283F // c0 (high single) .long 0xB268B966 // c0 (low single) .long 0x3F000000 // c1 (high 1 bit) .long 0xBDFE6529 // c1 (low single) .long 0xBD7B1953 // c2 .long 0x3D18E109 // c3 .long 0xBC4570B0 // c4 .long 0x3F3C7EDD // B' = pi/2 - B (high single) .long 0xB1000ADD // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE827420 // c0 (high single) .long 0x320B8B4D // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DFB9428 // c1 (low single) .long 0xBD7002B4 // c2 .long 0x3D142A6C // c3 .long 0xBC3A47FF // c4 .long 0x3F36365E // B' = pi/2 - B (high single) .long 0x320BB91C // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE7B9282 // c0 (high single) .long 0xB13383D2 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DF5D211 // c1 (low single) .long 0xBD6542B3 // c2 .long 0x3D0FE5E5 // c3 .long 0xBC31FB14 // c4 .long 0x3F2FEDDF // B' = pi/2 - B (high single) .long 0x329BBA77 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE724E73 // c0 (high single) .long 0x3120C3E2 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DF05283 // c1 (low single) .long 0xBD5AD45E // c2 .long 0x3D0BAFBF // c3 .long 0xBC27B8BB // c4 .long 0x3F29A560 // B' = pi/2 - B (high single) .long 0x32F19861 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE691B44 // c0 (high single) .long 0x31F18936 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DEB138B // c1 (low single) .long 0xBD50B2F7 // c2 .long 0x3D07BE3A // c3 .long 0xBC1E46A7 // c4 .long 0x3F235CE2 // B' = pi/2 - B (high single) .long 0xB2B889B6 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE5FF82C // c0 (high single) .long 0xB170723A // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DE61354 // c1 (low single) .long 0xBD46DA06 // c2 .long 0x3D0401F8 // c3 .long 0xBC14E013 // c4 .long 0x3F1D1463 // B' = pi/2 - B (high single) .long 0xB2455799 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE56E46B // c0 (high single) .long 0x31E3F001 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DE15025 // c1 (low single) .long 0xBD3D4550 // c2 .long 0x3D00462D // c3 .long 0xBC092C98 // c4 .long 0x3F16CBE4 // B' = pi/2 - B (high single) .long 0xB0CCDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE4DDF41 // c0 (high single) .long 0xB1AEA094 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DDCC85C // c1 (low single) .long 0xBD33F0BE // c2 .long 0x3CFA23B0 // c3 .long 0xBC01FCF7 // c4 .long 0x3F108365 // B' = pi/2 - B (high single) .long 0x3212200D // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE44E7F8 // c0 (high single) .long 0xB1CAA3CB // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DD87A74 // c1 (low single) .long 0xBD2AD885 // c2 .long 0x3CF3C785 // c3 .long 0xBBF1E348 // c4 .long 0x3F0A3AE6 // B' = pi/2 - B (high single) .long 0x329EEDF0 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE3BFDDC // c0 (high single) .long 0xB132521A // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DD464FC // c1 (low single) .long 0xBD21F8F1 // c2 .long 0x3CEE3076 // c3 .long 0xBBE6D263 // c4 .long 0x3F03F267 // B' = pi/2 - B (high single) .long 0x32F4CBD9 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE33203E // c0 (high single) .long 0x31FEF5BE // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DD0869C // c1 (low single) .long 0xBD194E8C // c2 .long 0x3CE8DCA9 // c3 .long 0xBBDADA55 // c4 .long 0x3EFB53D1 // B' = pi/2 - B (high single) .long 0x32155386 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE2A4E71 // c0 (high single) .long 0xB19CFCEC // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DCCDE11 // c1 (low single) .long 0xBD10D605 // c2 .long 0x3CE382A7 // c3 .long 0xBBC8BD97 // c4 .long 0x3EEEC2D4 // B' = pi/2 - B (high single) .long 0xB23EF0A7 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE2187D0 // c0 (high single) .long 0xB1B7C7F7 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC96A2B // c1 (low single) .long 0xBD088C22 // c2 .long 0x3CDE950E // c3 .long 0xBBB89AD1 // c4 .long 0x3EE231D6 // B' = pi/2 - B (high single) .long 0xB099A6A2 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE18CBB7 // c0 (high single) .long 0xAFE28430 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC629CE // c1 (low single) .long 0xBD006DCD // c2 .long 0x3CDA5A2C // c3 .long 0xBBB0B3D2 // c4 .long 0x3ED5A0D8 // B' = pi/2 - B (high single) .long 0x321886FF // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE101985 // c0 (high single) .long 0xB02FB2B8 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC31BF3 // c1 (low single) .long 0xBCF0F04D // c2 .long 0x3CD60BC7 // c3 .long 0xBBA138BA // c4 .long 0x3EC90FDB // B' = pi/2 - B (high single) .long 0xB23BBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBE07709D // c0 (high single) .long 0xB18A2A83 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC03FA2 // c1 (low single) .long 0xBCE15096 // c2 .long 0x3CD26472 // c3 .long 0xBB9A1270 // c4 .long 0x3EBC7EDD // B' = pi/2 - B (high single) .long 0xB0800ADD // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBDFDA0CB // c0 (high single) .long 0x2F14FCA0 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DBD93F7 // c1 (low single) .long 0xBCD1F71B // c2 .long 0x3CCEDD2B // c3 .long 0xBB905946 // c4 .long 0x3EAFEDDF // B' = pi/2 - B (high single) .long 0x321BBA77 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBDEC708C // c0 (high single) .long 0xB14895C4 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DBB181E // c1 (low single) .long 0xBCC2DEA6 // c2 .long 0x3CCB5027 // c3 .long 0xBB7F3969 // c4 .long 0x3EA35CE2 // B' = pi/2 - B (high single) .long 0xB23889B6 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBDDB4F55 // c0 (high single) .long 0x30F6437E // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB8CB52 // c1 (low single) .long 0xBCB40210 // c2 .long 0x3CC82D45 // c3 .long 0xBB643075 // c4 .long 0x3E96CBE4 // B' = pi/2 - B (high single) .long 0xB04CDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBDCA3BFF // c0 (high single) .long 0x311C95EA // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB6ACDE // c1 (low single) .long 0xBCA55C5B // c2 .long 0x3CC5BC04 // c3 .long 0xBB63A969 // c4 .long 0x3E8A3AE6 // B' = pi/2 - B (high single) .long 0x321EEDF0 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBDB93569 // c0 (high single) .long 0xAFB9ED00 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB4BC1F // c1 (low single) .long 0xBC96E905 // c2 .long 0x3CC2E6F5 // c3 .long 0xBB3E10A6 // c4 .long 0x3E7B53D1 // B' = pi/2 - B (high single) .long 0x31955386 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBDA83A77 // c0 (high single) .long 0x316D967A // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB2F87C // c1 (low single) .long 0xBC88A31F // c2 .long 0x3CC0E763 // c3 .long 0xBB3F1666 // c4 .long 0x3E6231D6 // B' = pi/2 - B (high single) .long 0xB019A6A2 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBD974A0D // c0 (high single) .long 0xB14F365B // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB1616F // c1 (low single) .long 0xBC750CD8 // c2 .long 0x3CBEB595 // c3 .long 0xBB22B883 // c4 .long 0x3E490FDB // B' = pi/2 - B (high single) .long 0xB1BBBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBD866317 // c0 (high single) .long 0xAFF02140 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAFF67D // c1 (low single) .long 0xBC591CD0 // c2 .long 0x3CBCBEAD // c3 .long 0xBB04BBEC // c4 .long 0x3E2FEDDF // B' = pi/2 - B (high single) .long 0x319BBA77 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBD6B08FF // c0 (high single) .long 0xB0EED236 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAEB739 // c1 (low single) .long 0xBC3D6D51 // c2 .long 0x3CBB485D // c3 .long 0xBAFFF5BA // c4 .long 0x3E16CBE4 // B' = pi/2 - B (high single) .long 0xAFCCDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBD495A6C // c0 (high single) .long 0xB0A427BD // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DADA345 // c1 (low single) .long 0xBC21F648 // c2 .long 0x3CB9D1B4 // c3 .long 0xBACB5567 // c4 .long 0x3DFB53D1 // B' = pi/2 - B (high single) .long 0x31155386 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBD27B856 // c0 (high single) .long 0xB0F7EE91 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DACBA4E // c1 (low single) .long 0xBC06AEE3 // c2 .long 0x3CB8E5DC // c3 .long 0xBAEC00EE // c4 .long 0x3DC90FDB // B' = pi/2 - B (high single) .long 0xB13BBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBD0620A3 // c0 (high single) .long 0xB0ECAB40 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DABFC11 // c1 (low single) .long 0xBBD7200F // c2 .long 0x3CB79475 // c3 .long 0xBA2B0ADC // c4 .long 0x3D96CBE4 // B' = pi/2 - B (high single) .long 0xAF4CDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBCC92278 // c0 (high single) .long 0x302F2E68 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAB6854 // c1 (low single) .long 0xBBA1214F // c2 .long 0x3CB6C1E9 // c3 .long 0x3843C2F3 // c4 .long 0x3D490FDB // B' = pi/2 - B (high single) .long 0xB0BBBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBC861015 // c0 (high single) .long 0xAFD68E2E // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAAFEEB // c1 (low single) .long 0xBB569F3F // c2 .long 0x3CB6A84E // c3 .long 0xBAC64194 // c4 .long 0x3CC90FDB // B' = pi/2 - B (high single) .long 0xB03BBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0xBC060BF3 // c0 (high single) .long 0x2FE251AE // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAABFB9 // c1 (low single) .long 0xBAD67C60 // c2 .long 0x3CB64CA5 // c3 .long 0xBACDE881 // c4 .long 0x00000000 // B' = pi/2 - B (high single) .long 0x00000000 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x00000000 // c0 (high single) .long 0x00000000 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAAAAAB // c1 (low single) .long 0x00000000 // c2 .long 0x3CB5E28B // c3 .long 0x00000000 // c4 .long 0xBCC90FDB // B' = pi/2 - B (high single) .long 0x303BBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3C060BF3 // c0 (high single) .long 0xAFE251AE // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAABFB9 // c1 (low single) .long 0x3AD67C60 // c2 .long 0x3CB64CA5 // c3 .long 0x3ACDE881 // c4 .long 0xBD490FDB // B' = pi/2 - B (high single) .long 0x30BBBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3C861015 // c0 (high single) .long 0x2FD68E2E // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAAFEEB // c1 (low single) .long 0x3B569F3F // c2 .long 0x3CB6A84E // c3 .long 0x3AC64194 // c4 .long 0xBD96CBE4 // B' = pi/2 - B (high single) .long 0x2F4CDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3CC92278 // c0 (high single) .long 0xB02F2E68 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAB6854 // c1 (low single) .long 0x3BA1214F // c2 .long 0x3CB6C1E9 // c3 .long 0xB843C2F2 // c4 .long 0xBDC90FDB // B' = pi/2 - B (high single) .long 0x313BBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3D0620A3 // c0 (high single) .long 0x30ECAB40 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DABFC11 // c1 (low single) .long 0x3BD7200F // c2 .long 0x3CB79475 // c3 .long 0x3A2B0ADC // c4 .long 0xBDFB53D1 // B' = pi/2 - B (high single) .long 0xB1155386 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3D27B856 // c0 (high single) .long 0x30F7EE91 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DACBA4E // c1 (low single) .long 0x3C06AEE3 // c2 .long 0x3CB8E5DC // c3 .long 0x3AEC00EE // c4 .long 0xBE16CBE4 // B' = pi/2 - B (high single) .long 0x2FCCDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3D495A6C // c0 (high single) .long 0x30A427BD // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DADA345 // c1 (low single) .long 0x3C21F648 // c2 .long 0x3CB9D1B4 // c3 .long 0x3ACB5567 // c4 .long 0xBE2FEDDF // B' = pi/2 - B (high single) .long 0xB19BBA77 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3D6B08FF // c0 (high single) .long 0x30EED236 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAEB739 // c1 (low single) .long 0x3C3D6D51 // c2 .long 0x3CBB485D // c3 .long 0x3AFFF5BA // c4 .long 0xBE490FDB // B' = pi/2 - B (high single) .long 0x31BBBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3D866317 // c0 (high single) .long 0x2FF02140 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DAFF67D // c1 (low single) .long 0x3C591CD0 // c2 .long 0x3CBCBEAD // c3 .long 0x3B04BBEC // c4 .long 0xBE6231D6 // B' = pi/2 - B (high single) .long 0x3019A6A2 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3D974A0D // c0 (high single) .long 0x314F365B // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB1616F // c1 (low single) .long 0x3C750CD8 // c2 .long 0x3CBEB595 // c3 .long 0x3B22B883 // c4 .long 0xBE7B53D1 // B' = pi/2 - B (high single) .long 0xB1955386 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3DA83A77 // c0 (high single) .long 0xB16D967A // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB2F87C // c1 (low single) .long 0x3C88A31F // c2 .long 0x3CC0E763 // c3 .long 0x3B3F1666 // c4 .long 0xBE8A3AE6 // B' = pi/2 - B (high single) .long 0xB21EEDF0 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3DB93569 // c0 (high single) .long 0x2FB9ED00 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB4BC1F // c1 (low single) .long 0x3C96E905 // c2 .long 0x3CC2E6F5 // c3 .long 0x3B3E10A6 // c4 .long 0xBE96CBE4 // B' = pi/2 - B (high single) .long 0x304CDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3DCA3BFF // c0 (high single) .long 0xB11C95EA // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB6ACDE // c1 (low single) .long 0x3CA55C5B // c2 .long 0x3CC5BC04 // c3 .long 0x3B63A969 // c4 .long 0xBEA35CE2 // B' = pi/2 - B (high single) .long 0x323889B6 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3DDB4F55 // c0 (high single) .long 0xB0F6437E // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DB8CB52 // c1 (low single) .long 0x3CB40210 // c2 .long 0x3CC82D45 // c3 .long 0x3B643075 // c4 .long 0xBEAFEDDF // B' = pi/2 - B (high single) .long 0xB21BBA77 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3DEC708C // c0 (high single) .long 0x314895C4 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DBB181E // c1 (low single) .long 0x3CC2DEA6 // c2 .long 0x3CCB5027 // c3 .long 0x3B7F3969 // c4 .long 0xBEBC7EDD // B' = pi/2 - B (high single) .long 0x30800ADD // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3DFDA0CB // c0 (high single) .long 0xAF14FCA0 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DBD93F7 // c1 (low single) .long 0x3CD1F71B // c2 .long 0x3CCEDD2B // c3 .long 0x3B905946 // c4 .long 0xBEC90FDB // B' = pi/2 - B (high single) .long 0x323BBD2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E07709D // c0 (high single) .long 0x318A2A83 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC03FA2 // c1 (low single) .long 0x3CE15096 // c2 .long 0x3CD26472 // c3 .long 0x3B9A1270 // c4 .long 0xBED5A0D8 // B' = pi/2 - B (high single) .long 0xB21886FF // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E101985 // c0 (high single) .long 0x302FB2B8 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC31BF3 // c1 (low single) .long 0x3CF0F04D // c2 .long 0x3CD60BC7 // c3 .long 0x3BA138BA // c4 .long 0xBEE231D6 // B' = pi/2 - B (high single) .long 0x3099A6A2 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E18CBB7 // c0 (high single) .long 0x2FE28430 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC629CE // c1 (low single) .long 0x3D006DCD // c2 .long 0x3CDA5A2C // c3 .long 0x3BB0B3D2 // c4 .long 0xBEEEC2D4 // B' = pi/2 - B (high single) .long 0x323EF0A7 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E2187D0 // c0 (high single) .long 0x31B7C7F7 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DC96A2B // c1 (low single) .long 0x3D088C22 // c2 .long 0x3CDE950E // c3 .long 0x3BB89AD1 // c4 .long 0xBEFB53D1 // B' = pi/2 - B (high single) .long 0xB2155386 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E2A4E71 // c0 (high single) .long 0x319CFCEC // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DCCDE11 // c1 (low single) .long 0x3D10D605 // c2 .long 0x3CE382A7 // c3 .long 0x3BC8BD97 // c4 .long 0xBF03F267 // B' = pi/2 - B (high single) .long 0xB2F4CBD9 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E33203E // c0 (high single) .long 0xB1FEF5BE // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DD0869C // c1 (low single) .long 0x3D194E8C // c2 .long 0x3CE8DCA9 // c3 .long 0x3BDADA55 // c4 .long 0xBF0A3AE6 // B' = pi/2 - B (high single) .long 0xB29EEDF0 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E3BFDDC // c0 (high single) .long 0x3132521A // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DD464FC // c1 (low single) .long 0x3D21F8F1 // c2 .long 0x3CEE3076 // c3 .long 0x3BE6D263 // c4 .long 0xBF108365 // B' = pi/2 - B (high single) .long 0xB212200D // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E44E7F8 // c0 (high single) .long 0x31CAA3CB // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DD87A74 // c1 (low single) .long 0x3D2AD885 // c2 .long 0x3CF3C785 // c3 .long 0x3BF1E348 // c4 .long 0xBF16CBE4 // B' = pi/2 - B (high single) .long 0x30CCDE2E // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E4DDF41 // c0 (high single) .long 0x31AEA094 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DDCC85C // c1 (low single) .long 0x3D33F0BE // c2 .long 0x3CFA23B0 // c3 .long 0x3C01FCF7 // c4 .long 0xBF1D1463 // B' = pi/2 - B (high single) .long 0x32455799 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E56E46B // c0 (high single) .long 0xB1E3F001 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DE15025 // c1 (low single) .long 0x3D3D4550 // c2 .long 0x3D00462D // c3 .long 0x3C092C98 // c4 .long 0xBF235CE2 // B' = pi/2 - B (high single) .long 0x32B889B6 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E5FF82C // c0 (high single) .long 0x3170723A // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DE61354 // c1 (low single) .long 0x3D46DA06 // c2 .long 0x3D0401F8 // c3 .long 0x3C14E013 // c4 .long 0xBF29A560 // B' = pi/2 - B (high single) .long 0xB2F19861 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E691B44 // c0 (high single) .long 0xB1F18936 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DEB138B // c1 (low single) .long 0x3D50B2F7 // c2 .long 0x3D07BE3A // c3 .long 0x3C1E46A7 // c4 .long 0xBF2FEDDF // B' = pi/2 - B (high single) .long 0xB29BBA77 // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E724E73 // c0 (high single) .long 0xB120C3E2 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DF05283 // c1 (low single) .long 0x3D5AD45E // c2 .long 0x3D0BAFBF // c3 .long 0x3C27B8BB // c4 .long 0xBF36365E // B' = pi/2 - B (high single) .long 0xB20BB91C // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E7B9282 // c0 (high single) .long 0x313383D2 // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DF5D211 // c1 (low single) .long 0x3D6542B3 // c2 .long 0x3D0FE5E5 // c3 .long 0x3C31FB14 // c4 .long 0xBF3C7EDD // B' = pi/2 - B (high single) .long 0x31000ADD // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E827420 // c0 (high single) .long 0xB20B8B4D // c0 (low single) .long 0x3E800000 // c1 (high 1 bit) .long 0x3DFB9428 // c1 (low single) .long 0x3D7002B4 // c2 .long 0x3D142A6C // c3 .long 0x3C3A47FF // c4 .long 0xBF42C75C // B' = pi/2 - B (high single) .long 0x324BBE8A // B' = pi/2 - B (low single) .long 0x3F800000 // tau (1 for cot path) .long 0x3E87283F // c0 (high single) .long 0x3268B966 // c0 (low single) .long 0x3F000000 // c1 (high 1 bit) .long 0xBDFE6529 // c1 (low single) .long 0x3D7B1953 // c2 .long 0x3D18E109 // c3 .long 0x3C4570B0 // c4 .long 0xBF490FDB // B' = pi/2 - B (high single) .long 0x32BBBD2E // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF800000 // c0 (high single) .long 0x2B410000 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xB3000000 // c1 (low single) .long 0xC0000000 // c2 .long 0x402AB7C8 // c3 .long 0xC05561DB // c4 .long 0xBF4F5859 // B' = pi/2 - B (high single) .long 0xB2EE64E8 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF73BB75 // c0 (high single) .long 0xB2FC908D // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBDBF94B0 // c1 (low single) .long 0xBFE8550F // c2 .long 0x40174F67 // c3 .long 0xC036C608 // c4 .long 0xBF55A0D8 // B' = pi/2 - B (high single) .long 0xB29886FF // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF68065E // c0 (high single) .long 0xB2670D1A // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBE36D1D6 // c1 (low single) .long 0xBFD35007 // c2 .long 0x4006A861 // c3 .long 0xC01D4BDA // c4 .long 0xBF5BE957 // B' = pi/2 - B (high single) .long 0xB205522A // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF5CD3BE // c0 (high single) .long 0xB1460308 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBE8306C5 // c1 (low single) .long 0xBFC09232 // c2 .long 0x3FF09632 // c3 .long 0xC007DB00 // c4 .long 0xBF6231D6 // B' = pi/2 - B (high single) .long 0x3119A6A2 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF521801 // c0 (high single) .long 0xB2AE4178 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBEA72938 // c1 (low single) .long 0xBFAFCC22 // c2 .long 0x3FD7BD4A // c3 .long 0xBFEBB01B // c4 .long 0xBF687A55 // B' = pi/2 - B (high single) .long 0x3252257B // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF47C8CC // c0 (high single) .long 0x3200F51A // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBEC82C6C // c1 (low single) .long 0xBFA0BAE9 // c2 .long 0x3FC2252F // c3 .long 0xBFCD24C7 // c4 .long 0xBF6EC2D4 // B' = pi/2 - B (high single) .long 0x32BEF0A7 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF3DDCCF // c0 (high single) .long 0xB2D29606 // c0 (low single) .long 0x40000000 // c1 (high 1 bit) .long 0xBEE6606F // c1 (low single) .long 0xBF9325D6 // c2 .long 0x3FAF4E69 // c3 .long 0xBFB3080C // c4 .long 0xBF750B52 // B' = pi/2 - B (high single) .long 0xB2EB316F // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF344BA9 // c0 (high single) .long 0x32B8B0EA // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EFDF4F7 // c1 (low single) .long 0xBF86DCA8 // c2 .long 0x3F9ED53B // c3 .long 0xBF9CBEDE // c4 .long 0xBF7B53D1 // B' = pi/2 - B (high single) .long 0xB2955386 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF2B0DC1 // c0 (high single) .long 0xB2AB7EBA // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EE496C2 // c1 (low single) .long 0xBF776C40 // c2 .long 0x3F9065C1 // c3 .long 0xBF89AFB6 // c4 .long 0xBF80CE28 // B' = pi/2 - B (high single) .long 0xB1FDD672 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF221C37 // c0 (high single) .long 0x320C61DC // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3ECD4F71 // c1 (low single) .long 0xBF631DAA // c2 .long 0x3F83B471 // c3 .long 0xBF7281EA // c4 .long 0xBF83F267 // B' = pi/2 - B (high single) .long 0xB374CBD9 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF1970C4 // c0 (high single) .long 0xB2904848 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EB7EFF8 // c1 (low single) .long 0xBF50907C // c2 .long 0x3F710FEA // c3 .long 0xBF561FED // c4 .long 0xBF8716A7 // B' = pi/2 - B (high single) .long 0x32588C6D // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF1105AF // c0 (high single) .long 0xB2F045B0 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3EA44EE2 // c1 (low single) .long 0xBF3F8FDB // c2 .long 0x3F5D3FD0 // c3 .long 0xBF3D0A23 // c4 .long 0xBF8A3AE6 // B' = pi/2 - B (high single) .long 0xB31EEDF0 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF08D5B9 // c0 (high single) .long 0x325EF98E // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E92478D // c1 (low single) .long 0xBF2FEDC9 // c2 .long 0x3F4BCD58 // c3 .long 0xBF27AE9E // c4 .long 0xBF8D5F26 // B' = pi/2 - B (high single) .long 0x330C0105 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBF00DC0D // c0 (high single) .long 0x3214AF72 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E81B994 // c1 (low single) .long 0xBF218233 // c2 .long 0x3F3C4531 // c3 .long 0xBF149688 // c4 .long 0xBF908365 // B' = pi/2 - B (high single) .long 0xB292200D // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBEF22870 // c0 (high single) .long 0xB25271F4 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E65107A // c1 (low single) .long 0xBF1429F0 // c2 .long 0x3F2E8AFC // c3 .long 0xBF040498 // c4 .long 0xBF93A7A5 // B' = pi/2 - B (high single) .long 0x3361DEEE // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBEE2F439 // c0 (high single) .long 0x31F4399E // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E49341C // c1 (low single) .long 0xBF07C61A // c2 .long 0x3F22560F // c3 .long 0xBEEAA81E // c4 .long 0xBF96CBE4 // B' = pi/2 - B (high single) .long 0x314CDE2E // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBED413CD // c0 (high single) .long 0x31C06152 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E2FB0CC // c1 (low single) .long 0xBEF876CB // c2 .long 0x3F177807 // c3 .long 0xBED08437 // c4 .long 0xBF99F023 // B' = pi/2 - B (high single) .long 0xB3484328 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBEC5800D // c0 (high single) .long 0x3214C3C1 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E185E54 // c1 (low single) .long 0xBEE2E342 // c2 .long 0x3F0DCA73 // c3 .long 0xBEB8CC21 // c4 .long 0xBF9D1463 // B' = pi/2 - B (high single) .long 0x32C55799 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBEB73250 // c0 (high single) .long 0x32028823 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3E0318F8 // c1 (low single) .long 0xBECEA678 // c2 .long 0x3F053C67 // c3 .long 0xBEA41E53 // c4 .long 0xBFA038A2 // B' = pi/2 - B (high single) .long 0xB2E4CA7E // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBEA92457 // c0 (high single) .long 0xB0B80830 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3DDF8200 // c1 (low single) .long 0xBEBB99E9 // c2 .long 0x3EFB4AA8 // c3 .long 0xBE9182BE // c4 .long 0xBFA35CE2 // B' = pi/2 - B (high single) .long 0x333889B6 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE9B5042 // c0 (high single) .long 0x322A3AEE // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3DBC7490 // c1 (low single) .long 0xBEA99AF5 // c2 .long 0x3EEDE107 // c3 .long 0xBE80E9AA // c4 .long 0xBFA68121 // B' = pi/2 - B (high single) .long 0xB1E43AAC // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE8DB082 // c0 (high single) .long 0x3132A234 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D9CD7D0 // c1 (low single) .long 0xBE988A60 // c2 .long 0x3EE203E3 // c3 .long 0xBE63582C // c4 .long 0xBFA9A560 // B' = pi/2 - B (high single) .long 0xB3719861 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE803FD4 // c0 (high single) .long 0x32279E66 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D807FC8 // c1 (low single) .long 0xBE884BD4 // c2 .long 0x3ED7812D // c3 .long 0xBE4636EB // c4 .long 0xBFACC9A0 // B' = pi/2 - B (high single) .long 0x32655A50 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE65F267 // c0 (high single) .long 0xB1B4B1DF // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D4E8B90 // c1 (low single) .long 0xBE718ACA // c2 .long 0x3ECE7164 // c3 .long 0xBE2DC161 // c4 .long 0xBFAFEDDF // B' = pi/2 - B (high single) .long 0xB31BBA77 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE4BAFAF // c0 (high single) .long 0xAF2A29E0 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3D221018 // c1 (low single) .long 0xBE53BED0 // c2 .long 0x3EC67E26 // c3 .long 0xBE1568E2 // c4 .long 0xBFB3121F // B' = pi/2 - B (high single) .long 0x330F347D // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE31AE4D // c0 (high single) .long 0x31F32251 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3CF6A500 // c1 (low single) .long 0xBE3707DA // c2 .long 0x3EBFA489 // c3 .long 0xBDFBD9C7 // c4 .long 0xBFB6365E // B' = pi/2 - B (high single) .long 0xB28BB91C // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBE17E564 // c0 (high single) .long 0x31C5A2E4 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3CB440D0 // c1 (low single) .long 0xBE1B3D00 // c2 .long 0x3EB9F664 // c3 .long 0xBDD647C0 // c4 .long 0xBFB95A9E // B' = pi/2 - B (high single) .long 0x33651267 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBDFC98C2 // c0 (high single) .long 0x30AE525C // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3C793D20 // c1 (low single) .long 0xBE003845 // c2 .long 0x3EB5271F // c3 .long 0xBDAC669E // c4 .long 0xBFBC7EDD // B' = pi/2 - B (high single) .long 0x31800ADD // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBDC9B5DC // c0 (high single) .long 0xB145AD86 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3C1EEF20 // c1 (low single) .long 0xBDCBAAEA // c2 .long 0x3EB14E5E // c3 .long 0xBD858BB2 // c4 .long 0xBFBFA31C // B' = pi/2 - B (high single) .long 0xB3450FB0 // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBD9711CE // c0 (high single) .long 0xB14FEB28 // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3BB24C00 // c1 (low single) .long 0xBD97E43A // c2 .long 0x3EAE6A89 // c3 .long 0xBD4D07E0 // c4 .long 0xBFC2C75C // B' = pi/2 - B (high single) .long 0x32CBBE8A // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBD49393C // c0 (high single) .long 0xB0A39F5B // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3B1E2B00 // c1 (low single) .long 0xBD49B5D4 // c2 .long 0x3EAC4F10 // c3 .long 0xBCFD9425 // c4 .long 0xBFC5EB9B // B' = pi/2 - B (high single) .long 0xB2DE638C // B' = pi/2 - B (low single) .long 0x00000000 // tau (1 for cot path) .long 0xBCC91A31 // c0 (high single) .long 0xAF8E8D1A // c0 (low single) .long 0x3F800000 // c1 (high 1 bit) .long 0x3A1DFA00 // c1 (low single) .long 0xBCC9392D // c2 .long 0x3EAB1889 // c3 .long 0xBC885D3B // c4 .align 16 .type __svml_stan_data_internal, @object .size __svml_stan_data_internal, .-__svml_stan_data_internal .space 16, 0x00 .align 16 #ifdef __svml_stan_reduction_data_internal_typedef typedef unsigned int VUINT32; typedef struct { __declspec(align(16)) VUINT32 _sPtable[256][3][1]; } __svml_stan_reduction_data_internal; #endif __svml_stan_reduction_data_internal: /* P_hi P_med P_lo */ .long 0x00000000, 0x00000000, 0x00000000 /* 0 */ .long 0x00000000, 0x00000000, 0x00000000 /* 1 */ .long 0x00000000, 0x00000000, 0x00000000 /* 2 */ .long 0x00000000, 0x00000000, 0x00000000 /* 3 */ .long 0x00000000, 0x00000000, 0x00000000 /* 4 */ .long 0x00000000, 0x00000000, 0x00000000 /* 5 */ .long 0x00000000, 0x00000000, 0x00000000 /* 6 */ .long 0x00000000, 0x00000000, 0x00000000 /* 7 */ .long 0x00000000, 0x00000000, 0x00000000 /* 8 */ .long 0x00000000, 0x00000000, 0x00000000 /* 9 */ .long 0x00000000, 0x00000000, 0x00000000 /* 10 */ .long 0x00000000, 0x00000000, 0x00000000 /* 11 */ .long 0x00000000, 0x00000000, 0x00000000 /* 12 */ .long 0x00000000, 0x00000000, 0x00000000 /* 13 */ .long 0x00000000, 0x00000000, 0x00000000 /* 14 */ .long 0x00000000, 0x00000000, 0x00000000 /* 15 */ .long 0x00000000, 0x00000000, 0x00000000 /* 16 */ .long 0x00000000, 0x00000000, 0x00000000 /* 17 */ .long 0x00000000, 0x00000000, 0x00000000 /* 18 */ .long 0x00000000, 0x00000000, 0x00000000 /* 19 */ .long 0x00000000, 0x00000000, 0x00000000 /* 20 */ .long 0x00000000, 0x00000000, 0x00000000 /* 21 */ .long 0x00000000, 0x00000000, 0x00000000 /* 22 */ .long 0x00000000, 0x00000000, 0x00000000 /* 23 */ .long 0x00000000, 0x00000000, 0x00000000 /* 24 */ .long 0x00000000, 0x00000000, 0x00000000 /* 25 */ .long 0x00000000, 0x00000000, 0x00000000 /* 26 */ .long 0x00000000, 0x00000000, 0x00000000 /* 27 */ .long 0x00000000, 0x00000000, 0x00000000 /* 28 */ .long 0x00000000, 0x00000000, 0x00000000 /* 29 */ .long 0x00000000, 0x00000000, 0x00000000 /* 30 */ .long 0x00000000, 0x00000000, 0x00000000 /* 31 */ .long 0x00000000, 0x00000000, 0x00000000 /* 32 */ .long 0x00000000, 0x00000000, 0x00000000 /* 33 */ .long 0x00000000, 0x00000000, 0x00000000 /* 34 */ .long 0x00000000, 0x00000000, 0x00000000 /* 35 */ .long 0x00000000, 0x00000000, 0x00000000 /* 36 */ .long 0x00000000, 0x00000000, 0x00000000 /* 37 */ .long 0x00000000, 0x00000000, 0x00000000 /* 38 */ .long 0x00000000, 0x00000000, 0x00000000 /* 39 */ .long 0x00000000, 0x00000000, 0x00000000 /* 40 */ .long 0x00000000, 0x00000000, 0x00000000 /* 41 */ .long 0x00000000, 0x00000000, 0x00000000 /* 42 */ .long 0x00000000, 0x00000000, 0x00000000 /* 43 */ .long 0x00000000, 0x00000000, 0x00000000 /* 44 */ .long 0x00000000, 0x00000000, 0x00000000 /* 45 */ .long 0x00000000, 0x00000000, 0x00000000 /* 46 */ .long 0x00000000, 0x00000000, 0x00000000 /* 47 */ .long 0x00000000, 0x00000000, 0x00000000 /* 48 */ .long 0x00000000, 0x00000000, 0x00000000 /* 49 */ .long 0x00000000, 0x00000000, 0x00000000 /* 50 */ .long 0x00000000, 0x00000000, 0x00000000 /* 51 */ .long 0x00000000, 0x00000000, 0x00000000 /* 52 */ .long 0x00000000, 0x00000000, 0x00000000 /* 53 */ .long 0x00000000, 0x00000000, 0x00000000 /* 54 */ .long 0x00000000, 0x00000000, 0x00000000 /* 55 */ .long 0x00000000, 0x00000000, 0x00000000 /* 56 */ .long 0x00000000, 0x00000000, 0x00000001 /* 57 */ .long 0x00000000, 0x00000000, 0x00000002 /* 58 */ .long 0x00000000, 0x00000000, 0x00000005 /* 59 */ .long 0x00000000, 0x00000000, 0x0000000A /* 60 */ .long 0x00000000, 0x00000000, 0x00000014 /* 61 */ .long 0x00000000, 0x00000000, 0x00000028 /* 62 */ .long 0x00000000, 0x00000000, 0x00000051 /* 63 */ .long 0x00000000, 0x00000000, 0x000000A2 /* 64 */ .long 0x00000000, 0x00000000, 0x00000145 /* 65 */ .long 0x00000000, 0x00000000, 0x0000028B /* 66 */ .long 0x00000000, 0x00000000, 0x00000517 /* 67 */ .long 0x00000000, 0x00000000, 0x00000A2F /* 68 */ .long 0x00000000, 0x00000000, 0x0000145F /* 69 */ .long 0x00000000, 0x00000000, 0x000028BE /* 70 */ .long 0x00000000, 0x00000000, 0x0000517C /* 71 */ .long 0x00000000, 0x00000000, 0x0000A2F9 /* 72 */ .long 0x00000000, 0x00000000, 0x000145F3 /* 73 */ .long 0x00000000, 0x00000000, 0x00028BE6 /* 74 */ .long 0x00000000, 0x00000000, 0x000517CC /* 75 */ .long 0x00000000, 0x00000000, 0x000A2F98 /* 76 */ .long 0x00000000, 0x00000000, 0x00145F30 /* 77 */ .long 0x00000000, 0x00000000, 0x0028BE60 /* 78 */ .long 0x00000000, 0x00000000, 0x00517CC1 /* 79 */ .long 0x00000000, 0x00000000, 0x00A2F983 /* 80 */ .long 0x00000000, 0x00000000, 0x0145F306 /* 81 */ .long 0x00000000, 0x00000000, 0x028BE60D /* 82 */ .long 0x00000000, 0x00000000, 0x0517CC1B /* 83 */ .long 0x00000000, 0x00000000, 0x0A2F9836 /* 84 */ .long 0x00000000, 0x00000000, 0x145F306D /* 85 */ .long 0x00000000, 0x00000000, 0x28BE60DB /* 86 */ .long 0x00000000, 0x00000000, 0x517CC1B7 /* 87 */ .long 0x00000000, 0x00000000, 0xA2F9836E /* 88 */ .long 0x00000000, 0x00000001, 0x45F306DC /* 89 */ .long 0x00000000, 0x00000002, 0x8BE60DB9 /* 90 */ .long 0x00000000, 0x00000005, 0x17CC1B72 /* 91 */ .long 0x00000000, 0x0000000A, 0x2F9836E4 /* 92 */ .long 0x00000000, 0x00000014, 0x5F306DC9 /* 93 */ .long 0x00000000, 0x00000028, 0xBE60DB93 /* 94 */ .long 0x00000000, 0x00000051, 0x7CC1B727 /* 95 */ .long 0x00000000, 0x000000A2, 0xF9836E4E /* 96 */ .long 0x00000000, 0x00000145, 0xF306DC9C /* 97 */ .long 0x00000000, 0x0000028B, 0xE60DB939 /* 98 */ .long 0x00000000, 0x00000517, 0xCC1B7272 /* 99 */ .long 0x00000000, 0x00000A2F, 0x9836E4E4 /* 100 */ .long 0x00000000, 0x0000145F, 0x306DC9C8 /* 101 */ .long 0x00000000, 0x000028BE, 0x60DB9391 /* 102 */ .long 0x00000000, 0x0000517C, 0xC1B72722 /* 103 */ .long 0x00000000, 0x0000A2F9, 0x836E4E44 /* 104 */ .long 0x00000000, 0x000145F3, 0x06DC9C88 /* 105 */ .long 0x00000000, 0x00028BE6, 0x0DB93910 /* 106 */ .long 0x00000000, 0x000517CC, 0x1B727220 /* 107 */ .long 0x00000000, 0x000A2F98, 0x36E4E441 /* 108 */ .long 0x00000000, 0x00145F30, 0x6DC9C882 /* 109 */ .long 0x00000000, 0x0028BE60, 0xDB939105 /* 110 */ .long 0x00000000, 0x00517CC1, 0xB727220A /* 111 */ .long 0x00000000, 0x00A2F983, 0x6E4E4415 /* 112 */ .long 0x00000000, 0x0145F306, 0xDC9C882A /* 113 */ .long 0x00000000, 0x028BE60D, 0xB9391054 /* 114 */ .long 0x00000000, 0x0517CC1B, 0x727220A9 /* 115 */ .long 0x00000000, 0x0A2F9836, 0xE4E44152 /* 116 */ .long 0x00000000, 0x145F306D, 0xC9C882A5 /* 117 */ .long 0x00000000, 0x28BE60DB, 0x9391054A /* 118 */ .long 0x00000000, 0x517CC1B7, 0x27220A94 /* 119 */ .long 0x00000000, 0xA2F9836E, 0x4E441529 /* 120 */ .long 0x00000001, 0x45F306DC, 0x9C882A53 /* 121 */ .long 0x00000002, 0x8BE60DB9, 0x391054A7 /* 122 */ .long 0x00000005, 0x17CC1B72, 0x7220A94F /* 123 */ .long 0x0000000A, 0x2F9836E4, 0xE441529F /* 124 */ .long 0x00000014, 0x5F306DC9, 0xC882A53F /* 125 */ .long 0x00000028, 0xBE60DB93, 0x91054A7F /* 126 */ .long 0x00000051, 0x7CC1B727, 0x220A94FE /* 127 */ .long 0x000000A2, 0xF9836E4E, 0x441529FC /* 128 */ .long 0x00000145, 0xF306DC9C, 0x882A53F8 /* 129 */ .long 0x0000028B, 0xE60DB939, 0x1054A7F0 /* 130 */ .long 0x00000517, 0xCC1B7272, 0x20A94FE1 /* 131 */ .long 0x00000A2F, 0x9836E4E4, 0x41529FC2 /* 132 */ .long 0x0000145F, 0x306DC9C8, 0x82A53F84 /* 133 */ .long 0x000028BE, 0x60DB9391, 0x054A7F09 /* 134 */ .long 0x0000517C, 0xC1B72722, 0x0A94FE13 /* 135 */ .long 0x0000A2F9, 0x836E4E44, 0x1529FC27 /* 136 */ .long 0x000145F3, 0x06DC9C88, 0x2A53F84E /* 137 */ .long 0x00028BE6, 0x0DB93910, 0x54A7F09D /* 138 */ .long 0x000517CC, 0x1B727220, 0xA94FE13A /* 139 */ .long 0x000A2F98, 0x36E4E441, 0x529FC275 /* 140 */ .long 0x00145F30, 0x6DC9C882, 0xA53F84EA /* 141 */ .long 0x0028BE60, 0xDB939105, 0x4A7F09D5 /* 142 */ .long 0x00517CC1, 0xB727220A, 0x94FE13AB /* 143 */ .long 0x00A2F983, 0x6E4E4415, 0x29FC2757 /* 144 */ .long 0x0145F306, 0xDC9C882A, 0x53F84EAF /* 145 */ .long 0x028BE60D, 0xB9391054, 0xA7F09D5F /* 146 */ .long 0x0517CC1B, 0x727220A9, 0x4FE13ABE /* 147 */ .long 0x0A2F9836, 0xE4E44152, 0x9FC2757D /* 148 */ .long 0x145F306D, 0xC9C882A5, 0x3F84EAFA /* 149 */ .long 0x28BE60DB, 0x9391054A, 0x7F09D5F4 /* 150 */ .long 0x517CC1B7, 0x27220A94, 0xFE13ABE8 /* 151 */ .long 0xA2F9836E, 0x4E441529, 0xFC2757D1 /* 152 */ .long 0x45F306DC, 0x9C882A53, 0xF84EAFA3 /* 153 */ .long 0x8BE60DB9, 0x391054A7, 0xF09D5F47 /* 154 */ .long 0x17CC1B72, 0x7220A94F, 0xE13ABE8F /* 155 */ .long 0x2F9836E4, 0xE441529F, 0xC2757D1F /* 156 */ .long 0x5F306DC9, 0xC882A53F, 0x84EAFA3E /* 157 */ .long 0xBE60DB93, 0x91054A7F, 0x09D5F47D /* 158 */ .long 0x7CC1B727, 0x220A94FE, 0x13ABE8FA /* 159 */ .long 0xF9836E4E, 0x441529FC, 0x2757D1F5 /* 160 */ .long 0xF306DC9C, 0x882A53F8, 0x4EAFA3EA /* 161 */ .long 0xE60DB939, 0x1054A7F0, 0x9D5F47D4 /* 162 */ .long 0xCC1B7272, 0x20A94FE1, 0x3ABE8FA9 /* 163 */ .long 0x9836E4E4, 0x41529FC2, 0x757D1F53 /* 164 */ .long 0x306DC9C8, 0x82A53F84, 0xEAFA3EA6 /* 165 */ .long 0x60DB9391, 0x054A7F09, 0xD5F47D4D /* 166 */ .long 0xC1B72722, 0x0A94FE13, 0xABE8FA9A /* 167 */ .long 0x836E4E44, 0x1529FC27, 0x57D1F534 /* 168 */ .long 0x06DC9C88, 0x2A53F84E, 0xAFA3EA69 /* 169 */ .long 0x0DB93910, 0x54A7F09D, 0x5F47D4D3 /* 170 */ .long 0x1B727220, 0xA94FE13A, 0xBE8FA9A6 /* 171 */ .long 0x36E4E441, 0x529FC275, 0x7D1F534D /* 172 */ .long 0x6DC9C882, 0xA53F84EA, 0xFA3EA69B /* 173 */ .long 0xDB939105, 0x4A7F09D5, 0xF47D4D37 /* 174 */ .long 0xB727220A, 0x94FE13AB, 0xE8FA9A6E /* 175 */ .long 0x6E4E4415, 0x29FC2757, 0xD1F534DD /* 176 */ .long 0xDC9C882A, 0x53F84EAF, 0xA3EA69BB /* 177 */ .long 0xB9391054, 0xA7F09D5F, 0x47D4D377 /* 178 */ .long 0x727220A9, 0x4FE13ABE, 0x8FA9A6EE /* 179 */ .long 0xE4E44152, 0x9FC2757D, 0x1F534DDC /* 180 */ .long 0xC9C882A5, 0x3F84EAFA, 0x3EA69BB8 /* 181 */ .long 0x9391054A, 0x7F09D5F4, 0x7D4D3770 /* 182 */ .long 0x27220A94, 0xFE13ABE8, 0xFA9A6EE0 /* 183 */ .long 0x4E441529, 0xFC2757D1, 0xF534DDC0 /* 184 */ .long 0x9C882A53, 0xF84EAFA3, 0xEA69BB81 /* 185 */ .long 0x391054A7, 0xF09D5F47, 0xD4D37703 /* 186 */ .long 0x7220A94F, 0xE13ABE8F, 0xA9A6EE06 /* 187 */ .long 0xE441529F, 0xC2757D1F, 0x534DDC0D /* 188 */ .long 0xC882A53F, 0x84EAFA3E, 0xA69BB81B /* 189 */ .long 0x91054A7F, 0x09D5F47D, 0x4D377036 /* 190 */ .long 0x220A94FE, 0x13ABE8FA, 0x9A6EE06D /* 191 */ .long 0x441529FC, 0x2757D1F5, 0x34DDC0DB /* 192 */ .long 0x882A53F8, 0x4EAFA3EA, 0x69BB81B6 /* 193 */ .long 0x1054A7F0, 0x9D5F47D4, 0xD377036D /* 194 */ .long 0x20A94FE1, 0x3ABE8FA9, 0xA6EE06DB /* 195 */ .long 0x41529FC2, 0x757D1F53, 0x4DDC0DB6 /* 196 */ .long 0x82A53F84, 0xEAFA3EA6, 0x9BB81B6C /* 197 */ .long 0x054A7F09, 0xD5F47D4D, 0x377036D8 /* 198 */ .long 0x0A94FE13, 0xABE8FA9A, 0x6EE06DB1 /* 199 */ .long 0x1529FC27, 0x57D1F534, 0xDDC0DB62 /* 200 */ .long 0x2A53F84E, 0xAFA3EA69, 0xBB81B6C5 /* 201 */ .long 0x54A7F09D, 0x5F47D4D3, 0x77036D8A /* 202 */ .long 0xA94FE13A, 0xBE8FA9A6, 0xEE06DB14 /* 203 */ .long 0x529FC275, 0x7D1F534D, 0xDC0DB629 /* 204 */ .long 0xA53F84EA, 0xFA3EA69B, 0xB81B6C52 /* 205 */ .long 0x4A7F09D5, 0xF47D4D37, 0x7036D8A5 /* 206 */ .long 0x94FE13AB, 0xE8FA9A6E, 0xE06DB14A /* 207 */ .long 0x29FC2757, 0xD1F534DD, 0xC0DB6295 /* 208 */ .long 0x53F84EAF, 0xA3EA69BB, 0x81B6C52B /* 209 */ .long 0xA7F09D5F, 0x47D4D377, 0x036D8A56 /* 210 */ .long 0x4FE13ABE, 0x8FA9A6EE, 0x06DB14AC /* 211 */ .long 0x9FC2757D, 0x1F534DDC, 0x0DB62959 /* 212 */ .long 0x3F84EAFA, 0x3EA69BB8, 0x1B6C52B3 /* 213 */ .long 0x7F09D5F4, 0x7D4D3770, 0x36D8A566 /* 214 */ .long 0xFE13ABE8, 0xFA9A6EE0, 0x6DB14ACC /* 215 */ .long 0xFC2757D1, 0xF534DDC0, 0xDB629599 /* 216 */ .long 0xF84EAFA3, 0xEA69BB81, 0xB6C52B32 /* 217 */ .long 0xF09D5F47, 0xD4D37703, 0x6D8A5664 /* 218 */ .long 0xE13ABE8F, 0xA9A6EE06, 0xDB14ACC9 /* 219 */ .long 0xC2757D1F, 0x534DDC0D, 0xB6295993 /* 220 */ .long 0x84EAFA3E, 0xA69BB81B, 0x6C52B327 /* 221 */ .long 0x09D5F47D, 0x4D377036, 0xD8A5664F /* 222 */ .long 0x13ABE8FA, 0x9A6EE06D, 0xB14ACC9E /* 223 */ .long 0x2757D1F5, 0x34DDC0DB, 0x6295993C /* 224 */ .long 0x4EAFA3EA, 0x69BB81B6, 0xC52B3278 /* 225 */ .long 0x9D5F47D4, 0xD377036D, 0x8A5664F1 /* 226 */ .long 0x3ABE8FA9, 0xA6EE06DB, 0x14ACC9E2 /* 227 */ .long 0x757D1F53, 0x4DDC0DB6, 0x295993C4 /* 228 */ .long 0xEAFA3EA6, 0x9BB81B6C, 0x52B32788 /* 229 */ .long 0xD5F47D4D, 0x377036D8, 0xA5664F10 /* 230 */ .long 0xABE8FA9A, 0x6EE06DB1, 0x4ACC9E21 /* 231 */ .long 0x57D1F534, 0xDDC0DB62, 0x95993C43 /* 232 */ .long 0xAFA3EA69, 0xBB81B6C5, 0x2B327887 /* 233 */ .long 0x5F47D4D3, 0x77036D8A, 0x5664F10E /* 234 */ .long 0xBE8FA9A6, 0xEE06DB14, 0xACC9E21C /* 235 */ .long 0x7D1F534D, 0xDC0DB629, 0x5993C439 /* 236 */ .long 0xFA3EA69B, 0xB81B6C52, 0xB3278872 /* 237 */ .long 0xF47D4D37, 0x7036D8A5, 0x664F10E4 /* 238 */ .long 0xE8FA9A6E, 0xE06DB14A, 0xCC9E21C8 /* 239 */ .long 0xD1F534DD, 0xC0DB6295, 0x993C4390 /* 240 */ .long 0xA3EA69BB, 0x81B6C52B, 0x32788720 /* 241 */ .long 0x47D4D377, 0x036D8A56, 0x64F10E41 /* 242 */ .long 0x8FA9A6EE, 0x06DB14AC, 0xC9E21C82 /* 243 */ .long 0x1F534DDC, 0x0DB62959, 0x93C43904 /* 244 */ .long 0x3EA69BB8, 0x1B6C52B3, 0x27887208 /* 245 */ .long 0x7D4D3770, 0x36D8A566, 0x4F10E410 /* 246 */ .long 0xFA9A6EE0, 0x6DB14ACC, 0x9E21C820 /* 247 */ .long 0xF534DDC0, 0xDB629599, 0x3C439041 /* 248 */ .long 0xEA69BB81, 0xB6C52B32, 0x78872083 /* 249 */ .long 0xD4D37703, 0x6D8A5664, 0xF10E4107 /* 250 */ .long 0xA9A6EE06, 0xDB14ACC9, 0xE21C820F /* 251 */ .long 0x534DDC0D, 0xB6295993, 0xC439041F /* 252 */ .long 0xA69BB81B, 0x6C52B327, 0x8872083F /* 253 */ .long 0x4D377036, 0xD8A5664F, 0x10E4107F /* 254 */ .long 0x9A6EE06D, 0xB14ACC9E, 0x21C820FF /* 255 */ .align 16 .type __svml_stan_reduction_data_internal, @object .size __svml_stan_reduction_data_internal, .-__svml_stan_reduction_data_internal .align 16 .FLT_16: .long 0xffffffff, 0x00000000, 0xffffffff, 0x00000000 .type .FLT_16, @object .size .FLT_16, 16