1/* _Float128 overrides for building ldbl-128 as _Float128.
2 Copyright (C) 2017-2018 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 <http://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 <ieee754_float128.h>
33#define ieee854_long_double_shape_type ieee854_float128_shape_type
34#define ieee854_long_double ieee854_float128
35
36#undef GET_LDOUBLE_LSW64
37#undef GET_LDOUBLE_MSW64
38#undef GET_LDOUBLE_WORDS64
39#undef SET_LDOUBLE_LSW64
40#undef SET_LDOUBLE_MSW64
41#undef SET_LDOUBLE_WORDS64
42#define GET_LDOUBLE_LSW64(x,y) GET_FLOAT128_LSW64 (x, y)
43#define GET_LDOUBLE_MSW64(x,y) GET_FLOAT128_MSW64 (x, y)
44#define GET_LDOUBLE_WORDS64(x,y,z) GET_FLOAT128_WORDS64 (x, y, z)
45#define SET_LDOUBLE_LSW64(x,y) SET_FLOAT128_LSW64 (x, y)
46#define SET_LDOUBLE_MSW64(x,y) SET_FLOAT128_MSW64 (x, y)
47#define SET_LDOUBLE_WORDS64(x,y,z) SET_FLOAT128_WORDS64 (x, y, z)
48
49#undef IEEE854_LONG_DOUBLE_BIAS
50#define IEEE854_LONG_DOUBLE_BIAS IEEE854_FLOAT128_BIAS
51
52#ifdef SET_RESTORE_ROUNDF128
53# undef SET_RESTORE_ROUNDL
54# define SET_RESTORE_ROUNDL(RM) SET_RESTORE_ROUNDF128 (RM)
55#endif
56
57
58/* misc macros from the header below. */
59#include <fix-fp-int-convert-overflow.h>
60#undef FIX_LDBL_LONG_CONVERT_OVERFLOW
61#undef FIX_LDBL_LLONG_CONVERT_OVERFLOW
62#define FIX_LDBL_LONG_CONVERT_OVERFLOW FIX_FLT128_LONG_CONVERT_OVERFLOW
63#define FIX_LDBL_LLONG_CONVERT_OVERFLOW FIX_FLT128_LLONG_CONVERT_OVERFLOW
64
65
66/* float.h constants. */
67#include <float.h>
68#undef LDBL_DIG
69#undef LDBL_EPSILON
70#undef LDBL_MANT_DIG
71#undef LDBL_MAX
72#undef LDBL_MAX_10_EXP
73#undef LDBL_MAX_EXP
74#undef LDBL_MIN
75#undef LDBL_MIN_10_EXP
76#undef LDBL_MIN_EXP
77#undef LDBL_TRUE_MIN
78#define LDBL_DIG FLT128_DIG
79#define LDBL_EPSILON FLT128_EPSILON
80#define LDBL_MANT_DIG FLT128_MANT_DIG
81#define LDBL_MAX FLT128_MAX
82#define LDBL_MAX_10_EXP FLT128_MAX_10_EXP
83#define LDBL_MAX_EXP FLT128_MAX_EXP
84#define LDBL_MIN FLT128_MIN
85#define LDBL_MIN_10_EXP FLT128_MIN_10_EXP
86#define LDBL_MIN_EXP FLT128_MIN_EXP
87#define LDBL_TRUE_MIN FLT128_TRUE_MIN
88
89
90/* math.h GNU constants. */
91#undef M_El
92#undef M_LOG2El
93#undef M_LOG10El
94#undef M_LN2l
95#undef M_LN10l
96#undef M_PIl
97#undef M_PI_2l
98#undef M_PI_4l
99#undef M_1_PIl
100#undef M_2_PIl
101#undef M_2_SQRTPIl
102#undef M_SQRT2l
103#undef M_SQRT1_2l
104#define M_El M_Ef128
105#define M_LOG2El M_LOG2Ef128
106#define M_LOG10El M_LOG10Ef128
107#define M_LN2l M_LN2f128
108#define M_LN10l M_LN10f128
109#define M_PIl M_PIf128
110#define M_PI_2l M_PI_2f128
111#define M_PI_4l M_PI_4f128
112#define M_1_PIl M_1_PIf128
113#define M_2_PIl M_2_PIf128
114#define M_2_SQRTPIl M_2_SQRTPIf128
115#define M_SQRT2l M_SQRT2f128
116#define M_SQRT1_2l M_SQRT1_2f128
117
118
119#include <libm-alias-ldouble.h>
120#include <libm-alias-float128.h>
121#undef libm_alias_ldouble_r
122#define libm_alias_ldouble_r(from, to, r) libm_alias_float128_r (from, to, r)
123
124
125/* IEEE function renames. */
126#define __ieee754_acoshl __ieee754_acoshf128
127#define __ieee754_acosl __ieee754_acosf128
128#define __ieee754_asinhl __ieee754_asinhf128
129#define __ieee754_asinl __ieee754_asinf128
130#define __ieee754_atan2l __ieee754_atan2f128
131#define __ieee754_atanhl __ieee754_atanhf128
132#define __ieee754_coshl __ieee754_coshf128
133#define __ieee754_cosl __ieee754_cosf128
134#define __ieee754_exp10l __ieee754_exp10f128
135#define __ieee754_exp2l __ieee754_exp2f128
136#define __ieee754_expl __ieee754_expf128
137#define __ieee754_fmodl __ieee754_fmodf128
138#define __ieee754_gammal_r __ieee754_gammaf128_r
139#define __ieee754_hypotl __ieee754_hypotf128
140#define __ieee754_ilogbl __ieee754_ilogbf128
141#define __ieee754_j0l __ieee754_j0f128
142#define __ieee754_j1l __ieee754_j1f128
143#define __ieee754_jnl __ieee754_jnf128
144#define __ieee754_lgammal_r __ieee754_lgammaf128_r
145#define __ieee754_log10l __ieee754_log10f128
146#define __ieee754_log2l __ieee754_log2f128
147#define __ieee754_logl __ieee754_logf128
148#define __ieee754_powl __ieee754_powf128
149#define __ieee754_rem_pio2l __ieee754_rem_pio2f128
150#define __ieee754_remainderl __ieee754_remainderf128
151#define __ieee754_sinhl __ieee754_sinhf128
152#define __ieee754_sqrtl __ieee754_sqrtf128
153#define __ieee754_y0l __ieee754_y0f128
154#define __ieee754_y1l __ieee754_y1f128
155#define __ieee754_ynl __ieee754_ynf128
156
157
158/* finite math entry points. */
159#define __acoshl_finite __acoshf128_finite
160#define __acosl_finite __acosf128_finite
161#define __asinl_finite __asinf128_finite
162#define __atan2l_finite __atan2f128_finite
163#define __atanhl_finite __atanhf128_finite
164#define __coshl_finite __coshf128_finite
165#define __cosl_finite __cosf128_finite
166#define __exp10l_finite __exp10f128_finite
167#define __exp2l_finite __exp2f128_finite
168#define __expl_finite __expf128_finite
169#define __fmodl_finite __fmodf128_finite
170#define __hypotl_finite __hypotf128_finite
171#define __ilogbl_finite __ilogbf128_finite
172#define __j0l_finite __j0f128_finite
173#define __j1l_finite __j1f128_finite
174#define __jnl_finite __jnf128_finite
175#define __lgammal_r_finite __lgammaf128_r_finite
176#define __log10l_finite __log10f128_finite
177#define __log2l_finite __log2f128_finite
178#define __logl_finite __logf128_finite
179#define __powl_finite __powf128_finite
180#define __remainderl_finite __remainderf128_finite
181#define __sinhl_finite __sinhf128_finite
182#define __y0l_finite __y0f128_finite
183#define __y1l_finite __y1f128_finite
184#define __ynl_finite __ynf128_finite
185
186
187/* internal function names. */
188#define __asinhl __asinhf128
189#define __atanl __atanf128
190#define __cbrtl __cbrtf128
191#define __ceill __ceilf128
192#define __copysignl __copysignf128
193#define __cosl __cosf128
194#define __erfcl __erfcf128
195#define __erfl __erff128
196#define __expl __expf128
197#define __expm1l __expm1f128
198#define __fabsl __fabsf128
199#define __fdiml __fdimf128
200#define __finitel __finitef128
201#define __floorl __floorf128
202#define __fmal __fmaf128
203#define __fmaxl __fmaxf128
204#define __fminl __fminf128
205#define __fpclassifyl __fpclassifyf128
206#define __frexpl __frexpf128
207#define __gammal_r_finite __gammaf128_r_finite
208#define __getpayloadl __getpayloadf128
209#define __isinfl __isinff128
210#define __isnanl __isnanf128
211#define __issignalingl __issignalingf128
212#define __ldexpl __ldexpf128
213#define __llrintl __llrintf128
214#define __llroundl __llroundf128
215#define __log1pl __log1pf128
216#define __logbl __logbf128
217#define __logl __logf128
218#define __lrintl __lrintf128
219#define __lroundl __lroundf128
220#define __modfl __modff128
221#define __nearbyintl __nearbyintf128
222#define __nextafterl __nextafterf128
223#define __nextdownl __nextdownf128
224#define __nextupl __nextupf128
225#define __remquol __remquof128
226#define __rintl __rintf128
227#define __roundevenl __roundevenf128
228#define __roundl __roundf128
229#define __scalblnl __scalblnf128
230#define __scalbnl __scalbnf128
231#define __signbitl __signbitf128
232#define __sincosl __sincosf128
233#define __sinl __sinf128
234#define __sqrtl __sqrtf128
235#define __tanhl __tanhf128
236#define __tanl __tanf128
237#define __totalorderl __totalorderf128
238#define __totalordermagl __totalordermagf128
239#define __truncl __truncf128
240#define __x2y2m1l __x2y2m1f128
241
242/* __nexttowardf128 is not _Float128 API. */
243#define __nexttowardl __nexttowardf128_do_not_use
244#define nexttowardl nexttowardf128_do_not_use
245
246
247/* public entry points. */
248#define asinhl asinhf128
249#define atanl atanf128
250#define cbrtl cbrtf128
251#define ceill ceilf128
252#define copysignl copysignf128
253#define cosl cosf128
254#define erfcl erfcf128
255#define erfl erff128
256#define expl expf128
257#define expm1l expm1f128
258#define fabsl fabsf128
259#define fdiml fdimf128
260#define finitel finitef128_do_not_use
261#define floorl floorf128
262#define fmal fmaf128
263#define fmaxl fmaxf128
264#define fminl fminf128
265#define frexpl frexpf128
266#define getpayloadl getpayloadf128
267#define isinfl isinff128_do_not_use
268#define isnanl isnanf128_do_not_use
269#define ldexpl ldexpf128
270#define llrintl llrintf128
271#define llroundl llroundf128
272#define log1pl log1pf128
273#define logbl logbf128
274#define logl logf128
275#define lrintl lrintf128
276#define lroundl lroundf128
277#define modfl modff128
278#define nanl nanf128
279#define nearbyintl nearbyintf128
280#define nextafterl nextafterf128
281#define nextdownl nextdownf128
282#define nextupl nextupf128
283#define remquol remquof128
284#define rintl rintf128
285#define roundevenl roundevenf128
286#define roundl roundf128
287#define scalbnl scalbnf128
288#define sincosl sincosf128
289#define sinl sinf128
290#define sqrtl sqrtf128
291#define tanhl tanhf128
292#define tanl tanf128
293#define totalorderl totalorderf128
294#define totalordermagl totalordermagf128
295#define truncl truncf128
296
297
298/* misc internal renames. */
299#define __builtin_fmal __builtin_fmaf128
300#define __expl_table __expf128_table
301#define __gamma_productl __gamma_productf128
302#define __kernel_cosl __kernel_cosf128
303#define __kernel_rem_pio2l __kernel_rem_pio2f128
304#define __kernel_sincosl __kernel_sincosf128
305#define __kernel_sinl __kernel_sinf128
306#define __kernel_tanl __kernel_tanf128
307#define __lgamma_negl __lgamma_negf128
308#define __lgamma_productl __lgamma_productf128
309#define __mpn_extract_long_double __mpn_extract_float128
310#define __sincosl_table __sincosf128_table
311#define mul_splitl mul_splitf128
312
313/* Builtin renames. */
314#define __builtin_copysignl __builtin_copysignf128
315#define __builtin_signbitl __builtin_signbit
316
317/* Get the constant suffix from bits/floatn-compat.h. */
318#define L(x) __f128 (x)
319
320static inline void
321mul_splitf128 (_Float128 *hi, _Float128 *lo, _Float128 x, _Float128 y)
322{
323#ifdef __FP_FAST_FMAF128
324 /* Fast built-in fused multiply-add. */
325 *hi = x * y;
326 *lo = __builtin_fmal (x, y, -*hi);
327#else
328 /* Apply Dekker's algorithm. */
329 *hi = x * y;
330# define C ((1LL << (FLT128_MANT_DIG + 1) / 2) + 1)
331 _Float128 x1 = x * C;
332 _Float128 y1 = y * C;
333# undef C
334 x1 = (x - x1) + x1;
335 y1 = (y - y1) + y1;
336 _Float128 x2 = x - x1;
337 _Float128 y2 = y - y1;
338 *lo = (((x1 * y1 - *hi) + x1 * y2) + x2 * y1) + x2 * y2;
339#endif
340}
341