diff options
Diffstat (limited to 'sysdeps/ieee754/float128/float128_private.h')
| -rw-r--r-- | sysdeps/ieee754/float128/float128_private.h | 366 |
1 files changed, 366 insertions, 0 deletions
diff --git a/sysdeps/ieee754/float128/float128_private.h b/sysdeps/ieee754/float128/float128_private.h new file mode 100644 index 0000000000..9dd15601e6 --- /dev/null +++ b/sysdeps/ieee754/float128/float128_private.h @@ -0,0 +1,366 @@ +/* _Float128 overrides for building ldbl-128 as _Float128. + Copyright (C) 2017-2018 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + <http://www.gnu.org/licenses/>. */ + +/* This must be included before the function renames below. */ +#include <gmp.h> +#include <math.h> +#undef HUGE_VALL +#define HUGE_VALL HUGE_VAL_F128 +#include <math/mul_splitl.h> + +/* This must be included before the renames of types and macros from + it. */ +#include <ieee754.h> + +/* Renames derived from math_private.h. */ +#include <math_private.h> +#include <ieee754_float128.h> +#define ieee854_long_double_shape_type ieee854_float128_shape_type +#define ieee854_long_double ieee854_float128 + +#undef GET_LDOUBLE_LSW64 +#undef GET_LDOUBLE_MSW64 +#undef GET_LDOUBLE_WORDS64 +#undef SET_LDOUBLE_LSW64 +#undef SET_LDOUBLE_MSW64 +#undef SET_LDOUBLE_WORDS64 +#define GET_LDOUBLE_LSW64(x,y) GET_FLOAT128_LSW64 (x, y) +#define GET_LDOUBLE_MSW64(x,y) GET_FLOAT128_MSW64 (x, y) +#define GET_LDOUBLE_WORDS64(x,y,z) GET_FLOAT128_WORDS64 (x, y, z) +#define SET_LDOUBLE_LSW64(x,y) SET_FLOAT128_LSW64 (x, y) +#define SET_LDOUBLE_MSW64(x,y) SET_FLOAT128_MSW64 (x, y) +#define SET_LDOUBLE_WORDS64(x,y,z) SET_FLOAT128_WORDS64 (x, y, z) + +#undef IEEE854_LONG_DOUBLE_BIAS +#define IEEE854_LONG_DOUBLE_BIAS IEEE854_FLOAT128_BIAS + +#ifdef SET_RESTORE_ROUNDF128 +# undef SET_RESTORE_ROUNDL +# define SET_RESTORE_ROUNDL(RM) SET_RESTORE_ROUNDF128 (RM) +#endif + +#ifdef libc_feholdexcept_setroundf128 +# undef libc_feholdexcept_setroundl +# define libc_feholdexcept_setroundl(ENV, RM) \ + libc_feholdexcept_setroundf128 (ENV, RM) +#endif + +#ifdef libc_feupdateenv_testf128 +# undef libc_feupdateenv_testl +# define libc_feupdateenv_testl(ENV, EX) libc_feupdateenv_testf128 (ENV, EX) +#endif + +/* misc macros from the header below. */ +#include <fix-fp-int-convert-overflow.h> +#undef FIX_LDBL_LONG_CONVERT_OVERFLOW +#undef FIX_LDBL_LLONG_CONVERT_OVERFLOW +#define FIX_LDBL_LONG_CONVERT_OVERFLOW FIX_FLT128_LONG_CONVERT_OVERFLOW +#define FIX_LDBL_LLONG_CONVERT_OVERFLOW FIX_FLT128_LLONG_CONVERT_OVERFLOW + + +/* float.h constants. */ +#include <float.h> +#undef LDBL_DIG +#undef LDBL_EPSILON +#undef LDBL_MANT_DIG +#undef LDBL_MAX +#undef LDBL_MAX_10_EXP +#undef LDBL_MAX_EXP +#undef LDBL_MIN +#undef LDBL_MIN_10_EXP +#undef LDBL_MIN_EXP +#undef LDBL_TRUE_MIN +#define LDBL_DIG FLT128_DIG +#define LDBL_EPSILON FLT128_EPSILON +#define LDBL_MANT_DIG FLT128_MANT_DIG +#define LDBL_MAX FLT128_MAX +#define LDBL_MAX_10_EXP FLT128_MAX_10_EXP +#define LDBL_MAX_EXP FLT128_MAX_EXP +#define LDBL_MIN FLT128_MIN +#define LDBL_MIN_10_EXP FLT128_MIN_10_EXP +#define LDBL_MIN_EXP FLT128_MIN_EXP +#define LDBL_TRUE_MIN FLT128_TRUE_MIN + + +/* math.h GNU constants. */ +#undef M_El +#undef M_LOG2El +#undef M_LOG10El +#undef M_LN2l +#undef M_LN10l +#undef M_PIl +#undef M_PI_2l +#undef M_PI_4l +#undef M_1_PIl +#undef M_2_PIl +#undef M_2_SQRTPIl +#undef M_SQRT2l +#undef M_SQRT1_2l +#define M_El M_Ef128 +#define M_LOG2El M_LOG2Ef128 +#define M_LOG10El M_LOG10Ef128 +#define M_LN2l M_LN2f128 +#define M_LN10l M_LN10f128 +#define M_PIl M_PIf128 +#define M_PI_2l M_PI_2f128 +#define M_PI_4l M_PI_4f128 +#define M_1_PIl M_1_PIf128 +#define M_2_PIl M_2_PIf128 +#define M_2_SQRTPIl M_2_SQRTPIf128 +#define M_SQRT2l M_SQRT2f128 +#define M_SQRT1_2l M_SQRT1_2f128 + + +#include <libm-alias-ldouble.h> +#include <libm-alias-float128.h> +#undef libm_alias_ldouble_r +#define libm_alias_ldouble_r(from, to, r) libm_alias_float128_r (from, to, r) + + +#include <math/math-narrow.h> +#undef libm_alias_float_ldouble +#define libm_alias_float_ldouble(func) libm_alias_float32_float128 (func) +#undef libm_alias_double_ldouble +#define libm_alias_double_ldouble(func) libm_alias_float64_float128 (func) + + +/* IEEE function renames. */ +#define __ieee754_acoshl __ieee754_acoshf128 +#define __ieee754_acosl __ieee754_acosf128 +#define __ieee754_asinhl __ieee754_asinhf128 +#define __ieee754_asinl __ieee754_asinf128 +#define __ieee754_atan2l __ieee754_atan2f128 +#define __ieee754_atanhl __ieee754_atanhf128 +#define __ieee754_coshl __ieee754_coshf128 +#define __ieee754_cosl __ieee754_cosf128 +#define __ieee754_exp10l __ieee754_exp10f128 +#define __ieee754_exp2l __ieee754_exp2f128 +#define __ieee754_expl __ieee754_expf128 +#define __ieee754_fmodl __ieee754_fmodf128 +#define __ieee754_gammal_r __ieee754_gammaf128_r +#define __ieee754_hypotl __ieee754_hypotf128 +#define __ieee754_ilogbl __ieee754_ilogbf128 +#define __ieee754_j0l __ieee754_j0f128 +#define __ieee754_j1l __ieee754_j1f128 +#define __ieee754_jnl __ieee754_jnf128 +#define __ieee754_lgammal_r __ieee754_lgammaf128_r +#define __ieee754_log10l __ieee754_log10f128 +#define __ieee754_log2l __ieee754_log2f128 +#define __ieee754_logl __ieee754_logf128 +#define __ieee754_powl __ieee754_powf128 +#define __ieee754_rem_pio2l __ieee754_rem_pio2f128 +#define __ieee754_remainderl __ieee754_remainderf128 +#define __ieee754_sinhl __ieee754_sinhf128 +#define __ieee754_sqrtl __ieee754_sqrtf128 +#define __ieee754_y0l __ieee754_y0f128 +#define __ieee754_y1l __ieee754_y1f128 +#define __ieee754_ynl __ieee754_ynf128 + + +/* finite math entry points. */ +#define __acoshl_finite __acoshf128_finite +#define __acosl_finite __acosf128_finite +#define __asinl_finite __asinf128_finite +#define __atan2l_finite __atan2f128_finite +#define __atanhl_finite __atanhf128_finite +#define __coshl_finite __coshf128_finite +#define __cosl_finite __cosf128_finite +#define __exp10l_finite __exp10f128_finite +#define __exp2l_finite __exp2f128_finite +#define __expl_finite __expf128_finite +#define __fmodl_finite __fmodf128_finite +#define __hypotl_finite __hypotf128_finite +#define __ilogbl_finite __ilogbf128_finite +#define __j0l_finite __j0f128_finite +#define __j1l_finite __j1f128_finite +#define __jnl_finite __jnf128_finite +#define __lgammal_r_finite __lgammaf128_r_finite +#define __log10l_finite __log10f128_finite +#define __log2l_finite __log2f128_finite +#define __logl_finite __logf128_finite +#define __powl_finite __powf128_finite +#define __remainderl_finite __remainderf128_finite +#define __sinhl_finite __sinhf128_finite +#define __y0l_finite __y0f128_finite +#define __y1l_finite __y1f128_finite +#define __ynl_finite __ynf128_finite + + +/* internal function names. */ +#define __asinhl __asinhf128 +#define __atanl __atanf128 +#define __cbrtl __cbrtf128 +#define __ceill __ceilf128 +#define __copysignl __copysignf128 +#define __cosl __cosf128 +#define __erfcl __erfcf128 +#define __erfl __erff128 +#define __expl __expf128 +#define __expm1l __expm1f128 +#define __fabsl __fabsf128 +#define __fdiml __fdimf128 +#define __finitel __finitef128 +#define __floorl __floorf128 +#define __fmal __fmaf128 +#define __fmaxl __fmaxf128 +#define __fminl __fminf128 +#define __fpclassifyl __fpclassifyf128 +#define __frexpl __frexpf128 +#define __gammal_r_finite __gammaf128_r_finite +#define __getpayloadl __getpayloadf128 +#define __isinfl __isinff128 +#define __isnanl __isnanf128 +#define __issignalingl __issignalingf128 +#define __ldexpl __ldexpf128 +#define __llrintl __llrintf128 +#define __llroundl __llroundf128 +#define __log1pl __log1pf128 +#define __logbl __logbf128 +#define __logl __logf128 +#define __lrintl __lrintf128 +#define __lroundl __lroundf128 +#define __modfl __modff128 +#define __nearbyintl __nearbyintf128 +#define __nextafterl __nextafterf128 +#define __nextdownl __nextdownf128 +#define __nextupl __nextupf128 +#define __remquol __remquof128 +#define __rintl __rintf128 +#define __roundevenl __roundevenf128 +#define __roundl __roundf128 +#define __scalblnl __scalblnf128 +#define __scalbnl __scalbnf128 +#define __signbitl __signbitf128 +#define __sincosl __sincosf128 +#define __sinl __sinf128 +#define __sqrtl __sqrtf128 +#define __tanhl __tanhf128 +#define __tanl __tanf128 +#define __totalorderl __totalorderf128 +#define __totalordermagl __totalordermagf128 +#define __truncl __truncf128 +#define __x2y2m1l __x2y2m1f128 + +#define __faddl __f32addf128 +#define __daddl __f64addf128 +#define __fdivl __f32divf128 +#define __ddivl __f64divf128 +#define __fmull __f32mulf128 +#define __dmull __f64mulf128 +#define __fsubl __f32subf128 +#define __dsubl __f64subf128 + +/* __nexttowardf128 is not _Float128 API. */ +#define __nexttowardl __nexttowardf128_do_not_use +#define nexttowardl nexttowardf128_do_not_use + + +/* public entry points. */ +#define asinhl asinhf128 +#define atanl atanf128 +#define cbrtl cbrtf128 +#define ceill ceilf128 +#define copysignl copysignf128 +#define cosl cosf128 +#define erfcl erfcf128 +#define erfl erff128 +#define expl expf128 +#define expm1l expm1f128 +#define fabsl fabsf128 +#define fdiml fdimf128 +#define finitel finitef128_do_not_use +#define floorl floorf128 +#define fmal fmaf128 +#define fmaxl fmaxf128 +#define fminl fminf128 +#define frexpl frexpf128 +#define getpayloadl getpayloadf128 +#define isinfl isinff128_do_not_use +#define isnanl isnanf128_do_not_use +#define ldexpl ldexpf128 +#define llrintl llrintf128 +#define llroundl llroundf128 +#define log1pl log1pf128 +#define logbl logbf128 +#define logl logf128 +#define lrintl lrintf128 +#define lroundl lroundf128 +#define modfl modff128 +#define nanl nanf128 +#define nearbyintl nearbyintf128 +#define nextafterl nextafterf128 +#define nextdownl nextdownf128 +#define nextupl nextupf128 +#define remquol remquof128 +#define rintl rintf128 +#define roundevenl roundevenf128 +#define roundl roundf128 +#define scalbnl scalbnf128 +#define sincosl sincosf128 +#define sinl sinf128 +#define sqrtl sqrtf128 +#define tanhl tanhf128 +#define tanl tanf128 +#define totalorderl totalorderf128 +#define totalordermagl totalordermagf128 +#define truncl truncf128 + + +/* misc internal renames. */ +#define __builtin_fmal __builtin_fmaf128 +#define __expl_table __expf128_table +#define __gamma_productl __gamma_productf128 +#define __kernel_cosl __kernel_cosf128 +#define __kernel_rem_pio2l __kernel_rem_pio2f128 +#define __kernel_sincosl __kernel_sincosf128 +#define __kernel_sinl __kernel_sinf128 +#define __kernel_tanl __kernel_tanf128 +#define __lgamma_negl __lgamma_negf128 +#define __lgamma_productl __lgamma_productf128 +#define __mpn_extract_long_double __mpn_extract_float128 +#define __sincosl_table __sincosf128_table +#define mul_splitl mul_splitf128 + +/* Builtin renames. */ +#define __builtin_copysignl __builtin_copysignf128 +#define __builtin_signbitl __builtin_signbit + +/* Get the constant suffix from bits/floatn-compat.h. */ +#define L(x) __f128 (x) + +static inline void +mul_splitf128 (_Float128 *hi, _Float128 *lo, _Float128 x, _Float128 y) +{ +#ifdef __FP_FAST_FMAF128 + /* Fast built-in fused multiply-add. */ + *hi = x * y; + *lo = __builtin_fmal (x, y, -*hi); +#else + /* Apply Dekker's algorithm. */ + *hi = x * y; +# define C ((1LL << (FLT128_MANT_DIG + 1) / 2) + 1) + _Float128 x1 = x * C; + _Float128 y1 = y * C; +# undef C + x1 = (x - x1) + x1; + y1 = (y - y1) + y1; + _Float128 x2 = x - x1; + _Float128 y2 = y - y1; + *lo = (((x1 * y1 - *hi) + x1 * y2) + x2 * y1) + x2 * y2; +#endif +} |