diff options
author | Samuel Thibault <samuel.thibault@ens-lyon.org> | 2018-12-27 17:58:54 +0000 |
---|---|---|
committer | Samuel Thibault <samuel.thibault@ens-lyon.org> | 2018-12-27 17:58:54 +0000 |
commit | 1ba06384fb715fb7732b20911ffc56a61d16cb0d (patch) | |
tree | 0fcfaf14f40169b289044495b1984a5fa3e0e965 /sysdeps/ieee754/flt-32/s_setpayloadf_main.c | |
parent | 1af57d5793424a5a6ec542d0a5ef0b7c4f312623 (diff) | |
parent | 963c37d5c0eb62b38f8764b23931c0dcdd497a13 (diff) |
Merge commit 'refs/top-bases/t/sysvshm' into t/sysvshm
Diffstat (limited to 'sysdeps/ieee754/flt-32/s_setpayloadf_main.c')
-rw-r--r-- | sysdeps/ieee754/flt-32/s_setpayloadf_main.c | 54 |
1 files changed, 54 insertions, 0 deletions
diff --git a/sysdeps/ieee754/flt-32/s_setpayloadf_main.c b/sysdeps/ieee754/flt-32/s_setpayloadf_main.c new file mode 100644 index 0000000000..a01c8677a3 --- /dev/null +++ b/sysdeps/ieee754/flt-32/s_setpayloadf_main.c @@ -0,0 +1,54 @@ +/* Set NaN payload. flt-32 version. + Copyright (C) 2016-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/>. */ + +#include <math.h> +#include <math_private.h> +#include <libm-alias-float.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 0x7f +#define PAYLOAD_DIG 22 +#define EXPLICIT_MANT_DIG 23 + +int +FUNC (float *x, float payload) +{ + uint32_t ix; + GET_FLOAT_WORD (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 & ((1U << (BIAS + EXPLICIT_MANT_DIG - exponent)) - 1)) != 0) + { + SET_FLOAT_WORD (*x, 0); + return 1; + } + if (ix != 0) + { + ix &= (1U << EXPLICIT_MANT_DIG) - 1; + ix |= 1U << EXPLICIT_MANT_DIG; + ix >>= BIAS + EXPLICIT_MANT_DIG - exponent; + } + ix |= 0x7f800000 | (SET_HIGH_BIT ? 0x400000 : 0); + SET_FLOAT_WORD (*x, ix); + return 0; +} |