1/*
2 * Polynomial coefficients and other constants for tgamma128.c.
3 *
4 * Copyright (c) 2006-2023, Arm Limited.
5 * SPDX-License-Identifier: MIT OR Apache-2.0 WITH LLVM-exception
6 */
7
8/* The largest positive value for which 128-bit tgamma does not overflow. */
9static const long double max_x =  0x1.b6e3180cd66a5c4206f128ba77f4p+10L;
10
11/* Coefficients of the polynomial used in the tgamma_large() subroutine */
12static const long double coeffs_large[] = {
13     0x1.8535745aa79569579b9eec0f3bbcp+0L,
14     0x1.0378f83c6fb8f0e51269f2b4a973p-3L,
15     0x1.59f6a05094f69686c3380f4e2783p-8L,
16    -0x1.0b291dee952a82764a4859b081a6p-8L,
17    -0x1.6dd301b2205bf936b5a3eaad0dbbp-12L,
18     0x1.387a8b5f38dd77e7f139b1021e86p-10L,
19     0x1.bca46637f65b13750c728cc29e40p-14L,
20    -0x1.d80401c00aef998c9e303151a51cp-11L,
21    -0x1.49cb6bb09f935a2053ccc2cf3711p-14L,
22     0x1.4e950204437dcaf2be77f73a6f45p-10L,
23     0x1.cb711a2d65f188bf60110934d6bep-14L,
24    -0x1.7d7ff4bc95dc7faefc5e767f70f1p-9L,
25    -0x1.0305ab9760cddb0d833e73766836p-12L,
26     0x1.3ef6c84bf1cd5c3f65ac2693bb5bp-7L,
27     0x1.bb4144740ad9290123fdcea684aap-11L,
28    -0x1.72ab4e88272a229bfafd192450f0p-5L,
29     0x1.80c70ac6eb3b7a698983d25a62b8p-12L,
30     0x1.e222791c6743ce3e3cae220fb236p-3L,
31     0x1.1a2dca1c82a9326c52b465f7cb7ap-2L,
32    -0x1.9d204fa235a42cd901b123d2ad47p+1L,
33     0x1.55b56d1158f77ddb1c95fc44ab02p+0L,
34     0x1.37f900a11dbd892abd7dde533e2dp+5L,
35    -0x1.2da49f4188dd89cb958369ef2401p+7L,
36     0x1.fdae5ec3ec6eb7dffc09edbe6c14p+7L,
37    -0x1.61433cebe649098c9611c4c7774ap+7L,
38};
39
40/* Coefficients of the polynomial used in the tgamma_tiny() subroutine */
41static const long double coeffs_tiny[] = {
42     0x1.0000000000000000000000000000p+0L,
43     0x1.2788cfc6fb618f49a37c7f0201fep-1L,
44    -0x1.4fcf4026afa2dceb8490ade22796p-1L,
45    -0x1.5815e8fa27047c8f42b5d9217244p-5L,
46     0x1.5512320b43fbe5dfa771333518f7p-3L,
47    -0x1.59af103c340927bffdd44f954bfcp-5L,
48    -0x1.3b4af28483e210479657e5543366p-7L,
49     0x1.d919c527f6070bfce9b29c2ace9cp-8L,
50    -0x1.317112ce35337def3556a18aa178p-10L,
51    -0x1.c364fe77a6f27677b985b1fa2e1dp-13L,
52     0x1.0c8a7a19a3fd40fe1f7e867efe7bp-13L,
53    -0x1.51cf9f090b5dc398ba86305e3634p-16L,
54    -0x1.4e80f64c04a339740de06ca9fa4ap-20L,
55     0x1.241ddc2aef2ec20e58b08f2fda17p-20L,
56};
57
58/* The location within the interval [1,2] where gamma has a minimum.
59 * Specified as the sum of two 128-bit values, for extra precision. */
60static const long double min_x_hi =  0x1.762d86356be3f6e1a9c8865e0a4fp+0L;
61static const long double min_x_lo =  0x1.ac54d7d218de21303a7c60f08840p-118L;
62
63/* The actual minimum value that gamma takes at that location.
64 * Again specified as the sum of two 128-bit values. */
65static const long double min_y_hi =  0x1.c56dc82a74aee8d8851566d40f32p-1L;
66static const long double min_y_lo =  0x1.8ed98685742c353ce55e5794686fp-114L;
67
68/* Coefficients of the polynomial used in the tgamma_central() subroutine
69 * for computing gamma on the interval [1,min_x] */
70static const long double coeffs_central_neg[] = {
71     0x1.b6c53f7377b83839c8a292e43b69p-2L,
72     0x1.0bae9f40c7d09ed76e732045850ap-3L,
73     0x1.4981175e14d04c3530e51d01c5fep-3L,
74     0x1.79f77aaf032c948af3a9edbd2061p-4L,
75     0x1.1e97bd10821095a5b79fbfdfa1a3p-4L,
76     0x1.8071ce0935e4dcf0b33b0fbec7c1p-5L,
77     0x1.0b44c2f92982f887b55ec36dfdb0p-5L,
78     0x1.6df1de1e178ef72ca7bd63d40870p-6L,
79     0x1.f63f502bde27e81c0f5e13479b43p-7L,
80     0x1.57fd67d901f40ea011353ad89a0ap-7L,
81     0x1.d7151376eed187eb753e2273cafcp-8L,
82     0x1.427162b5c6ff1d904c71ef53e37cp-8L,
83     0x1.b954b8c3a56cf93e49ef6538928ap-9L,
84     0x1.2dff2ec26a3ae5cd3aaccae7a09ep-9L,
85     0x1.9d35250d9b9378d9b59df734537ap-10L,
86     0x1.1b2c0c48b9855a28f6dbd6fdff3cp-10L,
87     0x1.7e0db39bb99cdb52b028d9359380p-11L,
88     0x1.2164b5e1d364a0b5eaf97c436aa7p-11L,
89     0x1.27521cf5fd24dcdf43524e6add11p-13L,
90     0x1.06461d62243bf9a826b42349672fp-10L,
91    -0x1.2b852abead28209b4e0c756dc46ep-9L,
92     0x1.be673c11a72c826115ec6d286c14p-8L,
93    -0x1.fd9ce330c215c31fcd3cb53c42ebp-7L,
94     0x1.fa362bd2dc68f41abef2d8600acdp-6L,
95    -0x1.a21585b2f52f8b23855de8e452edp-5L,
96     0x1.1f234431ed032052fc92e64e0493p-4L,
97    -0x1.40d332476ca0199c60cdae3f9132p-4L,
98     0x1.1d45dc665d86012eba2eea199cefp-4L,
99    -0x1.8491016cdd08dc9be7ade9b5fef3p-5L,
100     0x1.7e7e2fbc6d49ad484300d6add324p-6L,
101    -0x1.e63fe3f874a37276a8d7d8b705ecp-8L,
102     0x1.30a2a73944f8c84998314d69c23fp-10L,
103};
104
105/* Coefficients of the polynomial used in the tgamma_central() subroutine
106 * for computing gamma on the interval [min_x,2] */
107static const long double coeffs_central_pos[] = {
108     0x1.b6c53f7377b83839c8a292e22aa2p-2L,
109    -0x1.0bae9f40c7d09ed76e72e1c955dep-3L,
110     0x1.4981175e14d04c3530ee5e1ecebcp-3L,
111    -0x1.79f77aaf032c948ac983d77f3e07p-4L,
112     0x1.1e97bd10821095ab7dc94936cc11p-4L,
113    -0x1.8071ce0935e4d7edef8cbf2a1cf1p-5L,
114     0x1.0b44c2f929837fafef7b5d9e80f1p-5L,
115    -0x1.6df1de1e175fe2a51faa25cddbb4p-6L,
116     0x1.f63f502be57d11aed2cfe90843ffp-7L,
117    -0x1.57fd67d852f230015b9f64770273p-7L,
118     0x1.d715138adc07e5fce81077070357p-8L,
119    -0x1.4271618e9fda8992a667adb15f4fp-8L,
120     0x1.b954d15d9eb772e80fdd760672d7p-9L,
121    -0x1.2dfe391241d3cb79c8c15182843dp-9L,
122     0x1.9d44396fcd48451c3ba924cee814p-10L,
123    -0x1.1ac195fb99739e341589e39803e6p-10L,
124     0x1.82e46127b68f002770826e25f146p-11L,
125    -0x1.089dacd90d9f41493119ac178359p-11L,
126     0x1.6993c007b20394a057d21f3d37f8p-12L,
127    -0x1.ec43a709f4446560c099dec8e31bp-13L,
128     0x1.4ba36322f4074e9add9450f003cap-13L,
129    -0x1.b3f83a977965ca1b7937bf5b34cap-14L,
130     0x1.10af346abc09cb25a6d9fe810b6ep-14L,
131    -0x1.38d8ea1188f242f50203edc395bdp-15L,
132     0x1.39add987a948ec56f62b721a4475p-16L,
133    -0x1.02a4e141f286c8a967e2df9bc9adp-17L,
134     0x1.433b50af22425f546e87113062d7p-19L,
135    -0x1.0c7b73cb0013f00aafc103e8e382p-21L,
136     0x1.b852de313ec38da2297f6deaa6b4p-25L,
137};
138
139/* 128-bit float value of pi, used by the sin_pi_x_over_pi subroutine
140 */
141static const long double pi =  0x1.921fb54442d18469898cc51701b8p+1L;
142