annotate lib/signbitl.c @ 17363:5a51fb7777a9

sys_select, sys_time: port 2013-01-30 Solaris 2.6 fix to Cygwin Problem reported by Marco Atzeri in <http://lists.gnu.org/archive/html/bug-gnulib/2013-03/msg00000.html>. * lib/sys_select.in.h [HAVE_SYS_SELECT_H && _CYGWIN_SYS_TIME_H]: Simply delegate to the system <sys/select.h> in this case too. Also, pay attention to _GL_SYS_SELECT_H_REDIRECT_FROM_SYS_TIME_H only if OSF/1, since otherwise Cygwin breaks, and it doesn't seem to be needed on Solaris either. * lib/sys_time.in.h [_CYGWIN_SYS_TIME_H]: Simply delgate to the system <sys/time.h> in this case.
author Paul Eggert <eggert@cs.ucla.edu>
date Tue, 19 Mar 2013 09:08:47 -0700
parents e542fd46ad6f
children 344018b6e5d7
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
1 /* signbit() macro: Determine the sign bit of a floating-point number.
17249
e542fd46ad6f maint: update all copyright year number ranges
Eric Blake <eblake@redhat.com>
parents: 16201
diff changeset
2 Copyright (C) 2007, 2009-2013 Free Software Foundation, Inc.
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
3
9309
bbbbbf4cd1c5 Change copyright notice from GPLv2+ to GPLv3+.
Bruno Haible <bruno@clisp.org>
parents: 8659
diff changeset
4 This program is free software: you can redistribute it and/or modify
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
5 it under the terms of the GNU General Public License as published by
9309
bbbbbf4cd1c5 Change copyright notice from GPLv2+ to GPLv3+.
Bruno Haible <bruno@clisp.org>
parents: 8659
diff changeset
6 the Free Software Foundation; either version 3 of the License, or
bbbbbf4cd1c5 Change copyright notice from GPLv2+ to GPLv3+.
Bruno Haible <bruno@clisp.org>
parents: 8659
diff changeset
7 (at your option) any later version.
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
8
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
9 This program is distributed in the hope that it will be useful,
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
12 GNU General Public License for more details.
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
13
9309
bbbbbf4cd1c5 Change copyright notice from GPLv2+ to GPLv3+.
Bruno Haible <bruno@clisp.org>
parents: 8659
diff changeset
14 You should have received a copy of the GNU General Public License
bbbbbf4cd1c5 Change copyright notice from GPLv2+ to GPLv3+.
Bruno Haible <bruno@clisp.org>
parents: 8659
diff changeset
15 along with this program. If not, see <http://www.gnu.org/licenses/>. */
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
16
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
17 #include <config.h>
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
18
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
19 /* Specification. */
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
20 #include <math.h>
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
21
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
22 #include <string.h>
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
23 #include "isnanl-nolibm.h"
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
24 #include "float+.h"
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
25
11525
0b09102d516e Fix compilation error introduced on 2009-04-25.
Bruno Haible <bruno@clisp.org>
parents: 11506
diff changeset
26 #ifdef gl_signbitl_OPTIMIZED_MACRO
0b09102d516e Fix compilation error introduced on 2009-04-25.
Bruno Haible <bruno@clisp.org>
parents: 11506
diff changeset
27 # undef gl_signbitl
0b09102d516e Fix compilation error introduced on 2009-04-25.
Bruno Haible <bruno@clisp.org>
parents: 11506
diff changeset
28 #endif
0b09102d516e Fix compilation error introduced on 2009-04-25.
Bruno Haible <bruno@clisp.org>
parents: 11506
diff changeset
29
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
30 int
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
31 gl_signbitl (long double arg)
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
32 {
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
33 #if defined LDBL_SIGNBIT_WORD && defined LDBL_SIGNBIT_BIT
8659
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
34 /* The use of a union to extract the bits of the representation of a
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
35 'long double' is safe in practice, despite of the "aliasing rules" of
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
36 C99, because the GCC docs say
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
37 "Even with '-fstrict-aliasing', type-punning is allowed, provided the
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
38 memory is accessed through the union type."
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
39 and similarly for other compilers. */
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
40 # define NWORDS \
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
41 ((sizeof (long double) + sizeof (unsigned int) - 1) / sizeof (unsigned int))
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
42 union { long double value; unsigned int word[NWORDS]; } m;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
43 m.value = arg;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
44 return (m.word[LDBL_SIGNBIT_WORD] >> LDBL_SIGNBIT_BIT) & 1;
8659
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
45 #elif HAVE_COPYSIGNL_IN_LIBC
d6d05661430a Use copysign when implementing signbit, if the libc has it already.
Bruno Haible <bruno@clisp.org>
parents: 8652
diff changeset
46 return copysignl (1.0L, arg) < 0;
8652
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
47 #else
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
48 /* This does not do the right thing for NaN, but this is irrelevant for
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
49 most use cases. */
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
50 if (isnanl (arg))
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
51 return 0;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
52 if (arg < 0.0L)
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
53 return 1;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
54 else if (arg == 0.0L)
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
55 {
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
56 /* Distinguish 0.0L and -0.0L. */
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
57 static long double plus_zero = 0.0L;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
58 long double arg_mem = arg;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
59 return (memcmp (&plus_zero, &arg_mem, SIZEOF_LDBL) != 0);
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
60 }
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
61 else
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
62 return 0;
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
63 #endif
9aabfa2b3200 New module 'signbit'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
64 }