diff options
91 files changed, 765 insertions, 0 deletions
diff --git a/ChangeLog b/ChangeLog index bd66603025..c085c9b637 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,103 @@ +2017-05-25 Paul E. Murphy <murphyp@linux.vnet.ibm.com> + Gabriel F. T. Gomes <gftg@linux.vnet.ibm.com> + Tulio Magno Quites Machado Filho <tuliom@linux.vnet.ibm.com> + + * math/math.h [__HAVE_FLOAT128 && __GLIBC_USE (IEC_60559_TYPES_EXT)] + (SNANF128): New macro. + * math/e_sqrtf128.c: New file. + * math/s_fmaxmag_template.c: Include math_private.h in order to + make inline expansion of fabs128(). + * math/s_fminmag_template.c: Likewise. + + * sysdeps/ieee754/float128/e_acosf128.c: New file. + * sysdeps/ieee754/float128/e_acoshf128.c: New file. + * sysdeps/ieee754/float128/e_asinf128.c: New file. + * sysdeps/ieee754/float128/e_atan2f128.c: New file. + * sysdeps/ieee754/float128/e_atanhf128.c: New file. + * sysdeps/ieee754/float128/e_coshf128.c: New file. + * sysdeps/ieee754/float128/e_exp10f128.c: New file. + * sysdeps/ieee754/float128/e_expf128.c: New file. + * sysdeps/ieee754/float128/e_fmodf128.c: New file. + * sysdeps/ieee754/float128/e_gammaf128_r.c: New file. + * sysdeps/ieee754/float128/e_hypotf128.c: New file. + * sysdeps/ieee754/float128/e_ilogbf128.c: New file. + * sysdeps/ieee754/float128/e_j0f128.c: New file. + * sysdeps/ieee754/float128/e_j1f128.c: New file. + * sysdeps/ieee754/float128/e_jnf128.c: New file. + * sysdeps/ieee754/float128/e_lgammaf128_r.c: New file. + * sysdeps/ieee754/float128/e_log10f128.c: New file. + * sysdeps/ieee754/float128/e_log2f128.c: New file. + * sysdeps/ieee754/float128/e_logf128.c: New file. + * sysdeps/ieee754/float128/e_powf128.c: New file. + * sysdeps/ieee754/float128/e_rem_pio2f128.c: New file. + * sysdeps/ieee754/float128/e_remainderf128.c: New file. + * sysdeps/ieee754/float128/e_scalbf128.c: New file. + * sysdeps/ieee754/float128/e_sinhf128.c: New file. + * sysdeps/ieee754/float128/float128_private.h: New file. + * sysdeps/ieee754/float128/gamma_productf128.c: New file. + * sysdeps/ieee754/float128/ieee754_float128.h: New file. + * sysdeps/ieee754/float128/k_cosf128.c: New file. + * sysdeps/ieee754/float128/k_sincosf128.c: New file. + * sysdeps/ieee754/float128/k_sinf128.c: New file. + * sysdeps/ieee754/float128/k_tanf128.c: New file. + * sysdeps/ieee754/float128/lgamma_negf128.c: New file. + * sysdeps/ieee754/float128/lgamma_productf128.c: New file. + * sysdeps/ieee754/float128/s_asinhf128.c: New file. + * sysdeps/ieee754/float128/s_atanf128.c: New file. + * sysdeps/ieee754/float128/s_cbrtf128.c: New file. + * sysdeps/ieee754/float128/s_ceilf128.c: New file. + * sysdeps/ieee754/float128/s_copysignf128.c: New file. + * sysdeps/ieee754/float128/s_cosf128.c: New file. + * sysdeps/ieee754/float128/s_erff128.c: New file. + * sysdeps/ieee754/float128/s_expm1f128.c: New file. + * sysdeps/ieee754/float128/s_fabsf128.c: New file. + * sysdeps/ieee754/float128/s_finitef128.c: New file. + * sysdeps/ieee754/float128/s_floorf128.c: New file. + * sysdeps/ieee754/float128/s_fmaf128.c: New file. + * sysdeps/ieee754/float128/s_fpclassifyf128.c: New file. + * sysdeps/ieee754/float128/s_frexpf128.c: New file. + * sysdeps/ieee754/float128/s_fromfpf128.c: New file. + * sysdeps/ieee754/float128/s_fromfpxf128.c: New file. + * sysdeps/ieee754/float128/s_getpayloadf128.c: New file. + * sysdeps/ieee754/float128/s_isinff128.c: New file. + * sysdeps/ieee754/float128/s_isnanf128.c: New file. + * sysdeps/ieee754/float128/s_issignalingf128.c: New file. + * sysdeps/ieee754/float128/s_llrintf128.c: New file. + * sysdeps/ieee754/float128/s_llroundf128.c: New file. + * sysdeps/ieee754/float128/s_log1pf128.c: New file. + * sysdeps/ieee754/float128/s_logbf128.c: New file. + * sysdeps/ieee754/float128/s_lrintf128.c: New file. + * sysdeps/ieee754/float128/s_lroundf128.c: New file. + * sysdeps/ieee754/float128/s_modff128.c: New file. + * sysdeps/ieee754/float128/s_nearbyintf128.c: New file. + * sysdeps/ieee754/float128/s_nextafterf128.c: New file. + * sysdeps/ieee754/float128/s_nexttowardf128.c: New file. + * sysdeps/ieee754/float128/s_nextupf128.c: New file. + * sysdeps/ieee754/float128/s_remquof128.c: New file. + * sysdeps/ieee754/float128/s_rintf128.c: New file. + * sysdeps/ieee754/float128/s_roundevenf128.c: New file. + * sysdeps/ieee754/float128/s_roundf128.c: New file. + * sysdeps/ieee754/float128/s_scalblnf128.c: New file. + * sysdeps/ieee754/float128/s_scalbnf128.c: New file. + * sysdeps/ieee754/float128/s_setpayloadf128.c: New file. + * sysdeps/ieee754/float128/s_setpayloadsigf128.c: New file. + * sysdeps/ieee754/float128/s_signbitf128.c: New file. + * sysdeps/ieee754/float128/s_significandf128.c: New file. + * sysdeps/ieee754/float128/s_sincosf128.c: New file. + * sysdeps/ieee754/float128/s_sinf128.c: New file. + * sysdeps/ieee754/float128/s_tanf128.c: New file. + * sysdeps/ieee754/float128/s_tanhf128.c: New file. + * sysdeps/ieee754/float128/s_totalorderf128.c: New file. + * sysdeps/ieee754/float128/s_totalordermagf128.c: New file. + * sysdeps/ieee754/float128/s_truncf128.c: New file. + * sysdeps/ieee754/float128/s_ufromfpf128.c: New file. + * sysdeps/ieee754/float128/s_ufromfpxf128.c: New file. + * sysdeps/ieee754/float128/t_sincosf128.c: New file. + * sysdeps/ieee754/float128/x2y2m1f128.c: New file. + + * sysdeps/ieee754/ldbl-128ibm/bits/iscanonical.h: + (__iscanonicalf128): Define as a macro. + 2017-05-24 Steve Ellcey <sellcey@caviumnetworks.com> * sysdeps/aarch64/memcpy.S (MEMMOVE, MEMCPY): New macros. diff --git a/math/e_sqrtf128.c b/math/e_sqrtf128.c new file mode 100644 index 0000000000..acf8ba45e3 --- /dev/null +++ b/math/e_sqrtf128.c @@ -0,0 +1,14 @@ +#include <math.h> +#include <stdio.h> +#include <errno.h> + +_Float128 +__ieee754_sqrtf128 (_Float128 x) +{ + fputs ("__ieee754_sqrtf128 not implemented\n", stderr); + __set_errno (ENOSYS); + return 0.0; +} +strong_alias (__ieee754_sqrtf128, __sqrtf128_finite) + +stub_warning (sqrtf128) diff --git a/math/math.h b/math/math.h index 5f006a8358..e21708045a 100644 --- a/math/math.h +++ b/math/math.h @@ -64,6 +64,9 @@ __BEGIN_DECLS # define SNANL (__builtin_nansl ("")) # endif #endif +#if __HAVE_FLOAT128 && __GLIBC_USE (IEC_60559_TYPES_EXT) +# define SNANF128 (__builtin_nansf128 ("")) +#endif /* Get __GLIBC_FLT_EVAL_METHOD. */ #include <bits/flt-eval-method.h> diff --git a/math/s_fmaxmag_template.c b/math/s_fmaxmag_template.c index 6beef03292..f7c824886c 100644 --- a/math/s_fmaxmag_template.c +++ b/math/s_fmaxmag_template.c @@ -17,6 +17,7 @@ <http://www.gnu.org/licenses/>. */ #include <math.h> +#include <math_private.h> FLOAT M_DECL_FUNC (__fmaxmag) (FLOAT x, FLOAT y) diff --git a/math/s_fminmag_template.c b/math/s_fminmag_template.c index a1e4c2b135..1054c38411 100644 --- a/math/s_fminmag_template.c +++ b/math/s_fminmag_template.c @@ -17,6 +17,7 @@ <http://www.gnu.org/licenses/>. */ #include <math.h> +#include <math_private.h> FLOAT M_DECL_FUNC (__fminmag) (FLOAT x, FLOAT y) diff --git a/sysdeps/ieee754/float128/e_acosf128.c b/sysdeps/ieee754/float128/e_acosf128.c new file mode 100644 index 0000000000..7ddf7dcdf1 --- /dev/null +++ b/sysdeps/ieee754/float128/e_acosf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_acosl.c" diff --git a/sysdeps/ieee754/float128/e_acoshf128.c b/sysdeps/ieee754/float128/e_acoshf128.c new file mode 100644 index 0000000000..f6dd40cd88 --- /dev/null +++ b/sysdeps/ieee754/float128/e_acoshf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_acoshl.c" diff --git a/sysdeps/ieee754/float128/e_asinf128.c b/sysdeps/ieee754/float128/e_asinf128.c new file mode 100644 index 0000000000..133ab8d875 --- /dev/null +++ b/sysdeps/ieee754/float128/e_asinf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_asinl.c" diff --git a/sysdeps/ieee754/float128/e_atan2f128.c b/sysdeps/ieee754/float128/e_atan2f128.c new file mode 100644 index 0000000000..9aa740f770 --- /dev/null +++ b/sysdeps/ieee754/float128/e_atan2f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_atan2l.c" diff --git a/sysdeps/ieee754/float128/e_atanhf128.c b/sysdeps/ieee754/float128/e_atanhf128.c new file mode 100644 index 0000000000..f26c8d54d2 --- /dev/null +++ b/sysdeps/ieee754/float128/e_atanhf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_atanhl.c" diff --git a/sysdeps/ieee754/float128/e_coshf128.c b/sysdeps/ieee754/float128/e_coshf128.c new file mode 100644 index 0000000000..2abf067f64 --- /dev/null +++ b/sysdeps/ieee754/float128/e_coshf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_coshl.c" diff --git a/sysdeps/ieee754/float128/e_exp10f128.c b/sysdeps/ieee754/float128/e_exp10f128.c new file mode 100644 index 0000000000..b3468d256b --- /dev/null +++ b/sysdeps/ieee754/float128/e_exp10f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_exp10l.c" diff --git a/sysdeps/ieee754/float128/e_expf128.c b/sysdeps/ieee754/float128/e_expf128.c new file mode 100644 index 0000000000..b727b17cc9 --- /dev/null +++ b/sysdeps/ieee754/float128/e_expf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_expl.c" diff --git a/sysdeps/ieee754/float128/e_fmodf128.c b/sysdeps/ieee754/float128/e_fmodf128.c new file mode 100644 index 0000000000..ed8a7491ed --- /dev/null +++ b/sysdeps/ieee754/float128/e_fmodf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_fmodl.c" diff --git a/sysdeps/ieee754/float128/e_gammaf128_r.c b/sysdeps/ieee754/float128/e_gammaf128_r.c new file mode 100644 index 0000000000..895ac6374d --- /dev/null +++ b/sysdeps/ieee754/float128/e_gammaf128_r.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_gammal_r.c" diff --git a/sysdeps/ieee754/float128/e_hypotf128.c b/sysdeps/ieee754/float128/e_hypotf128.c new file mode 100644 index 0000000000..1f06555505 --- /dev/null +++ b/sysdeps/ieee754/float128/e_hypotf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_hypotl.c" diff --git a/sysdeps/ieee754/float128/e_ilogbf128.c b/sysdeps/ieee754/float128/e_ilogbf128.c new file mode 100644 index 0000000000..2861801854 --- /dev/null +++ b/sysdeps/ieee754/float128/e_ilogbf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_ilogbl.c" diff --git a/sysdeps/ieee754/float128/e_j0f128.c b/sysdeps/ieee754/float128/e_j0f128.c new file mode 100644 index 0000000000..b624b5c596 --- /dev/null +++ b/sysdeps/ieee754/float128/e_j0f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_j0l.c" diff --git a/sysdeps/ieee754/float128/e_j1f128.c b/sysdeps/ieee754/float128/e_j1f128.c new file mode 100644 index 0000000000..445428e742 --- /dev/null +++ b/sysdeps/ieee754/float128/e_j1f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_j1l.c" diff --git a/sysdeps/ieee754/float128/e_jnf128.c b/sysdeps/ieee754/float128/e_jnf128.c new file mode 100644 index 0000000000..7854e11a1a --- /dev/null +++ b/sysdeps/ieee754/float128/e_jnf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_jnl.c" diff --git a/sysdeps/ieee754/float128/e_lgammaf128_r.c b/sysdeps/ieee754/float128/e_lgammaf128_r.c new file mode 100644 index 0000000000..3517ac32af --- /dev/null +++ b/sysdeps/ieee754/float128/e_lgammaf128_r.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_lgammal_r.c" diff --git a/sysdeps/ieee754/float128/e_log10f128.c b/sysdeps/ieee754/float128/e_log10f128.c new file mode 100644 index 0000000000..1c3341e412 --- /dev/null +++ b/sysdeps/ieee754/float128/e_log10f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_log10l.c" diff --git a/sysdeps/ieee754/float128/e_log2f128.c b/sysdeps/ieee754/float128/e_log2f128.c new file mode 100644 index 0000000000..36becaadc4 --- /dev/null +++ b/sysdeps/ieee754/float128/e_log2f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_log2l.c" diff --git a/sysdeps/ieee754/float128/e_logf128.c b/sysdeps/ieee754/float128/e_logf128.c new file mode 100644 index 0000000000..b0c9975caf --- /dev/null +++ b/sysdeps/ieee754/float128/e_logf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_logl.c" diff --git a/sysdeps/ieee754/float128/e_powf128.c b/sysdeps/ieee754/float128/e_powf128.c new file mode 100644 index 0000000000..3afaf7f6a6 --- /dev/null +++ b/sysdeps/ieee754/float128/e_powf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_powl.c" diff --git a/sysdeps/ieee754/float128/e_rem_pio2f128.c b/sysdeps/ieee754/float128/e_rem_pio2f128.c new file mode 100644 index 0000000000..86c2ca1024 --- /dev/null +++ b/sysdeps/ieee754/float128/e_rem_pio2f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_rem_pio2l.c" diff --git a/sysdeps/ieee754/float128/e_remainderf128.c b/sysdeps/ieee754/float128/e_remainderf128.c new file mode 100644 index 0000000000..90c18f8493 --- /dev/null +++ b/sysdeps/ieee754/float128/e_remainderf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_remainderl.c" diff --git a/sysdeps/ieee754/float128/e_scalbf128.c b/sysdeps/ieee754/float128/e_scalbf128.c new file mode 100644 index 0000000000..067b724164 --- /dev/null +++ b/sysdeps/ieee754/float128/e_scalbf128.c @@ -0,0 +1 @@ +/* Not defined for _FloatN types. */ diff --git a/sysdeps/ieee754/float128/e_sinhf128.c b/sysdeps/ieee754/float128/e_sinhf128.c new file mode 100644 index 0000000000..42a54e0015 --- /dev/null +++ b/sysdeps/ieee754/float128/e_sinhf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/e_sinhl.c" diff --git a/sysdeps/ieee754/float128/float128_private.h b/sysdeps/ieee754/float128/float128_private.h new file mode 100644 index 0000000000..894e6ac8cb --- /dev/null +++ b/sysdeps/ieee754/float128/float128_private.h @@ -0,0 +1,324 @@ +/* _Float128 overrides for building ldbl-128 as _Float128. + Copyright (C) 2017 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 <math.h> +#undef HUGE_VALL +#define HUGE_VALL HUGE_VAL_F128 +#include <math/mul_splitl.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() SET_RESTORE_ROUNDF128() +#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 + + +/* 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 __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 __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 __truncl __truncf128 +#define __x2y2m1l __x2y2m1f128 + +/* __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 __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 +} diff --git a/sysdeps/ieee754/float128/gamma_productf128.c b/sysdeps/ieee754/float128/gamma_productf128.c new file mode 100644 index 0000000000..be2271f12f --- /dev/null +++ b/sysdeps/ieee754/float128/gamma_productf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/gamma_productl.c" diff --git a/sysdeps/ieee754/float128/ieee754_float128.h b/sysdeps/ieee754/float128/ieee754_float128.h new file mode 100644 index 0000000000..e8e7211d88 --- /dev/null +++ b/sysdeps/ieee754/float128/ieee754_float128.h @@ -0,0 +1,140 @@ +/* _Float128 IEEE like macros. + Copyright (C) 2017 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/>. */ +#ifndef _IEEE754_FLOAT128_H +#define _IEEE754_FLOAT128_H + +#include <endian.h> +#include <stdint.h> + +# if __FLOAT_WORD_ORDER == BIG_ENDIAN +# define __FLT_EORDER2(t, a, b) t a; t b; +# define __FLT_EORDER4(t, a, b, c, d) \ + t a; t b; t c; t d; +# define __FLT_EORDER6(t, a, b, c, d, e, f) \ + t a; t b; t c; t d; t e; t f; +# define __FLT_EORDER7(t, a, b, c, d, e, f, g) \ + t a; t b; t c; t d; t e; t f; t g; +# else +# define __FLT_EORDER2(t, a, b) \ + t b; t a; +# define __FLT_EORDER4(t, a, b, c, d) \ + t d; t c; t b; t a; +# define __FLT_EORDER6(t, a, b, c, d, e, f) \ + t f; t e; t d; t c; t b; t a; +# define __FLT_EORDER7(t, a, b, c, d, e, f, g) \ + t g; t f; t e; t d; t c; t b; t a; +# endif + +/* A union which permits us to convert between _Float128 and + four 32 bit ints or two 64 bit ints. */ + +typedef union +{ + _Float128 value; + struct + { + __FLT_EORDER2 (uint64_t, msw, lsw); + } parts64; + struct + { + __FLT_EORDER4 (uint32_t, w0, w1, w2, w3); + } parts32; +} ieee854_float128_shape_type; + +/* Get two 64 bit ints from a _Float128. */ + +# define GET_FLOAT128_WORDS64(ix0,ix1,d) \ +do { \ + ieee854_float128_shape_type qw_u; \ + qw_u.value = (d); \ + (ix0) = qw_u.parts64.msw; \ + (ix1) = qw_u.parts64.lsw; \ +} while (0) + +/* Set a _Float128 from two 64 bit ints. */ + +# define SET_FLOAT128_WORDS64(d,ix0,ix1) \ +do { \ + ieee854_float128_shape_type qw_u; \ + qw_u.parts64.msw = (ix0); \ + qw_u.parts64.lsw = (ix1); \ + (d) = qw_u.value; \ +} while (0) + +/* Get the more significant 64 bits of a _Float128 mantissa. */ + +# define GET_FLOAT128_MSW64(v,d) \ +do { \ + ieee854_float128_shape_type sh_u; \ + sh_u.value = (d); \ + (v) = sh_u.parts64.msw; \ +} while (0) + +/* Set the more significant 64 bits of a _Float128 mantissa from an int. */ + +# define SET_FLOAT128_MSW64(d,v) \ +do { \ + ieee854_float128_shape_type sh_u; \ + sh_u.value = (d); \ + sh_u.parts64.msw = (v); \ + (d) = sh_u.value; \ +} while (0) + +/* Get the least significant 64 bits of a _Float128 mantissa. */ + +# define GET_FLOAT128_LSW64(v,d) \ +do { \ + ieee854_float128_shape_type sh_u; \ + sh_u.value = (d); \ + (v) = sh_u.parts64.lsw; \ +} while (0) + +/* Likewise, some helper macros which are exposed via ieee754.h for + C99 real types, but not _Float128. */ + +union ieee854_float128 + { + _Float128 d; + + /* This is the IEEE 854 quad-precision format. */ + struct + { + __FLT_EORDER6 (unsigned int, negative:1, + exponent:15, + mantissa0:16, + mantissa1:32, + mantissa2:32, + mantissa3:32) + } ieee; + + /* This format makes it easier to see if a NaN is a signalling NaN. */ + struct + { + __FLT_EORDER7 (unsigned int, negative:1, + exponent:15, + quiet_nan:1, + mantissa0:15, + mantissa1:32, + mantissa2:32, + mantissa3:32) + } ieee_nan; + }; + +#define IEEE854_FLOAT128_BIAS 0x3fff /* Added to exponent. */ + +#endif diff --git a/sysdeps/ieee754/float128/k_cosf128.c b/sysdeps/ieee754/float128/k_cosf128.c new file mode 100644 index 0000000000..9db0906e9a --- /dev/null +++ b/sysdeps/ieee754/float128/k_cosf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/k_cosl.c" diff --git a/sysdeps/ieee754/float128/k_sincosf128.c b/sysdeps/ieee754/float128/k_sincosf128.c new file mode 100644 index 0000000000..14c0f1eccf --- /dev/null +++ b/sysdeps/ieee754/float128/k_sincosf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/k_sincosl.c" diff --git a/sysdeps/ieee754/float128/k_sinf128.c b/sysdeps/ieee754/float128/k_sinf128.c new file mode 100644 index 0000000000..f3acf1cfe1 --- /dev/null +++ b/sysdeps/ieee754/float128/k_sinf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/k_sinl.c" diff --git a/sysdeps/ieee754/float128/k_tanf128.c b/sysdeps/ieee754/float128/k_tanf128.c new file mode 100644 index 0000000000..ca6be539f7 --- /dev/null +++ b/sysdeps/ieee754/float128/k_tanf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/k_tanl.c" diff --git a/sysdeps/ieee754/float128/lgamma_negf128.c b/sysdeps/ieee754/float128/lgamma_negf128.c new file mode 100644 index 0000000000..9c16f93025 --- /dev/null +++ b/sysdeps/ieee754/float128/lgamma_negf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/lgamma_negl.c" diff --git a/sysdeps/ieee754/float128/lgamma_productf128.c b/sysdeps/ieee754/float128/lgamma_productf128.c new file mode 100644 index 0000000000..5efe5dd576 --- /dev/null +++ b/sysdeps/ieee754/float128/lgamma_productf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/lgamma_productl.c" diff --git a/sysdeps/ieee754/float128/s_asinhf128.c b/sysdeps/ieee754/float128/s_asinhf128.c new file mode 100644 index 0000000000..7b93d8cf3a --- /dev/null +++ b/sysdeps/ieee754/float128/s_asinhf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_asinhl.c" diff --git a/sysdeps/ieee754/float128/s_atanf128.c b/sysdeps/ieee754/float128/s_atanf128.c new file mode 100644 index 0000000000..9b4d7ecec3 --- /dev/null +++ b/sysdeps/ieee754/float128/s_atanf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_atanl.c" diff --git a/sysdeps/ieee754/float128/s_cbrtf128.c b/sysdeps/ieee754/float128/s_cbrtf128.c new file mode 100644 index 0000000000..3bd5797373 --- /dev/null +++ b/sysdeps/ieee754/float128/s_cbrtf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_cbrtl.c" diff --git a/sysdeps/ieee754/float128/s_ceilf128.c b/sysdeps/ieee754/float128/s_ceilf128.c new file mode 100644 index 0000000000..0af15f5f03 --- /dev/null +++ b/sysdeps/ieee754/float128/s_ceilf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_ceill.c" diff --git a/sysdeps/ieee754/float128/s_copysignf128.c b/sysdeps/ieee754/float128/s_copysignf128.c new file mode 100644 index 0000000000..808f7abbc0 --- /dev/null +++ b/sysdeps/ieee754/float128/s_copysignf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_copysignl.c" diff --git a/sysdeps/ieee754/float128/s_cosf128.c b/sysdeps/ieee754/float128/s_cosf128.c new file mode 100644 index 0000000000..8ba552695c --- /dev/null +++ b/sysdeps/ieee754/float128/s_cosf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_cosl.c" diff --git a/sysdeps/ieee754/float128/s_erff128.c b/sysdeps/ieee754/float128/s_erff128.c new file mode 100644 index 0000000000..ac16ad6665 --- /dev/null +++ b/sysdeps/ieee754/float128/s_erff128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_erfl.c" diff --git a/sysdeps/ieee754/float128/s_expm1f128.c b/sysdeps/ieee754/float128/s_expm1f128.c new file mode 100644 index 0000000000..ea28d89db4 --- /dev/null +++ b/sysdeps/ieee754/float128/s_expm1f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_expm1l.c" diff --git a/sysdeps/ieee754/float128/s_fabsf128.c b/sysdeps/ieee754/float128/s_fabsf128.c new file mode 100644 index 0000000000..79ba47c3fd --- /dev/null +++ b/sysdeps/ieee754/float128/s_fabsf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_fabsl.c" diff --git a/sysdeps/ieee754/float128/s_finitef128.c b/sysdeps/ieee754/float128/s_finitef128.c new file mode 100644 index 0000000000..801de88e32 --- /dev/null +++ b/sysdeps/ieee754/float128/s_finitef128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_finitel.c" diff --git a/sysdeps/ieee754/float128/s_floorf128.c b/sysdeps/ieee754/float128/s_floorf128.c new file mode 100644 index 0000000000..18298436a1 --- /dev/null +++ b/sysdeps/ieee754/float128/s_floorf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_floorl.c" diff --git a/sysdeps/ieee754/float128/s_fmaf128.c b/sysdeps/ieee754/float128/s_fmaf128.c new file mode 100644 index 0000000000..6497895c8d --- /dev/null +++ b/sysdeps/ieee754/float128/s_fmaf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_fmal.c" diff --git a/sysdeps/ieee754/float128/s_fpclassifyf128.c b/sysdeps/ieee754/float128/s_fpclassifyf128.c new file mode 100644 index 0000000000..15131dc4a2 --- /dev/null +++ b/sysdeps/ieee754/float128/s_fpclassifyf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_fpclassifyl.c" diff --git a/sysdeps/ieee754/float128/s_frexpf128.c b/sysdeps/ieee754/float128/s_frexpf128.c new file mode 100644 index 0000000000..7b040b3e33 --- /dev/null +++ b/sysdeps/ieee754/float128/s_frexpf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_frexpl.c" diff --git a/sysdeps/ieee754/float128/s_fromfpf128.c b/sysdeps/ieee754/float128/s_fromfpf128.c new file mode 100644 index 0000000000..891de3d448 --- /dev/null +++ b/sysdeps/ieee754/float128/s_fromfpf128.c @@ -0,0 +1,5 @@ +#define UNSIGNED 0 +#define INEXACT 0 +#define FUNC fromfpf128 +#include <float128_private.h> +#include "../ldbl-128/s_fromfpl_main.c" diff --git a/sysdeps/ieee754/float128/s_fromfpxf128.c b/sysdeps/ieee754/float128/s_fromfpxf128.c new file mode 100644 index 0000000000..21676fab03 --- /dev/null +++ b/sysdeps/ieee754/float128/s_fromfpxf128.c @@ -0,0 +1,5 @@ +#define UNSIGNED 0 +#define INEXACT 1 +#define FUNC fromfpxf128 +#include <float128_private.h> +#include "../ldbl-128/s_fromfpl_main.c" diff --git a/sysdeps/ieee754/float128/s_getpayloadf128.c b/sysdeps/ieee754/float128/s_getpayloadf128.c new file mode 100644 index 0000000000..2e2607a81b --- /dev/null +++ b/sysdeps/ieee754/float128/s_getpayloadf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_getpayloadl.c" diff --git a/sysdeps/ieee754/float128/s_isinff128.c b/sysdeps/ieee754/float128/s_isinff128.c new file mode 100644 index 0000000000..62cc424b16 --- /dev/null +++ b/sysdeps/ieee754/float128/s_isinff128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_isinfl.c" diff --git a/sysdeps/ieee754/float128/s_isnanf128.c b/sysdeps/ieee754/float128/s_isnanf128.c new file mode 100644 index 0000000000..efba24059a --- /dev/null +++ b/sysdeps/ieee754/float128/s_isnanf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_isnanl.c" diff --git a/sysdeps/ieee754/float128/s_issignalingf128.c b/sysdeps/ieee754/float128/s_issignalingf128.c new file mode 100644 index 0000000000..1d4599544f --- /dev/null +++ b/sysdeps/ieee754/float128/s_issignalingf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_issignalingl.c" diff --git a/sysdeps/ieee754/float128/s_llrintf128.c b/sysdeps/ieee754/float128/s_llrintf128.c new file mode 100644 index 0000000000..bb9ca580cd --- /dev/null +++ b/sysdeps/ieee754/float128/s_llrintf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_llrintl.c" diff --git a/sysdeps/ieee754/float128/s_llroundf128.c b/sysdeps/ieee754/float128/s_llroundf128.c new file mode 100644 index 0000000000..be54a90608 --- /dev/null +++ b/sysdeps/ieee754/float128/s_llroundf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_llroundl.c" diff --git a/sysdeps/ieee754/float128/s_log1pf128.c b/sysdeps/ieee754/float128/s_log1pf128.c new file mode 100644 index 0000000000..48bb84f987 --- /dev/null +++ b/sysdeps/ieee754/float128/s_log1pf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_log1pl.c" diff --git a/sysdeps/ieee754/float128/s_logbf128.c b/sysdeps/ieee754/float128/s_logbf128.c new file mode 100644 index 0000000000..167384a8d4 --- /dev/null +++ b/sysdeps/ieee754/float128/s_logbf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_logbl.c" diff --git a/sysdeps/ieee754/float128/s_lrintf128.c b/sysdeps/ieee754/float128/s_lrintf128.c new file mode 100644 index 0000000000..1cfa9d7ca4 --- /dev/null +++ b/sysdeps/ieee754/float128/s_lrintf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_lrintl.c" diff --git a/sysdeps/ieee754/float128/s_lroundf128.c b/sysdeps/ieee754/float128/s_lroundf128.c new file mode 100644 index 0000000000..13ba9f220d --- /dev/null +++ b/sysdeps/ieee754/float128/s_lroundf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_lroundl.c" diff --git a/sysdeps/ieee754/float128/s_modff128.c b/sysdeps/ieee754/float128/s_modff128.c new file mode 100644 index 0000000000..4618c6c380 --- /dev/null +++ b/sysdeps/ieee754/float128/s_modff128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_modfl.c" diff --git a/sysdeps/ieee754/float128/s_nearbyintf128.c b/sysdeps/ieee754/float128/s_nearbyintf128.c new file mode 100644 index 0000000000..e61a3b3bfb --- /dev/null +++ b/sysdeps/ieee754/float128/s_nearbyintf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_nearbyintl.c" diff --git a/sysdeps/ieee754/float128/s_nextafterf128.c b/sysdeps/ieee754/float128/s_nextafterf128.c new file mode 100644 index 0000000000..2c43a00384 --- /dev/null +++ b/sysdeps/ieee754/float128/s_nextafterf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_nextafterl.c" diff --git a/sysdeps/ieee754/float128/s_nexttowardf128.c b/sysdeps/ieee754/float128/s_nexttowardf128.c new file mode 100644 index 0000000000..006e4c98b1 --- /dev/null +++ b/sysdeps/ieee754/float128/s_nexttowardf128.c @@ -0,0 +1 @@ +/* This function does not exist for _FloatN types. */ diff --git a/sysdeps/ieee754/float128/s_nextupf128.c b/sysdeps/ieee754/float128/s_nextupf128.c new file mode 100644 index 0000000000..7d5d0b8c72 --- /dev/null +++ b/sysdeps/ieee754/float128/s_nextupf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_nextupl.c" diff --git a/sysdeps/ieee754/float128/s_remquof128.c b/sysdeps/ieee754/float128/s_remquof128.c new file mode 100644 index 0000000000..1cef61ab1b --- /dev/null +++ b/sysdeps/ieee754/float128/s_remquof128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_remquol.c" diff --git a/sysdeps/ieee754/float128/s_rintf128.c b/sysdeps/ieee754/float128/s_rintf128.c new file mode 100644 index 0000000000..2adb95f360 --- /dev/null +++ b/sysdeps/ieee754/float128/s_rintf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_rintl.c" diff --git a/sysdeps/ieee754/float128/s_roundevenf128.c b/sysdeps/ieee754/float128/s_roundevenf128.c new file mode 100644 index 0000000000..5a9b3f395f --- /dev/null +++ b/sysdeps/ieee754/float128/s_roundevenf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_roundevenl.c" diff --git a/sysdeps/ieee754/float128/s_roundf128.c b/sysdeps/ieee754/float128/s_roundf128.c new file mode 100644 index 0000000000..1eb36f2a5e --- /dev/null +++ b/sysdeps/ieee754/float128/s_roundf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_roundl.c" diff --git a/sysdeps/ieee754/float128/s_scalblnf128.c b/sysdeps/ieee754/float128/s_scalblnf128.c new file mode 100644 index 0000000000..999223c517 --- /dev/null +++ b/sysdeps/ieee754/float128/s_scalblnf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_scalblnl.c" diff --git a/sysdeps/ieee754/float128/s_scalbnf128.c b/sysdeps/ieee754/float128/s_scalbnf128.c new file mode 100644 index 0000000000..0e7ab2663b --- /dev/null +++ b/sysdeps/ieee754/float128/s_scalbnf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_scalbnl.c" diff --git a/sysdeps/ieee754/float128/s_setpayloadf128.c b/sysdeps/ieee754/float128/s_setpayloadf128.c new file mode 100644 index 0000000000..63e046a269 --- /dev/null +++ b/sysdeps/ieee754/float128/s_setpayloadf128.c @@ -0,0 +1,4 @@ +#include <float128_private.h> +#define SIG 0 +#define FUNC setpayloadf128 +#include "../ldbl-128/s_setpayloadl_main.c" diff --git a/sysdeps/ieee754/float128/s_setpayloadsigf128.c b/sysdeps/ieee754/float128/s_setpayloadsigf128.c new file mode 100644 index 0000000000..85b2c4a1a8 --- /dev/null +++ b/sysdeps/ieee754/float128/s_setpayloadsigf128.c @@ -0,0 +1,4 @@ +#include <float128_private.h> +#define SIG 1 +#define FUNC setpayloadsigf128 +#include "../ldbl-128/s_setpayloadl_main.c" diff --git a/sysdeps/ieee754/float128/s_signbitf128.c b/sysdeps/ieee754/float128/s_signbitf128.c new file mode 100644 index 0000000000..71c1ca3a34 --- /dev/null +++ b/sysdeps/ieee754/float128/s_signbitf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_signbitl.c" diff --git a/sysdeps/ieee754/float128/s_significandf128.c b/sysdeps/ieee754/float128/s_significandf128.c new file mode 100644 index 0000000000..067b724164 --- /dev/null +++ b/sysdeps/ieee754/float128/s_significandf128.c @@ -0,0 +1 @@ +/* Not defined for _FloatN types. */ diff --git a/sysdeps/ieee754/float128/s_sincosf128.c b/sysdeps/ieee754/float128/s_sincosf128.c new file mode 100644 index 0000000000..472adde17f --- /dev/null +++ b/sysdeps/ieee754/float128/s_sincosf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_sincosl.c" diff --git a/sysdeps/ieee754/float128/s_sinf128.c b/sysdeps/ieee754/float128/s_sinf128.c new file mode 100644 index 0000000000..d79a1163a5 --- /dev/null +++ b/sysdeps/ieee754/float128/s_sinf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_sinl.c" diff --git a/sysdeps/ieee754/float128/s_tanf128.c b/sysdeps/ieee754/float128/s_tanf128.c new file mode 100644 index 0000000000..382961aada --- /dev/null +++ b/sysdeps/ieee754/float128/s_tanf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_tanl.c" diff --git a/sysdeps/ieee754/float128/s_tanhf128.c b/sysdeps/ieee754/float128/s_tanhf128.c new file mode 100644 index 0000000000..e02c9a6005 --- /dev/null +++ b/sysdeps/ieee754/float128/s_tanhf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_tanhl.c" diff --git a/sysdeps/ieee754/float128/s_totalorderf128.c b/sysdeps/ieee754/float128/s_totalorderf128.c new file mode 100644 index 0000000000..1b115d8307 --- /dev/null +++ b/sysdeps/ieee754/float128/s_totalorderf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_totalorderl.c" diff --git a/sysdeps/ieee754/float128/s_totalordermagf128.c b/sysdeps/ieee754/float128/s_totalordermagf128.c new file mode 100644 index 0000000000..e44c657275 --- /dev/null +++ b/sysdeps/ieee754/float128/s_totalordermagf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_totalordermagl.c" diff --git a/sysdeps/ieee754/float128/s_truncf128.c b/sysdeps/ieee754/float128/s_truncf128.c new file mode 100644 index 0000000000..474d9dc77e --- /dev/null +++ b/sysdeps/ieee754/float128/s_truncf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/s_truncl.c" diff --git a/sysdeps/ieee754/float128/s_ufromfpf128.c b/sysdeps/ieee754/float128/s_ufromfpf128.c new file mode 100644 index 0000000000..0cd2281035 --- /dev/null +++ b/sysdeps/ieee754/float128/s_ufromfpf128.c @@ -0,0 +1,5 @@ +#define UNSIGNED 1 +#define INEXACT 0 +#define FUNC ufromfpf128 +#include <float128_private.h> +#include "../ldbl-128/s_fromfpl_main.c" diff --git a/sysdeps/ieee754/float128/s_ufromfpxf128.c b/sysdeps/ieee754/float128/s_ufromfpxf128.c new file mode 100644 index 0000000000..c0cd7e3bc8 --- /dev/null +++ b/sysdeps/ieee754/float128/s_ufromfpxf128.c @@ -0,0 +1,5 @@ +#define UNSIGNED 1 +#define INEXACT 1 +#define FUNC ufromfpxf128 +#include <float128_private.h> +#include "../ldbl-128/s_fromfpl_main.c" diff --git a/sysdeps/ieee754/float128/t_sincosf128.c b/sysdeps/ieee754/float128/t_sincosf128.c new file mode 100644 index 0000000000..7e699d3c8e --- /dev/null +++ b/sysdeps/ieee754/float128/t_sincosf128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/t_sincosl.c" diff --git a/sysdeps/ieee754/float128/x2y2m1f128.c b/sysdeps/ieee754/float128/x2y2m1f128.c new file mode 100644 index 0000000000..68880792e6 --- /dev/null +++ b/sysdeps/ieee754/float128/x2y2m1f128.c @@ -0,0 +1,2 @@ +#include <float128_private.h> +#include "../ldbl-128/x2y2m1l.c" diff --git a/sysdeps/ieee754/ldbl-128ibm/bits/iscanonical.h b/sysdeps/ieee754/ldbl-128ibm/bits/iscanonical.h index 198234c9ef..7ddb368d26 100644 --- a/sysdeps/ieee754/ldbl-128ibm/bits/iscanonical.h +++ b/sysdeps/ieee754/ldbl-128ibm/bits/iscanonical.h @@ -27,6 +27,9 @@ extern int __iscanonicall (long double __x) __THROW __attribute__ ((__const__)); # define __iscanonicalf(x) ((void) (__typeof (x)) (x), 1) # define __iscanonical(x) ((void) (__typeof (x)) (x), 1) +# if __HAVE_DISTINCT_FLOAT128 +# define __iscanonicalf128(x) ((void) (__typeof (x)) (x), 1) +# endif /* Return nonzero value if X is canonical. In IEEE interchange binary formats, all values are canonical, but the argument must still be |