?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
Run DCE after a LoopFlatten test to reduce spurious output [nfc]
[llvm-project.git] / libc / src / math / generic / log2.cpp
blobd72b0931c14b81b678184c0a4e9fcb7e83712b49
1 //===-- Double-precision log2(x) function ---------------------------------===//
2 //
3 // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4 // See https://llvm.org/LICENSE.txt for license information.
5 // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6 //
7 //===----------------------------------------------------------------------===//
8
9 #include "src/math/log2.h"
10 #include "src/__support/FPUtil/FEnvImpl.h"
11 #include "src/__support/FPUtil/FPBits.h"
12 #include "src/__support/FPUtil/PolyEval.h"
13 #include "src/__support/FPUtil/double_double.h"
14 #include "src/__support/FPUtil/dyadic_float.h"
15 #include "src/__support/FPUtil/multiply_add.h"
16 #include "src/__support/common.h"
17 #include "src/__support/macros/optimization.h" // LIBC_UNLIKELY
18
19 #include "common_constants.h"
20 #include "log_range_reduction.h"
21
22 namespace LIBC_NAMESPACE {
23
24 // 128-bit precision dyadic floating point numbers.
25 using Float128 = typename fputil::DyadicFloat<128>;
26 using MType = typename Float128::MantissaType;
27
28 namespace {
29
30 constexpr fputil::DoubleDouble LOG2_E = {0x1.777d0ffda0d24p-56,
31 0x1.71547652b82fep0};
32
33 // Extra errors from P is from using x^2 to reduce evaluation latency.
34 constexpr double P_ERR = 0x1.0p-49;
35
36 const fputil::DoubleDouble LOG_R1[128] = {
37 {0.0, 0.0},
38 {0x1.46662d417cedp-62, 0x1.010157588de71p-7},
39 {0x1.27c8e8416e71fp-60, 0x1.0205658935847p-6},
40 {-0x1.d192d0619fa67p-60, 0x1.8492528c8cabfp-6},
41 {0x1.c05cf1d753622p-59, 0x1.0415d89e74444p-5},
42 {-0x1.cdd6f7f4a137ep-59, 0x1.466aed42de3eap-5},
43 {0x1.a8be97660a23dp-60, 0x1.894aa149fb343p-5},
44 {-0x1.e48fb0500efd4p-59, 0x1.ccb73cdddb2ccp-5},
45 {-0x1.dd7009902bf32p-58, 0x1.08598b59e3a07p-4},
46 {-0x1.7558367a6acf6p-59, 0x1.1973bd1465567p-4},
47 {0x1.7a976d3b5b45fp-59, 0x1.3bdf5a7d1ee64p-4},
48 {0x1.f38745c5c450ap-58, 0x1.5e95a4d9791cbp-4},
49 {-0x1.72566212cdd05p-61, 0x1.700d30aeac0e1p-4},
50 {-0x1.478a85704ccb7p-58, 0x1.9335e5d594989p-4},
51 {-0x1.0057eed1ca59fp-59, 0x1.b6ac88dad5b1cp-4},
52 {0x1.a38cb559a6706p-58, 0x1.c885801bc4b23p-4},
53 {-0x1.a2bf991780d3fp-59, 0x1.ec739830a112p-4},
54 {-0x1.ac9f4215f9393p-58, 0x1.fe89139dbd566p-4},
55 {-0x1.0e63a5f01c691p-58, 0x1.1178e8227e47cp-3},
56 {-0x1.c6ef1d9b2ef7ep-59, 0x1.1aa2b7e23f72ap-3},
57 {-0x1.499a3f25af95fp-58, 0x1.2d1610c86813ap-3},
58 {0x1.7d411a5b944adp-58, 0x1.365fcb0159016p-3},
59 {-0x1.0d5604930f135p-58, 0x1.4913d8333b561p-3},
60 {-0x1.71a9682395bfdp-61, 0x1.527e5e4a1b58dp-3},
61 {-0x1.d34f0f4621bedp-60, 0x1.6574ebe8c133ap-3},
62 {-0x1.8de59c21e166cp-57, 0x1.6f0128b756abcp-3},
63 {-0x1.1232ce70be781p-57, 0x1.823c16551a3c2p-3},
64 {0x1.55aa8b6997a4p-58, 0x1.8beafeb38fe8cp-3},
65 {0x1.142c507fb7a3dp-58, 0x1.95a5adcf7017fp-3},
66 {0x1.bcafa9de97203p-57, 0x1.a93ed3c8ad9e3p-3},
67 {-0x1.6353ab386a94dp-57, 0x1.b31d8575bce3dp-3},
68 {0x1.dd355f6a516d7p-60, 0x1.bd087383bd8adp-3},
69 {0x1.60629242471a2p-57, 0x1.d1037f2655e7bp-3},
70 {0x1.aa11d49f96cb9p-58, 0x1.db13db0d4894p-3},
71 {0x1.2276041f43042p-59, 0x1.e530effe71012p-3},
72 {-0x1.08ab2ddc708ap-58, 0x1.ef5ade4dcffe6p-3},
73 {0x1.f665066f980a2p-57, 0x1.f991c6cb3b379p-3},
74 {0x1.cdb16ed4e9138p-56, 0x1.07138604d5862p-2},
75 {0x1.162c79d5d11eep-58, 0x1.0c42d676162e3p-2},
76 {-0x1.0e63a5f01c691p-57, 0x1.1178e8227e47cp-2},
77 {0x1.66fbd28b40935p-56, 0x1.16b5ccbacfb73p-2},
78 {-0x1.12aeb84249223p-57, 0x1.1bf99635a6b95p-2},
79 {0x1.e0efadd9db02bp-56, 0x1.269621134db92p-2},
80 {-0x1.82dad7fd86088p-56, 0x1.2bef07cdc9354p-2},
81 {-0x1.3d69909e5c3dcp-56, 0x1.314f1e1d35ce4p-2},
82 {-0x1.324f0e883858ep-58, 0x1.36b6776be1117p-2},
83 {-0x1.2ad27e50a8ec6p-56, 0x1.3c25277333184p-2},
84 {0x1.0dbb243827392p-57, 0x1.419b423d5e8c7p-2},
85 {0x1.8fb4c14c56eefp-60, 0x1.4718dc271c41bp-2},
86 {-0x1.123615b147a5dp-58, 0x1.4c9e09e172c3cp-2},
87 {-0x1.8f7e9b38a6979p-57, 0x1.522ae0738a3d8p-2},
88 {-0x1.0908d15f88b63p-57, 0x1.57bf753c8d1fbp-2},
89 {-0x1.6541148cbb8a2p-56, 0x1.5d5bddf595f3p-2},
90 {0x1.dc18ce51fff99p-57, 0x1.630030b3aac49p-2},
91 {0x1.a64eadd740178p-58, 0x1.68ac83e9c6a14p-2},
92 {0x1.657c222d868cdp-58, 0x1.6e60ee6af1972p-2},
93 {0x1.84a4ee3059583p-56, 0x1.741d876c67bb1p-2},
94 {-0x1.c168817443f22p-56, 0x1.79e26687cfb3ep-2},
95 {-0x1.219024acd3b77p-58, 0x1.7fafa3bd8151cp-2},
96 {-0x1.486666443b153p-56, 0x1.85855776dcbfbp-2},
97 {-0x1.70f2f38238303p-56, 0x1.8b639a88b2df5p-2},
98 {-0x1.ad4bb98c1f2c5p-56, 0x1.914a8635bf68ap-2},
99 {-0x1.89d2816cf838fp-57, 0x1.973a3431356aep-2},
100 {0x1.87bcbcfd3e187p-59, 0x1.9d32bea15ed3bp-2},
101 {-0x1.ba8062860ae23p-57, 0x1.a33440224fa79p-2},
102 {-0x1.ba8062860ae23p-57, 0x1.a33440224fa79p-2},
103 {0x1.bcafa9de97203p-56, 0x1.a93ed3c8ad9e3p-2},
104 {0x1.9d56c45dd3e86p-56, 0x1.af5295248cddp-2},
105 {0x1.494b610665378p-56, 0x1.b56fa04462909p-2},
106 {0x1.6fd02999b21e1p-59, 0x1.bb9611b80e2fbp-2},
107 {-0x1.bfc00b8f3feaap-56, 0x1.c1c60693fa39ep-2},
108 {-0x1.bfc00b8f3feaap-56, 0x1.c1c60693fa39ep-2},
109 {0x1.223eadb651b4ap-57, 0x1.c7ff9c74554c9p-2},
110 {0x1.0798270b29f39p-56, 0x1.ce42f18064743p-2},
111 {0x1.d7f4d3b3d406bp-56, 0x1.d490246defa6bp-2},
112 {-0x1.0b5837185a661p-56, 0x1.dae75484c9616p-2},
113 {-0x1.ac81cc8a4dfb8p-56, 0x1.e148a1a2726cep-2},
114 {-0x1.ac81cc8a4dfb8p-56, 0x1.e148a1a2726cep-2},
115 {0x1.57d646a17bc6ap-56, 0x1.e7b42c3ddad73p-2},
116 {-0x1.74b71fb5e57e3p-62, 0x1.ee2a156b413e5p-2},
117 {-0x1.0d487f5aba5e5p-57, 0x1.f4aa7ee03192dp-2},
118 {-0x1.0d487f5aba5e5p-57, 0x1.f4aa7ee03192dp-2},
119 {0x1.7e8f05924d259p-57, 0x1.fb358af7a4884p-2},
120 {0x1.1713a36138e19p-57, 0x1.00e5ae5b207abp-1},
121 {-0x1.17f9e54e78104p-57, 0x1.04360be7603adp-1},
122 {-0x1.17f9e54e78104p-57, 0x1.04360be7603adp-1},
123 {0x1.2241edf5fd1f7p-57, 0x1.078bf0533c568p-1},
124 {0x1.0d710fcfc4e0dp-55, 0x1.0ae76e2d054fap-1},
125 {0x1.0d710fcfc4e0dp-55, 0x1.0ae76e2d054fap-1},
126 {0x1.3300f002e836ep-55, 0x1.0e4898611cce1p-1},
127 {-0x1.91eee7772c7c2p-55, 0x1.11af823c75aa8p-1},
128 {-0x1.91eee7772c7c2p-55, 0x1.11af823c75aa8p-1},
129 {0x1.342eb628dba17p-56, 0x1.151c3f6f29612p-1},
130 {0x1.89df1568ca0bp-55, 0x1.188ee40f23ca6p-1},
131 {0x1.89df1568ca0bp-55, 0x1.188ee40f23ca6p-1},
132 {0x1.59bddae1ccce2p-56, 0x1.1c07849ae6007p-1},
133 {-0x1.2164ff40e9817p-56, 0x1.1f8635fc61659p-1},
134 {-0x1.2164ff40e9817p-56, 0x1.1f8635fc61659p-1},
135 {-0x1.fcc8dbccc25cbp-57, 0x1.230b0d8bebc98p-1},
136 {0x1.e0efadd9db02bp-55, 0x1.269621134db92p-1},
137 {0x1.e0efadd9db02bp-55, 0x1.269621134db92p-1},
138 {-0x1.6a0c343be95dcp-56, 0x1.2a2786d0ec107p-1},
139 {-0x1.b941ee770436bp-56, 0x1.2dbf557b0df43p-1},
140 {-0x1.b941ee770436bp-56, 0x1.2dbf557b0df43p-1},
141 {0x1.6c3a5f12642c9p-57, 0x1.315da4434068bp-1},
142 {0x1.6c3a5f12642c9p-57, 0x1.315da4434068bp-1},
143 {-0x1.f01ab6065515cp-56, 0x1.35028ad9d8c86p-1},
144 {0x1.21512aa596ea3p-55, 0x1.38ae2171976e7p-1},
145 {0x1.21512aa596ea3p-55, 0x1.38ae2171976e7p-1},
146 {0x1.1930603d87b6ep-56, 0x1.3c6080c36bfb5p-1},
147 {0x1.1930603d87b6ep-56, 0x1.3c6080c36bfb5p-1},
148 {0x1.86cf0f38b461ap-57, 0x1.4019c2125ca93p-1},
149 {-0x1.84f481051f71ap-56, 0x1.43d9ff2f923c5p-1},
150 {-0x1.84f481051f71ap-56, 0x1.43d9ff2f923c5p-1},
151 {0x1.2541aca7d5844p-55, 0x1.47a1527e8a2d3p-1},
152 {0x1.2541aca7d5844p-55, 0x1.47a1527e8a2d3p-1},
153 {0x1.c457b531506f6p-55, 0x1.4b6fd6f970c1fp-1},
154 {0x1.c457b531506f6p-55, 0x1.4b6fd6f970c1fp-1},
155 {0x1.d749362382a77p-56, 0x1.4f45a835a4e19p-1},
156 {0x1.d749362382a77p-56, 0x1.4f45a835a4e19p-1},
157 {0x1.988ba4aea614dp-56, 0x1.5322e26867857p-1},
158 {0x1.988ba4aea614dp-56, 0x1.5322e26867857p-1},
159 {0x1.80bff3303dd48p-55, 0x1.5707a26bb8c66p-1},
160 {0x1.80bff3303dd48p-55, 0x1.5707a26bb8c66p-1},
161 {-0x1.6714fbcd8135bp-55, 0x1.5af405c3649ep-1},
162 {-0x1.6714fbcd8135bp-55, 0x1.5af405c3649ep-1},
163 {0x1.1c066d235ee63p-56, 0x1.5ee82aa24192p-1},
164 {0.0, 0.0},
165 };
166
167 const LogRR LOG2_TABLE = {
168 // -log2(r) with 128-bit precision generated by SageMath with:
169 //
170 // for i in range(1, 127):
171 // r = 2^-8 * ceil( 2^8 * (1 - 2^(-8)) / (1 + i*2^(-7)) );
172 // s, m, e = RealField(128)(r).log2().sign_mantissa_exponent();
173 // print("{false,", e, ", MType({", hex(m % 2^64), ",", hex((m >> 64) %
174 // 2^64),
175 // "})},");
176 /* .step_1 = */ {
177 {false, 0, MType(0)},
178 {false, -134, MType({0xe8c251630adb856a, 0xb963dd107b993ada})},
179 {false, -133, MType({0xa41b08fbe05f82d0, 0xba1f7430f9aab1b2})},
180 {false, -132, MType({0x1f06c085bc1b865d, 0x8c25c7262b57c149})},
181 {false, -132, MType({0x2e1c07f0438ebac0, 0xbb9ca64ecac6aaef})},
182 {false, -132, MType({0xaacc0e21d6541224, 0xeb75e8f8ff5ff022})},
183 {false, -131, MType({0x31514aef39ce6303, 0x8dd9953002a4e866})},
184 {false, -131, MType({0x50799beaaab2940c, 0xa62b07f3457c4070})},
185 {false, -131, MType({0xda288fc615a727dc, 0xbeb024b67dda6339})},
186 {false, -131, MType({0x22dbbaced44516ce, 0xcb0657cd5dbe4f6f})},
187 {false, -131, MType({0xd939dceecdd9ce05, 0xe3da945b878e27d0})},
188 {false, -131, MType({0x9596a8e2e84c8f45, 0xfce4aee0e88b2749})},
189 {false, -130, MType({0x243efd9325954cfe, 0x84bf1c673032495d})},
190 {false, -130, MType({0x91d79938e7226384, 0x916d6e1559a4b696})},
191 {false, -130, MType({0x22563c9ed9462091, 0x9e37db2866f2850b})},
192 {false, -130, MType({0x3a53ca1181015ada, 0xa4a7c31dc6f9a5d5})},
193 {false, -130, MType({0x3eb8023eed65d601, 0xb19d45fa1be70855})},
194 {false, -130, MType({0xce5cabbd2d753d9b, 0xb823018e3cfc25f0})},
195 {false, -130, MType({0x54dbf16fb0695ee3, 0xc544c055fde99333})},
196 {false, -130, MType({0x5196a85a067c6739, 0xcbe0e589e3f6042d})},
197 {false, -130, MType({0xf349845e48955078, 0xd930124bea9a2c66})},
198 {false, -130, MType({0x815ef705cfaef035, 0xdfe33d3fffa66037})},
199 {false, -130, MType({0x2ba704dcaa76f41d, 0xed61169f220e97f2})},
200 {false, -130, MType({0x2062f36bc14d0d93, 0xf42be9e9b09b3def})},
201 {false, -129, MType({0x132880194144b02b, 0x80ecdde7d30ea2ed})},
202 {false, -129, MType({0x54880de63812fd49, 0x845e706cafd1bf61})},
203 {false, -129, MType({0xa87c02eaf36e2c29, 0x8b4e029b1f8ac391})},
204 {false, -129, MType({0x9804237ec8d9431d, 0x8ecc164ea93841ae})},
205 {false, -129, MType({0x20f81ca95d9e7968, 0x924e69589e6b6268})},
206 {false, -129, MType({0x124bc6f1acf95dc4, 0x995ff71b8773432d})},
207 {false, -129, MType({0x5a5e8e21bff3336b, 0x9cef470aacfb7bf9})},
208 {false, -129, MType({0x4e53fa3329f65894, 0xa08300be1f651473})},
209 {false, -129, MType({0x2742d7296a39eed6, 0xa7b7dd96762cc3c7})},
210 {false, -129, MType({0xf359c5544bc5e134, 0xab591735abc724e4})},
211 {false, -129, MType({0x6b6c874dd96e1d75, 0xaefee78f75707221})},
212 {false, -129, MType({0x21006678c0a5c390, 0xb2a95a4cc313bb59})},
213 {false, -129, MType({0x6d40900b25024b32, 0xb6587b432e47501b})},
214 {false, -129, MType({0x89e2eb553b279b3d, 0xbdc4f8167955698f})},
215 {false, -129, MType({0xd58525aad392ca50, 0xc1826c8608fe9951})},
216 {false, -129, MType({0x54dbf16fb0695ee3, 0xc544c055fde99333})},
217 {false, -129, MType({0x88d5eae3326327bb, 0xc90c004926e9dbfb})},
218 {false, -129, MType({0x46dfa05bddfded8c, 0xccd83954b6359379})},
219 {false, -129, MType({0xbfe9dbebf2e8a45e, 0xd47fcb8c0852f0c0})},
220 {false, -129, MType({0x7b11f1c5160c515c, 0xd85b3fa7a3407fa8})},
221 {false, -129, MType({0x1339e5677ec44dd0, 0xdc3be2bd8d837f7f})},
222 {false, -129, MType({0xea2b8c7bb0ee9c8b, 0xe021c2cf17ed9bdb})},
223 {false, -129, MType({0xaec562332791fe38, 0xe40cee16a2ff21c4})},
224 {false, -129, MType({0x71682ebacca79cfa, 0xe7fd7308d6895b14})},
225 {false, -129, MType({0xa5ad5ce9fb5a7bb6, 0xebf36055e1abc61e})},
226 {false, -129, MType({0x3225190531a852c5, 0xefeec4eac371584e})},
227 {false, -129, MType({0xda8ad649da21eab0, 0xf3efaff29c559a77})},
228 {false, -129, MType({0x4c3e2ea7c15c3d1e, 0xf7f630d808fc2ada})},
229 {false, -129, MType({0xbcb9bfa9852e0d35, 0xfc02574686680cc6})},
230 {false, -128, MType({0xce032f41d1e774e8, 0x800a1995f0019518})},
231 {false, -128, MType({0x9b39ffeebc29372a, 0x8215ea5cd3e4c4c7})},
232 {false, -128, MType({0x87f95f1befb6f806, 0x8424a6335c777e0b})},
233 {false, -128, MType({0xb987b42e3bb332a1, 0x8636557862acb7ce})},
234 {false, -128, MType({0x139a7ba83bf2d136, 0x884b00aef726cec5})},
235 {false, -128, MType({0x50799beaaab2941, 0x8a62b07f3457c407})},
236 {false, -128, MType({0x8bd744617e9b7d52, 0x8c7d6db7169e0cda})},
237 {false, -128, MType({0x46ad444333ceb10, 0x8e9b414b5a92a606})},
238 {false, -128, MType({0xef4c737fba4f5d66, 0x90bc345861bf3d52})},
239 {false, -128, MType({0xae441c09d761c549, 0x92e050231df57d6f})},
240 {false, -128, MType({0x6e36aa9ce90a3879, 0x95079e1a0382dc79})},
241 {false, -128, MType({0xefca1a184e93809, 0x973227d6027ebd8a})},
242 {false, -128, MType({0xefca1a184e93809, 0x973227d6027ebd8a})},
243 {false, -128, MType({0x124bc6f1acf95dc4, 0x995ff71b8773432d})},
244 {false, -128, MType({0x352bea51e58ea9e8, 0x9b9115db83a3dd2d})},
245 {false, -128, MType({0x266d6cdc959153bc, 0x9dc58e347d37696d})},
246 {false, -128, MType({0x4527d82c8214ddca, 0x9ffd6a73a78eaf35})},
247 {false, -128, MType({0x404cabb76d600e3c, 0xa238b5160413106e})},
248 {false, -128, MType({0x404cabb76d600e3c, 0xa238b5160413106e})},
249 {false, -128, MType({0xcab7d2ec23f0eef3, 0xa47778c98bcc86a1})},
250 {false, -128, MType({0x761c48dd859de2d3, 0xa6b9c06e6211646b})},
251 {false, -128, MType({0x7fd3b7d7e5d148bb, 0xa8ff971810a5e181})},
252 {false, -128, MType({0xc27c6780d92b4d11, 0xab49080ecda53208})},
253 {false, -128, MType({0xdb502402c94092cd, 0xad961ed0cb91d406})},
254 {false, -128, MType({0xdb502402c94092cd, 0xad961ed0cb91d406})},
255 {false, -128, MType({0x3432ef6b732b6843, 0xafe6e71393eeda29})},
256 {false, -128, MType({0xbb324da7e046e792, 0xb23b6cc56cc84c99})},
257 {false, -128, MType({0xb21709ce430c8e24, 0xb493bc0ec9954243})},
258 {false, -128, MType({0xb21709ce430c8e24, 0xb493bc0ec9954243})},
259 {false, -128, MType({0xe91ad16ecff10111, 0xb6efe153c7e319f6})},
260 {false, -128, MType({0xce31e481cd797e79, 0xb94fe935b83e3eb5})},
261 {false, -128, MType({0xda3e961a96c580fa, 0xbbb3e094b3d228d3})},
262 {false, -128, MType({0xda3e961a96c580fa, 0xbbb3e094b3d228d3})},
263 {false, -128, MType({0xf396598aae91499a, 0xbe1bd4913f3fda43})},
264 {false, -128, MType({0xae4cceb0f621941b, 0xc087d28dfb2febb8})},
265 {false, -128, MType({0xae4cceb0f621941b, 0xc087d28dfb2febb8})},
266 {false, -128, MType({0x6c1855c42078f81b, 0xc2f7e831632b6670})},
267 {false, -128, MType({0x169535fb8bf577c8, 0xc56c23679b4d206e})},
268 {false, -128, MType({0x169535fb8bf577c8, 0xc56c23679b4d206e})},
269 {false, -128, MType({0x3b24cecc60217942, 0xc7e492644d64237e})},
270 {false, -128, MType({0x3dc2687fcf939696, 0xca6143a49626d820})},
271 {false, -128, MType({0x3dc2687fcf939696, 0xca6143a49626d820})},
272 {false, -128, MType({0xa62e6add1a901a0, 0xcce245f1031e41fa})},
273 {false, -128, MType({0x5bb6e23138ad51e1, 0xcf67a85fa1f89a04})},
274 {false, -128, MType({0x5bb6e23138ad51e1, 0xcf67a85fa1f89a04})},
275 {false, -128, MType({0x7fc60a5103092bae, 0xd1f17a5621fb01ac})},
276 {false, -128, MType({0xbfe9dbebf2e8a45e, 0xd47fcb8c0852f0c0})},
277 {false, -128, MType({0xbfe9dbebf2e8a45e, 0xd47fcb8c0852f0c0})},
278 {false, -128, MType({0x8e2d7d378127d823, 0xd712ac0cf811659d})},
279 {false, -128, MType({0x5c1a7f14b168b365, 0xd9aa2c3b0ea3cbc1})},
280 {false, -128, MType({0x5c1a7f14b168b365, 0xd9aa2c3b0ea3cbc1})},
281 {false, -128, MType({0xb7579f0f8d3d514b, 0xdc465cd155a90942})},
282 {false, -128, MType({0xb7579f0f8d3d514b, 0xdc465cd155a90942})},
283 {false, -128, MType({0xb087205eb55aea85, 0xdee74ee64b0c38d3})},
284 {false, -128, MType({0x424a2623d60dfb16, 0xe18d13ee805a4de3})},
285 {false, -128, MType({0x424a2623d60dfb16, 0xe18d13ee805a4de3})},
286 {false, -128, MType({0x4d3a591ae6854787, 0xe437bdbf5254459c})},
287 {false, -128, MType({0x4d3a591ae6854787, 0xe437bdbf5254459c})},
288 {false, -128, MType({0x8dcdb6b24c5c5cdf, 0xe6e75e91b9cca551})},
289 {false, -128, MType({0x33ac7d9ebba8a53c, 0xe99c090536ece983})},
290 {false, -128, MType({0x33ac7d9ebba8a53c, 0xe99c090536ece983})},
291 {false, -128, MType({0xfb2eede4b59d8959, 0xec55d022d80e3d27})},
292 {false, -128, MType({0xfb2eede4b59d8959, 0xec55d022d80e3d27})},
293 {false, -128, MType({0x308b454666de8f99, 0xef14c7605d60654c})},
294 {false, -128, MType({0x308b454666de8f99, 0xef14c7605d60654c})},
295 {false, -128, MType({0x8383cb0ce23bebd4, 0xf1d902a37aaa5085})},
296 {false, -128, MType({0x8383cb0ce23bebd4, 0xf1d902a37aaa5085})},
297 {false, -128, MType({0x64fc87b4a41f7b70, 0xf4a2964538813c67})},
298 {false, -128, MType({0x64fc87b4a41f7b70, 0xf4a2964538813c67})},
299 {false, -128, MType({0x3f5d7d82b65c5686, 0xf77197157665f689})},
300 {false, -128, MType({0x3f5d7d82b65c5686, 0xf77197157665f689})},
301 {false, -128, MType({0x6476077b9fbd41ae, 0xfa461a5e8f4b759d})},
302 {false, -128, MType({0x6476077b9fbd41ae, 0xfa461a5e8f4b759d})},
303 {false, -128, MType({0xe3909ffd0d61778, 0xfd2035e9221ef5d0})},
304 {false, 0, MType(0)},
305 },
306 // -log2(r) for the second step, generated by SageMath with:
307 //
308 // for i in range(-2^6, 2^7 + 1):
309 // r = 2^-16 * round( 2^16 / (1 + i*2^(-14)) );
310 // s, m, e = RealField(128)(r).log2().sign_mantissa_exponent();
311 // print("{true," if s == 1 else "{false,", e, ",
312 // MType({", hex(m % 2^64), ",", hex((m >> 64) % 2^64), "})},");
313 /* .step_2 = */
314 {
315 {true, -135, MType({0xb5cfed58337e848a, 0xb906155918954401})},
316 {true, -135, MType({0xffaf2ac1b1d20910, 0xb6264958a3c7fa2b})},
317 {true, -135, MType({0x52521a3950ea2ed8, 0xb34671e439aa448e})},
318 {true, -135, MType({0xf87e1abdee10fd95, 0xb0668efb7ef48ab7})},
319 {true, -135, MType({0xfbd43bbcc24c5e43, 0xad86a09e185af0e8})},
320 {true, -135, MType({0x2f4f5d48f9796742, 0xaaa6a6cbaa8d57ce})},
321 {true, -135, MType({0x3477fd67c1cab6b3, 0xa7c6a183da375c3d})},
322 {true, -135, MType({0x7b4d33eb381fe558, 0xa4e690c64c0056f0})},
323 {true, -135, MType({0x3ce25e48cb498dea, 0xa2067492a48b5c43})},
324 {true, -135, MType({0x70b0fcc9e4330983, 0x9f264ce888773bed})},
325 {true, -135, MType({0xbc9e4267d3189b22, 0x9c4619c79c5e80bf})},
326 {true, -135, MType({0x5fb3d896326615c4, 0x9965db2f84d7705f})},
327 {true, -135, MType({0x178b58311e96d323, 0x9685911fe6740b02})},
328 {true, -135, MType({0x6bf8b6cf73d847, 0x93a53b9865c20b2a})},
329 {true, -135, MType({0x7019f6e64a580a02, 0x90c4da98a74ae561})},
330 {true, -135, MType({0xcb5733cf0eb4191d, 0x8de46e204f93c7f6})},
331 {true, -135, MType({0x56148d4fc5e415b6, 0x8b03f62f031d9ab8})},
332 {true, -135, MType({0xfe5370f425872623, 0x882372c46664feaf})},
333 {true, -135, MType({0x21b72a1457ee70d6, 0x8542e3e01de24ddf})},
334 {true, -135, MType({0xabff4f89968bed0b, 0x81aa211f1e332fcf})},
335 {true, -136, MType({0x86410a676480a5a7, 0xfd92f0cf88d75f24})},
336 {true, -136, MType({0x44280889021970e4, 0xf7d1886b2a876289})},
337 {true, -136, MType({0x32eb139d9812090d, 0xf21009106a42bc14})},
338 {true, -136, MType({0xbef9dd41e8e42810, 0xec4e72be90cd2d2d})},
339 {true, -136, MType({0x689d08ca6c7c3eb1, 0xe68cc574e6e1e5d7})},
340 {true, -136, MType({0x1ef259a7f69821d, 0xe0cb0132b5338423})},
341 {true, -136, MType({0xe22cea71b7bb8467, 0xdb0925f7446c13a9})},
342 {true, -136, MType({0xe5bb27303f542fe, 0xd54733c1dd2d0d04})},
343 {true, -136, MType({0x57453c8d5dc64ce1, 0xcf852a91c80f553f})},
344 {true, -136, MType({0x6cc7add1fc09ef92, 0xc9c30a664da33d56})},
345 {true, -136, MType({0xe678d7280de1c07f, 0xc400d33eb67081a7})},
346 {true, -136, MType({0x419bbeb2239bdc39, 0xbe3e851a4af6496d})},
347 {true, -136, MType({0xd4676d1d81755809, 0xb87c1ff853ab2631})},
348 {true, -136, MType({0xb69dfef7ac2e2890, 0xb2b9a3d818fd1349})},
349 {true, -136, MType({0x9f72fa0a8fccabc0, 0xacf710b8e3517548})},
350 {true, -136, MType({0xb8bfe6a3addb988e, 0xa7346699fb051978})},
351 {true, -136, MType({0x67862c8ec9dcd60d, 0xa171a57aa86c3551})},
352 {true, -136, MType({0x9bd3370909e28a6, 0x9baecd5a33d265ee})},
353 {true, -136, MType({0xa96bc611b991419b, 0x95ebde37e57aaf84})},
354 {true, -136, MType({0xa50bb80f203f0d62, 0x9028d813059f7cdc})},
355 {true, -136, MType({0x4d36cd474f65a317, 0x8a65baeadc729ec5})},
356 {true, -136, MType({0x779be241ef4874a3, 0x84a286beb21d4b8c})},
357 {true, -137, MType({0xe76a962fa65ace3, 0xfdbe771b9d803cea})},
358 {true, -137, MType({0xd3d35627464a5267, 0xf237b2aef4e62e5a})},
359 {true, -137, MType({0x162ef4b0e838c363, 0xe6b0c035fa8b328c})},
360 {true, -137, MType({0x77bb10b976b3b9ca, 0xdb299faf3e7cd74f})},
361 {true, -137, MType({0x209853cee70bc58b, 0xcfa2511950b77014})},
362 {true, -137, MType({0x63f9b57cbaf2e58d, 0xc41ad472c12614d3})},
363 {true, -137, MType({0x4fca1c931bd6e6d6, 0xb89329ba1fa2a0fd})},
364 {true, -137, MType({0x26d26e434a53490a, 0xad0b50edfbf5b265})},
365 {true, -137, MType({0xc55e079078dc86a0, 0xa1834a0ce5d6a82d})},
366 {true, -137, MType({0xf05b9d5bd28f540b, 0x95fb15156ceba1b5})},
367 {true, -137, MType({0x8ef87f1a11cdb727, 0x8a72b20620c97d84})},
368 {true, -138, MType({0x9d6870114c1183cf, 0xfdd441bb21e7b069})},
369 {true, -138, MType({0x63d514fff97e86f3, 0xe6c2c33499ba16c4})},
370 {true, -138, MType({0x11a381901eadd883, 0xcfb0e875c7cc5929})},
371 {true, -138, MType({0xa9d69d37bc0a5bac, 0xb89eb17bcabe1857})},
372 {true, -138, MType({0x2dc97c9ffefd2497, 0xa18c1e43c10c6898})},
373 {true, -138, MType({0xdcdc8afcb2ac09a, 0x8a792ecac911cf92})},
374 {true, -139, MType({0xdd454eb3a1489470, 0xe6cbc61c020c8446})},
375 {true, -139, MType({0x878035864d84b319, 0xb8a476150dfe4470})},
376 {true, -139, MType({0x7ce595cc53b8342c, 0x8a7c6d7af1de7942})},
377 {true, -140, MType({0x4710b59049899141, 0xb8a7588fd29b1baa})},
378 {true, -141, MType({0x5957f633309d74e3, 0xb8a8c9d8be9ae994})},
379 {false, 0, MType({0x0, 0x0})},
380 {false, -141, MType({0x8268aba030b1adf6, 0xb8abac81ab576f3b})},
381 {false, -140, MType({0x1511cba2fb213a10, 0xb8ad1de1ac9ea6a5})},
382 {false, -139, MType({0x6379fb9fd9bc6235, 0x8a82eb7708262500})},
383 {false, -139, MType({0xb6fe1bf601ee27d5, 0xb8b000b8c65957cc})},
384 {false, -139, MType({0x8c6e60693a14e6d0, 0xe6ddcebbd72d3f7f})},
385 {false, -138, MType({0xe9bcfd0c62eaa2ca, 0x8a862ac30095c084})},
386 {false, -138, MType({0x73b214209a5234a7, 0xa19dca8e85918b6d})},
387 {false, -138, MType({0x347d4ca3109fe4db, 0xb8b5c6c35e142a9b})},
388 {false, -138, MType({0x37a62c48783bb066, 0xcfce1f646dca7745})},
389 {false, -138, MType({0x794b6437fb56344, 0xe6e6d4749883fbe3})},
390 {false, -138, MType({0x1cb9a45ed90318e6, 0xfdffe5f6c232f658})},
391 {false, -137, MType({0xbc118e5dbbef7dbc, 0x8a8ca9f6e7762d0f})},
392 {false, -137, MType({0xb4c0fb9535907cf8, 0x96198f2e5173e93b})},
393 {false, -137, MType({0xc051d2c5f00a9bb9, 0xa1a6a2a3113fe246})},
394 {false, -137, MType({0x553269878c1e5110, 0xad33e4569918a8d5})},
395 {false, -137, MType({0xbc906750b0ce372c, 0xb8c1544a5b4e2caf})},
396 {false, -137, MType({0x4c50eaa63be294b6, 0xc44ef27fca41bdd8})},
397 {false, -137, MType({0xb6cb28db8c065b44, 0xcfdcbef858660da1})},
398 {false, -137, MType({0x70479336830ceb05, 0xdb6ab9b5783f2fc5})},
399 {false, -137, MType({0x2a458c831f6aeb49, 0xe6f8e2b89c629b7a})},
400 {false, -137, MType({0x6489ba5bd391e206, 0xf2873a0337772c8a})},
401 {false, -137, MType({0x13f6fda510aeec3b, 0xfe15bf96bc35246b})},
402 {false, -136, MType({0x2f9a0ef9e8250836, 0x84d239ba4eb315a9})},
403 {false, -136, MType({0x389019e822b70f1e, 0x8a99aacf26f2a8a7})},
404 {false, -136, MType({0x308beeffa12cf669, 0x9061330aa04f87ae})},
405 {false, -136, MType({0x9886a71b25a2085d, 0x9628d26d7448a43f})},
406 {false, -136, MType({0x70ba9cebe0b969c3, 0x9bf088f85c65a56b})},
407 {false, -136, MType({0xcd855dc705ea2bea, 0xa1b856ac1236e85b})},
408 {false, -136, MType({0x7736196b11afb331, 0xa7803b894f5580e0})},
409 {false, -136, MType({0x94c99761b8eab3d8, 0xad483790cd6339fa})},
410 {false, -136, MType({0x6194b8c040814736, 0xb3104ac3460a9668})},
411 {false, -136, MType({0xedde8d24c7a999cc, 0xb8d8752172fed130})},
412 {false, -136, MType({0xea6b01ebde42f1d0, 0xbea0b6ac0dfbde2f})},
413 {false, -136, MType({0x7ef732b69334cf50, 0xc4690f63d0c66aa1})},
414 {false, -136, MType({0x2ba86275fcfc2d72, 0xca317f49752bddae})},
415 {false, -136, MType({0xb56ea44e185bf99f, 0xcffa065db50258f6})},
416 {false, -136, MType({0x1d5c3bbeb6902bfe, 0xd5c2a4a14a28b920})},
417 {false, -136, MType({0xa2f2bb9e156b0f37, 0xdb8b5a14ee86965f})},
418 {false, -136, MType({0xd166eb8da06ab5ef, 0xe15426b95c0c4506})},
419 {false, -136, MType({0x97dc7bae4219de0f, 0xe71d0a8f4cb2d60f})},
420 {false, -136, MType({0x6c9a8e7698f416c4, 0xece605977a7c17a8})},
421 {false, -136, MType({0x7b3a20aa5289695e, 0xf2af17d29f7295c0})},
422 {false, -136, MType({0xddcf578ee2c2897b, 0xf878414175a99a93})},
423 {false, -136, MType({0xe10ebd96c3ec30ec, 0xfe4181e4b73d2f37})},
424 {false, -135, MType({0xa9b7baecb34ba577, 0x82056cde8f290e13})},
425 {false, -135, MType({0x2da910dc61c182da, 0x8430f56d5e1edfd1})},
426 {false, -135, MType({0xfaca09dc7e0ba8b5, 0x8715b5a8f27bed90})},
427 {false, -135, MType({0xd723876173c0947, 0x89fa818019a2cace})},
428 {false, -135, MType({0x4e6651df154e8f8c, 0x8cdf58f330b64515})},
429 {false, -135, MType({0xee54b77d3bc34b6d, 0x8fc43c0294dd8af3})},
430 {false, -135, MType({0xad07dde9b5f92cce, 0x92a92aaea3442c3d})},
431 {false, -135, MType({0x261aacf944b638f0, 0x958e24f7b91a1a53})},
432 {false, -135, MType({0x232f5d64a85b219d, 0x98732ade3393a868})},
433 {false, -135, MType({0xf3a958bb706093fc, 0x9b583c626fe98bc9})},
434 {false, -135, MType({0xc9eaa059e7b0333a, 0x9e3d5984cb58dc25})},
435 {false, -135, MType({0x1e154029663243c0, 0xa1228245a32313cf})},
436 {false, -135, MType({0x16515200e283d006, 0xa407b6a5548e1006})},
437 {false, -135, MType({0xf498168a3337ca4f, 0xa6ecf6a43ce4113d})},
438 {false, -135, MType({0x8a04a89f0548a10f, 0xa9d24242b973bb63})},
439 {false, -135, MType({0xafaad01f25772805, 0xacb7998127901623})},
440 {false, -135, MType({0xc4f47950543fe0b8, 0xaf9cfc5fe4908d31})},
441 {false, -135, MType({0x338655e677d0d3ec, 0xb2826adf4dd0f08e})},
442 {false, -135, MType({0xf8ac2ce19d009541, 0xb567e4ffc0b174cc})},
443 {false, -135, MType({0x344d5e7dd7b2f465, 0xb84d6ac19a96b35c})},
444 {false, -135, MType({0xbd6a217fb4598ec7, 0xbb32fc2538e9aaca})},
445 {false, -135, MType({0xbc21ff368f562b75, 0xbe18992af917bf0e})},
446 {false, -135, MType({0x4944139ccbf2cb9a, 0xc0fe41d33892b9cc})},
447 {false, -135, MType({0x1369970c8b67e6b5, 0xc3e3f61e54d0ca9c})},
448 {false, -135, MType({0x99b370e2d04a530, 0xc6c9b60cab4c8752})},
449 {false, -135, MType({0xb81c3d48aff589f, 0xc9af819e9984ec44})},
450 {false, -135, MType({0x9f22b80993be311b, 0xcc9558d47cfd5c90})},
451 {false, -135, MType({0xac29209c8d8985ae, 0xcf7b3baeb33da265})},
452 {false, -135, MType({0x3cbb6a520292351d, 0xd2612a2d99d1ef47})},
453 {false, -135, MType({0x43de9ae40507ef24, 0xd54724518e4adc56})},
454 {false, -135, MType({0x69677b902ea4df3a, 0xd82d2a1aee3d6a97})},
455 {false, -135, MType({0xdb7a3aff74967bd5, 0xdb133b8a17430339})},
456 {false, -135, MType({0x25990c82a0066ac6, 0xddf9589f66f977de})},
457 {false, -135, MType({0xd424aacf4babf55, 0xe0df815b3b0302dd})},
458 {false, -135, MType({0xf8e3e7eb5a7bdebb, 0xe30c278d9936c595})},
459 {false, -135, MType({0x5ef8bf5adf5deebe, 0xe5f264adb62d5810})},
460 {false, -135, MType({0x331d19965368fc82, 0xe8d8ad75590bdf92})},
461 {false, -135, MType({0x901c30c427e358b8, 0xebbf01e4df85219e})},
462 {false, -135, MType({0xaeac7e9857253b06, 0xeea561fca7504dc1})},
463 {false, -135, MType({0xe2113e5893ab5b40, 0xf18bcdbd0e28fdd7})},
464 {false, -135, MType({0x9a4efc80ae977826, 0xf472452671cf3654})},
465 {false, -135, MType({0x6bf3ba8319332c9f, 0xf758c83930076689})},
466 {false, -135, MType({0x1d732d302e75018b, 0xfa3f56f5a69a68ed})},
467 {false, -135, MType({0xba179c5dbcceec01, 0xfd25f15c33558362})},
468 {false, -134, MType({0x5543f53b8ad85039, 0x80064bb69a0533c0})},
469 {false, -134, MType({0xe971a5565b93cb67, 0x8179a4948347996b})},
470 {false, -134, MType({0x5b399644ba714691, 0x82ed0348045f379d})},
471 {false, -134, MType({0x5079f1e0ec4b8496, 0x846067d14c3b8982})},
472 {false, -134, MType({0x6aba4990a32e8873, 0x85d3d23089ce40b0})},
473 {false, -134, MType({0xe16770c3a404291c, 0x87474265ec0b4548})},
474 {false, -134, MType({0x1edb7ffb1d6b3eab, 0x88bab871a1e8b61c})},
475 {false, -134, MType({0x603243e1ba7c7865, 0x8a2e3453da5ee8cd})},
476 {false, -134, MType({0x57ea5c03ea4621dd, 0x8ba1b60cc46869f6})},
477 {false, -134, MType({0xd3534cbf43bd7fd8, 0x8d153d9c8f01fd4a})},
478 {false, -134, MType({0x62c8c8075dc91cd5, 0x8e88cb03692a9dbc})},
479 {false, -134, MType({0x4bb70a5e3db7b85, 0x8ffc5e4181e37d9e})},
480 {false, -134, MType({0xd3875ba32159547a, 0x916ff757083006c7})},
481 {false, -134, MType({0x5c94c80e7a8f66b1, 0x9286adfca91ba28d})},
482 {false, -134, MType({0x52d313c47b4f91db, 0x93fa514ba0517623})},
483 {false, -134, MType({0x80829e9f3957a4c3, 0x956dfa72866fc57d})},
484 {false, -134, MType({0x1cd4917972015ae7, 0x96e1a9718a824be5})},
485 {false, -134, MType({0x1af23c29ef3032da, 0x98555e48db96fcd2})},
486 {false, -134, MType({0xe7f7bf240be67b80, 0x99c918f8a8be040e})},
487 {false, -134, MType({0x2bbe3cd4f7d868fa, 0x9b3cd9812109c5dc})},
488 {false, -134, MType({0x8c75d6a4c5ae460d, 0x9cb09fe2738edf14})},
489 {false, -134, MType({0x750fb989c9a06186, 0x9e246c1ccf642550})},
490 {false, -134, MType({0xde787e244901bdf9, 0x9f983e3063a2a709})},
491 {false, -134, MType({0x1ba3205ff729efa4, 0xa10c161d5f65abc0})},
492 {false, -134, MType({0xa864d2a038fb19cd, 0xa27ff3e3f1cab41b})},
493 {false, -134, MType({0xfb21f083a5fec56d, 0xa3f3d78449f17a11})},
494 {false, -134, MType({0x594c5552bcc377f5, 0xa567c0fe96fbf109})},
495 {false, -134, MType({0xaeb35a353fc5a503, 0xa6dbb053080e45fc})},
496 {false, -134, MType({0x67a5c05130c0f330, 0xa84fa581cc4edf9f})},
497 {false, -134, MType({0x4de5cafde1caf46f, 0xa9c3a08b12e65e81})},
498 {false, -134, MType({0x686fce3d160e88fd, 0xab37a16f0aff9d32})},
499 {false, -134, MType({0xde1375b3af6749a6, 0xacaba82de3c7b066})},
500 {false, -134, MType({0x243569048ac4affe, 0xadc2b114c632da56})},
501 {false, -134, MType({0xd6796227dcd39551, 0xaf36c21319b80ea2})},
502 {false, -134, MType({0xabc9265386172074, 0xb0aad8eccfb38d51})},
503 {false, -134, MType({0xcaac9f17896f2ce, 0xb21ef5a2175ac65e})},
504 {false, -134, MType({0x1c65a3c7f828972b, 0xb39318331fe56492})},
505 {false, -134, MType({0xabdc66446a4286d9, 0xb50740a0188d4daa})},
506 {false, -134, MType({0x2f3bbe8e8d72abec, 0xb67b6ee9308ea27b})},
507 {false, -134, MType({0xb67dbdd7f03d168c, 0xb7efa30e9727bf11})},
508 },
509 // -log2(r) for the third step, generated by SageMath with:
510 //
511 // for i in range(-80, 81):
512 // r = 2^-21 * round( 2^21 / (1 + i*2^(-21)) );
513 // s, m, e = RealField(128)(r).log2().sign_mantissa_exponent();
514 // print("{true," if (s == 1) else "{false,", e, ",
515 // MType({", hex(m % 2^64), ",", hex((m >> 64) % 2^64), "})},");
516 /* .step_3 = */
517 {
518 {true, -142, MType({0x26f2c63c0827ccbb, 0xe6d3a96b978fc16e})},
519 {true, -142, MType({0x4b56fe667c8ec091, 0xe3f107a9fbfc50ca})},
520 {true, -142, MType({0x647d76181aec10fc, 0xe10e65d14b937265})},
521 {true, -142, MType({0x99e8f4d5379eca79, 0xde2bc3e18653b4f5})},
522 {true, -142, MType({0xf07da89990c20623, 0xdb4921daac3ba730})},
523 {true, -142, MType({0x4a8121848531851a, 0xd8667fbcbd49d7cd})},
524 {true, -142, MType({0x679a4d854ae13619, 0xd583dd87b97cd580})},
525 {true, -142, MType({0xe4d174072487a514, 0xd2a13b3ba0d32eff})},
526 {true, -142, MType({0x3c90319d969b54be, 0xcfbe98d8734b7301})},
527 {true, -142, MType({0xc6a173b09ba301e6, 0xccdbf65e30e43039})},
528 {true, -142, MType({0xb8317428d7d8d06b, 0xc9f953ccd99bf55e})},
529 {true, -142, MType({0x23cdb51bcc2061cd, 0xc716b1246d715125})},
530 {true, -142, MType({0xf964fc78084fd515, 0xc4340e64ec62d241})},
531 {true, -142, MType({0x6474fb15ccbb015, 0xc1516b8e566f076a})},
532 {true, -142, MType({0xf525ef6d0b75b1c3, 0xbe6ec8a0ab947f51})},
533 {true, -142, MType({0x4e13532df7ee8da7, 0xbb8c259bebd1c8ae})},
534 {true, -142, MType({0x76832500d72a9027, 0xb8a9828017257233})},
535 {true, -142, MType({0xb14a3d285e592ba0, 0xb5c6df4d2d8e0a95})},
536 {true, -142, MType({0x1e9e9dc9711f6e20, 0xb2e43c032f0a2089})},
537 {true, -142, MType({0xbc176e974f255fac, 0xb00198a21b9842c1})},
538 {true, -142, MType({0x64acf87fc0f648e6, 0xad1ef529f336fff3})},
539 {true, -142, MType({0xd0b8a1574433e1f8, 0xaa3c519ab5e4e6d1})},
540 {true, -142, MType({0x95f4e785371c69a9, 0xa759adf463a08610})},
541 {true, -142, MType({0x277d5db00363a46f, 0xa4770a36fc686c63})},
542 {true, -142, MType({0xd5cea669485ec36c, 0xa1946662803b287c})},
543 {true, -142, MType({0xcec66fda04833322, 0x9eb1c276ef174910})},
544 {true, -142, MType({0x1da36f6ebe3851db, 0x9bcf1e7448fb5cd2})},
545 {true, -142, MType({0xab055d83abfc0d82, 0x98ec7a5a8de5f273})},
546 {true, -142, MType({0x3cecf110dbda68e9, 0x9609d629bdd598a8})},
547 {true, -142, MType({0x76bbdb565a37e84b, 0x932731e1d8c8de22})},
548 {true, -142, MType({0xd934c38857eee4f3, 0x90448d82debe5194})},
549 {true, -142, MType({0xc27b427b4fbfc7db, 0x8d61e90ccfb481b1})},
550 {true, -142, MType({0x6e13de502b142b39, 0x8a7f447faba9fd2b})},
551 {true, -142, MType({0xf4e406206614e2ba, 0x879c9fdb729d52b3})},
552 {true, -142, MType({0x4d320daa3312ea6c, 0x84b9fb20248d10fd})},
553 {true, -142, MType({0x4aa528fc9d433c1a, 0x81d7564dc177c6b9})},
554 {true, -143, MType({0x3c8ad047559b1622, 0xfde962c892b80533})},
555 {true, -143, MType({0xacf765a8fc5bcc31, 0xf82418c77870a69f})},
556 {true, -143, MType({0xbe238832edd27f20, 0xf25ece9834168f1a})},
557 {true, -143, MType({0x2644bfca329b708, 0xec99843ac5a6dc07})},
558 {true, -143, MType({0xc6d05a788e614744, 0xe6d439af2d1eaac6})},
559 {true, -143, MType({0x133fe9cc57a8c1d0, 0xe10eeef56a7b18bc})},
560 {true, -143, MType({0xaa4cb429195fb5dd, 0xdb49a40d7db94348})},
561 {true, -143, MType({0x951ef239abbb959, 0xd58458f766d647ce})},
562 {true, -143, MType({0x686c430c89143d35, 0xcfbf0db325cf43ad})},
563 {true, -143, MType({0xba79c248afd42c12, 0xc9f9c240baa15447})},
564 {true, -143, MType({0xad19e0a92f115327, 0xc43476a0254996fd})},
565 {true, -143, MType({0xa8ad6ac3b0c99520, 0xbe6f2ad165c5292f})},
566 {true, -143, MType({0xd0567d4a9cc5e6a1, 0xb8a9ded47c11283d})},
567 {true, -143, MType({0x1f87c654b231443, 0xb2e492a9682ab188})},
568 {true, -143, MType({0xd6380b08358051bc, 0xad1f46502a0ee26d})},
569 {true, -143, MType({0xa07b024d26d391f6, 0xa759f9c8c1bad84e})},
570 {true, -143, MType({0x6ee868cb69e3a7d8, 0xa194ad132f2bb089})},
571 {true, -143, MType({0xa6869eff6682f73, 0x9bcf602f725e887d})},
572 {true, -143, MType({0xf6a44d559ccf3f61, 0x960a131d8b507d87})},
573 {true, -143, MType({0x72066e1d30a8e210, 0x9044c5dd79fead08})},
574 {true, -143, MType({0x75ba3245b1b856af, 0x8a7f786f3e66345c})},
575 {true, -143, MType({0xb5ac020473ab198f, 0x84ba2ad2d88430e1})},
576 {true, -144, MType({0x41127e3a88eb6741, 0xfde9ba1090ab7feb})},
577 {true, -144, MType({0xbf80787522aca1c4, 0xf25f1e1f1baffdea})},
578 {true, -144, MType({0xaf00688b14fa3adc, 0xe6d481d15210167b})},
579 {true, -144, MType({0x4d72837c8ab4d1e5, 0xdb49e52733c60457})},
580 {true, -144, MType({0x4e38ac27bb252090, 0xcfbf4820c0cc0236})},
581 {true, -144, MType({0xda3661f9292f59e8, 0xc434aabdf91c4ad0})},
582 {true, -144, MType({0x8fd0af9bdfd21488, 0xb8aa0cfedcb118de})},
583 {true, -144, MType({0x82ee19a9abf0bfa5, 0xad1f6ee36b84a716})},
584 {true, -144, MType({0x3cf68d5b5369a251, 0xa194d06ba591302f})},
585 {true, -144, MType({0xbcd34f38c977647e, 0x960a31978ad0eede})},
586 {true, -144, MType({0x76eee9c9605e2143, 0x8a7f92671b3e1dda})},
587 {true, -145, MType({0xaa6a3887f0c803ab, 0xfde9e5b4ada5efae})},
588 {true, -145, MType({0x6e25927e582ac191, 0xe6d4a5e27b136f13})},
589 {true, -145, MType({0xe2ebcac2f3a8e9eb, 0xcfbf65579eb92f4a})},
590 {true, -145, MType({0x9d9acc22d5690751, 0xb8aa2414188ba5bb})},
591 {true, -145, MType({0x1e12604b6d4132ef, 0xa194e217e87f47cb})},
592 {true, -145, MType({0xcf340d2acb9b92a9, 0x8a7f9f630e888add})},
593 {true, -146, MType({0xdc5e49fbde3c520, 0xe6d4b7eb1537c8ae})},
594 {true, -146, MType({0xc074c9557c01188, 0xb8aa2f9eb95b9332})},
595 {true, -146, MType({0xf0f82818ff9b654f, 0x8a7fa5e109656009})},
596 {true, -147, MType({0xd4cd612078bbe9b0, 0xb8aa35640a7c33eb})},
597 {true, -148, MType({0xf08cf68f42e09fa0, 0xb8aa3846b33aaecf})},
598 {false, 0, MType({0x0, 0x0})},
599 {false, -148, MType({0x68bd0facdf0ddaaf, 0xb8aa3e0c0513f9b1})},
600 {false, -147, MType({0x192af653dd41575b, 0xb8aa40eeae2ec9b3})},
601 {false, -146, MType({0x3b5c89842e540a51, 0x8a7fb2dd018e4892})},
602 {false, -146, MType({0x34ad8ebdd8b2750c, 0xb8aa46b400c0bee3})},
603 {false, -146, MType({0x70b12bd698e5be74, 0xe6d4dbfc54c5dd1b})},
604 {false, -145, MType({0x8c7e424efbd90e1, 0x8a7fb95afeda5c46})},
605 {false, -145, MType({0x31b8eba774a1de77, 0xa19505707dd23344})},
606 {false, -145, MType({0xee400e8c68838733, 0xb8aa523ea755fe32})},
607 {false, -145, MType({0xe71fa0b5603bc2f, 0xcfbf9fc57b7147be})},
608 {false, -145, MType({0x7763c919d8ac65f1, 0xe6d4ee04fa2f9a92})},
609 {false, -145, MType({0x232b270bb6046ec1, 0xfdea3cfd239c815e})},
610 {false, -144, MType({0x106f39197e068972, 0x8a7fc656fbe1c368})},
611 {false, -144, MType({0x4a4a6f4012941bd9, 0x960a6e8bbb581acc})},
612 {false, -144, MType({0x5bb34c1120b3e54b, 0xa195171cd0370c34})},
613 {false, -144, MType({0x6bb6731392a3147a, 0xad1fc00a3a845cf9})},
614 {false, -144, MType({0x2be1268dcee3c8fc, 0xb8aa6953fa45d275})},
615 {false, -144, MType({0xd84158d5d50251a9, 0xc43512fa0f813201})},
616 {false, -144, MType({0x3765bda15d0ef0fa, 0xcfbfbcfc7a3c40fa})},
617 {false, -144, MType({0x9a5ddb55f9cc27d9, 0xdb4a675b3a7cc4b9})},
618 {false, -144, MType({0xdcba1c593d918775, 0xe6d512165048829b})},
619 {false, -144, MType({0x648be060e1e30a95, 0xf25fbd2dbba53ffd})},
620 {false, -144, MType({0x22658dc2f1bcf6e8, 0xfdea68a17c98c23b})},
621 {false, -143, MType({0x48ad5162fb4a236e, 0x84ba8a38c9946759})},
622 {false, -143, MType({0xdb7fe3789405ce3a, 0x8a7fe04effad9560})},
623 {false, -143, MType({0x91b56e2e4f2e5ed8, 0x90453693609acde3})},
624 {false, -143, MType({0xf8998880c3bb4d76, 0x960a8d05ec5ef390})},
625 {false, -143, MType({0xe2b878052f67efee, 0x9bcfe3a6a2fce918})},
626 {false, -143, MType({0x67df399193f707c0, 0xa1953a758477912b})},
627 {false, -143, MType({0xe51b89e4d5d095e1, 0xa75a917290d1ce78})},
628 {false, -143, MType({0xfcbbee4edbf9f47d, 0xad1fe89dc80e83b1})},
629 {false, -143, MType({0x964fbd58b168371b, 0xb2e53ff72a309387})},
630 {false, -143, MType({0xdea7276ca7acd135, 0xb8aa977eb73ae0aa})},
631 {false, -143, MType({0x47d33f7e7afc83a6, 0xbe6fef346f304dcd})},
632 {false, -143, MType({0x892603b377909123, 0xc43547185213bda0})},
633 {false, -143, MType({0x9f32660aa06239fb, 0xc9fa9f2a5fe812d6})},
634 {false, -143, MType({0xcbcc5504d7407f6c, 0xcfbff76a98b03021})},
635 {false, -143, MType({0x9608c44d06402ebe, 0xd5854fd8fc6ef834})},
636 {false, -143, MType({0xca3db5604a863477, 0xdb4aa8758b274dc1})},
637 {false, -143, MType({0x7a024036206c37d6, 0xe110014044dc137c})},
638 {false, -143, MType({0xfc2e9be890ff7ee3, 0xe6d55a3929902c17})},
639 {false, -143, MType({0xecdc275c60da1b53, 0xec9ab36039467a47})},
640 {false, -143, MType({0x2d6571e94056607f, 0xf2600cb57401e0c0})},
641 {false, -143, MType({0xe4664401fd1ca2a7, 0xf8256638d9c54234})},
642 {false, -143, MType({0x7dbba7dcb50b3fd7, 0xfdeabfea6a93815a})},
643 {false, -142, MType({0xd541f90d853c794b, 0x81d80ce51337c072})},
644 {false, -142, MType({0xb08f65392ce8b75b, 0x84bab9ec06ae11c5})},
645 {false, -142, MType({0x6e969a29f8462436, 0x879d670a0fae2600})},
646 {false, -142, MType({0xcfc8cbcaa2bf130c, 0x8a80143f2e396e7d})},
647 {false, -142, MType({0xb737e48c19421e68, 0x8d62c18b62515c98})},
648 {false, -142, MType({0x2a9689b997c50c0b, 0x90456eeeabf761ac})},
649 {false, -142, MType({0x52381fccc774d66b, 0x93281c690b2cef13})},
650 {false, -142, MType({0x7910cec1dd92dc10, 0x960ac9fa7ff37629})},
651 {false, -142, MType({0xcb5866bbaff34cb, 0x98ed77a30a4c684a})},
652 {false, -142, MType({0x9d5c02c80c702d11, 0x9bd02562aa3936d0})},
653 {false, -142, MType({0xdddad0536b56e775, 0x9eb2d3395fbb5318})},
654 {false, -142, MType({0xa3a9505d7f71247a, 0xa19581272ad42e7e})},
655 {false, -142, MType({0xe6dfbd5d210830d7, 0xa4782f2c0b853a5d})},
656 {false, -142, MType({0xc2372f447bdcfa45, 0xa75add4801cfe812})},
657 {false, -142, MType({0x73099fd532c14b05, 0xaa3d8b7b0db5a8f9})},
658 {false, -142, MType({0x5951eef483de2c37, 0xad2039c52f37ee6e})},
659 {false, -142, MType({0xf7abe6ff6da76f1e, 0xb002e826665829cd})},
660 {false, -142, MType({0xf354411ed47c5d7b, 0xb2e5969eb317cc74})},
661 {false, -142, MType({0x1428a99ba8f5911f, 0xb5c8452e157847c0})},
662 {false, -142, MType({0x44a7c4330edff2c8, 0xb8aaf3d48d7b0d0c})},
663 {false, -142, MType({0x91f1306a84e4e07b, 0xbb8da2921b218db6})},
664 {false, -142, MType({0x2bc58de40cdf7b6a, 0xbe705166be6d3b1c})},
665 {false, -142, MType({0x648680b254df1d99, 0xc1530052775f869a})},
666 {false, -142, MType({0xb136b5ace0d6f74d, 0xc435af5545f9e18e})},
667 {false, -142, MType({0xa979e6c434fad480, 0xc7185e6f2a3dbd56})},
668 {false, -142, MType({0x794df5600c90a5a, 0xc9fb0da0242c8b50})},
669 {false, -142, MType({0xa86d80814ac18cf1, 0xccddbce833c7bcd8})},
670 {false, -142, MType({0x8b8ac57a9cca2d56, 0xcfc06c475910c34e})},
671 {false, -142, MType({0xd314c7e03140001f, 0xd2a31bbd9409100f})},
672 {false, -142, MType({0xc3d4c40e20b5ec89, 0xd585cb4ae4b2147a})},
673 {false, -142, MType({0xc5351d729060644e, 0xd8687aef4b0d41ed})},
674 {false, -142, MType({0x614162e1e12e445d, 0xdb4b2aaac71c09c7})},
675 {false, -142, MType({0x44a652eadf8ede85, 0xde2dda7d58dfdd66})},
676 {false, -142, MType({0x3eb1e02af3e52c3c, 0xe1108a67005a2e29})},
677 {false, -142, MType({0x415335a253a82aa2, 0xe3f33a67bd8c6d6f})},
678 {false, -142, MType({0x611abb0833305fe1, 0xe6d5ea7f90780c97})},
679 },
680 // -log2(r) for the fourth step, generated by SageMath with:
681 //
682 // for i in range(-65, 65):
683 // r = 2^-28 * round( 2^28 / (1 + i*2^(-28)) );
684 // s, m, e = RealField(128)(r).log2().sign_mantissa_exponent();
685 // print("{true," if (s == 1) else "{false,", e, ",
686 // MType({", hex(m % 2^64), ",", hex((m >> 64) % 2^64), "})},");
687 /* .step_4 = */
688 {
689 {true, -149, MType({0xef1bffe565ce0a46, 0xbb8ce2990b5d0b90})},
690 {true, -149, MType({0xbea3244560ca3d99, 0xb8aa39b807a576e4})},
691 {true, -149, MType({0x8b91f71ceefa31a2, 0xb5c790d6d5c354df})},
692 {true, -149, MType({0x9096e3d684001c0e, 0xb2e4e7f575b6a57b})},
693 {true, -149, MType({0x86054c794367f36, 0xb0023f13e77f68b3})},
694 {true, -149, MType({0x2d9cb33094afe4de, 0xad1f96322b1d9e80})},
695 {true, -149, MType({0x3afa673cfb3698f3, 0xaa3ced50409146dd})},
696 {true, -149, MType({0x6b27d8033e4c6450, 0xa75a446e27da61c4})},
697 {true, -149, MType({0xf8d36b84d52a477b, 0xa4779b8be0f8ef2f})},
698 {true, -149, MType({0x1eab86ae37c03565, 0xa194f2a96becef1a})},
699 {true, -149, MType({0x175e8d56deb4ce2c, 0x9eb249c6c8b6617d})},
700 {true, -149, MType({0x1d9ae241436519da, 0x9bcfa0e3f7554653})},
701 {true, -149, MType({0x6c0ee71adfe44325, 0x98ecf800f7c99d96})},
702 {true, -149, MType({0x3d68fc7c2efb522f, 0x960a4f1dca136741})},
703 {true, -149, MType({0xcc5781e8ac28e749, 0x9327a63a6e32a34d})},
704 {true, -149, MType({0x5388d5ced3a0f5af, 0x9044fd56e42751b6})},
705 {true, -149, MType({0xdab5588224c7e4a, 0x8d6254732bf17275})},
706 {true, -149, MType({0x356d5d5915c94a70, 0x8a7fab8f45910584})},
707 {true, -149, MType({0x57d48712c69a6a7, 0x879d02ab31060ade})},
708 {true, -149, MType({0xb88970eae5341d60, 0x84ba59c6ee50827c})},
709 {true, -149, MType({0x89402fcbbfe331bb, 0x81d7b0e27d706c5a})},
710 {true, -150, MType({0x649fba0879ca348b, 0xfdea0ffbbccb90e3})},
711 {true, -150, MType({0xdccd9edfbab6f777, 0xf824be3222612d78})},
712 {true, -150, MType({0xf066b9aa4636478e, 0xf25f6c682ba1ae69})},
713 {true, -150, MType({0x14c7b3cb21578781, 0xec9a1a9dd88d13ab})},
714 {true, -150, MType({0xbf4d347b528f56e1, 0xe6d4c8d329235d30})},
715 {true, -150, MType({0x6553e0c9e1b70799, 0xe10f77081d648aef})},
716 {true, -150, MType({0x7c385b9bd80c1375, 0xdb4a253cb5509cdb})},
717 {true, -150, MType({0x795745ac402f919d, 0xd584d370f0e792e9})},
718 {true, -150, MType({0xd20d3d8c2625ac1b, 0xcfbf81a4d0296d0d})},
719 {true, -150, MType({0xfbb6dfa297551554, 0xc9fa2fd853162b3c})},
720 {true, -150, MType({0x6bb0c62ca2867d91, 0xc434de0b79adcd6b})},
721 {true, -150, MType({0x9757893d57e40877, 0xbe6f8c3e43f0538d})},
722 {true, -150, MType({0xf407bebdc8f8c28e, 0xb8aa3a70b1ddbd97})},
723 {true, -150, MType({0xf71dfa6d08b016be, 0xb2e4e8a2c3760b7e})},
724 {true, -150, MType({0x15f6cde02b5543ce, 0xad1f96d478b93d37})},
725 {true, -150, MType({0xc5eec8824692d1e9, 0xa75a4505d1a752b4})},
726 {true, -150, MType({0x7c6277947172081a, 0xa194f336ce404bec})},
727 {true, -150, MType({0xaeae662dc45a61ce, 0x9bcfa1676e8428d2})},
728 {true, -150, MType({0xd22f1d3b59110455, 0x960a4f97b272e95b})},
729 {true, -150, MType({0x5c4123804ab83462, 0x9044fdc79a0c8d7c})},
730 {true, -150, MType({0xc240fd95b5cecb89, 0x8a7fabf725511528})},
731 {true, -150, MType({0x798b2deab82fadc4, 0x84ba5a2654408055})},
732 {true, -151, MType({0xeef86988e2227ddb, 0xfdea10aa4db59ded})},
733 {true, -151, MType({0x62e1207c0209b090, 0xf25f6d073a400203})},
734 {true, -151, MType({0x3989789113ec7bee, 0xe6d4c9636e202cd4})},
735 {true, -151, MType({0x5daa65565e562909, 0xdb4a25bee9561e49})},
736 {true, -151, MType({0xb9fcd6062a84acbd, 0xcfbf8219abe1d64b})},
737 {true, -151, MType({0x3939b586c46792b3, 0xc434de73b5c354c4})},
738 {true, -151, MType({0xc619ea6a7a9ee85e, 0xb8aa3acd06fa999b})},
739 {true, -151, MType({0x4b5656ef9e7a27fd, 0xad1f97259f87a4bb})},
740 {true, -151, MType({0xb3a7d90083f7239c, 0xa194f37d7f6a760b})},
741 {true, -151, MType({0xe9c74a3381c0f016, 0x960a4fd4a6a30d75})},
742 {true, -151, MType({0xd86d7fcaf12ed012, 0x8a7fac2b15316ae2})},
743 {true, -152, MType({0xd4a6956a5c863e0f, 0xfdea1101962b1c76})},
744 {true, -152, MType({0x1462ef192f547877, 0xe6d4c9ab909eeed1})},
745 {true, -152, MType({0x45819d2f1d72eb8b, 0xcfbf825419be4ca6})},
746 {true, -152, MType({0x3d742790eedbe719, 0xb8aa3afb318935c8})},
747 {true, -152, MType({0xd1ac0d7b70d74492, 0xa194f3a0d7ffaa08})},
748 {true, -152, MType({0xd79ac58375f83d0c, 0x8a7fac450d21a939})},
749 {true, -153, MType({0x49637b2bac367e87, 0xe6d4c9cfa1de665a})},
750 {true, -153, MType({0x1cc4b5eedcc78b35, 0xb8aa3b1246d08f69})},
751 {true, -153, MType({0xd43bf48a42745836, 0x8a7fac520919cd43})},
752 {true, -154, MType({0x3557bdcf592619eb, 0xb8aa3b1dd1743f1c})},
753 {true, -155, MType({0x6bdc2e83d3ebb0c4, 0xb8aa3b2396c617ae})},
754 {false, 0, MType({0x0, 0x0})},
755 {false, -155, MType({0x2d5b40050e44e8ab, 0xb8aa3b2f2169ca44})},
756 {false, -154, MType({0xb8560371b8f04afe, 0xb8aa3b34e6bba447})},
757 {false, -153, MType({0xc79a43ccc70459cc, 0x8a7fac6c010a1f14})},
758 {false, -153, MType({0x22c25632f519f77f, 0xb8aa3b40715f59c0})},
759 {false, -153, MType({0x42c10a314e35fb9e, 0xe6d4ca17c45d8282})},
760 {false, -152, MType({0xbe5a212ed7b949e4, 0x8a7fac78fd024cdb})},
761 {false, -152, MType({0x12dcf94ef5c5b918, 0xa194f3e7892a4fde})},
762 {false, -152, MType({0x49781013e57110ce, 0xb8aa3b5786a6ca76})},
763 {false, -152, MType({0x8cba70c085c12cb3, 0xcfbf82c8f577bcd2})},
764 {false, -152, MType({0x7332f3fb09328b8, 0xe6d4ca3bd59d2721})},
765 {false, -152, MType({0xe37168243a9d8b14, 0xfdea11b02717098f})},
766 {false, -151, MType({0xa602205479b93722, 0x8a7fac92f4f2b226})},
767 {false, -151, MType({0xb5bd735852c0d583, 0x960a504e8f041bc3})},
768 {false, -151, MType({0x363248630b0d812d, 0xa194f40ae1bfc1b6})},
769 {false, -151, MType({0x3ca83f0e02b823c0, 0xad1f97c7ed25a415})},
770 {false, -151, MType({0xde66fb46974bc4fd, 0xb8aa3b85b135c2f7})},
771 {false, -151, MType({0x30b6254e23c69fc2, 0xc434df442df01e75})},
772 {false, -151, MType({0x48dd69ba009b370c, 0xcfbf83036354b6a4})},
773 {false, -151, MType({0x3c24797383b16af5, 0xdb4a26c351638b9c})},
774 {false, -151, MType({0x1fd309b800678db7, 0xe6d4ca83f81c9d74})},
775 {false, -151, MType({0x930d418c79378a3, 0xf25f6e45577fec43})},
776 {false, -151, MType({0xd85967b2783a12c, 0xfdea12076f8d7820})},
777 {false, -150, MType({0x210c898c360016ed, 0x84ba5ae52022a091})},
778 {false, -150, MType({0x5e19883eef2605ab, 0x8a7facc6e4d3a3b0})},
779 {false, -150, MType({0x488dacc6629300ae, 0x9044fea905d9c579})},
780 {false, -150, MType({0x6b0cdebd3264e3e3, 0x960a508b833505f7})},
781 {false, -150, MType({0x503b07e7ff788dc2, 0x9bcfa26e5ce56536})},
782 {false, -150, MType({0x82bc1435696a69d1, 0xa194f45192eae341})},
783 {false, -150, MType({0x8d33f1be0e96fb1f, 0xa75a463525458024})},
784 {false, -150, MType({0xfa4690c48c1b66c9, 0xad1f981913f53bea})},
785 {false, -150, MType({0x5497e3b57dd5fe75, 0xb2e4e9fd5efa16a0})},
786 {false, -150, MType({0x26cbdf277e66cad5, 0xb8aa3be206541050})},
787 {false, -150, MType({0xfb8679db27301625, 0xbe6f8dc70a032905})},
788 {false, -150, MType({0x5d6bacbb1056f6aa, 0xc434dfac6a0760cd})},
789 {false, -150, MType({0xd71f72dbd0c3d936, 0xc9fa31922660b7b1})},
790 {false, -150, MType({0xf345c97bfe230ba2, 0xcfbf83783f0f2dbe})},
791 {false, -150, MType({0x3c82b0042ce54751, 0xd584d55eb412c300})},
792 {false, -150, MType({0x3d7a2806f0403bae, 0xdb4a2745856b7781})},
793 {false, -150, MType({0x80d03540da2f18ae, 0xe10f792cb3194b4d})},
794 {false, -150, MType({0x9128dd987b73194f, 0xe6d4cb143d1c3e70})},
795 {false, -150, MType({0xf928291e63940e14, 0xec9a1cfc237450f5})},
796 {false, -150, MType({0x4372220d20e0e78a, 0xf25f6ee4662182e9})},
797 {false, -150, MType({0xfaaad4c9407040c7, 0xf824c0cd0523d455})},
798 {false, -150, MType({0xa9764fe14e20e9e4, 0xfdea12b6007b4547})},
799 {false, -149, MType({0xed3c5206ea4d3942, 0x81d7b24fac13eae4})},
800 {false, -149, MType({0xc2af218aea6da27, 0x84ba5b448614c2f4})},
801 {false, -149, MType({0xf6d912ac383aaeba, 0x879d04398e402ad6})},
802 {false, -149, MType({0x7298bf5cca8b3d95, 0x8a7fad2ec4962293})},
803 {false, -149, MType({0x44bc04daa8808214, 0x8d6256242916aa2f})},
804 {false, -149, MType({0x3294f0eb14683198, 0x9044ff19bbc1c1b0})},
805 {false, -149, MType({0x17592684ff600c3, 0x9327a80f7c97691c})},
806 {false, -149, MType({0x76aff9419c43e8b9, 0x960a51056b97a078})},
807 {false, -149, MType({0x5796367b39d26c63, 0x98ecf9fb88c267cb})},
808 {false, -149, MType({0x697a5c2e6888ddaa, 0x9bcfa2f1d417bf1a})},
809 {false, -149, MType({0x71ae7d8967b5a2b7, 0x9eb24be84d97a66b})},
810 {false, -149, MType({0x3584aecf760e7b39, 0xa194f4def5421dc4})},
811 {false, -149, MType({0x7a4f0558d1b0c59e, 0xa4779dd5cb17252a})},
812 {false, -149, MType({0x55f9792b821c455, 0xa75a46cccf16bca4})},
813 {false, -149, MType({0x9c087cff664ee311, 0xaa3cefc40140e436})},
814 {false, -149, MType({0x39bce36188dfc04, 0xad1f98bb61959be8})},
815 {false, -149, MType({0x16ba4e30a9d9d21, 0xb00241b2f014e3be})},
816 {false, -149, MType({0x5aca1bc777a54d5e, 0xb2e4eaaaacbebbbe})},
817 {false, -149, MType({0xd5094eb99a35d1f0, 0xb5c793a2979323ee})},
818 {false, -149, MType({0x357b5aa4ac49738d, 0xb8aa3c9ab0921c55})},
819 }};
820
821 // > P = fpminimax(log2(1 + x)/x, 3, [|128...|], [-0x1.0002143p-29 , 0x1p-29]);
822 // > P;
823 // > dirtyinfnorm(log2(1 + x)/x - P, [-0x1.0002143p-29 , 0x1p-29]);
824 // 0x1.27ad5...p-121
825 const Float128 BIG_COEFFS[4]{
826 {true, -129, MType({0x3eccf6940d66bbcc, 0xb8aa3b295c2b21e3})},
827 {false, -129, MType({0xee39a6d649394bb1, 0xf6384ee1d01febc9})},
828 {true, -128, MType({0xbe87fed067ea2ad5, 0xb8aa3b295c17f0bb})},
829 {false, -127, MType({0xbe87fed0691d3e3f, 0xb8aa3b295c17f0bb})},
830 };
831
832 // Reuse the output of the fast pass range reduction.
833 // -2^-8 <= m_x < 2^-7
834 double log2_accurate(int e_x, int index, double m_x) {
835
836 Float128 sum(static_cast<float>(e_x));
837 sum = fputil::quick_add(sum, LOG2_TABLE.step_1[index]);
838
839 Float128 v_f128 = log_range_reduction(m_x, LOG2_TABLE, sum);
840
841 // Polynomial approximation
842 Float128 p = fputil::quick_mul(v_f128, BIG_COEFFS[0]);
843 p = fputil::quick_mul(v_f128, fputil::quick_add(p, BIG_COEFFS[1]));
844 p = fputil::quick_mul(v_f128, fputil::quick_add(p, BIG_COEFFS[2]));
845 p = fputil::quick_mul(v_f128, fputil::quick_add(p, BIG_COEFFS[3]));
846
847 Float128 r = fputil::quick_add(sum, p);
848
849 return static_cast<double>(r);
850 }
851
852 } // namespace
853
854 LLVM_LIBC_FUNCTION(double, log2, (double x)) {
855 using FPBits_t = typename fputil::FPBits<double>;
856 FPBits_t xbits(x);
857 uint64_t x_u = xbits.uintval();
858
859 int x_e = -FPBits_t::EXPONENT_BIAS;
860
861 if (LIBC_UNLIKELY(x_u == 0x3FF0'0000'0000'0000ULL)) {
862 // log2(1.0) = +0.0
863 return 0.0;
864 }
865
866 if (LIBC_UNLIKELY(xbits.uintval() < FPBits_t::MIN_NORMAL ||
867 xbits.uintval() > FPBits_t::MAX_NORMAL)) {
868 if (xbits.is_zero()) {
869 // return -Inf and raise FE_DIVBYZERO.
870 fputil::set_errno_if_required(ERANGE);
871 fputil::raise_except_if_required(FE_DIVBYZERO);
872 return static_cast<double>(FPBits_t::neg_inf());
873 }
874 if (xbits.get_sign() && !xbits.is_nan()) {
875 fputil::set_errno_if_required(EDOM);
876 fputil::raise_except_if_required(FE_INVALID);
877 return FPBits_t::build_quiet_nan(0);
878 }
879 if (xbits.is_inf_or_nan()) {
880 return x;
881 }
882 // Normalize denormal inputs.
883 xbits.set_val(x * 0x1.0p52);
884 x_e -= 52;
885 x_u = xbits.uintval();
886 }
887
888 // log2(x) = log2(2^x_e * x_m)
889 // = x_e + log2(x_m)
890 // Range reduction for log2(x_m):
891 // For each x_m, we would like to find r such that:
892 // -2^-8 <= r * x_m - 1 < 2^-7
893 int shifted = x_u >> 45;
894 int index = shifted & 0x7F;
895 double r = RD[index];
896
897 // Add unbiased exponent. Add an extra 1 if the 8 leading fractional bits are
898 // all 1's.
899 x_e += static_cast<int>((x_u + (1ULL << 45)) >> 52);
900 double e_x = static_cast<double>(x_e);
901
902 // Set m = 1.mantissa.
903 uint64_t x_m = (x_u & 0x000F'FFFF'FFFF'FFFFULL) | 0x3FF0'0000'0000'0000ULL;
904 double m = FPBits_t(x_m).get_val();
905
906 double u, u_sq, err;
907 fputil::DoubleDouble r1;
908
909 // Perform exact range reduction
910 #ifdef LIBC_TARGET_CPU_HAS_FMA
911 u = fputil::multiply_add(r, m, -1.0); // exact
912 #else
913 uint64_t c_m = x_m & 0x3FFF'E000'0000'0000ULL;
914 double c = FPBits_t(c_m).get_val();
915 u = fputil::multiply_add(r, m - c, CD[index]); // exact
916 #endif // LIBC_TARGET_CPU_HAS_FMA
917
918 // Exact sum:
919 // r1.hi + r1.lo = e_x * log(2)_hi - log(r)_hi + u
920 r1 = fputil::exact_add(LOG_R1[index].hi, u);
921
922 // Error of u_sq = ulp(u^2);
923 u_sq = u * u;
924 // Total error is bounded by ~ C * ulp(u^2).
925 err = u_sq * P_ERR;
926 // Degree-7 minimax polynomial
927 double p0 = fputil::multiply_add(u, LOG_COEFFS[1], LOG_COEFFS[0]);
928 double p1 = fputil::multiply_add(u, LOG_COEFFS[3], LOG_COEFFS[2]);
929 double p2 = fputil::multiply_add(u, LOG_COEFFS[5], LOG_COEFFS[4]);
930 double p = fputil::polyeval(u_sq, LOG_R1[index].lo, p0, p1, p2);
931
932 r1.lo += p;
933
934 // Quick double-double multiplication:
935 // r2.hi + r2.lo ~ r1 * log2(e),
936 // with error bounded by:
937 // 4*ulp( ulp(r2.hi) )
938 fputil::DoubleDouble r2 = fputil::quick_mult(r1, LOG2_E);
939 fputil::DoubleDouble r3 = fputil::exact_add(e_x, r2.hi);
940 r3.lo += r2.lo;
941
942 // Overall, if we choose sufficiently large constant C, the total error is
943 // bounded by (C * ulp(u^2)).
944
945 // Lower bound from the result
946 double left = r3.hi + (r3.lo - err);
947 // Upper bound from the result
948 double right = r3.hi + (r3.lo + err);
949
950 // Ziv's test if fast pass is accurate enough.
951 if (left == right)
952 return left;
953
954 return log2_accurate(x_e, index, u);
955 }
956
957 } // namespace LIBC_NAMESPACE
LLVM monorepo