Build dlopenable modules with `-module'.
[bpt/guile.git] / lib / math.in.h
CommitLineData
2e65b52f
LC
1/* A GNU-like <math.h>.
2
3 Copyright (C) 2002-2003, 2007-2010 Free Software Foundation, Inc.
4
5 This program is free software: you can redistribute it and/or modify
6 it under the terms of the GNU Lesser General Public License as published by
7 the Free Software Foundation; either version 3 of the License, or
8 (at your option) any later version.
9
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 Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public License
16 along with this program. If not, see <http://www.gnu.org/licenses/>. */
17
18#ifndef _GL_MATH_H
19
20#if __GNUC__ >= 3
21@PRAGMA_SYSTEM_HEADER@
22#endif
23
24/* The include_next requires a split double-inclusion guard. */
25#@INCLUDE_NEXT_AS_FIRST_DIRECTIVE@ @NEXT_AS_FIRST_DIRECTIVE_MATH_H@
26
27#ifndef _GL_MATH_H
28#define _GL_MATH_H
29
30
31/* The definitions of _GL_FUNCDECL_RPL etc. are copied here. */
32
33/* The definition of _GL_ARG_NONNULL is copied here. */
34
35/* The definition of _GL_WARN_ON_USE is copied here. */
36
37/* Helper macros to define a portability warning for the
38 classification macro FUNC called with VALUE. POSIX declares the
39 classification macros with an argument of real-floating (that is,
40 one of float, double, or long double). */
41#define _GL_WARN_REAL_FLOATING_DECL(func) \
42static inline int \
43rpl_ ## func ## f (float f) \
44{ \
45 return func (f); \
46} \
47static inline int \
48rpl_ ## func ## d (double d) \
49{ \
50 return func (d); \
51} \
52static inline int \
53rpl_ ## func ## l (long double l) \
54{ \
55 return func (l); \
56} \
57_GL_WARN_ON_USE (rpl_ ## func ## f, #func " is unportable - " \
58 "use gnulib module " #func " for portability"); \
59_GL_WARN_ON_USE (rpl_ ## func ## d, #func " is unportable - " \
60 "use gnulib module " #func " for portability"); \
61_GL_WARN_ON_USE (rpl_ ## func ## l, #func " is unportable - " \
62 "use gnulib module " #func " for portability")
63#define _GL_WARN_REAL_FLOATING_IMPL(func, value) \
64 (sizeof (value) == sizeof (float) ? rpl_ ## func ## f (value) \
65 : sizeof (value) == sizeof (double) ? rpl_ ## func ## d (value) \
66 : rpl_ ## func ## l (value))
67
68
69/* POSIX allows platforms that don't support NAN. But all major
70 machines in the past 15 years have supported something close to
71 IEEE NaN, so we define this unconditionally. We also must define
72 it on platforms like Solaris 10, where NAN is present but defined
73 as a function pointer rather than a floating point constant. */
74#if !defined NAN || @REPLACE_NAN@
75# undef NAN
76 /* The Compaq (ex-DEC) C 6.4 compiler chokes on the expression 0.0 / 0.0. */
77# ifdef __DECC
78static float
79_NaN ()
80{
81 static float zero = 0.0f;
82 return zero / zero;
83}
84# define NAN (_NaN())
85# else
86# define NAN (0.0f / 0.0f)
87# endif
88#endif
89
90/* Solaris 10 defines HUGE_VAL, but as a function pointer rather
91 than a floating point constant. */
92#if @REPLACE_HUGE_VAL@
93# undef HUGE_VAL
94# define HUGE_VAL (1.0 / 0.0)
95#endif
96
97
98/* Write x as
99 x = mantissa * 2^exp
100 where
101 If x finite and nonzero: 0.5 <= |mantissa| < 1.0.
102 If x is zero: mantissa = x, exp = 0.
103 If x is infinite or NaN: mantissa = x, exp unspecified.
104 Store exp in *EXPPTR and return mantissa. */
105#if @GNULIB_FREXP@
106# if @REPLACE_FREXP@
107# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
108# define frexp rpl_frexp
109# endif
110_GL_FUNCDECL_RPL (frexp, double, (double x, int *expptr) _GL_ARG_NONNULL ((2)));
111_GL_CXXALIAS_RPL (frexp, double, (double x, int *expptr));
112# else
113_GL_CXXALIAS_SYS (frexp, double, (double x, int *expptr));
114# endif
115_GL_CXXALIASWARN (frexp);
116#elif defined GNULIB_POSIXCHECK
117# undef frexp
118/* Assume frexp is always declared. */
119_GL_WARN_ON_USE (frexp, "frexp is unportable - "
120 "use gnulib module frexp for portability");
121#endif
122
123
124#if @GNULIB_LOGB@
125# if !@HAVE_DECL_LOGB@
126_GL_EXTERN_C double logb (double x);
127# endif
128#elif defined GNULIB_POSIXCHECK
129# undef logb
130# if HAVE_RAW_DECL_LOGB
131_GL_WARN_ON_USE (logb, "logb is unportable - "
132 "use gnulib module logb for portability");
133# endif
134#endif
135
136
137#if @GNULIB_ACOSL@
138# if !@HAVE_ACOSL@ || !@HAVE_DECL_ACOSL@
139_GL_FUNCDECL_SYS (acosl, long double, (long double x));
140# endif
141_GL_CXXALIAS_SYS (acosl, long double, (long double x));
142_GL_CXXALIASWARN (acosl);
143#elif defined GNULIB_POSIXCHECK
144# undef acosl
145# if HAVE_RAW_DECL_ACOSL
146_GL_WARN_ON_USE (acosl, "acosl is unportable - "
147 "use gnulib module mathl for portability");
148# endif
149#endif
150
151
152#if @GNULIB_ASINL@
153# if !@HAVE_ASINL@ || !@HAVE_DECL_ASINL@
154_GL_FUNCDECL_SYS (asinl, long double, (long double x));
155# endif
156_GL_CXXALIAS_SYS (asinl, long double, (long double x));
157_GL_CXXALIASWARN (asinl);
158#elif defined GNULIB_POSIXCHECK
159# undef asinl
160# if HAVE_RAW_DECL_ASINL
161_GL_WARN_ON_USE (asinl, "asinl is unportable - "
162 "use gnulib module mathl for portability");
163# endif
164#endif
165
166
167#if @GNULIB_ATANL@
168# if !@HAVE_ATANL@ || !@HAVE_DECL_ATANL@
169_GL_FUNCDECL_SYS (atanl, long double, (long double x));
170# endif
171_GL_CXXALIAS_SYS (atanl, long double, (long double x));
172_GL_CXXALIASWARN (atanl);
173#elif defined GNULIB_POSIXCHECK
174# undef atanl
175# if HAVE_RAW_DECL_ATANL
176_GL_WARN_ON_USE (atanl, "atanl is unportable - "
177 "use gnulib module mathl for portability");
178# endif
179#endif
180
181
182#if @GNULIB_CEILF@
183# if @REPLACE_CEILF@
184# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
185# define ceilf rpl_ceilf
186# endif
187_GL_FUNCDECL_RPL (ceilf, float, (float x));
188_GL_CXXALIAS_RPL (ceilf, float, (float x));
189# else
190# if !@HAVE_DECL_CEILF@
191_GL_FUNCDECL_SYS (ceilf, float, (float x));
192# endif
193_GL_CXXALIAS_SYS (ceilf, float, (float x));
194# endif
195_GL_CXXALIASWARN (ceilf);
196#elif defined GNULIB_POSIXCHECK
197# undef ceilf
198# if HAVE_RAW_DECL_CEILF
199_GL_WARN_ON_USE (ceilf, "ceilf is unportable - "
200 "use gnulib module ceilf for portability");
201# endif
202#endif
203
204#if @GNULIB_CEILL@
205# if @REPLACE_CEILL@
206# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
207# define ceill rpl_ceill
208# endif
209_GL_FUNCDECL_RPL (ceill, long double, (long double x));
210_GL_CXXALIAS_RPL (ceill, long double, (long double x));
211# else
212# if !@HAVE_DECL_CEILL@
213_GL_FUNCDECL_SYS (ceill, long double, (long double x));
214# endif
215_GL_CXXALIAS_SYS (ceill, long double, (long double x));
216# endif
217_GL_CXXALIASWARN (ceill);
218#elif defined GNULIB_POSIXCHECK
219# undef ceill
220# if HAVE_RAW_DECL_CEILL
221_GL_WARN_ON_USE (ceill, "ceill is unportable - "
222 "use gnulib module ceill for portability");
223# endif
224#endif
225
226
227#if @GNULIB_COSL@
228# if !@HAVE_COSL@ || !@HAVE_DECL_COSL@
229_GL_FUNCDECL_SYS (cosl, long double, (long double x));
230# endif
231_GL_CXXALIAS_SYS (cosl, long double, (long double x));
232_GL_CXXALIASWARN (cosl);
233#elif defined GNULIB_POSIXCHECK
234# undef cosl
235# if HAVE_RAW_DECL_COSL
236_GL_WARN_ON_USE (cosl, "cosl is unportable - "
237 "use gnulib module mathl for portability");
238# endif
239#endif
240
241
242#if @GNULIB_EXPL@
243# if !@HAVE_EXPL@ || !@HAVE_DECL_EXPL@
244_GL_FUNCDECL_SYS (expl, long double, (long double x));
245# endif
246_GL_CXXALIAS_SYS (expl, long double, (long double x));
247_GL_CXXALIASWARN (expl);
248#elif defined GNULIB_POSIXCHECK
249# undef expl
250# if HAVE_RAW_DECL_EXPL
251_GL_WARN_ON_USE (expl, "expl is unportable - "
252 "use gnulib module mathl for portability");
253# endif
254#endif
255
256
257#if @GNULIB_FLOORF@
258# if @REPLACE_FLOORF@
259# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
260# define floorf rpl_floorf
261# endif
262_GL_FUNCDECL_RPL (floorf, float, (float x));
263_GL_CXXALIAS_RPL (floorf, float, (float x));
264#else
265# if !@HAVE_DECL_FLOORF@
266_GL_FUNCDECL_SYS (floorf, float, (float x));
267# endif
268_GL_CXXALIAS_SYS (floorf, float, (float x));
269# endif
270_GL_CXXALIASWARN (floorf);
271#elif defined GNULIB_POSIXCHECK
272# undef floorf
273# if HAVE_RAW_DECL_FLOORF
274_GL_WARN_ON_USE (floorf, "floorf is unportable - "
275 "use gnulib module floorf for portability");
276# endif
277#endif
278
279#if @GNULIB_FLOORL@
280# if @REPLACE_FLOORL@
281# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
282# define floorl rpl_floorl
283# endif
284_GL_FUNCDECL_RPL (floorl, long double, (long double x));
285_GL_CXXALIAS_RPL (floorl, long double, (long double x));
286# else
287# if !@HAVE_DECL_FLOORL@
288_GL_FUNCDECL_SYS (floorl, long double, (long double x));
289# endif
290_GL_CXXALIAS_SYS (floorl, long double, (long double x));
291# endif
292_GL_CXXALIASWARN (floorl);
293#elif defined GNULIB_POSIXCHECK
294# undef floorl
295# if HAVE_RAW_DECL_FLOORL
296_GL_WARN_ON_USE (floorl, "floorl is unportable - "
297 "use gnulib module floorl for portability");
298# endif
299#endif
300
301
302/* Write x as
303 x = mantissa * 2^exp
304 where
305 If x finite and nonzero: 0.5 <= |mantissa| < 1.0.
306 If x is zero: mantissa = x, exp = 0.
307 If x is infinite or NaN: mantissa = x, exp unspecified.
308 Store exp in *EXPPTR and return mantissa. */
309#if @GNULIB_FREXPL@ && @REPLACE_FREXPL@
310# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
311# define frexpl rpl_frexpl
312# endif
313_GL_FUNCDECL_RPL (frexpl, long double,
314 (long double x, int *expptr) _GL_ARG_NONNULL ((2)));
315_GL_CXXALIAS_RPL (frexpl, long double, (long double x, int *expptr));
316#else
317# if !@HAVE_DECL_FREXPL@
318_GL_FUNCDECL_SYS (frexpl, long double,
319 (long double x, int *expptr) _GL_ARG_NONNULL ((2)));
320# endif
321# if @GNULIB_FREXPL@
322_GL_CXXALIAS_SYS (frexpl, long double, (long double x, int *expptr));
323# endif
324#endif
325#if @GNULIB_FREXPL@ && !(@REPLACE_FREXPL@ && !@HAVE_DECL_FREXPL@)
326_GL_CXXALIASWARN (frexpl);
327#endif
328#if !@GNULIB_FREXPL@ && defined GNULIB_POSIXCHECK
329# undef frexpl
330# if HAVE_RAW_DECL_FREXPL
331_GL_WARN_ON_USE (frexpl, "frexpl is unportable - "
332 "use gnulib module frexpl for portability");
333# endif
334#endif
335
336
337/* Return x * 2^exp. */
338#if @GNULIB_LDEXPL@ && @REPLACE_LDEXPL@
339# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
340# define ldexpl rpl_ldexpl
341# endif
342_GL_FUNCDECL_RPL (ldexpl, long double, (long double x, int exp));
343_GL_CXXALIAS_RPL (ldexpl, long double, (long double x, int exp));
344#else
345# if !@HAVE_DECL_LDEXPL@
346_GL_FUNCDECL_SYS (ldexpl, long double, (long double x, int exp));
347# endif
348# if @GNULIB_LDEXPL@
349_GL_CXXALIAS_SYS (ldexpl, long double, (long double x, int exp));
350# endif
351#endif
352#if @GNULIB_LDEXPL@
353_GL_CXXALIASWARN (ldexpl);
354#endif
355#if !@GNULIB_LDEXPL@ && defined GNULIB_POSIXCHECK
356# undef ldexpl
357# if HAVE_RAW_DECL_LDEXPL
358_GL_WARN_ON_USE (ldexpl, "ldexpl is unportable - "
359 "use gnulib module ldexpl for portability");
360# endif
361#endif
362
363
364#if @GNULIB_LOGL@
365# if !@HAVE_LOGL@ || !@HAVE_DECL_LOGL@
366_GL_FUNCDECL_SYS (logl, long double, (long double x));
367# endif
368_GL_CXXALIAS_SYS (logl, long double, (long double x));
369_GL_CXXALIASWARN (logl);
370#elif defined GNULIB_POSIXCHECK
371# undef logl
372# if HAVE_RAW_DECL_LOGL
373_GL_WARN_ON_USE (logl, "logl is unportable - "
374 "use gnulib module mathl for portability");
375# endif
376#endif
377
378
379#if @GNULIB_ROUNDF@
380# if @REPLACE_ROUNDF@
381# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
382# undef roundf
383# define roundf rpl_roundf
384# endif
385_GL_FUNCDECL_RPL (roundf, float, (float x));
386_GL_CXXALIAS_RPL (roundf, float, (float x));
387# else
388# if !@HAVE_DECL_ROUNDF@
389_GL_FUNCDECL_SYS (roundf, float, (float x));
390# endif
391_GL_CXXALIAS_SYS (roundf, float, (float x));
392# endif
393_GL_CXXALIASWARN (roundf);
394#elif defined GNULIB_POSIXCHECK
395# undef roundf
396# if HAVE_RAW_DECL_ROUNDF
397_GL_WARN_ON_USE (roundf, "roundf is unportable - "
398 "use gnulib module roundf for portability");
399# endif
400#endif
401
402#if @GNULIB_ROUND@
403# if @REPLACE_ROUND@
404# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
405# undef round
406# define round rpl_round
407# endif
408_GL_FUNCDECL_RPL (round, double, (double x));
409_GL_CXXALIAS_RPL (round, double, (double x));
410# else
411# if !@HAVE_DECL_ROUND@
412_GL_FUNCDECL_SYS (round, double, (double x));
413# endif
414_GL_CXXALIAS_SYS (round, double, (double x));
415# endif
416_GL_CXXALIASWARN (round);
417#elif defined GNULIB_POSIXCHECK
418# undef round
419# if HAVE_RAW_DECL_ROUND
420_GL_WARN_ON_USE (round, "round is unportable - "
421 "use gnulib module round for portability");
422# endif
423#endif
424
425#if @GNULIB_ROUNDL@
426# if @REPLACE_ROUNDL@
427# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
428# undef roundl
429# define roundl rpl_roundl
430# endif
431_GL_FUNCDECL_RPL (roundl, long double, (long double x));
432_GL_CXXALIAS_RPL (roundl, long double, (long double x));
433# else
434# if !@HAVE_DECL_ROUNDL@
435_GL_FUNCDECL_SYS (roundl, long double, (long double x));
436# endif
437_GL_CXXALIAS_SYS (roundl, long double, (long double x));
438# endif
439_GL_CXXALIASWARN (roundl);
440#elif defined GNULIB_POSIXCHECK
441# undef roundl
442# if HAVE_RAW_DECL_ROUNDL
443_GL_WARN_ON_USE (roundl, "roundl is unportable - "
444 "use gnulib module roundl for portability");
445# endif
446#endif
447
448
449#if @GNULIB_SINL@
450# if !@HAVE_SINL@ || !@HAVE_DECL_SINL@
451_GL_FUNCDECL_SYS (sinl, long double, (long double x));
452# endif
453_GL_CXXALIAS_SYS (sinl, long double, (long double x));
454_GL_CXXALIASWARN (sinl);
455#elif defined GNULIB_POSIXCHECK
456# undef sinl
457# if HAVE_RAW_DECL_SINL
458_GL_WARN_ON_USE (sinl, "sinl is unportable - "
459 "use gnulib module mathl for portability");
460# endif
461#endif
462
463
464#if @GNULIB_SQRTL@
465# if !@HAVE_SQRTL@ || !@HAVE_DECL_SQRTL@
466_GL_FUNCDECL_SYS (sqrtl, long double, (long double x));
467# endif
468_GL_CXXALIAS_SYS (sqrtl, long double, (long double x));
469_GL_CXXALIASWARN (sqrtl);
470#elif defined GNULIB_POSIXCHECK
471# undef sqrtl
472# if HAVE_RAW_DECL_SQRTL
473_GL_WARN_ON_USE (sqrtl, "sqrtl is unportable - "
474 "use gnulib module mathl for portability");
475# endif
476#endif
477
478
479#if @GNULIB_TANL@
480# if !@HAVE_TANL@ || !@HAVE_DECL_TANL@
481_GL_FUNCDECL_SYS (tanl, long double, (long double x));
482# endif
483_GL_CXXALIAS_SYS (tanl, long double, (long double x));
484_GL_CXXALIASWARN (tanl);
485#elif defined GNULIB_POSIXCHECK
486# undef tanl
487# if HAVE_RAW_DECL_TANL
488_GL_WARN_ON_USE (tanl, "tanl is unportable - "
489 "use gnulib module mathl for portability");
490# endif
491#endif
492
493
494#if @GNULIB_TRUNCF@
495# if !@HAVE_DECL_TRUNCF@
496_GL_FUNCDECL_SYS (truncf, float, (float x));
497# endif
498_GL_CXXALIAS_SYS (truncf, float, (float x));
499_GL_CXXALIASWARN (truncf);
500#elif defined GNULIB_POSIXCHECK
501# undef truncf
502# if HAVE_RAW_DECL_TRUNCF
503_GL_WARN_ON_USE (truncf, "truncf is unportable - "
504 "use gnulib module truncf for portability");
505# endif
506#endif
507
508#if @GNULIB_TRUNC@
509# if !@HAVE_DECL_TRUNC@
510_GL_FUNCDECL_SYS (trunc, double, (double x));
511# endif
512_GL_CXXALIAS_SYS (trunc, double, (double x));
513_GL_CXXALIASWARN (trunc);
514#elif defined GNULIB_POSIXCHECK
515# undef trunc
516# if HAVE_RAW_DECL_TRUNC
517_GL_WARN_ON_USE (trunc, "trunc is unportable - "
518 "use gnulib module trunc for portability");
519# endif
520#endif
521
522#if @GNULIB_TRUNCL@
523# if @REPLACE_TRUNCL@
524# if !(defined __cplusplus && defined GNULIB_NAMESPACE)
525# undef truncl
526# define truncl rpl_truncl
527# endif
528_GL_FUNCDECL_RPL (truncl, long double, (long double x));
529_GL_CXXALIAS_RPL (truncl, long double, (long double x));
530# else
531# if !@HAVE_DECL_TRUNCL@
532_GL_FUNCDECL_SYS (truncl, long double, (long double x));
533# endif
534_GL_CXXALIAS_SYS (truncl, long double, (long double x));
535# endif
536_GL_CXXALIASWARN (truncl);
537#elif defined GNULIB_POSIXCHECK
538# undef truncl
539# if HAVE_RAW_DECL_TRUNCL
540_GL_WARN_ON_USE (truncl, "truncl is unportable - "
541 "use gnulib module truncl for portability");
542# endif
543#endif
544
545
546#if @GNULIB_ISFINITE@
547# if @REPLACE_ISFINITE@
548_GL_EXTERN_C int gl_isfinitef (float x);
549_GL_EXTERN_C int gl_isfinited (double x);
550_GL_EXTERN_C int gl_isfinitel (long double x);
551# undef isfinite
552# define isfinite(x) \
553 (sizeof (x) == sizeof (long double) ? gl_isfinitel (x) : \
554 sizeof (x) == sizeof (double) ? gl_isfinited (x) : \
555 gl_isfinitef (x))
556# endif
557#elif defined GNULIB_POSIXCHECK
558# if defined isfinite
559_GL_WARN_REAL_FLOATING_DECL (isfinite);
560# undef isfinite
561# define isfinite(x) _GL_WARN_REAL_FLOATING_IMPL (isfinite, x)
562# endif
563#endif
564
565
566#if @GNULIB_ISINF@
567# if @REPLACE_ISINF@
568_GL_EXTERN_C int gl_isinff (float x);
569_GL_EXTERN_C int gl_isinfd (double x);
570_GL_EXTERN_C int gl_isinfl (long double x);
571# undef isinf
572# define isinf(x) \
573 (sizeof (x) == sizeof (long double) ? gl_isinfl (x) : \
574 sizeof (x) == sizeof (double) ? gl_isinfd (x) : \
575 gl_isinff (x))
576# endif
577#elif defined GNULIB_POSIXCHECK
578# if defined isinf
579_GL_WARN_REAL_FLOATING_DECL (isinf);
580# undef isinf
581# define isinf(x) _GL_WARN_REAL_FLOATING_IMPL (isinf, x)
582# endif
583#endif
584
585
586#if @GNULIB_ISNANF@
587/* Test for NaN for 'float' numbers. */
588# if @HAVE_ISNANF@
589/* The original <math.h> included above provides a declaration of isnan macro
590 or (older) isnanf function. */
591# if __GNUC__ >= 4
592 /* GCC 4.0 and newer provides three built-ins for isnan. */
593# undef isnanf
594# define isnanf(x) __builtin_isnanf ((float)(x))
595# elif defined isnan
596# undef isnanf
597# define isnanf(x) isnan ((float)(x))
598# endif
599# else
600/* Test whether X is a NaN. */
601# undef isnanf
602# define isnanf rpl_isnanf
603_GL_EXTERN_C int isnanf (float x);
604# endif
605#endif
606
607#if @GNULIB_ISNAND@
608/* Test for NaN for 'double' numbers.
609 This function is a gnulib extension, unlike isnan() which applied only
610 to 'double' numbers earlier but now is a type-generic macro. */
611# if @HAVE_ISNAND@
612/* The original <math.h> included above provides a declaration of isnan macro. */
613# if __GNUC__ >= 4
614 /* GCC 4.0 and newer provides three built-ins for isnan. */
615# undef isnand
616# define isnand(x) __builtin_isnan ((double)(x))
617# else
618# undef isnand
619# define isnand(x) isnan ((double)(x))
620# endif
621# else
622/* Test whether X is a NaN. */
623# undef isnand
624# define isnand rpl_isnand
625_GL_EXTERN_C int isnand (double x);
626# endif
627#endif
628
629#if @GNULIB_ISNANL@
630/* Test for NaN for 'long double' numbers. */
631# if @HAVE_ISNANL@
632/* The original <math.h> included above provides a declaration of isnan macro or (older) isnanl function. */
633# if __GNUC__ >= 4
634 /* GCC 4.0 and newer provides three built-ins for isnan. */
635# undef isnanl
636# define isnanl(x) __builtin_isnanl ((long double)(x))
637# elif defined isnan
638# undef isnanl
639# define isnanl(x) isnan ((long double)(x))
640# endif
641# else
642/* Test whether X is a NaN. */
643# undef isnanl
644# define isnanl rpl_isnanl
645_GL_EXTERN_C int isnanl (long double x);
646# endif
647#endif
648
649/* This must come *after* the snippets for GNULIB_ISNANF and GNULIB_ISNANL! */
650#if @GNULIB_ISNAN@
651# if @REPLACE_ISNAN@
652/* We can't just use the isnanf macro (e.g.) as exposed by
653 isnanf.h (e.g.) here, because those may end up being macros
654 that recursively expand back to isnan. So use the gnulib
655 replacements for them directly. */
656# if @HAVE_ISNANF@ && __GNUC__ >= 4
657# define gl_isnan_f(x) __builtin_isnan ((float)(x))
658# else
659_GL_EXTERN_C int rpl_isnanf (float x);
660# define gl_isnan_f(x) rpl_isnanf (x)
661# endif
662# if @HAVE_ISNAND@ && __GNUC__ >= 4
663# define gl_isnan_d(x) __builtin_isnan ((double)(x))
664# else
665_GL_EXTERN_C int rpl_isnand (double x);
666# define gl_isnan_d(x) rpl_isnand (x)
667# endif
668# if @HAVE_ISNANL@ && __GNUC__ >= 4
669# define gl_isnan_l(x) __builtin_isnan ((long double)(x))
670# else
671_GL_EXTERN_C int rpl_isnanl (long double x);
672# define gl_isnan_l(x) rpl_isnanl (x)
673# endif
674# undef isnan
675# define isnan(x) \
676 (sizeof (x) == sizeof (long double) ? gl_isnan_l (x) : \
677 sizeof (x) == sizeof (double) ? gl_isnan_d (x) : \
678 gl_isnan_f (x))
679# endif
680#elif defined GNULIB_POSIXCHECK
681# if defined isnan
682_GL_WARN_REAL_FLOATING_DECL (isnan);
683# undef isnan
684# define isnan(x) _GL_WARN_REAL_FLOATING_IMPL (isnan, x)
685# endif
686#endif
687
688
689#if @GNULIB_SIGNBIT@
690# if @REPLACE_SIGNBIT_USING_GCC@
691# undef signbit
692 /* GCC 4.0 and newer provides three built-ins for signbit. */
693# define signbit(x) \
694 (sizeof (x) == sizeof (long double) ? __builtin_signbitl (x) : \
695 sizeof (x) == sizeof (double) ? __builtin_signbit (x) : \
696 __builtin_signbitf (x))
697# endif
698# if @REPLACE_SIGNBIT@
699# undef signbit
700_GL_EXTERN_C int gl_signbitf (float arg);
701_GL_EXTERN_C int gl_signbitd (double arg);
702_GL_EXTERN_C int gl_signbitl (long double arg);
703# if __GNUC__ >= 2 && !__STRICT_ANSI__
704# if defined FLT_SIGNBIT_WORD && defined FLT_SIGNBIT_BIT && !defined gl_signbitf
705# define gl_signbitf_OPTIMIZED_MACRO
706# define gl_signbitf(arg) \
707 ({ union { float _value; \
708 unsigned int _word[(sizeof (float) + sizeof (unsigned int) - 1) / sizeof (unsigned int)]; \
709 } _m; \
710 _m._value = (arg); \
711 (_m._word[FLT_SIGNBIT_WORD] >> FLT_SIGNBIT_BIT) & 1; \
712 })
713# endif
714# if defined DBL_SIGNBIT_WORD && defined DBL_SIGNBIT_BIT && !defined gl_signbitd
715# define gl_signbitd_OPTIMIZED_MACRO
716# define gl_signbitd(arg) \
717 ({ union { double _value; \
718 unsigned int _word[(sizeof (double) + sizeof (unsigned int) - 1) / sizeof (unsigned int)]; \
719 } _m; \
720 _m._value = (arg); \
721 (_m._word[DBL_SIGNBIT_WORD] >> DBL_SIGNBIT_BIT) & 1; \
722 })
723# endif
724# if defined LDBL_SIGNBIT_WORD && defined LDBL_SIGNBIT_BIT && !defined gl_signbitl
725# define gl_signbitl_OPTIMIZED_MACRO
726# define gl_signbitl(arg) \
727 ({ union { long double _value; \
728 unsigned int _word[(sizeof (long double) + sizeof (unsigned int) - 1) / sizeof (unsigned int)]; \
729 } _m; \
730 _m._value = (arg); \
731 (_m._word[LDBL_SIGNBIT_WORD] >> LDBL_SIGNBIT_BIT) & 1; \
732 })
733# endif
734# endif
735# define signbit(x) \
736 (sizeof (x) == sizeof (long double) ? gl_signbitl (x) : \
737 sizeof (x) == sizeof (double) ? gl_signbitd (x) : \
738 gl_signbitf (x))
739# endif
740#elif defined GNULIB_POSIXCHECK
741# if defined signbit
742_GL_WARN_REAL_FLOATING_DECL (signbit);
743# undef signbit
744# define signbit(x) _GL_WARN_REAL_FLOATING_IMPL (signbit, x)
745# endif
746#endif
747
748
749#endif /* _GL_MATH_H */
750#endif /* _GL_MATH_H */