xref: /freebsd/contrib/arm-optimized-routines/pl/math/expf_data.c (revision 072a4ba82a01476eaee33781ccd241033eefcf0b)
1*072a4ba8SAndrew Turner /*
2*072a4ba8SAndrew Turner  * Coeffs and table entries for single-precision exp. Copied from
3*072a4ba8SAndrew Turner  * math/exp2f_data.c, with EXP2F_TABLE_BITS == 32.
4*072a4ba8SAndrew Turner  *
5*072a4ba8SAndrew Turner  * Copyright (c) 2017-2023, Arm Limited.
6*072a4ba8SAndrew Turner  * SPDX-License-Identifier: MIT OR Apache-2.0 WITH LLVM-exception
7*072a4ba8SAndrew Turner  */
8*072a4ba8SAndrew Turner 
9*072a4ba8SAndrew Turner #include "math_config.h"
10*072a4ba8SAndrew Turner 
11*072a4ba8SAndrew Turner #define N (1 << EXPF_TABLE_BITS)
12*072a4ba8SAndrew Turner 
13*072a4ba8SAndrew Turner const struct expf_data __expf_data = {
14*072a4ba8SAndrew Turner   /* tab[i] = uint(2^(i/N)) - (i << 52-BITS)
15*072a4ba8SAndrew Turner      used for computing 2^(k/N) for an int |k| < 150 N as
16*072a4ba8SAndrew Turner      double(tab[k%N] + (k << 52-BITS)) */
17*072a4ba8SAndrew Turner   .tab = {
18*072a4ba8SAndrew Turner 0x3ff0000000000000, 0x3fefd9b0d3158574, 0x3fefb5586cf9890f, 0x3fef9301d0125b51,
19*072a4ba8SAndrew Turner 0x3fef72b83c7d517b, 0x3fef54873168b9aa, 0x3fef387a6e756238, 0x3fef1e9df51fdee1,
20*072a4ba8SAndrew Turner 0x3fef06fe0a31b715, 0x3feef1a7373aa9cb, 0x3feedea64c123422, 0x3feece086061892d,
21*072a4ba8SAndrew Turner 0x3feebfdad5362a27, 0x3feeb42b569d4f82, 0x3feeab07dd485429, 0x3feea47eb03a5585,
22*072a4ba8SAndrew Turner 0x3feea09e667f3bcd, 0x3fee9f75e8ec5f74, 0x3feea11473eb0187, 0x3feea589994cce13,
23*072a4ba8SAndrew Turner 0x3feeace5422aa0db, 0x3feeb737b0cdc5e5, 0x3feec49182a3f090, 0x3feed503b23e255d,
24*072a4ba8SAndrew Turner 0x3feee89f995ad3ad, 0x3feeff76f2fb5e47, 0x3fef199bdd85529c, 0x3fef3720dcef9069,
25*072a4ba8SAndrew Turner 0x3fef5818dcfba487, 0x3fef7c97337b9b5f, 0x3fefa4afa2a490da, 0x3fefd0765b6e4540,
26*072a4ba8SAndrew Turner   },
27*072a4ba8SAndrew Turner   .invln2_scaled = 0x1.71547652b82fep+0 * N,
28*072a4ba8SAndrew Turner   .poly_scaled = {
29*072a4ba8SAndrew Turner   0x1.c6af84b912394p-5/N/N/N, 0x1.ebfce50fac4f3p-3/N/N, 0x1.62e42ff0c52d6p-1/N,
30*072a4ba8SAndrew Turner   },
31*072a4ba8SAndrew Turner };
32