1 | /* Function sinh 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 | * Compute sinh(x) as (exp(x)-exp(-x))/2, |
23 | * where exp is calculated as |
24 | * exp(M*ln2 + ln2*(j/2^k) + r) = 2^M * 2^(j/2^k) * exp(r) |
25 | * |
26 | * Special cases: |
27 | * |
28 | * sinh(NaN) = quiet NaN, and raise invalid exception |
29 | * sinh(INF) = that INF |
30 | * sinh(x) = x for subnormals |
31 | * sinh(x) overflows for big x and returns MAXLOG+log(2) |
32 | * |
33 | */ |
34 | |
35 | /* Offsets for data table __svml_dsinh_data_internal |
36 | */ |
37 | #define _dbInvLn2 0 |
38 | #define _dbLn2hi 16 |
39 | #define _dbLn2lo 32 |
40 | #define _dSign 48 |
41 | #define _dbT 64 |
42 | #define _dbShifter 2112 |
43 | #define _iDomainRange 2128 |
44 | #define _dPC2 2144 |
45 | #define _dPC3 2160 |
46 | #define _dPC4 2176 |
47 | #define _dPC5 2192 |
48 | #define _lIndexMask 2208 |
49 | |
50 | #include <sysdep.h> |
51 | |
52 | .section .text.sse4, "ax" , @progbits |
53 | ENTRY(_ZGVbN2v_sinh_sse4) |
54 | subq $72, %rsp |
55 | cfi_def_cfa_offset(80) |
56 | movaps %xmm0, %xmm2 |
57 | |
58 | /* Abs argument */ |
59 | movups _dSign+__svml_dsinh_data_internal(%rip), %xmm0 |
60 | lea _dbT+8+__svml_dsinh_data_internal(%rip), %rsi |
61 | andps %xmm2, %xmm0 |
62 | movaps %xmm0, %xmm1 |
63 | |
64 | /* |
65 | * Load argument |
66 | * dM = x*2^K/log(2) + RShifter |
67 | */ |
68 | movups _dbInvLn2+__svml_dsinh_data_internal(%rip), %xmm10 |
69 | pxor %xmm2, %xmm1 |
70 | mulpd %xmm1, %xmm10 |
71 | movups _dbShifter+__svml_dsinh_data_internal(%rip), %xmm5 |
72 | addpd %xmm5, %xmm10 |
73 | |
74 | /* |
75 | * R |
76 | * dN = dM - RShifter |
77 | */ |
78 | movaps %xmm10, %xmm7 |
79 | subpd %xmm5, %xmm7 |
80 | |
81 | /* dR = dX - dN*Log2_hi/2^K */ |
82 | movups _dbLn2hi+__svml_dsinh_data_internal(%rip), %xmm6 |
83 | mulpd %xmm7, %xmm6 |
84 | |
85 | /* dR = (dX - dN*Log2_hi/2^K) - dN*Log2_lo/2^K */ |
86 | movups _dbLn2lo+__svml_dsinh_data_internal(%rip), %xmm8 |
87 | mulpd %xmm7, %xmm8 |
88 | |
89 | /* |
90 | * Check for overflow\underflow |
91 | * |
92 | */ |
93 | pshufd $221, %xmm1, %xmm4 |
94 | subpd %xmm6, %xmm1 |
95 | subpd %xmm8, %xmm1 |
96 | |
97 | /* VLOAD_CONST( D, dPC[0], TAB._dPC1 ); */ |
98 | movq _iDomainRange+__svml_dsinh_data_internal(%rip), %xmm3 |
99 | pcmpgtd %xmm3, %xmm4 |
100 | |
101 | /* dR2 = dR^2 */ |
102 | movaps %xmm1, %xmm3 |
103 | mulpd %xmm1, %xmm3 |
104 | movmskps %xmm4, %edx |
105 | |
106 | /* |
107 | * sinh(r) = r*((a1=1)+r^2*(a3+r^2*a5)) = r + r*(r^2*(a3+r^2*a5)) .... |
108 | * dSinh_r = (a3+r^2*a5) |
109 | */ |
110 | movups _dPC5+__svml_dsinh_data_internal(%rip), %xmm12 |
111 | |
112 | /* |
113 | * poly(r) = (dG2+dG1)+dG3*sinh(dR)+dG1*sinh(dR)+(dG1+dG2)*dR2*(a2 +a4*dR2) |
114 | * dOut = (a2 +a4*dR2) |
115 | */ |
116 | movups _dPC4+__svml_dsinh_data_internal(%rip), %xmm13 |
117 | mulpd %xmm3, %xmm12 |
118 | mulpd %xmm3, %xmm13 |
119 | addpd _dPC3+__svml_dsinh_data_internal(%rip), %xmm12 |
120 | addpd _dPC2+__svml_dsinh_data_internal(%rip), %xmm13 |
121 | |
122 | /* dSinh_r = r^2*(a3+r^2*a5) */ |
123 | mulpd %xmm3, %xmm12 |
124 | |
125 | /* dOut = dR2*(a2 +a4*dR2) */ |
126 | mulpd %xmm13, %xmm3 |
127 | |
128 | /* dSinh_r = r + r*(r^2*(a3+r^2*a5)) */ |
129 | mulpd %xmm1, %xmm12 |
130 | |
131 | /* |
132 | * Index and lookup |
133 | * j |
134 | */ |
135 | movups _lIndexMask+__svml_dsinh_data_internal(%rip), %xmm9 |
136 | andps %xmm10, %xmm9 |
137 | movd %xmm9, %eax |
138 | |
139 | /* split j and N */ |
140 | pxor %xmm9, %xmm10 |
141 | |
142 | /* |
143 | * G1, G2, G3: dTdif, dTn * 2^N, 2^(-N) |
144 | * lM now is an EXP(2^N) |
145 | */ |
146 | psllq $45, %xmm10 |
147 | |
148 | /* */ |
149 | movaps %xmm10, %xmm4 |
150 | pextrw $4, %xmm9, %ecx |
151 | addpd %xmm12, %xmm1 |
152 | shll $4, %eax |
153 | shll $4, %ecx |
154 | movq (%rax, %rsi), %xmm11 |
155 | movhpd (%rcx, %rsi), %xmm11 |
156 | paddq %xmm11, %xmm4 |
157 | |
158 | /* */ |
159 | psubq %xmm10, %xmm11 |
160 | |
161 | /* dG3 = dTn*2^N + dTn*2^-N */ |
162 | movdqa %xmm4, %xmm14 |
163 | addpd %xmm11, %xmm14 |
164 | |
165 | /* dG2 = dTn*2^N - dTn*2^-N */ |
166 | subpd %xmm11, %xmm4 |
167 | movq -8(%rax, %rsi), %xmm15 |
168 | movhpd -8(%rcx, %rsi), %xmm15 |
169 | paddq %xmm10, %xmm15 |
170 | |
171 | /* dG2 += dG1 */ |
172 | addpd %xmm15, %xmm4 |
173 | |
174 | /* dG1 += dG3 */ |
175 | addpd %xmm14, %xmm15 |
176 | |
177 | /* dOut = dG2*dR2*(a2 +a4*dR2) */ |
178 | mulpd %xmm4, %xmm3 |
179 | |
180 | /* dOut = dG1*sinh(dR)+dG2*dR2*(a2 +a4*dR2) */ |
181 | mulpd %xmm15, %xmm1 |
182 | addpd %xmm1, %xmm3 |
183 | |
184 | /* dOut = dG2 + dG1*sinh(dR)+dG2*dR2*(a2 +a4*dR2) */ |
185 | addpd %xmm3, %xmm4 |
186 | |
187 | /* Ret H */ |
188 | orps %xmm4, %xmm0 |
189 | andl $3, %edx |
190 | |
191 | /* Go to special inputs processing branch */ |
192 | jne L(SPECIAL_VALUES_BRANCH) |
193 | # LOE rbx rbp r12 r13 r14 r15 edx xmm0 xmm2 |
194 | |
195 | /* Restore registers |
196 | * and exit the function |
197 | */ |
198 | |
199 | L(EXIT): |
200 | addq $72, %rsp |
201 | cfi_def_cfa_offset(8) |
202 | ret |
203 | cfi_def_cfa_offset(80) |
204 | |
205 | /* Branch to process |
206 | * special inputs |
207 | */ |
208 | |
209 | L(SPECIAL_VALUES_BRANCH): |
210 | movups %xmm2, 32(%rsp) |
211 | movups %xmm0, 48(%rsp) |
212 | # LOE rbx rbp r12 r13 r14 r15 edx xmm0 |
213 | |
214 | xorl %eax, %eax |
215 | movq %r12, 16(%rsp) |
216 | cfi_offset(12, -64) |
217 | movl %eax, %r12d |
218 | movq %r13, 8(%rsp) |
219 | cfi_offset(13, -72) |
220 | movl %edx, %r13d |
221 | movq %r14, (%rsp) |
222 | cfi_offset(14, -80) |
223 | # LOE rbx rbp r15 r12d r13d |
224 | |
225 | /* Range mask |
226 | * bits check |
227 | */ |
228 | |
229 | L(RANGEMASK_CHECK): |
230 | btl %r12d, %r13d |
231 | |
232 | /* Call scalar math function */ |
233 | jc L(SCALAR_MATH_CALL) |
234 | # LOE rbx rbp r15 r12d r13d |
235 | |
236 | /* Special inputs |
237 | * processing loop |
238 | */ |
239 | |
240 | L(SPECIAL_VALUES_LOOP): |
241 | incl %r12d |
242 | cmpl $2, %r12d |
243 | |
244 | /* Check bits in range mask */ |
245 | jl L(RANGEMASK_CHECK) |
246 | # LOE rbx rbp r15 r12d r13d |
247 | |
248 | movq 16(%rsp), %r12 |
249 | cfi_restore(12) |
250 | movq 8(%rsp), %r13 |
251 | cfi_restore(13) |
252 | movq (%rsp), %r14 |
253 | cfi_restore(14) |
254 | movups 48(%rsp), %xmm0 |
255 | |
256 | /* Go to exit */ |
257 | jmp L(EXIT) |
258 | cfi_offset(12, -64) |
259 | cfi_offset(13, -72) |
260 | cfi_offset(14, -80) |
261 | # LOE rbx rbp r12 r13 r14 r15 xmm0 |
262 | |
263 | /* Scalar math function call |
264 | * to process special input |
265 | */ |
266 | |
267 | L(SCALAR_MATH_CALL): |
268 | movl %r12d, %r14d |
269 | movsd 32(%rsp, %r14, 8), %xmm0 |
270 | call sinh@PLT |
271 | # LOE rbx rbp r14 r15 r12d r13d xmm0 |
272 | |
273 | movsd %xmm0, 48(%rsp, %r14, 8) |
274 | |
275 | /* Process special inputs in loop */ |
276 | jmp L(SPECIAL_VALUES_LOOP) |
277 | # LOE rbx rbp r15 r12d r13d |
278 | END(_ZGVbN2v_sinh_sse4) |
279 | |
280 | .section .rodata, "a" |
281 | .align 16 |
282 | |
283 | #ifdef __svml_dsinh_data_internal_typedef |
284 | typedef unsigned int VUINT32; |
285 | typedef struct { |
286 | __declspec(align(16)) VUINT32 _dbInvLn2[2][2]; |
287 | __declspec(align(16)) VUINT32 _dbLn2hi[2][2]; |
288 | __declspec(align(16)) VUINT32 _dbLn2lo[2][2]; |
289 | __declspec(align(16)) VUINT32 _dSign[2][2]; // 0x8000000000000000 |
290 | __declspec(align(16)) VUINT32 _dbT[(1<<7)][2][2]; // precalc poly coeff |
291 | __declspec(align(16)) VUINT32 _dbShifter[2][2]; |
292 | __declspec(align(16)) VUINT32 _iDomainRange[4][1]; |
293 | __declspec(align(16)) VUINT32 _dPC2[2][2]; |
294 | __declspec(align(16)) VUINT32 _dPC3[2][2]; |
295 | __declspec(align(16)) VUINT32 _dPC4[2][2]; |
296 | __declspec(align(16)) VUINT32 _dPC5[2][2]; |
297 | __declspec(align(16)) VUINT32 _lIndexMask[2][2]; |
298 | } __svml_dsinh_data_internal; |
299 | #endif |
300 | __svml_dsinh_data_internal: |
301 | .quad 0x3FF71547652B82FE, 0x3FF71547652B82FE /* _dbInvLn2 = 1/log(2) */ |
302 | .align 16 |
303 | .quad 0x3FE62E42FEFA0000, 0x3FE62E42FEFA0000 /* _dbLn2hi = log(2) hi */ |
304 | .align 16 |
305 | .quad 0x3D7CF79ABC9E3B3A, 0x3D7CF79ABC9E3B3A /* _dbLn2lo = log(2) lo */ |
306 | .align 16 |
307 | .quad 0x8000000000000000, 0x8000000000000000 /* _dSign */ |
308 | //_dbT |
309 | .align 16 |
310 | .quad 0x0000000000000000, 0x3FE0000000000000 // 2^( 0 /128-1) - 2^(- 0 /128-1), 2^(- 0 /128-1) |
311 | .quad 0x3F762E4A19BD1E74, 0x3FDFD3C22B8F71F1 // 2^( 1 /128-1) - 2^(- 1 /128-1), 2^(- 1 /128-1) |
312 | .quad 0x3F862E5F6A0DFD36, 0x3FDFA7C1819E90D8 // 2^( 2 /128-1) - 2^(- 2 /128-1), 2^(- 2 /128-1) |
313 | .quad 0x3F90A2E234040F5F, 0x3FDF7BFDAD9CBE14 // 2^( 3 /128-1) - 2^(- 3 /128-1), 2^(- 3 /128-1) |
314 | .quad 0x3F962EB4ABCC5A81, 0x3FDF50765B6E4540 // 2^( 4 /128-1) - 2^(- 4 /128-1), 2^(- 4 /128-1) |
315 | .quad 0x3F9BBAB1C5033244, 0x3FDF252B376BBA97 // 2^( 5 /128-1) - 2^(- 5 /128-1), 2^(- 5 /128-1) |
316 | .quad 0x3FA0A372144EEB45, 0x3FDEFA1BEE615A27 // 2^( 6 /128-1) - 2^(- 6 /128-1), 2^(- 6 /128-1) |
317 | .quad 0x3FA369AB3FFBF8B0, 0x3FDECF482D8E67F1 // 2^( 7 /128-1) - 2^(- 7 /128-1), 2^(- 7 /128-1) |
318 | .quad 0x3FA63009BA740A2A, 0x3FDEA4AFA2A490DA // 2^( 8 /128-1) - 2^(- 8 /128-1), 2^(- 8 /128-1) |
319 | .quad 0x3FA8F692D8EA1B5A, 0x3FDE7A51FBC74C83 // 2^( 9 /128-1) - 2^(- 9 /128-1), 2^(- 9 /128-1) |
320 | .quad 0x3FABBD4BF0E31A6F, 0x3FDE502EE78B3FF6 // 2^( 10 /128-1) - 2^(- 10 /128-1), 2^(- 10 /128-1) |
321 | .quad 0x3FAE843A5840286A, 0x3FDE264614F5A129 // 2^( 11 /128-1) - 2^(- 11 /128-1), 2^(- 11 /128-1) |
322 | .quad 0x3FB0A5B1B2A46D0A, 0x3FDDFC97337B9B5F // 2^( 12 /128-1) - 2^(- 12 /128-1), 2^(- 12 /128-1) |
323 | .quad 0x3FB20966375ABCDF, 0x3FDDD321F301B460 // 2^( 13 /128-1) - 2^(- 13 /128-1), 2^(- 13 /128-1) |
324 | .quad 0x3FB36D3D65DCA4E8, 0x3FDDA9E603DB3285 // 2^( 14 /128-1) - 2^(- 14 /128-1), 2^(- 14 /128-1) |
325 | .quad 0x3FB4D139EA06642A, 0x3FDD80E316C98398 // 2^( 15 /128-1) - 2^(- 15 /128-1), 2^(- 15 /128-1) |
326 | .quad 0x3FB6355E6FFBF9BA, 0x3FDD5818DCFBA487 // 2^( 16 /128-1) - 2^(- 16 /128-1), 2^(- 16 /128-1) |
327 | .quad 0x3FB799ADA42E4788, 0x3FDD2F87080D89F2 // 2^( 17 /128-1) - 2^(- 17 /128-1), 2^(- 17 /128-1) |
328 | .quad 0x3FB8FE2A336035BC, 0x3FDD072D4A07897C // 2^( 18 /128-1) - 2^(- 18 /128-1), 2^(- 18 /128-1) |
329 | .quad 0x3FBA62D6CAABD6B6, 0x3FDCDF0B555DC3FA // 2^( 19 /128-1) - 2^(- 19 /128-1), 2^(- 19 /128-1) |
330 | .quad 0x3FBBC7B617878BAF, 0x3FDCB720DCEF9069 // 2^( 20 /128-1) - 2^(- 20 /128-1), 2^(- 20 /128-1) |
331 | .quad 0x3FBD2CCAC7CB2A11, 0x3FDC8F6D9406E7B5 // 2^( 21 /128-1) - 2^(- 21 /128-1), 2^(- 21 /128-1) |
332 | .quad 0x3FBE921789B52185, 0x3FDC67F12E57D14B // 2^( 22 /128-1) - 2^(- 22 /128-1), 2^(- 22 /128-1) |
333 | .quad 0x3FBFF79F0BEFA2C7, 0x3FDC40AB5FFFD07A // 2^( 23 /128-1) - 2^(- 23 /128-1), 2^(- 23 /128-1) |
334 | .quad 0x3FC0AEB1FECAE3A9, 0x3FDC199BDD85529C // 2^( 24 /128-1) - 2^(- 24 /128-1), 2^(- 24 /128-1) |
335 | .quad 0x3FC161B4871C5CEC, 0x3FDBF2C25BD71E09 // 2^( 25 /128-1) - 2^(- 25 /128-1), 2^(- 25 /128-1) |
336 | .quad 0x3FC214D876F26FD0, 0x3FDBCC1E904BC1D2 // 2^( 26 /128-1) - 2^(- 26 /128-1), 2^(- 26 /128-1) |
337 | .quad 0x3FC2C81F2693816F, 0x3FDBA5B030A1064A // 2^( 27 /128-1) - 2^(- 27 /128-1), 2^(- 27 /128-1) |
338 | .quad 0x3FC37B89EE88BEF7, 0x3FDB7F76F2FB5E47 // 2^( 28 /128-1) - 2^(- 28 /128-1), 2^(- 28 /128-1) |
339 | .quad 0x3FC42F1A27A0B3CD, 0x3FDB59728DE5593A // 2^( 29 /128-1) - 2^(- 29 /128-1), 2^(- 29 /128-1) |
340 | .quad 0x3FC4E2D12AF1E037, 0x3FDB33A2B84F15FB // 2^( 30 /128-1) - 2^(- 30 /128-1), 2^(- 30 /128-1) |
341 | .quad 0x3FC596B051DD508D, 0x3FDB0E07298DB666 // 2^( 31 /128-1) - 2^(- 31 /128-1), 2^(- 31 /128-1) |
342 | .quad 0x3FC64AB8F61134FA, 0x3FDAE89F995AD3AD // 2^( 32 /128-1) - 2^(- 32 /128-1), 2^(- 32 /128-1) |
343 | .quad 0x3FC6FEEC718B79D1, 0x3FDAC36BBFD3F37A // 2^( 33 /128-1) - 2^(- 33 /128-1), 2^(- 33 /128-1) |
344 | .quad 0x3FC7B34C1E9C607F, 0x3FDA9E6B5579FDBF // 2^( 34 /128-1) - 2^(- 34 /128-1), 2^(- 34 /128-1) |
345 | .quad 0x3FC867D957E91912, 0x3FDA799E1330B358 // 2^( 35 /128-1) - 2^(- 35 /128-1), 2^(- 35 /128-1) |
346 | .quad 0x3FC91C95786E5C72, 0x3FDA5503B23E255D // 2^( 36 /128-1) - 2^(- 36 /128-1), 2^(- 36 /128-1) |
347 | .quad 0x3FC9D181DB83072F, 0x3FDA309BEC4A2D33 // 2^( 37 /128-1) - 2^(- 37 /128-1), 2^(- 37 /128-1) |
348 | .quad 0x3FCA869FDCDAB512, 0x3FDA0C667B5DE565 // 2^( 38 /128-1) - 2^(- 38 /128-1), 2^(- 38 /128-1) |
349 | .quad 0x3FCB3BF0D8885D4C, 0x3FD9E86319E32323 // 2^( 39 /128-1) - 2^(- 39 /128-1), 2^(- 39 /128-1) |
350 | .quad 0x3FCBF1762B00EF69, 0x3FD9C49182A3F090 // 2^( 40 /128-1) - 2^(- 40 /128-1), 2^(- 40 /128-1) |
351 | .quad 0x3FCCA731311DF0FB, 0x3FD9A0F170CA07BA // 2^( 41 /128-1) - 2^(- 41 /128-1), 2^(- 41 /128-1) |
352 | .quad 0x3FCD5D2348201C09, 0x3FD97D829FDE4E50 // 2^( 42 /128-1) - 2^(- 42 /128-1), 2^(- 42 /128-1) |
353 | .quad 0x3FCE134DCDB1FE3E, 0x3FD95A44CBC8520F // 2^( 43 /128-1) - 2^(- 43 /128-1), 2^(- 43 /128-1) |
354 | .quad 0x3FCEC9B21FEA98EA, 0x3FD93737B0CDC5E5 // 2^( 44 /128-1) - 2^(- 44 /128-1), 2^(- 44 /128-1) |
355 | .quad 0x3FCF80519D5001D3, 0x3FD9145B0B91FFC6 // 2^( 45 /128-1) - 2^(- 45 /128-1), 2^(- 45 /128-1) |
356 | .quad 0x3FD01B96D26D026A, 0x3FD8F1AE99157736 // 2^( 46 /128-1) - 2^(- 46 /128-1), 2^(- 46 /128-1) |
357 | .quad 0x3FD07723CAFA6331, 0x3FD8CF3216B5448C // 2^( 47 /128-1) - 2^(- 47 /128-1), 2^(- 47 /128-1) |
358 | .quad 0x3FD0D2D06841B373, 0x3FD8ACE5422AA0DB // 2^( 48 /128-1) - 2^(- 48 /128-1), 2^(- 48 /128-1) |
359 | .quad 0x3FD12E9D5A715381, 0x3FD88AC7D98A6699 // 2^( 49 /128-1) - 2^(- 49 /128-1), 2^(- 49 /128-1) |
360 | .quad 0x3FD18A8B51F5C661, 0x3FD868D99B4492ED // 2^( 50 /128-1) - 2^(- 50 /128-1), 2^(- 50 /128-1) |
361 | .quad 0x3FD1E69AFF7B04D7, 0x3FD8471A4623C7AD // 2^( 51 /128-1) - 2^(- 51 /128-1), 2^(- 51 /128-1) |
362 | .quad 0x3FD242CD13EDD0F1, 0x3FD82589994CCE13 // 2^( 52 /128-1) - 2^(- 52 /128-1), 2^(- 52 /128-1) |
363 | .quad 0x3FD29F22407D0A0C, 0x3FD80427543E1A12 // 2^( 53 /128-1) - 2^(- 53 /128-1), 2^(- 53 /128-1) |
364 | .quad 0x3FD2FB9B369B0153, 0x3FD7E2F336CF4E62 // 2^( 54 /128-1) - 2^(- 54 /128-1), 2^(- 54 /128-1) |
365 | .quad 0x3FD35838A7FECEC8, 0x3FD7C1ED0130C132 // 2^( 55 /128-1) - 2^(- 55 /128-1), 2^(- 55 /128-1) |
366 | .quad 0x3FD3B4FB46A5A6CC, 0x3FD7A11473EB0187 // 2^( 56 /128-1) - 2^(- 56 /128-1), 2^(- 56 /128-1) |
367 | .quad 0x3FD411E3C4D4302F, 0x3FD780694FDE5D3F // 2^( 57 /128-1) - 2^(- 57 /128-1), 2^(- 57 /128-1) |
368 | .quad 0x3FD46EF2D517DAC8, 0x3FD75FEB564267C9 // 2^( 58 /128-1) - 2^(- 58 /128-1), 2^(- 58 /128-1) |
369 | .quad 0x3FD4CC292A48369E, 0x3FD73F9A48A58174 // 2^( 59 /128-1) - 2^(- 59 /128-1), 2^(- 59 /128-1) |
370 | .quad 0x3FD5298777884B96, 0x3FD71F75E8EC5F74 // 2^( 60 /128-1) - 2^(- 60 /128-1), 2^(- 60 /128-1) |
371 | .quad 0x3FD5870E7047F1BC, 0x3FD6FF7DF9519484 // 2^( 61 /128-1) - 2^(- 61 /128-1), 2^(- 61 /128-1) |
372 | .quad 0x3FD5E4BEC8452A1A, 0x3FD6DFB23C651A2F // 2^( 62 /128-1) - 2^(- 62 /128-1), 2^(- 62 /128-1) |
373 | .quad 0x3FD64299338D7827, 0x3FD6C012750BDABF // 2^( 63 /128-1) - 2^(- 63 /128-1), 2^(- 63 /128-1) |
374 | .quad 0x3FD6A09E667F3BCD, 0x3FD6A09E667F3BCD // 2^( 64 /128-1) - 2^(- 64 /128-1), 2^(- 64 /128-1) |
375 | .quad 0x3FD6FECF15CB0C0B, 0x3FD68155D44CA973 // 2^( 65 /128-1) - 2^(- 65 /128-1), 2^(- 65 /128-1) |
376 | .quad 0x3FD75D2BF6751239, 0x3FD6623882552225 // 2^( 66 /128-1) - 2^(- 66 /128-1), 2^(- 66 /128-1) |
377 | .quad 0x3FD7BBB5BDD665E8, 0x3FD6434634CCC320 // 2^( 67 /128-1) - 2^(- 67 /128-1), 2^(- 67 /128-1) |
378 | .quad 0x3FD81A6D219E6963, 0x3FD6247EB03A5585 // 2^( 68 /128-1) - 2^(- 68 /128-1), 2^(- 68 /128-1) |
379 | .quad 0x3FD87952D7D426DF, 0x3FD605E1B976DC09 // 2^( 69 /128-1) - 2^(- 69 /128-1), 2^(- 69 /128-1) |
380 | .quad 0x3FD8D86796D7AE49, 0x3FD5E76F15AD2148 // 2^( 70 /128-1) - 2^(- 70 /128-1), 2^(- 70 /128-1) |
381 | .quad 0x3FD937AC156373C8, 0x3FD5C9268A5946B7 // 2^( 71 /128-1) - 2^(- 71 /128-1), 2^(- 71 /128-1) |
382 | .quad 0x3FD997210A8DAEE4, 0x3FD5AB07DD485429 // 2^( 72 /128-1) - 2^(- 72 /128-1), 2^(- 72 /128-1) |
383 | .quad 0x3FD9F6C72DC9BA68, 0x3FD58D12D497C7FD // 2^( 73 /128-1) - 2^(- 73 /128-1), 2^(- 73 /128-1) |
384 | .quad 0x3FDA569F36E974EA, 0x3FD56F4736B527DA // 2^( 74 /128-1) - 2^(- 74 /128-1), 2^(- 74 /128-1) |
385 | .quad 0x3FDAB6A9DE1EA215, 0x3FD551A4CA5D920F // 2^( 75 /128-1) - 2^(- 75 /128-1), 2^(- 75 /128-1) |
386 | .quad 0x3FDB16E7DBFC4CA3, 0x3FD5342B569D4F82 // 2^( 76 /128-1) - 2^(- 76 /128-1), 2^(- 76 /128-1) |
387 | .quad 0x3FDB7759E9782918, 0x3FD516DAA2CF6642 // 2^( 77 /128-1) - 2^(- 77 /128-1), 2^(- 77 /128-1) |
388 | .quad 0x3FDBD800BFEBF932, 0x3FD4F9B2769D2CA7 // 2^( 78 /128-1) - 2^(- 78 /128-1), 2^(- 78 /128-1) |
389 | .quad 0x3FDC38DD1916F025, 0x3FD4DCB299FDDD0D // 2^( 79 /128-1) - 2^(- 79 /128-1), 2^(- 79 /128-1) |
390 | .quad 0x3FDC99EFAF1F1790, 0x3FD4BFDAD5362A27 // 2^( 80 /128-1) - 2^(- 80 /128-1), 2^(- 80 /128-1) |
391 | .quad 0x3FDCFB393C92B539, 0x3FD4A32AF0D7D3DE // 2^( 81 /128-1) - 2^(- 81 /128-1), 2^(- 81 /128-1) |
392 | .quad 0x3FDD5CBA7C69B19C, 0x3FD486A2B5C13CD0 // 2^( 82 /128-1) - 2^(- 82 /128-1), 2^(- 82 /128-1) |
393 | .quad 0x3FDDBE742A06FF34, 0x3FD46A41ED1D0057 // 2^( 83 /128-1) - 2^(- 83 /128-1), 2^(- 83 /128-1) |
394 | .quad 0x3FDE2067013A029D, 0x3FD44E086061892D // 2^( 84 /128-1) - 2^(- 84 /128-1), 2^(- 84 /128-1) |
395 | .quad 0x3FDE8293BE3FFB87, 0x3FD431F5D950A897 // 2^( 85 /128-1) - 2^(- 85 /128-1), 2^(- 85 /128-1) |
396 | .quad 0x3FDEE4FB1DC56E75, 0x3FD4160A21F72E2A // 2^( 86 /128-1) - 2^(- 86 /128-1), 2^(- 86 /128-1) |
397 | .quad 0x3FDF479DDCE78F58, 0x3FD3FA4504AC801C // 2^( 87 /128-1) - 2^(- 87 /128-1), 2^(- 87 /128-1) |
398 | .quad 0x3FDFAA7CB935ACFE, 0x3FD3DEA64C123422 // 2^( 88 /128-1) - 2^(- 88 /128-1), 2^(- 88 /128-1) |
399 | .quad 0x3FE006CC38594EB1, 0x3FD3C32DC313A8E5 // 2^( 89 /128-1) - 2^(- 89 /128-1), 2^(- 89 /128-1) |
400 | .quad 0x3FE03878E0EB1569, 0x3FD3A7DB34E59FF7 // 2^( 90 /128-1) - 2^(- 90 /128-1), 2^(- 90 /128-1) |
401 | .quad 0x3FE06A44B5C74101, 0x3FD38CAE6D05D866 // 2^( 91 /128-1) - 2^(- 91 /128-1), 2^(- 91 /128-1) |
402 | .quad 0x3FE09C3016A0D077, 0x3FD371A7373AA9CB // 2^( 92 /128-1) - 2^(- 92 /128-1), 2^(- 92 /128-1) |
403 | .quad 0x3FE0CE3B63676360, 0x3FD356C55F929FF1 // 2^( 93 /128-1) - 2^(- 93 /128-1), 2^(- 93 /128-1) |
404 | .quad 0x3FE10066FC47F240, 0x3FD33C08B26416FF // 2^( 94 /128-1) - 2^(- 94 /128-1), 2^(- 94 /128-1) |
405 | .quad 0x3FE132B341AD8761, 0x3FD32170FC4CD831 // 2^( 95 /128-1) - 2^(- 95 /128-1), 2^(- 95 /128-1) |
406 | .quad 0x3FE165209441F823, 0x3FD306FE0A31B715 // 2^( 96 /128-1) - 2^(- 96 /128-1), 2^(- 96 /128-1) |
407 | .quad 0x3FE197AF54EE9EBB, 0x3FD2ECAFA93E2F56 // 2^( 97 /128-1) - 2^(- 97 /128-1), 2^(- 97 /128-1) |
408 | .quad 0x3FE1CA5FE4DD1475, 0x3FD2D285A6E4030B // 2^( 98 /128-1) - 2^(- 98 /128-1), 2^(- 98 /128-1) |
409 | .quad 0x3FE1FD32A577EC72, 0x3FD2B87FD0DAD990 // 2^( 99 /128-1) - 2^(- 99 /128-1), 2^(- 99 /128-1) |
410 | .quad 0x3FE23027F86B6ED6, 0x3FD29E9DF51FDEE1 // 2^( 100 /128-1) - 2^(- 100 /128-1), 2^(- 100 /128-1) |
411 | .quad 0x3FE263403FA65489, 0x3FD284DFE1F56381 // 2^( 101 /128-1) - 2^(- 101 /128-1), 2^(- 101 /128-1) |
412 | .quad 0x3FE2967BDD5A8364, 0x3FD26B4565E27CDD // 2^( 102 /128-1) - 2^(- 102 /128-1), 2^(- 102 /128-1) |
413 | .quad 0x3FE2C9DB33FDCAE9, 0x3FD251CE4FB2A63F // 2^( 103 /128-1) - 2^(- 103 /128-1), 2^(- 103 /128-1) |
414 | .quad 0x3FE2FD5EA64AA180, 0x3FD2387A6E756238 // 2^( 104 /128-1) - 2^(- 104 /128-1), 2^(- 104 /128-1) |
415 | .quad 0x3FE331069740E22F, 0x3FD21F49917DDC96 // 2^( 105 /128-1) - 2^(- 105 /128-1), 2^(- 105 /128-1) |
416 | .quad 0x3FE364D36A268AE0, 0x3FD2063B88628CD6 // 2^( 106 /128-1) - 2^(- 106 /128-1), 2^(- 106 /128-1) |
417 | .quad 0x3FE398C582887B27, 0x3FD1ED5022FCD91D // 2^( 107 /128-1) - 2^(- 107 /128-1), 2^(- 107 /128-1) |
418 | .quad 0x3FE3CCDD443B3394, 0x3FD1D4873168B9AA // 2^( 108 /128-1) - 2^(- 108 /128-1), 2^(- 108 /128-1) |
419 | .quad 0x3FE4011B135B9590, 0x3FD1BBE084045CD4 // 2^( 109 /128-1) - 2^(- 109 /128-1), 2^(- 109 /128-1) |
420 | .quad 0x3FE4357F544FA3C1, 0x3FD1A35BEB6FCB75 // 2^( 110 /128-1) - 2^(- 110 /128-1), 2^(- 110 /128-1) |
421 | .quad 0x3FE46A0A6BC742FD, 0x3FD18AF9388C8DEA // 2^( 111 /128-1) - 2^(- 111 /128-1), 2^(- 111 /128-1) |
422 | .quad 0x3FE49EBCBEBCFBCA, 0x3FD172B83C7D517B // 2^( 112 /128-1) - 2^(- 112 /128-1), 2^(- 112 /128-1) |
423 | .quad 0x3FE4D396B276BC6F, 0x3FD15A98C8A58E51 // 2^( 113 /128-1) - 2^(- 113 /128-1), 2^(- 113 /128-1) |
424 | .quad 0x3FE50898AC869B96, 0x3FD1429AAEA92DE0 // 2^( 114 /128-1) - 2^(- 114 /128-1), 2^(- 114 /128-1) |
425 | .quad 0x3FE53DC312CB9B7A, 0x3FD12ABDC06C31CC // 2^( 115 /128-1) - 2^(- 115 /128-1), 2^(- 115 /128-1) |
426 | .quad 0x3FE573164B726DB6, 0x3FD11301D0125B51 // 2^( 116 /128-1) - 2^(- 116 /128-1), 2^(- 116 /128-1) |
427 | .quad 0x3FE5A892BCF6379B, 0x3FD0FB66AFFED31B // 2^( 117 /128-1) - 2^(- 117 /128-1), 2^(- 117 /128-1) |
428 | .quad 0x3FE5DE38CE215725, 0x3FD0E3EC32D3D1A2 // 2^( 118 /128-1) - 2^(- 118 /128-1), 2^(- 118 /128-1) |
429 | .quad 0x3FE61408E60E2888, 0x3FD0CC922B7247F7 // 2^( 119 /128-1) - 2^(- 119 /128-1), 2^(- 119 /128-1) |
430 | .quad 0x3FE64A036C27CC52, 0x3FD0B5586CF9890F // 2^( 120 /128-1) - 2^(- 120 /128-1), 2^(- 120 /128-1) |
431 | .quad 0x3FE68028C82AEE2F, 0x3FD09E3ECAC6F383 // 2^( 121 /128-1) - 2^(- 121 /128-1), 2^(- 121 /128-1) |
432 | .quad 0x3FE6B67962268C43, 0x3FD0874518759BC8 // 2^( 122 /128-1) - 2^(- 122 /128-1), 2^(- 122 /128-1) |
433 | .quad 0x3FE6ECF5A27CBF28, 0x3FD0706B29DDF6DE // 2^( 123 /128-1) - 2^(- 123 /128-1), 2^(- 123 /128-1) |
434 | .quad 0x3FE7239DF1E38286, 0x3FD059B0D3158574 // 2^( 124 /128-1) - 2^(- 124 /128-1), 2^(- 124 /128-1) |
435 | .quad 0x3FE75A72B9657E51, 0x3FD04315E86E7F85 // 2^( 125 /128-1) - 2^(- 125 /128-1), 2^(- 125 /128-1) |
436 | .quad 0x3FE791746262D0A8, 0x3FD02C9A3E778061 // 2^( 126 /128-1) - 2^(- 126 /128-1), 2^(- 126 /128-1) |
437 | .quad 0x3FE7C8A35691D856, 0x3FD0163DA9FB3335 // 2^( 127 /128-1) - 2^(- 127 /128-1), 2^(- 127 /128-1) |
438 | .align 16 |
439 | .quad 0x42C8000000000000, 0x42C8000000000000 /* _dbShifter = 1.5 * 2^(52-k) */ |
440 | .align 16 |
441 | .long 0x40861d99, 0x40861d99, 0x40861d99, 0x40861d99 /* _iDomainRange 0x40861d9ac12a3e85 =(1021*2^K-0.5)*log(2)/2^K -needed for quick exp */ |
442 | .align 16 |
443 | .quad 0x3FDFFFFFFFFFFDBD, 0x3FDFFFFFFFFFFDBD /* _dPC2 */ |
444 | .align 16 |
445 | .quad 0x3FC55555555554AD, 0x3FC55555555554AD /* _dPC3 */ |
446 | .align 16 |
447 | .quad 0x3FA55555CF16D299, 0x3FA55555CF16D299 /* _dPC4 */ |
448 | .align 16 |
449 | .quad 0x3F8111115712F425, 0x3F8111115712F425 /* _dPC5 */ |
450 | .align 16 |
451 | .quad 0x000000000000007f, 0x000000000000007f /* _lIndexMask */ |
452 | .align 16 |
453 | .type __svml_dsinh_data_internal, @object |
454 | .size __svml_dsinh_data_internal, .-__svml_dsinh_data_internal |
455 | |