From b7f88156174f9bbf9c9393651b2ad0cdd3893e83 Mon Sep 17 00:00:00 2001 From: Siddhesh Poyarekar Date: Wed, 23 Dec 2020 13:53:34 +0530 Subject: x86 long double: Support pseudo numbers in fpclassifyl Also move sysdeps/i386/fpu/s_fpclassifyl.c to sysdeps/x86/fpu/s_fpclassifyl.c and remove sysdeps/x86_64/fpu/s_fpclassifyl.c Reviewed-by: Adhemerval Zanella --- sysdeps/x86/fpu/s_fpclassifyl.c | 46 +++++++++++++++++++++++++++++++++++++++++ 1 file changed, 46 insertions(+) create mode 100644 sysdeps/x86/fpu/s_fpclassifyl.c (limited to 'sysdeps/x86/fpu') diff --git a/sysdeps/x86/fpu/s_fpclassifyl.c b/sysdeps/x86/fpu/s_fpclassifyl.c new file mode 100644 index 0000000000..a48b7685c4 --- /dev/null +++ b/sysdeps/x86/fpu/s_fpclassifyl.c @@ -0,0 +1,46 @@ +/* Return classification value corresponding to argument. + Copyright (C) 1997-2020 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Ulrich Drepper , 1997. + + 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 + . */ + +#include + +#include + + +int +__fpclassifyl (long double x) +{ + uint32_t ex, hx, lx; + int retval = FP_NORMAL; + + GET_LDOUBLE_WORDS (ex, hx, lx, x); + ex &= 0x7fff; + if ((ex | lx | hx) == 0) + retval = FP_ZERO; + else if (ex == 0 && (hx & 0x80000000) == 0) + retval = FP_SUBNORMAL; + /* Pseudo-normals, i.e. pseudo-zero, pseudo-infinity and un-normals. They + behave like NaNs, so categorize them as such. */ + else if ((hx & 0x80000000) == 0) + retval = FP_NAN; + else if (ex == 0x7fff) + retval = ((hx & 0x7fffffff) | lx) != 0 ? FP_NAN : FP_INFINITE; + + return retval; +} +libm_hidden_def (__fpclassifyl) -- cgit 1.4.1