lib/math.in.h

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