1 /* A GNU-like <math.h>.
3 Copyright (C) 2002-2003, 2007 Free Software Foundation, Inc.
5 This program is free software; you can redistribute it and/or modify
6 it under the terms of the GNU General Public License as published by
7 the Free Software Foundation; either version 2, or (at your option)
10 This program is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 GNU General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software Foundation,
17 Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */
19 #if defined __DECC && __DECC_VER >= 60000000
20 # include_next <math.h>
26 #if !(defined __DECC && __DECC_VER >= 60000000)
27 # include @ABSOLUTE_MATH_H@
31 /* The definition of GL_LINK_WARNING is copied here. */
42 If x finite and nonzero: 0.5 <= |mantissa| < 1.0.
43 If x is zero: mantissa = x, exp = 0.
44 If x is infinite or NaN: mantissa = x, exp unspecified.
45 Store exp and return mantissa. */
48 # define frexp rpl_frexp
49 extern double frexp (double x, int *exp);
51 #elif defined GNULIB_POSIXCHECK
54 (GL_LINK_WARNING ("frexp is unportable - " \
55 "use gnulib module frexp for portability"), \
60 #if @GNULIB_MATHL@ || !@HAVE_DECL_ACOSL@
61 extern long double acosl (long double x);
63 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
66 (GL_LINK_WARNING ("acosl is unportable - " \
67 "use gnulib module mathl for portability"), \
71 #if @GNULIB_MATHL@ || !@HAVE_DECL_ASINL@
72 extern long double asinl (long double x);
74 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
77 (GL_LINK_WARNING ("asinl is unportable - " \
78 "use gnulib module mathl for portability"), \
82 #if @GNULIB_MATHL@ || !@HAVE_DECL_ATANL@
83 extern long double atanl (long double x);
85 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
88 (GL_LINK_WARNING ("atanl is unportable - " \
89 "use gnulib module mathl for portability"), \
93 #if @GNULIB_MATHL@ || !@HAVE_DECL_CEILL@
94 extern long double ceill (long double x);
96 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
99 (GL_LINK_WARNING ("ceill is unportable - " \
100 "use gnulib module mathl for portability"), \
104 #if @GNULIB_MATHL@ || !@HAVE_DECL_COSL@
105 extern long double cosl (long double x);
107 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
110 (GL_LINK_WARNING ("cosl is unportable - " \
111 "use gnulib module mathl for portability"), \
115 #if @GNULIB_MATHL@ || !@HAVE_DECL_EXPL@
116 extern long double expl (long double x);
118 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
121 (GL_LINK_WARNING ("expl is unportable - " \
122 "use gnulib module mathl for portability"), \
126 #if @GNULIB_MATHL@ || !@HAVE_DECL_FLOORL@
127 extern long double floorl (long double x);
129 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
132 (GL_LINK_WARNING ("floorl is unportable - " \
133 "use gnulib module mathl for portability"), \
140 If x finite and nonzero: 0.5 <= |mantissa| < 1.0.
141 If x is zero: mantissa = x, exp = 0.
142 If x is infinite or NaN: mantissa = x, exp unspecified.
143 Store exp and return mantissa. */
144 #if @GNULIB_FREXPL@ && @REPLACE_FREXPL@
145 # define frexpl rpl_frexpl
147 #if (@GNULIB_FREXPL@ && @REPLACE_FREXPL@) || !@HAVE_DECL_FREXPL@
148 extern long double frexpl (long double x, int *exp);
150 #if !@GNULIB_FREXPL@ && defined GNULIB_POSIXCHECK
152 # define frexpl(x,e) \
153 (GL_LINK_WARNING ("frexpl is unportable - " \
154 "use gnulib module frexpl for portability"), \
158 /* Return x * 2^exp. */
159 #if @GNULIB_LDEXPL@ && @REPLACE_LDEXPL@
160 # define ldexpl rpl_ldexpl
162 #if (@GNULIB_LDEXPL@ && @REPLACE_LDEXPL@) || !@HAVE_DECL_LDEXPL@
163 extern long double ldexpl (long double x, int exp);
165 #if !@GNULIB_LDEXPL@ && defined GNULIB_POSIXCHECK
167 # define ldexpl(x,e) \
168 (GL_LINK_WARNING ("ldexpl is unportable - " \
169 "use gnulib module ldexpl for portability"), \
173 #if @GNULIB_MATHL@ || !@HAVE_DECL_LOGL@
174 extern long double logl (long double x);
176 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
179 (GL_LINK_WARNING ("logl is unportable - " \
180 "use gnulib module mathl for portability"), \
184 #if @GNULIB_MATHL@ || !@HAVE_DECL_SINL@
185 extern long double sinl (long double x);
187 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
190 (GL_LINK_WARNING ("sinl is unportable - " \
191 "use gnulib module mathl for portability"), \
195 #if @GNULIB_MATHL@ || !@HAVE_DECL_SQRTL@
196 extern long double sqrtl (long double x);
198 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
201 (GL_LINK_WARNING ("sqrtl is unportable - " \
202 "use gnulib module mathl for portability"), \
206 #if @GNULIB_MATHL@ || !@HAVE_DECL_TANL@
207 extern long double tanl (long double x);
209 #if !@GNULIB_MATHL@ && defined GNULIB_POSIXCHECK
212 (GL_LINK_WARNING ("tanl is unportable - " \
213 "use gnulib module mathl for portability"), \
219 # if @REPLACE_SIGNBIT@
221 extern int gl_signbitf (float arg);
222 extern int gl_signbitd (double arg);
223 extern int gl_signbitl (long double arg);
224 # if __GNUC__ >= 2 && !__STRICT_ANSI__
225 # if defined FLT_SIGNBIT_WORD && defined FLT_SIGNBIT_BIT
226 # define gl_signbitf(arg) \
227 ({ union { float _value; \
228 unsigned int _word[(sizeof (float) + sizeof (unsigned int) - 1) / sizeof (unsigned int)]; \
231 (_m._word[FLT_SIGNBIT_WORD] >> FLT_SIGNBIT_BIT) & 1; \
234 # if defined DBL_SIGNBIT_WORD && defined DBL_SIGNBIT_BIT
235 # define gl_signbitd(arg) \
236 ({ union { double _value; \
237 unsigned int _word[(sizeof (double) + sizeof (unsigned int) - 1) / sizeof (unsigned int)]; \
240 (_m._word[DBL_SIGNBIT_WORD] >> DBL_SIGNBIT_BIT) & 1; \
243 # if defined LDBL_SIGNBIT_WORD && defined LDBL_SIGNBIT_BIT
244 # define gl_signbitl(arg) \
245 ({ union { long double _value; \
246 unsigned int _word[(sizeof (long double) + sizeof (unsigned int) - 1) / sizeof (unsigned int)]; \
249 (_m._word[LDBL_SIGNBIT_WORD] >> LDBL_SIGNBIT_BIT) & 1; \
253 # define signbit(x) \
254 (sizeof (x) == sizeof (long double) ? gl_signbitl (x) : \
255 sizeof (x) == sizeof (double) ? gl_signbitd (x) : \
258 #elif defined GNULIB_POSIXCHECK
259 /* How to override a macro? */