From 302bb1a3c327d546b19447d7304bb2b27e9a34de Mon Sep 17 00:00:00 2001 From: "Paul E. Murphy" Date: Mon, 8 Aug 2016 15:49:31 -0500 Subject: float128: Add wrappers to override ldbl-128 as float128. This change defines float128_private.h which contains macros used to override long double naming conventions when building a ldbl file. * 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. --- sysdeps/ieee754/float128/ieee754_float128.h | 140 ++++++++++++++++++++++++++++ 1 file changed, 140 insertions(+) create mode 100644 sysdeps/ieee754/float128/ieee754_float128.h (limited to 'sysdeps/ieee754/float128/ieee754_float128.h') 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 + . */ +#ifndef _IEEE754_FLOAT128_H +#define _IEEE754_FLOAT128_H + +#include +#include + +# 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 -- cgit 1.4.1