From 3d165855f2de19f3f5dc51d193dbe9431121a5c5 Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Mon, 4 Dec 2000 01:26:09 +0000 Subject: Update. * sysdeps/i386/fpu/fpclassify.c: New file. --- ChangeLog | 2 ++ sysdeps/i386/fpu/s_fpclassifyl.c | 42 ++++++++++++++++++++++++++++++++++++++++ 2 files changed, 44 insertions(+) create mode 100644 sysdeps/i386/fpu/s_fpclassifyl.c diff --git a/ChangeLog b/ChangeLog index 92809c5..3070399 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,5 +1,7 @@ 2000-12-03 Ulrich Drepper + * sysdeps/i386/fpu/fpclassify.c: New file. + * sysdeps/i386/ldbl2mpn.c (__mpn_extract_long_double): Handle the pseudo denormal representation of LDBL_MIN. diff --git a/sysdeps/i386/fpu/s_fpclassifyl.c b/sysdeps/i386/fpu/s_fpclassifyl.c new file mode 100644 index 0000000..61481f0 --- /dev/null +++ b/sysdeps/i386/fpu/s_fpclassifyl.c @@ -0,0 +1,42 @@ +/* Return classification value corresponding to argument. + Copyright (C) 1997, 2000 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 Library General Public License as + published by the Free Software Foundation; either version 2 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 + Library General Public License for more details. + + You should have received a copy of the GNU Library General Public + License along with the GNU C Library; see the file COPYING.LIB. If not, + write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, + Boston, MA 02111-1307, USA. */ + +#include + +#include "math_private.h" + + +int +__fpclassifyl (long double x) +{ + u_int32_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; + else if (ex == 0x7fff) + retval = ((hx & 0x7fffffff) | lx) != 0 ? FP_NAN : FP_INFINITE; + + return retval; +} -- cgit v1.1