1 /* signbit() macro: Determine the sign bit of a floating-point number.
2 Copyright (C) 2007 Free Software Foundation, Inc.
4 This program is free software; you can redistribute it and/or modify
5 it under the terms of the GNU General Public License as published by
6 the Free Software Foundation; either version 2, or (at your option)
9 This program is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 GNU General Public License for more details.
14 You should have received a copy of the GNU General Public License along
15 with this program; if not, write to the Free Software Foundation,
16 Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */
24 #include "isnanl-nolibm.h"
30 gl_signbitl (long double arg)
32 #if defined LDBL_SIGNBIT_WORD && defined LDBL_SIGNBIT_BIT
34 ((sizeof (long double) + sizeof (unsigned int) - 1) / sizeof (unsigned int))
35 union { long double value; unsigned int word[NWORDS]; } m;
37 return (m.word[LDBL_SIGNBIT_WORD] >> LDBL_SIGNBIT_BIT) & 1;
39 /* This does not do the right thing for NaN, but this is irrelevant for
47 /* Distinguish 0.0L and -0.0L. */
48 static long double plus_zero = 0.0L;
49 long double arg_mem = arg;
50 return (memcmp (&plus_zero, &arg_mem, SIZEOF_LDBL) != 0);