1/* _Float128 overrides for building ldbl-128 as _Float128.
2 Copyright (C) 2017-2021 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/* This must be included before the function renames below. */
20#include <gmp.h>
21#include <math.h>
22#undef HUGE_VALL
23#define HUGE_VALL HUGE_VAL_F128
24#include <math/mul_splitl.h>
25
26/* This must be included before the renames of types and macros from
27 it. */
28#include <ieee754.h>
29
30/* Renames derived from math_private.h. */
31#include <math_private.h>
32#include <fenv_private.h>
33#include <ieee754_float128.h>
34#define ieee854_long_double_shape_type ieee854_float128_shape_type
35#define ieee854_long_double ieee854_float128
36
37#undef GET_LDOUBLE_LSW64
38#undef GET_LDOUBLE_MSW64
39#undef GET_LDOUBLE_WORDS64
40#undef SET_LDOUBLE_LSW64
41#undef SET_LDOUBLE_MSW64
42#undef SET_LDOUBLE_WORDS64
43#define GET_LDOUBLE_LSW64(x,y) GET_FLOAT128_LSW64 (x, y)
44#define GET_LDOUBLE_MSW64(x,y) GET_FLOAT128_MSW64 (x, y)
45#define GET_LDOUBLE_WORDS64(x,y,z) GET_FLOAT128_WORDS64 (x, y, z)
46#define SET_LDOUBLE_LSW64(x,y) SET_FLOAT128_LSW64 (x, y)
47#define SET_LDOUBLE_MSW64(x,y) SET_FLOAT128_MSW64 (x, y)
48#define SET_LDOUBLE_WORDS64(x,y,z) SET_FLOAT128_WORDS64 (x, y, z)
49
50#undef IEEE854_LONG_DOUBLE_BIAS
51#define IEEE854_LONG_DOUBLE_BIAS IEEE854_FLOAT128_BIAS
52
53#ifdef SET_RESTORE_ROUNDF128
54# undef SET_RESTORE_ROUNDL
55# define SET_RESTORE_ROUNDL(RM) SET_RESTORE_ROUNDF128 (RM)
56#endif
57
58#ifdef libc_feholdexcept_setroundf128
59# undef libc_feholdexcept_setroundl
60# define libc_feholdexcept_setroundl(ENV, RM) \
61 libc_feholdexcept_setroundf128 (ENV, RM)
62#endif
63
64#ifdef libc_feupdateenv_testf128
65# undef libc_feupdateenv_testl
66# define libc_feupdateenv_testl(ENV, EX) libc_feupdateenv_testf128 (ENV, EX)
67#endif
68
69/* misc macros from the header below. */
70#include <fix-fp-int-convert-overflow.h>
71#undef FIX_LDBL_LONG_CONVERT_OVERFLOW
72#undef FIX_LDBL_LLONG_CONVERT_OVERFLOW
73#define FIX_LDBL_LONG_CONVERT_OVERFLOW FIX_FLT128_LONG_CONVERT_OVERFLOW
74#define FIX_LDBL_LLONG_CONVERT_OVERFLOW FIX_FLT128_LLONG_CONVERT_OVERFLOW
75
76
77/* float.h constants. */
78#include <float.h>
79#undef LDBL_DIG
80#undef LDBL_EPSILON
81#undef LDBL_MANT_DIG
82#undef LDBL_MAX
83#undef LDBL_MAX_10_EXP
84#undef LDBL_MAX_EXP
85#undef LDBL_MIN
86#undef LDBL_MIN_10_EXP
87#undef LDBL_MIN_EXP
88#undef LDBL_TRUE_MIN
89#define LDBL_DIG FLT128_DIG
90#define LDBL_EPSILON FLT128_EPSILON
91#define LDBL_MANT_DIG FLT128_MANT_DIG
92#define LDBL_MAX FLT128_MAX
93#define LDBL_MAX_10_EXP FLT128_MAX_10_EXP
94#define LDBL_MAX_EXP FLT128_MAX_EXP
95#define LDBL_MIN FLT128_MIN
96#define LDBL_MIN_10_EXP FLT128_MIN_10_EXP
97#define LDBL_MIN_EXP FLT128_MIN_EXP
98#define LDBL_TRUE_MIN FLT128_TRUE_MIN
99
100
101/* math.h GNU constants. */
102#undef M_El
103#undef M_LOG2El
104#undef M_LOG10El
105#undef M_LN2l
106#undef M_LN10l
107#undef M_PIl
108#undef M_PI_2l
109#undef M_PI_4l
110#undef M_1_PIl
111#undef M_2_PIl
112#undef M_2_SQRTPIl
113#undef M_SQRT2l
114#undef M_SQRT1_2l
115#define M_El M_Ef128
116#define M_LOG2El M_LOG2Ef128
117#define M_LOG10El M_LOG10Ef128
118#define M_LN2l M_LN2f128
119#define M_LN10l M_LN10f128
120#define M_PIl M_PIf128
121#define M_PI_2l M_PI_2f128
122#define M_PI_4l M_PI_4f128
123#define M_1_PIl M_1_PIf128
124#define M_2_PIl M_2_PIf128
125#define M_2_SQRTPIl M_2_SQRTPIf128
126#define M_SQRT2l M_SQRT2f128
127#define M_SQRT1_2l M_SQRT1_2f128
128
129
130#include <libm-alias-ldouble.h>
131#include <libm-alias-float128.h>
132#undef libm_alias_ldouble_r
133#define libm_alias_ldouble_r(from, to, r) libm_alias_float128_r (from, to, r)
134
135
136#include <math/math-narrow.h>
137#undef libm_alias_float_ldouble
138#define libm_alias_float_ldouble(func) libm_alias_float32_float128 (func)
139#undef libm_alias_double_ldouble
140#define libm_alias_double_ldouble(func) libm_alias_float64_float128 (func)
141
142#include <math-use-builtins.h>
143#undef USE_NEARBYINTL_BUILTIN
144#define USE_NEARBYINTL_BUILTIN USE_NEARBYINTF128_BUILTIN
145#undef USE_RINTL_BUILTIN
146#define USE_RINTL_BUILTIN USE_RINTF128_BUILTIN
147#undef USE_FLOORL_BUILTIN
148#define USE_FLOORL_BUILTIN USE_FLOORF128_BUILTIN
149#undef USE_CEILL_BUILTIN
150#define USE_CEILL_BUILTIN USE_CEILF128_BUILTIN
151#undef USE_TRUNCL_BUILTIN
152#define USE_TRUNCL_BUILTIN USE_TRUNCF128_BUILTIN
153#undef USE_ROUNDL_BUILTIN
154#define USE_ROUNDL_BUILTIN USE_ROUNDF128_BUILTIN
155#undef USE_ROUNDEVENL_BUILTIN
156#define USE_ROUNDEVENL_BUILTIN USE_ROUNDEVENF128_BUILTIN
157#undef USE_COPYSIGNL_BUILTIN
158#define USE_COPYSIGNL_BUILTIN USE_COPYSIGNF128_BUILTIN
159#undef USE_FMAL_BUILTIN
160#define USE_FMAL_BUILTIN USE_FMAF128_BUILTIN
161
162/* IEEE function renames. */
163#define __ieee754_acoshl __ieee754_acoshf128
164#define __ieee754_acosl __ieee754_acosf128
165#define __ieee754_asinhl __ieee754_asinhf128
166#define __ieee754_asinl __ieee754_asinf128
167#define __ieee754_atan2l __ieee754_atan2f128
168#define __ieee754_atanhl __ieee754_atanhf128
169#define __ieee754_coshl __ieee754_coshf128
170#define __ieee754_cosl __ieee754_cosf128
171#define __ieee754_exp10l __ieee754_exp10f128
172#define __ieee754_exp2l __ieee754_exp2f128
173#define __ieee754_expl __ieee754_expf128
174#define __ieee754_fmodl __ieee754_fmodf128
175#define __ieee754_gammal_r __ieee754_gammaf128_r
176#define __ieee754_hypotl __ieee754_hypotf128
177#define __ieee754_ilogbl __ieee754_ilogbf128
178#define __ieee754_j0l __ieee754_j0f128
179#define __ieee754_j1l __ieee754_j1f128
180#define __ieee754_jnl __ieee754_jnf128
181#define __ieee754_lgammal_r __ieee754_lgammaf128_r
182#define __ieee754_log10l __ieee754_log10f128
183#define __ieee754_log2l __ieee754_log2f128
184#define __ieee754_logl __ieee754_logf128
185#define __ieee754_powl __ieee754_powf128
186#define __ieee754_rem_pio2l __ieee754_rem_pio2f128
187#define __ieee754_remainderl __ieee754_remainderf128
188#define __ieee754_sinhl __ieee754_sinhf128
189#define __ieee754_sqrtl __ieee754_sqrtf128
190#define __ieee754_y0l __ieee754_y0f128
191#define __ieee754_y1l __ieee754_y1f128
192#define __ieee754_ynl __ieee754_ynf128
193
194
195/* finite math entry points. */
196#define __acoshl_finite __acoshf128_finite
197#define __acosl_finite __acosf128_finite
198#define __asinl_finite __asinf128_finite
199#define __atan2l_finite __atan2f128_finite
200#define __atanhl_finite __atanhf128_finite
201#define __coshl_finite __coshf128_finite
202#define __cosl_finite __cosf128_finite
203#define __exp10l_finite __exp10f128_finite
204#define __exp2l_finite __exp2f128_finite
205#define __expl_finite __expf128_finite
206#define __fmodl_finite __fmodf128_finite
207#define __hypotl_finite __hypotf128_finite
208#define __ilogbl_finite __ilogbf128_finite
209#define __j0l_finite __j0f128_finite
210#define __j1l_finite __j1f128_finite
211#define __jnl_finite __jnf128_finite
212#define __lgammal_r_finite __lgammaf128_r_finite
213#define __log10l_finite __log10f128_finite
214#define __log2l_finite __log2f128_finite
215#define __logl_finite __logf128_finite
216#define __powl_finite __powf128_finite
217#define __remainderl_finite __remainderf128_finite
218#define __sinhl_finite __sinhf128_finite
219#define __y0l_finite __y0f128_finite
220#define __y1l_finite __y1f128_finite
221#define __ynl_finite __ynf128_finite
222
223
224/* internal function names. */
225#define __asinhl __asinhf128
226#define __atanl __atanf128
227#define __cbrtl __cbrtf128
228#define __ceill __ceilf128
229#define __copysignl __copysignf128
230#define __cosl __cosf128
231#define __erfcl __erfcf128
232#define __erfl __erff128
233#define __expl __expf128
234#define __expm1l __expm1f128
235#define __fabsl __fabsf128
236#define __fdiml __fdimf128
237#define __finitel __finitef128
238#define __floorl __floorf128
239#define __fmal __fmaf128
240#define __fmaxl __fmaxf128
241#define __fminl __fminf128
242#define __fpclassifyl __fpclassifyf128
243#define __frexpl __frexpf128
244#define __gammal_r_finite __gammaf128_r_finite
245#define __getpayloadl __getpayloadf128
246#define __isinfl __isinff128
247#define __isnanl __isnanf128
248#define __issignalingl __issignalingf128
249#define __ldexpl __ldexpf128
250#define __llrintl __llrintf128
251#define __llroundl __llroundf128
252#define __log1pl __log1pf128
253#define __logbl __logbf128
254#define __logl __logf128
255#define __lrintl __lrintf128
256#define __lroundl __lroundf128
257#define __modfl __modff128
258#define __nearbyintl __nearbyintf128
259#define __nextafterl __nextafterf128
260#define __nextdownl __nextdownf128
261#define __nextupl __nextupf128
262#define __remquol __remquof128
263#define __rintl __rintf128
264#define __roundevenl __roundevenf128
265#define __roundl __roundf128
266#define __scalblnl __scalblnf128
267#define __scalbnl __scalbnf128
268#define __signbitl __signbitf128
269#define __sincosl __sincosf128
270#define __sinl __sinf128
271#define __sqrtl __sqrtf128
272#define __tanhl __tanhf128
273#define __tanl __tanf128
274#define __totalorderl __totalorderf128
275#define __totalorder_compatl __totalorder_compatf128
276#define __totalordermagl __totalordermagf128
277#define __totalordermag_compatl __totalordermag_compatf128
278#define __truncl __truncf128
279#define __x2y2m1l __x2y2m1f128
280
281#define __faddl __f32addf128
282#define __daddl __f64addf128
283#define __fdivl __f32divf128
284#define __ddivl __f64divf128
285#define __fmull __f32mulf128
286#define __dmull __f64mulf128
287#define __fsubl __f32subf128
288#define __dsubl __f64subf128
289
290/* Used on __finite compat alias. */
291#define __acosl __acosf128
292#define __acoshl __acoshf128
293#define __asinl __asinf128
294#define __atan2l __atan2f128
295#define __atanhl __atanhf128
296#define __coshl __coshf128
297#define __exp10l __exp10f128
298#define __expl __expf128
299#define __fmodl __fmodf128
300#define __gammal_r __gammaf128_r
301#define __hypotl __hypotf128
302#define __j0l __j0f128
303#define __j1l __j1f128
304#define __jnl __jnf128
305#define __lgammal_r __lgammaf128_r
306#define __log10l __log10f128
307#define __log2l __log2f128
308#define __logl __logf128
309#define __powl __powf128
310#define __remainderl __remainderf128
311#define __sinhl __sinhf128
312#define __y0l __y0f128
313#define __y1l __y1f128
314#define __ynl __ynf128
315
316/* __nexttowardf128 is not _Float128 API. */
317#define __nexttowardl __nexttowardf128_do_not_use
318#define nexttowardl nexttowardf128_do_not_use
319
320
321/* public entry points. */
322#define asinhl asinhf128
323#define atanl atanf128
324#define cbrtl cbrtf128
325#define ceill ceilf128
326#define copysignl copysignf128
327#define cosl cosf128
328#define erfcl erfcf128
329#define erfl erff128
330#define expl expf128
331#define expm1l expm1f128
332#define fabsl fabsf128
333#define fdiml fdimf128
334#define finitel finitef128_do_not_use
335#define floorl floorf128
336#define fmal fmaf128
337#define fmaxl fmaxf128
338#define fminl fminf128
339#define frexpl frexpf128
340#define getpayloadl getpayloadf128
341#define isinfl isinff128_do_not_use
342#define isnanl isnanf128_do_not_use
343#define ldexpl ldexpf128
344#define llrintl llrintf128
345#define llroundl llroundf128
346#define log1pl log1pf128
347#define logbl logbf128
348#define logl logf128
349#define lrintl lrintf128
350#define lroundl lroundf128
351#define modfl modff128
352#define nanl nanf128
353#define nearbyintl nearbyintf128
354#define nextafterl nextafterf128
355#define nextdownl nextdownf128
356#define nextupl nextupf128
357#define remquol remquof128
358#define rintl rintf128
359#define roundevenl roundevenf128
360#define roundl roundf128
361#define scalbnl scalbnf128
362#define sincosl sincosf128
363#define sinl sinf128
364#define sqrtl sqrtf128
365#define tanhl tanhf128
366#define tanl tanf128
367#define totalorderl totalorderf128
368#define totalordermagl totalordermagf128
369#define truncl truncf128
370
371
372/* misc internal renames. */
373#define __builtin_fmal __builtin_fmaf128
374#define __expl_table __expf128_table
375#define __gamma_productl __gamma_productf128
376#define __kernel_cosl __kernel_cosf128
377#define __kernel_rem_pio2l __kernel_rem_pio2f128
378#define __kernel_sincosl __kernel_sincosf128
379#define __kernel_sinl __kernel_sinf128
380#define __kernel_tanl __kernel_tanf128
381#define __lgamma_negl __lgamma_negf128
382#define __lgamma_productl __lgamma_productf128
383#define __mpn_extract_long_double __mpn_extract_float128
384#define __sincosl_table __sincosf128_table
385#define mul_splitl mul_splitf128
386
387/* Builtin renames. */
388#define __builtin_copysignl __builtin_copysignf128
389#define __builtin_signbitl __builtin_signbit
390#define __builtin_nearbyintl __builtin_nearbyintf128
391#define __builtin_rintl __builtin_rintf128
392#define __builtin_floorl __builtin_floorf128
393#define __builtin_ceill __builtin_ceilf128
394#define __builtin_truncl __builtin_truncf128
395#define __builtin_roundl __builtin_roundf128
396#define __builtin_copysignl __builtin_copysignf128
397
398/* Get the constant suffix from bits/floatn-compat.h. */
399#define L(x) __f128 (x)
400
401static inline void
402mul_splitf128 (_Float128 *hi, _Float128 *lo, _Float128 x, _Float128 y)
403{
404#ifdef __FP_FAST_FMAF128
405 /* Fast built-in fused multiply-add. */
406 *hi = x * y;
407 *lo = __builtin_fmal (x, y, -*hi);
408#else
409 /* Apply Dekker's algorithm. */
410 *hi = x * y;
411# define C ((1LL << (FLT128_MANT_DIG + 1) / 2) + 1)
412 _Float128 x1 = x * C;
413 _Float128 y1 = y * C;
414# undef C
415 x1 = (x - x1) + x1;
416 y1 = (y - y1) + y1;
417 _Float128 x2 = x - x1;
418 _Float128 y2 = y - y1;
419 *lo = (((x1 * y1 - *hi) + x1 * y2) + x2 * y1) + x2 * y2;
420#endif
421}
422