diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/s_setpayload.c | 6 | ||||
| -rw-r--r-- | sysdeps/ieee754/dbl-64/s_setpayload_main.c | 69 | ||||
| -rw-r--r-- | sysdeps/ieee754/dbl-64/wordsize-64/s_setpayload_main.c | 53 |
3 files changed, 128 insertions, 0 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_setpayload.c b/sysdeps/ieee754/dbl-64/s_setpayload.c new file mode 100644 index 0000000000..5ab70dee73 --- /dev/null +++ b/sysdeps/ieee754/dbl-64/s_setpayload.c @@ -0,0 +1,6 @@ +#define SIG 0 +#define FUNC setpayload +#include <s_setpayload_main.c> +#ifdef NO_LONG_DOUBLE +weak_alias (setpayload, setpayloadl) +#endif diff --git a/sysdeps/ieee754/dbl-64/s_setpayload_main.c b/sysdeps/ieee754/dbl-64/s_setpayload_main.c new file mode 100644 index 0000000000..d7a3609c39 --- /dev/null +++ b/sysdeps/ieee754/dbl-64/s_setpayload_main.c @@ -0,0 +1,69 @@ +/* Set NaN payload. dbl-64 version. + Copyright (C) 2016 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/>. */ + +#include <math.h> +#include <math_private.h> +#include <nan-high-order-bit.h> +#include <stdint.h> + +#define SET_HIGH_BIT (HIGH_ORDER_BIT_IS_SET_FOR_SNAN ? SIG : !SIG) +#define BIAS 0x3ff +#define PAYLOAD_DIG 51 +#define EXPLICIT_MANT_DIG 52 + +int +FUNC (double *x, double payload) +{ + uint32_t hx, lx; + EXTRACT_WORDS (hx, lx, payload); + int exponent = hx >> (EXPLICIT_MANT_DIG - 32); + /* Test if argument is (a) negative or too large; (b) too small, + except for 0 when allowed; (c) not an integer. */ + if (exponent >= BIAS + PAYLOAD_DIG + || (exponent < BIAS && !(SET_HIGH_BIT && hx == 0 && lx == 0))) + { + INSERT_WORDS (*x, 0, 0); + return 1; + } + int shift = BIAS + EXPLICIT_MANT_DIG - exponent; + if (shift < 32 + ? (lx & ((1U << shift) - 1)) != 0 + : (lx != 0 || (hx & ((1U << (shift - 32)) - 1)) != 0)) + { + INSERT_WORDS (*x, 0, 0); + return 1; + } + if (exponent != 0) + { + hx &= (1U << (EXPLICIT_MANT_DIG - 32)) - 1; + hx |= 1U << (EXPLICIT_MANT_DIG - 32); + if (shift >= 32) + { + lx = hx >> (shift - 32); + hx = 0; + } + else if (shift != 0) + { + lx = (lx >> shift) | (hx << (32 - shift)); + hx >>= shift; + } + } + hx |= 0x7ff00000 | (SET_HIGH_BIT ? 0x80000 : 0); + INSERT_WORDS (*x, hx, lx); + return 0; +} diff --git a/sysdeps/ieee754/dbl-64/wordsize-64/s_setpayload_main.c b/sysdeps/ieee754/dbl-64/wordsize-64/s_setpayload_main.c new file mode 100644 index 0000000000..691579295c --- /dev/null +++ b/sysdeps/ieee754/dbl-64/wordsize-64/s_setpayload_main.c @@ -0,0 +1,53 @@ +/* Set NaN payload. dbl-64/wordsize-64 version. + Copyright (C) 2016 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/>. */ + +#include <math.h> +#include <math_private.h> +#include <nan-high-order-bit.h> +#include <stdint.h> + +#define SET_HIGH_BIT (HIGH_ORDER_BIT_IS_SET_FOR_SNAN ? SIG : !SIG) +#define BIAS 0x3ff +#define PAYLOAD_DIG 51 +#define EXPLICIT_MANT_DIG 52 + +int +FUNC (double *x, double payload) +{ + uint64_t ix; + EXTRACT_WORDS64 (ix, payload); + int exponent = ix >> EXPLICIT_MANT_DIG; + /* Test if argument is (a) negative or too large; (b) too small, + except for 0 when allowed; (c) not an integer. */ + if (exponent >= BIAS + PAYLOAD_DIG + || (exponent < BIAS && !(SET_HIGH_BIT && ix == 0)) + || (ix & ((1ULL << (BIAS + EXPLICIT_MANT_DIG - exponent)) - 1)) != 0) + { + INSERT_WORDS64 (*x, 0); + return 1; + } + if (ix != 0) + { + ix &= (1ULL << EXPLICIT_MANT_DIG) - 1; + ix |= 1ULL << EXPLICIT_MANT_DIG; + ix >>= BIAS + EXPLICIT_MANT_DIG - exponent; + } + ix |= 0x7ff0000000000000ULL | (SET_HIGH_BIT ? 0x8000000000000ULL : 0); + INSERT_WORDS64 (*x, ix); + return 0; +} |