xref: /aosp_15_r20/external/arm-optimized-routines/math/exp2f_data.c (revision 412f47f9e737e10ed5cc46ec6a8d7fa2264f8a14)
1*412f47f9SXin Li /*
2*412f47f9SXin Li  * Shared data between expf, exp2f and powf.
3*412f47f9SXin Li  *
4*412f47f9SXin Li  * Copyright (c) 2017-2018, Arm Limited.
5*412f47f9SXin Li  * SPDX-License-Identifier: MIT OR Apache-2.0 WITH LLVM-exception
6*412f47f9SXin Li  */
7*412f47f9SXin Li 
8*412f47f9SXin Li #include "math_config.h"
9*412f47f9SXin Li 
10*412f47f9SXin Li #define N (1 << EXP2F_TABLE_BITS)
11*412f47f9SXin Li 
12*412f47f9SXin Li const struct exp2f_data __exp2f_data = {
13*412f47f9SXin Li   /* tab[i] = uint(2^(i/N)) - (i << 52-BITS)
14*412f47f9SXin Li      used for computing 2^(k/N) for an int |k| < 150 N as
15*412f47f9SXin Li      double(tab[k%N] + (k << 52-BITS)) */
16*412f47f9SXin Li   .tab = {
17*412f47f9SXin Li #if N == 8
18*412f47f9SXin Li 0x3ff0000000000000, 0x3fef72b83c7d517b, 0x3fef06fe0a31b715, 0x3feebfdad5362a27,
19*412f47f9SXin Li 0x3feea09e667f3bcd, 0x3feeace5422aa0db, 0x3feee89f995ad3ad, 0x3fef5818dcfba487,
20*412f47f9SXin Li #elif N == 16
21*412f47f9SXin Li 0x3ff0000000000000, 0x3fefb5586cf9890f, 0x3fef72b83c7d517b, 0x3fef387a6e756238,
22*412f47f9SXin Li 0x3fef06fe0a31b715, 0x3feedea64c123422, 0x3feebfdad5362a27, 0x3feeab07dd485429,
23*412f47f9SXin Li 0x3feea09e667f3bcd, 0x3feea11473eb0187, 0x3feeace5422aa0db, 0x3feec49182a3f090,
24*412f47f9SXin Li 0x3feee89f995ad3ad, 0x3fef199bdd85529c, 0x3fef5818dcfba487, 0x3fefa4afa2a490da,
25*412f47f9SXin Li #elif N == 32
26*412f47f9SXin Li 0x3ff0000000000000, 0x3fefd9b0d3158574, 0x3fefb5586cf9890f, 0x3fef9301d0125b51,
27*412f47f9SXin Li 0x3fef72b83c7d517b, 0x3fef54873168b9aa, 0x3fef387a6e756238, 0x3fef1e9df51fdee1,
28*412f47f9SXin Li 0x3fef06fe0a31b715, 0x3feef1a7373aa9cb, 0x3feedea64c123422, 0x3feece086061892d,
29*412f47f9SXin Li 0x3feebfdad5362a27, 0x3feeb42b569d4f82, 0x3feeab07dd485429, 0x3feea47eb03a5585,
30*412f47f9SXin Li 0x3feea09e667f3bcd, 0x3fee9f75e8ec5f74, 0x3feea11473eb0187, 0x3feea589994cce13,
31*412f47f9SXin Li 0x3feeace5422aa0db, 0x3feeb737b0cdc5e5, 0x3feec49182a3f090, 0x3feed503b23e255d,
32*412f47f9SXin Li 0x3feee89f995ad3ad, 0x3feeff76f2fb5e47, 0x3fef199bdd85529c, 0x3fef3720dcef9069,
33*412f47f9SXin Li 0x3fef5818dcfba487, 0x3fef7c97337b9b5f, 0x3fefa4afa2a490da, 0x3fefd0765b6e4540,
34*412f47f9SXin Li #elif N == 64
35*412f47f9SXin Li 0x3ff0000000000000, 0x3fefec9a3e778061, 0x3fefd9b0d3158574, 0x3fefc74518759bc8,
36*412f47f9SXin Li 0x3fefb5586cf9890f, 0x3fefa3ec32d3d1a2, 0x3fef9301d0125b51, 0x3fef829aaea92de0,
37*412f47f9SXin Li 0x3fef72b83c7d517b, 0x3fef635beb6fcb75, 0x3fef54873168b9aa, 0x3fef463b88628cd6,
38*412f47f9SXin Li 0x3fef387a6e756238, 0x3fef2b4565e27cdd, 0x3fef1e9df51fdee1, 0x3fef1285a6e4030b,
39*412f47f9SXin Li 0x3fef06fe0a31b715, 0x3feefc08b26416ff, 0x3feef1a7373aa9cb, 0x3feee7db34e59ff7,
40*412f47f9SXin Li 0x3feedea64c123422, 0x3feed60a21f72e2a, 0x3feece086061892d, 0x3feec6a2b5c13cd0,
41*412f47f9SXin Li 0x3feebfdad5362a27, 0x3feeb9b2769d2ca7, 0x3feeb42b569d4f82, 0x3feeaf4736b527da,
42*412f47f9SXin Li 0x3feeab07dd485429, 0x3feea76f15ad2148, 0x3feea47eb03a5585, 0x3feea23882552225,
43*412f47f9SXin Li 0x3feea09e667f3bcd, 0x3fee9fb23c651a2f, 0x3fee9f75e8ec5f74, 0x3fee9feb564267c9,
44*412f47f9SXin Li 0x3feea11473eb0187, 0x3feea2f336cf4e62, 0x3feea589994cce13, 0x3feea8d99b4492ed,
45*412f47f9SXin Li 0x3feeace5422aa0db, 0x3feeb1ae99157736, 0x3feeb737b0cdc5e5, 0x3feebd829fde4e50,
46*412f47f9SXin Li 0x3feec49182a3f090, 0x3feecc667b5de565, 0x3feed503b23e255d, 0x3feede6b5579fdbf,
47*412f47f9SXin Li 0x3feee89f995ad3ad, 0x3feef3a2b84f15fb, 0x3feeff76f2fb5e47, 0x3fef0c1e904bc1d2,
48*412f47f9SXin Li 0x3fef199bdd85529c, 0x3fef27f12e57d14b, 0x3fef3720dcef9069, 0x3fef472d4a07897c,
49*412f47f9SXin Li 0x3fef5818dcfba487, 0x3fef69e603db3285, 0x3fef7c97337b9b5f, 0x3fef902ee78b3ff6,
50*412f47f9SXin Li 0x3fefa4afa2a490da, 0x3fefba1bee615a27, 0x3fefd0765b6e4540, 0x3fefe7c1819e90d8,
51*412f47f9SXin Li #endif
52*412f47f9SXin Li   },
53*412f47f9SXin Li   .shift_scaled = 0x1.8p+52 / N,
54*412f47f9SXin Li   .poly = {
55*412f47f9SXin Li #if N == 8
56*412f47f9SXin Li   0x1.c6a00335106e2p-5, 0x1.ec0c313449f55p-3, 0x1.62e431111f69fp-1,
57*412f47f9SXin Li #elif N == 16
58*412f47f9SXin Li   0x1.c6ac6aa313963p-5, 0x1.ebfff4532d9bap-3, 0x1.62e43001bc49fp-1,
59*412f47f9SXin Li #elif N == 32
60*412f47f9SXin Li   0x1.c6af84b912394p-5, 0x1.ebfce50fac4f3p-3, 0x1.62e42ff0c52d6p-1,
61*412f47f9SXin Li #elif N == 64
62*412f47f9SXin Li   0x1.c6b04b4221b2ap-5, 0x1.ebfc213e184d7p-3, 0x1.62e42fefb5b7fp-1,
63*412f47f9SXin Li #endif
64*412f47f9SXin Li   },
65*412f47f9SXin Li   .shift = 0x1.8p+52,
66*412f47f9SXin Li   .invln2_scaled = 0x1.71547652b82fep+0 * N,
67*412f47f9SXin Li   .poly_scaled = {
68*412f47f9SXin Li #if N == 8
69*412f47f9SXin Li   0x1.c6a00335106e2p-5/N/N/N, 0x1.ec0c313449f55p-3/N/N, 0x1.62e431111f69fp-1/N,
70*412f47f9SXin Li #elif N == 16
71*412f47f9SXin Li   0x1.c6ac6aa313963p-5/N/N/N, 0x1.ebfff4532d9bap-3/N/N, 0x1.62e43001bc49fp-1/N,
72*412f47f9SXin Li #elif N == 32
73*412f47f9SXin Li   0x1.c6af84b912394p-5/N/N/N, 0x1.ebfce50fac4f3p-3/N/N, 0x1.62e42ff0c52d6p-1/N,
74*412f47f9SXin Li #elif N == 64
75*412f47f9SXin Li   0x1.c6b04b4221b2ap-5/N/N/N, 0x1.ebfc213e184d7p-3/N/N, 0x1.62e42fefb5b7fp-1/N,
76*412f47f9SXin Li #endif
77*412f47f9SXin Li   },
78*412f47f9SXin Li };
79