contrib/libc++/include/cmath

   1 // -*- C++ -*-
   2 //===---------------------------- cmath -----------------------------------===//
   3 //
   4 //                     The LLVM Compiler Infrastructure
   5 //
   6 // This file is dual licensed under the MIT and the University of Illinois Open
   7 // Source Licenses. See LICENSE.TXT for details.
   8 //
   9 //===----------------------------------------------------------------------===//
  10
  11 #ifndef _LIBCPP_CMATH
  12 #define _LIBCPP_CMATH
  13
  14 /*
  15     cmath synopsis
  16
  17 Macros:
  18
  19     HUGE_VAL
  20     HUGE_VALF               // C99
  21     HUGE_VALL               // C99
  22     INFINITY                // C99
  23     NAN                     // C99
  24     FP_INFINITE             // C99
  25     FP_NAN                  // C99
  26     FP_NORMAL               // C99
  27     FP_SUBNORMAL            // C99
  28     FP_ZERO                 // C99
  29     FP_FAST_FMA             // C99
  30     FP_FAST_FMAF            // C99
  31     FP_FAST_FMAL            // C99
  32     FP_ILOGB0               // C99
  33     FP_ILOGBNAN             // C99
  34     MATH_ERRNO              // C99
  35     MATH_ERREXCEPT          // C99
  36     math_errhandling        // C99
  37
  38 namespace std
  39 {
  40
  41 Types:
  42
  43     float_t                 // C99
  44     double_t                // C99
  45
  46 // C90
  47
  48 floating_point abs(floating_point x);
  49
  50 floating_point acos (arithmetic x);
  51 float          acosf(float x);
  52 long double    acosl(long double x);
  53
  54 floating_point asin (arithmetic x);
  55 float          asinf(float x);
  56 long double    asinl(long double x);
  57
  58 floating_point atan (arithmetic x);
  59 float          atanf(float x);
  60 long double    atanl(long double x);
  61
  62 floating_point atan2 (arithmetic y, arithmetic x);
  63 float          atan2f(float y, float x);
  64 long double    atan2l(long double y, long double x);
  65
  66 floating_point ceil (arithmetic x);
  67 float          ceilf(float x);
  68 long double    ceill(long double x);
  69
  70 floating_point cos (arithmetic x);
  71 float          cosf(float x);
  72 long double    cosl(long double x);
  73
  74 floating_point cosh (arithmetic x);
  75 float          coshf(float x);
  76 long double    coshl(long double x);
  77
  78 floating_point exp (arithmetic x);
  79 float          expf(float x);
  80 long double    expl(long double x);
  81
  82 floating_point fabs (arithmetic x);
  83 float          fabsf(float x);
  84 long double    fabsl(long double x);
  85
  86 floating_point floor (arithmetic x);
  87 float          floorf(float x);
  88 long double    floorl(long double x);
  89
  90 floating_point fmod (arithmetic x, arithmetic y);
  91 float          fmodf(float x, float y);
  92 long double    fmodl(long double x, long double y);
  93
  94 floating_point frexp (arithmetic value, int* exp);
  95 float          frexpf(float value, int* exp);
  96 long double    frexpl(long double value, int* exp);
  97
  98 floating_point ldexp (arithmetic value, int exp);
  99 float          ldexpf(float value, int exp);
 100 long double    ldexpl(long double value, int exp);
 101
 102 floating_point log (arithmetic x);
 103 float          logf(float x);
 104 long double    logl(long double x);
 105
 106 floating_point log10 (arithmetic x);
 107 float          log10f(float x);
 108 long double    log10l(long double x);
 109
 110 floating_point modf (floating_point value, floating_point* iptr);
 111 float          modff(float value, float* iptr);
 112 long double    modfl(long double value, long double* iptr);
 113
 114 floating_point pow (arithmetic x, arithmetic y);
 115 float          powf(float x, float y);
 116 long double    powl(long double x, long double y);
 117
 118 floating_point sin (arithmetic x);
 119 float          sinf(float x);
 120 long double    sinl(long double x);
 121
 122 floating_point sinh (arithmetic x);
 123 float          sinhf(float x);
 124 long double    sinhl(long double x);
 125
 126 floating_point sqrt (arithmetic x);
 127 float          sqrtf(float x);
 128 long double    sqrtl(long double x);
 129
 130 floating_point tan (arithmetic x);
 131 float          tanf(float x);
 132 long double    tanl(long double x);
 133
 134 floating_point tanh (arithmetic x);
 135 float          tanhf(float x);
 136 long double    tanhl(long double x);
 137
 138 //  C99
 139
 140 bool signbit(floating_point x);
 141
 142 int fpclassify(floating_point x);
 143
 144 bool isfinite(floating_point x);
 145 bool isinf(floating_point x);
 146 bool isnan(floating_point x);
 147 bool isnormal(floating_point x);
 148
 149 bool isgreater(floating_point x, floating_point y);
 150 bool isgreaterequal(floating_point x, floating_point y);
 151 bool isless(floating_point x, floating_point y);
 152 bool islessequal(floating_point x, floating_point y);
 153 bool islessgreater(floating_point x, floating_point y);
 154 bool isunordered(floating_point x, floating_point y);
 155
 156 floating_point acosh (arithmetic x);
 157 float          acoshf(float x);
 158 long double    acoshl(long double x);
 159
 160 floating_point asinh (arithmetic x);
 161 float          asinhf(float x);
 162 long double    asinhl(long double x);
 163
 164 floating_point atanh (arithmetic x);
 165 float          atanhf(float x);
 166 long double    atanhl(long double x);
 167
 168 floating_point cbrt (arithmetic x);
 169 float          cbrtf(float x);
 170 long double    cbrtl(long double x);
 171
 172 floating_point copysign (arithmetic x, arithmetic y);
 173 float          copysignf(float x, float y);
 174 long double    copysignl(long double x, long double y);
 175
 176 floating_point erf (arithmetic x);
 177 float          erff(float x);
 178 long double    erfl(long double x);
 179
 180 floating_point erfc (arithmetic x);
 181 float          erfcf(float x);
 182 long double    erfcl(long double x);
 183
 184 floating_point exp2 (arithmetic x);
 185 float          exp2f(float x);
 186 long double    exp2l(long double x);
 187
 188 floating_point expm1 (arithmetic x);
 189 float          expm1f(float x);
 190 long double    expm1l(long double x);
 191
 192 floating_point fdim (arithmetic x, arithmetic y);
 193 float          fdimf(float x, float y);
 194 long double    fdiml(long double x, long double y);
 195
 196 floating_point fma (arithmetic x, arithmetic y, arithmetic z);
 197 float          fmaf(float x, float y, float z);
 198 long double    fmal(long double x, long double y, long double z);
 199
 200 floating_point fmax (arithmetic x, arithmetic y);
 201 float          fmaxf(float x, float y);
 202 long double    fmaxl(long double x, long double y);
 203
 204 floating_point fmin (arithmetic x, arithmetic y);
 205 float          fminf(float x, float y);
 206 long double    fminl(long double x, long double y);
 207
 208 floating_point hypot (arithmetic x, arithmetic y);
 209 float          hypotf(float x, float y);
 210 long double    hypotl(long double x, long double y);
 211
 212 int ilogb (arithmetic x);
 213 int ilogbf(float x);
 214 int ilogbl(long double x);
 215
 216 floating_point lgamma (arithmetic x);
 217 float          lgammaf(float x);
 218 long double    lgammal(long double x);
 219
 220 long long llrint (arithmetic x);
 221 long long llrintf(float x);
 222 long long llrintl(long double x);
 223
 224 long long llround (arithmetic x);
 225 long long llroundf(float x);
 226 long long llroundl(long double x);
 227
 228 floating_point log1p (arithmetic x);
 229 float          log1pf(float x);
 230 long double    log1pl(long double x);
 231
 232 floating_point log2 (arithmetic x);
 233 float          log2f(float x);
 234 long double    log2l(long double x);
 235
 236 floating_point logb (arithmetic x);
 237 float          logbf(float x);
 238 long double    logbl(long double x);
 239
 240 long lrint (arithmetic x);
 241 long lrintf(float x);
 242 long lrintl(long double x);
 243
 244 long lround (arithmetic x);
 245 long lroundf(float x);
 246 long lroundl(long double x);
 247
 248 double      nan (const char* str);
 249 float       nanf(const char* str);
 250 long double nanl(const char* str);
 251
 252 floating_point nearbyint (arithmetic x);
 253 float          nearbyintf(float x);
 254 long double    nearbyintl(long double x);
 255
 256 floating_point nextafter (arithmetic x, arithmetic y);
 257 float          nextafterf(float x, float y);
 258 long double    nextafterl(long double x, long double y);
 259
 260 floating_point nexttoward (arithmetic x, long double y);
 261 float          nexttowardf(float x, long double y);
 262 long double    nexttowardl(long double x, long double y);
 263
 264 floating_point remainder (arithmetic x, arithmetic y);
 265 float          remainderf(float x, float y);
 266 long double    remainderl(long double x, long double y);
 267
 268 floating_point remquo (arithmetic x, arithmetic y, int* pquo);
 269 float          remquof(float x, float y, int* pquo);
 270 long double    remquol(long double x, long double y, int* pquo);
 271
 272 floating_point rint (arithmetic x);
 273 float          rintf(float x);
 274 long double    rintl(long double x);
 275
 276 floating_point round (arithmetic x);
 277 float          roundf(float x);
 278 long double    roundl(long double x);
 279
 280 floating_point scalbln (arithmetic x, long ex);
 281 float          scalblnf(float x, long ex);
 282 long double    scalblnl(long double x, long ex);
 283
 284 floating_point scalbn (arithmetic x, int ex);
 285 float          scalbnf(float x, int ex);
 286 long double    scalbnl(long double x, int ex);
 287
 288 floating_point tgamma (arithmetic x);
 289 float          tgammaf(float x);
 290 long double    tgammal(long double x);
 291
 292 floating_point trunc (arithmetic x);
 293 float          truncf(float x);
 294 long double    truncl(long double x);
 295
 296 }  // std
 297
 298 */
 299
 300 #include <__config>
 301 #include <math.h>
 302 #include <type_traits>
 303
 304 #ifdef _MSC_VER
 305 #include "support/win32/math_win32.h"
 306 #endif
 307
 308 #if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER)
 309 #pragma GCC system_header
 310 #endif
 311
 312 // signbit
 313
 314 #ifdef signbit
 315
 316 template <class _A1>
 317 _LIBCPP_ALWAYS_INLINE
 318 bool
 319 __libcpp_signbit(_A1 __x)
 320 {
 321     return signbit(__x);
 322 }
 323
 324 #undef signbit
 325
 326 template <class _A1>
 327 inline _LIBCPP_INLINE_VISIBILITY
 328 typename std::enable_if<std::is_floating_point<_A1>::value, bool>::type
 329 signbit(_A1 __x)
 330 {
 331     return __libcpp_signbit(__x);
 332 }
 333
 334 #endif  // signbit
 335
 336 // fpclassify
 337
 338 #ifdef fpclassify
 339
 340 template <class _A1>
 341 _LIBCPP_ALWAYS_INLINE
 342 int
 343 __libcpp_fpclassify(_A1 __x)
 344 {
 345     return fpclassify(__x);
 346 }
 347
 348 #undef fpclassify
 349
 350 template <class _A1>
 351 inline _LIBCPP_INLINE_VISIBILITY
 352 typename std::enable_if<std::is_floating_point<_A1>::value, int>::type
 353 fpclassify(_A1 __x)
 354 {
 355     return __libcpp_fpclassify(__x);
 356 }
 357
 358 #endif  // fpclassify
 359
 360 // isfinite
 361
 362 #ifdef isfinite
 363
 364 template <class _A1>
 365 _LIBCPP_ALWAYS_INLINE
 366 bool
 367 __libcpp_isfinite(_A1 __x)
 368 {
 369     return isfinite(__x);
 370 }
 371
 372 #undef isfinite
 373
 374 template <class _A1>
 375 inline _LIBCPP_INLINE_VISIBILITY
 376 typename std::enable_if<std::is_floating_point<_A1>::value, bool>::type
 377 isfinite(_A1 __x)
 378 {
 379     return __libcpp_isfinite(__x);
 380 }
 381
 382 #endif  // isfinite
 383
 384 // isinf
 385
 386 #ifdef isinf
 387
 388 template <class _A1>
 389 _LIBCPP_ALWAYS_INLINE
 390 bool
 391 __libcpp_isinf(_A1 __x)
 392 {
 393     return isinf(__x);
 394 }
 395
 396 #undef isinf
 397
 398 template <class _A1>
 399 inline _LIBCPP_INLINE_VISIBILITY
 400 typename std::enable_if<std::is_floating_point<_A1>::value, bool>::type
 401 isinf(_A1 __x)
 402 {
 403     return __libcpp_isinf(__x);
 404 }
 405
 406 #endif  // isinf
 407
 408 // isnan
 409
 410 #ifdef isnan
 411
 412 template <class _A1>
 413 _LIBCPP_ALWAYS_INLINE
 414 bool
 415 __libcpp_isnan(_A1 __x)
 416 {
 417     return isnan(__x);
 418 }
 419
 420 #undef isnan
 421
 422 template <class _A1>
 423 inline _LIBCPP_INLINE_VISIBILITY
 424 typename std::enable_if<std::is_floating_point<_A1>::value, bool>::type
 425 isnan(_A1 __x)
 426 {
 427     return __libcpp_isnan(__x);
 428 }
 429
 430 #endif  // isnan
 431
 432 // isnormal
 433
 434 #ifdef isnormal
 435
 436 template <class _A1>
 437 _LIBCPP_ALWAYS_INLINE
 438 bool
 439 __libcpp_isnormal(_A1 __x)
 440 {
 441     return isnormal(__x);
 442 }
 443
 444 #undef isnormal
 445
 446 template <class _A1>
 447 inline _LIBCPP_INLINE_VISIBILITY
 448 typename std::enable_if<std::is_floating_point<_A1>::value, bool>::type
 449 isnormal(_A1 __x)
 450 {
 451     return __libcpp_isnormal(__x);
 452 }
 453
 454 #endif  // isnormal
 455
 456 // isgreater
 457
 458 #ifdef isgreater
 459
 460 template <class _A1, class _A2>
 461 _LIBCPP_ALWAYS_INLINE
 462 bool
 463 __libcpp_isgreater(_A1 __x, _A2 __y)
 464 {
 465     return isgreater(__x, __y);
 466 }
 467
 468 #undef isgreater
 469
 470 template <class _A1, class _A2>
 471 inline _LIBCPP_INLINE_VISIBILITY
 472 typename std::enable_if
 473 <
 474     std::is_floating_point<_A1>::value &&
 475     std::is_floating_point<_A2>::value,
 476     bool
 477 >::type
 478 isgreater(_A1 __x, _A2 __y)
 479 {
 480     return __libcpp_isgreater(__x, __y);
 481 }
 482
 483 #endif  // isgreater
 484
 485 // isgreaterequal
 486
 487 #ifdef isgreaterequal
 488
 489 template <class _A1, class _A2>
 490 _LIBCPP_ALWAYS_INLINE
 491 bool
 492 __libcpp_isgreaterequal(_A1 __x, _A2 __y)
 493 {
 494     return isgreaterequal(__x, __y);
 495 }
 496
 497 #undef isgreaterequal
 498
 499 template <class _A1, class _A2>
 500 inline _LIBCPP_INLINE_VISIBILITY
 501 typename std::enable_if
 502 <
 503     std::is_floating_point<_A1>::value &&
 504     std::is_floating_point<_A2>::value,
 505     bool
 506 >::type
 507 isgreaterequal(_A1 __x, _A2 __y)
 508 {
 509     return __libcpp_isgreaterequal(__x, __y);
 510 }
 511
 512 #endif  // isgreaterequal
 513
 514 // isless
 515
 516 #ifdef isless
 517
 518 template <class _A1, class _A2>
 519 _LIBCPP_ALWAYS_INLINE
 520 bool
 521 __libcpp_isless(_A1 __x, _A2 __y)
 522 {
 523     return isless(__x, __y);
 524 }
 525
 526 #undef isless
 527
 528 template <class _A1, class _A2>
 529 inline _LIBCPP_INLINE_VISIBILITY
 530 typename std::enable_if
 531 <
 532     std::is_floating_point<_A1>::value &&
 533     std::is_floating_point<_A2>::value,
 534     bool
 535 >::type
 536 isless(_A1 __x, _A2 __y)
 537 {
 538     return __libcpp_isless(__x, __y);
 539 }
 540
 541 #endif  // isless
 542
 543 // islessequal
 544
 545 #ifdef islessequal
 546
 547 template <class _A1, class _A2>
 548 _LIBCPP_ALWAYS_INLINE
 549 bool
 550 __libcpp_islessequal(_A1 __x, _A2 __y)
 551 {
 552     return islessequal(__x, __y);
 553 }
 554
 555 #undef islessequal
 556
 557 template <class _A1, class _A2>
 558 inline _LIBCPP_INLINE_VISIBILITY
 559 typename std::enable_if
 560 <
 561     std::is_floating_point<_A1>::value &&
 562     std::is_floating_point<_A2>::value,
 563     bool
 564 >::type
 565 islessequal(_A1 __x, _A2 __y)
 566 {
 567     return __libcpp_islessequal(__x, __y);
 568 }
 569
 570 #endif  // islessequal
 571
 572 // islessgreater
 573
 574 #ifdef islessgreater
 575
 576 template <class _A1, class _A2>
 577 _LIBCPP_ALWAYS_INLINE
 578 bool
 579 __libcpp_islessgreater(_A1 __x, _A2 __y)
 580 {
 581     return islessgreater(__x, __y);
 582 }
 583
 584 #undef islessgreater
 585
 586 template <class _A1, class _A2>
 587 inline _LIBCPP_INLINE_VISIBILITY
 588 typename std::enable_if
 589 <
 590     std::is_floating_point<_A1>::value &&
 591     std::is_floating_point<_A2>::value,
 592     bool
 593 >::type
 594 islessgreater(_A1 __x, _A2 __y)
 595 {
 596     return __libcpp_islessgreater(__x, __y);
 597 }
 598
 599 #endif  // islessgreater
 600
 601 // isunordered
 602
 603 #ifdef isunordered
 604
 605 template <class _A1, class _A2>
 606 _LIBCPP_ALWAYS_INLINE
 607 bool
 608 __libcpp_isunordered(_A1 __x, _A2 __y)
 609 {
 610     return isunordered(__x, __y);
 611 }
 612
 613 #undef isunordered
 614
 615 template <class _A1, class _A2>
 616 inline _LIBCPP_INLINE_VISIBILITY
 617 typename std::enable_if
 618 <
 619     std::is_floating_point<_A1>::value &&
 620     std::is_floating_point<_A2>::value,
 621     bool
 622 >::type
 623 isunordered(_A1 __x, _A2 __y)
 624 {
 625     return __libcpp_isunordered(__x, __y);
 626 }
 627
 628 #endif  // isunordered
 629
 630 _LIBCPP_BEGIN_NAMESPACE_STD
 631
 632 using ::signbit;
 633 using ::fpclassify;
 634 using ::isfinite;
 635 using ::isinf;
 636 using ::isnan;
 637 using ::isnormal;
 638 using ::isgreater;
 639 using ::isgreaterequal;
 640 using ::isless;
 641 using ::islessequal;
 642 using ::islessgreater;
 643 using ::isunordered;
 644 using ::isunordered;
 645
 646 using ::float_t;
 647 using ::double_t;
 648
 649 // abs
 650
 651 template <class _A1>
 652 inline _LIBCPP_INLINE_VISIBILITY
 653 typename enable_if<is_floating_point<_A1>::value, _A1>::type
 654 abs(_A1 __x) {return fabs(__x);}
 655
 656 #ifndef __sun__
 657
 658 // acos
 659
 660 using ::acos;
 661 using ::acosf;
 662
 663 #ifndef _MSC_VER
 664 inline _LIBCPP_INLINE_VISIBILITY float       acos(float __x)       {return acosf(__x);}
 665 inline _LIBCPP_INLINE_VISIBILITY long double acos(long double __x) {return acosl(__x);}
 666 #endif
 667
 668 template <class _A1>
 669 inline _LIBCPP_INLINE_VISIBILITY
 670 typename enable_if<is_integral<_A1>::value, double>::type
 671 acos(_A1 __x) {return acos((double)__x);}
 672
 673 // asin
 674
 675 using ::asin;
 676 using ::asinf;
 677
 678 #ifndef _MSC_VER
 679 inline _LIBCPP_INLINE_VISIBILITY float       asin(float __x)       {return asinf(__x);}
 680 inline _LIBCPP_INLINE_VISIBILITY long double asin(long double __x) {return asinl(__x);}
 681 #endif
 682
 683 template <class _A1>
 684 inline _LIBCPP_INLINE_VISIBILITY
 685 typename enable_if<is_integral<_A1>::value, double>::type
 686 asin(_A1 __x) {return asin((double)__x);}
 687
 688 // atan
 689
 690 using ::atan;
 691 using ::atanf;
 692
 693 #ifndef _MSC_VER
 694 inline _LIBCPP_INLINE_VISIBILITY float       atan(float __x)       {return atanf(__x);}
 695 inline _LIBCPP_INLINE_VISIBILITY long double atan(long double __x) {return atanl(__x);}
 696 #endif
 697
 698 template <class _A1>
 699 inline _LIBCPP_INLINE_VISIBILITY
 700 typename enable_if<is_integral<_A1>::value, double>::type
 701 atan(_A1 __x) {return atan((double)__x);}
 702
 703 // atan2
 704
 705 using ::atan2;
 706 using ::atan2f;
 707
 708 #ifndef _MSC_VER
 709 inline _LIBCPP_INLINE_VISIBILITY float       atan2(float __y, float __x)             {return atan2f(__y, __x);}
 710 inline _LIBCPP_INLINE_VISIBILITY long double atan2(long double __y, long double __x) {return atan2l(__y, __x);}
 711 #endif
 712
 713 template <class _A1, class _A2>
 714 inline _LIBCPP_INLINE_VISIBILITY
 715 typename enable_if
 716 <
 717     is_arithmetic<_A1>::value &&
 718     is_arithmetic<_A2>::value,
 719     typename __promote<_A1, _A2>::type
 720 >::type
 721 atan2(_A1 __y, _A2 __x)
 722 {
 723     typedef typename __promote<_A1, _A2>::type __result_type;
 724     static_assert((!(is_same<_A1, __result_type>::value &&
 725                       is_same<_A2, __result_type>::value)), "");
 726     return atan2((__result_type)__y, (__result_type)__x);
 727 }
 728
 729 // ceil
 730
 731 using ::ceil;
 732 using ::ceilf;
 733
 734 #ifndef _MSC_VER
 735 inline _LIBCPP_INLINE_VISIBILITY float       ceil(float __x)       {return ceilf(__x);}
 736 inline _LIBCPP_INLINE_VISIBILITY long double ceil(long double __x) {return ceill(__x);}
 737 #endif
 738
 739 template <class _A1>
 740 inline _LIBCPP_INLINE_VISIBILITY
 741 typename enable_if<is_integral<_A1>::value, double>::type
 742 ceil(_A1 __x) {return ceil((double)__x);}
 743
 744 // cos
 745
 746 using ::cos;
 747 using ::cosf;
 748
 749 #ifndef _MSC_VER
 750 inline _LIBCPP_INLINE_VISIBILITY float       cos(float __x)       {return cosf(__x);}
 751 inline _LIBCPP_INLINE_VISIBILITY long double cos(long double __x) {return cosl(__x);}
 752 #endif
 753
 754 template <class _A1>
 755 inline _LIBCPP_ALWAYS_INLINE _LIBCPP_INLINE_VISIBILITY
 756 typename enable_if<is_integral<_A1>::value, double>::type
 757 cos(_A1 __x) {return cos((double)__x);}
 758
 759 // cosh
 760
 761 using ::cosh;
 762 using ::coshf;
 763
 764 #ifndef _MSC_VER
 765 inline _LIBCPP_INLINE_VISIBILITY float       cosh(float __x)       {return coshf(__x);}
 766 inline _LIBCPP_INLINE_VISIBILITY long double cosh(long double __x) {return coshl(__x);}
 767 #endif
 768
 769 template <class _A1>
 770 inline _LIBCPP_INLINE_VISIBILITY
 771 typename enable_if<is_integral<_A1>::value, double>::type
 772 cosh(_A1 __x) {return cosh((double)__x);}
 773
 774 #endif // __sun__
 775 // exp
 776
 777 using ::exp;
 778 using ::expf;
 779
 780 #ifndef __sun__
 781
 782 #ifndef _MSC_VER
 783 inline _LIBCPP_INLINE_VISIBILITY float       exp(float __x)       {return expf(__x);}
 784 inline _LIBCPP_INLINE_VISIBILITY long double exp(long double __x) {return expl(__x);}
 785 #endif
 786
 787
 788 template <class _A1>
 789 inline _LIBCPP_INLINE_VISIBILITY
 790 typename enable_if<is_integral<_A1>::value, double>::type
 791 exp(_A1 __x) {return exp((double)__x);}
 792
 793 // fabs
 794
 795 using ::fabs;
 796 using ::fabsf;
 797
 798 #ifndef _MSC_VER
 799 inline _LIBCPP_INLINE_VISIBILITY float       fabs(float __x)       {return fabsf(__x);}
 800 inline _LIBCPP_INLINE_VISIBILITY long double fabs(long double __x) {return fabsl(__x);}
 801 #endif
 802
 803 template <class _A1>
 804 inline _LIBCPP_INLINE_VISIBILITY
 805 typename enable_if<is_integral<_A1>::value, double>::type
 806 fabs(_A1 __x) {return fabs((double)__x);}
 807
 808 // floor
 809
 810 using ::floor;
 811 using ::floorf;
 812
 813 #ifndef _MSC_VER
 814 inline _LIBCPP_INLINE_VISIBILITY float       floor(float __x)       {return floorf(__x);}
 815 inline _LIBCPP_INLINE_VISIBILITY long double floor(long double __x) {return floorl(__x);}
 816 #endif
 817
 818 template <class _A1>
 819 inline _LIBCPP_INLINE_VISIBILITY
 820 typename enable_if<is_integral<_A1>::value, double>::type
 821 floor(_A1 __x) {return floor((double)__x);}
 822
 823 // fmod
 824
 825 #endif //__sun__
 826 using ::fmod;
 827 using ::fmodf;
 828 #ifndef __sun__
 829
 830 #ifndef _MSC_VER
 831 inline _LIBCPP_INLINE_VISIBILITY float       fmod(float __x, float __y)             {return fmodf(__x, __y);}
 832 inline _LIBCPP_INLINE_VISIBILITY long double fmod(long double __x, long double __y) {return fmodl(__x, __y);}
 833 #endif
 834
 835 template <class _A1, class _A2>
 836 inline _LIBCPP_INLINE_VISIBILITY
 837 typename enable_if
 838 <
 839     is_arithmetic<_A1>::value &&
 840     is_arithmetic<_A2>::value,
 841     typename __promote<_A1, _A2>::type
 842 >::type
 843 fmod(_A1 __x, _A2 __y)
 844 {
 845     typedef typename __promote<_A1, _A2>::type __result_type;
 846     static_assert((!(is_same<_A1, __result_type>::value &&
 847                       is_same<_A2, __result_type>::value)), "");
 848     return fmod((__result_type)__x, (__result_type)__y);
 849 }
 850
 851
 852 // frexp
 853
 854 using ::frexp;
 855 using ::frexpf;
 856
 857 #ifndef _MSC_VER
 858 inline _LIBCPP_INLINE_VISIBILITY float       frexp(float __x, int* __e)       {return frexpf(__x, __e);}
 859 inline _LIBCPP_INLINE_VISIBILITY long double frexp(long double __x, int* __e) {return frexpl(__x, __e);}
 860 #endif
 861
 862 template <class _A1>
 863 inline _LIBCPP_INLINE_VISIBILITY
 864 typename enable_if<is_integral<_A1>::value, double>::type
 865 frexp(_A1 __x, int* __e) {return frexp((double)__x, __e);}
 866
 867 // ldexp
 868
 869 using ::ldexp;
 870 using ::ldexpf;
 871
 872 #ifndef _MSC_VER
 873 inline _LIBCPP_INLINE_VISIBILITY float       ldexp(float __x, int __e)       {return ldexpf(__x, __e);}
 874 inline _LIBCPP_INLINE_VISIBILITY long double ldexp(long double __x, int __e) {return ldexpl(__x, __e);}
 875 #endif
 876
 877 template <class _A1>
 878 inline _LIBCPP_INLINE_VISIBILITY
 879 typename enable_if<is_integral<_A1>::value, double>::type
 880 ldexp(_A1 __x, int __e) {return ldexp((double)__x, __e);}
 881
 882 // log
 883
 884 #endif // __sun__
 885 using ::log;
 886 using ::logf;
 887 #ifndef __sun__
 888
 889 #ifndef _MSC_VER
 890 inline _LIBCPP_INLINE_VISIBILITY float       log(float __x)       {return logf(__x);}
 891 inline _LIBCPP_INLINE_VISIBILITY long double log(long double __x) {return logl(__x);}
 892 #endif
 893
 894 template <class _A1>
 895 inline _LIBCPP_INLINE_VISIBILITY
 896 typename enable_if<is_integral<_A1>::value, double>::type
 897 log(_A1 __x) {return log((double)__x);}
 898
 899
 900 // log10
 901
 902 using ::log10;
 903 using ::log10f;
 904
 905 #ifndef _MSC_VER
 906 inline _LIBCPP_INLINE_VISIBILITY float       log10(float __x)       {return log10f(__x);}
 907 inline _LIBCPP_INLINE_VISIBILITY long double log10(long double __x) {return log10l(__x);}
 908 #endif
 909
 910 template <class _A1>
 911 inline _LIBCPP_INLINE_VISIBILITY
 912 typename enable_if<is_integral<_A1>::value, double>::type
 913 log10(_A1 __x) {return log10((double)__x);}
 914
 915 // modf
 916
 917 using ::modf;
 918 using ::modff;
 919
 920 #ifndef _MSC_VER
 921 inline _LIBCPP_INLINE_VISIBILITY float       modf(float __x, float* __y)             {return modff(__x, __y);}
 922 inline _LIBCPP_INLINE_VISIBILITY long double modf(long double __x, long double* __y) {return modfl(__x, __y);}
 923 #endif
 924
 925 // pow
 926
 927 #endif // __sun__
 928 using ::pow;
 929 using ::powf;
 930
 931 #ifndef __sun__
 932
 933 #ifndef _MSC_VER
 934 inline _LIBCPP_INLINE_VISIBILITY float       pow(float __x, float __y)             {return powf(__x, __y);}
 935 inline _LIBCPP_INLINE_VISIBILITY long double pow(long double __x, long double __y) {return powl(__x, __y);}
 936 #endif
 937
 938 template <class _A1, class _A2>
 939 inline _LIBCPP_INLINE_VISIBILITY
 940 typename enable_if
 941 <
 942     is_arithmetic<_A1>::value &&
 943     is_arithmetic<_A2>::value,
 944     typename __promote<_A1, _A2>::type
 945 >::type
 946 pow(_A1 __x, _A2 __y)
 947 {
 948     typedef typename __promote<_A1, _A2>::type __result_type;
 949     static_assert((!(is_same<_A1, __result_type>::value &&
 950                       is_same<_A2, __result_type>::value)), "");
 951     return pow((__result_type)__x, (__result_type)__y);
 952 }
 953
 954
 955 // sin
 956
 957 using ::sin;
 958 using ::sinf;
 959
 960 #ifndef _MSC_VER
 961 inline _LIBCPP_INLINE_VISIBILITY float       sin(float __x)       {return sinf(__x);}
 962 inline _LIBCPP_INLINE_VISIBILITY long double sin(long double __x) {return sinl(__x);}
 963 #endif
 964
 965 template <class _A1>
 966 inline _LIBCPP_INLINE_VISIBILITY
 967 typename enable_if<is_integral<_A1>::value, double>::type
 968 sin(_A1 __x) {return sin((double)__x);}
 969
 970 // sinh
 971
 972 using ::sinh;
 973 using ::sinhf;
 974
 975 #ifndef _MSC_VER
 976 inline _LIBCPP_INLINE_VISIBILITY float       sinh(float __x)       {return sinhf(__x);}
 977 inline _LIBCPP_INLINE_VISIBILITY long double sinh(long double __x) {return sinhl(__x);}
 978 #endif
 979
 980 template <class _A1>
 981 inline _LIBCPP_INLINE_VISIBILITY
 982 typename enable_if<is_integral<_A1>::value, double>::type
 983 sinh(_A1 __x) {return sinh((double)__x);}
 984
 985 // sqrt
 986
 987 #endif // __sun__
 988 using ::sqrt;
 989 using ::sqrtf;
 990
 991
 992 #if !(defined(_MSC_VER) || defined(__sun__))
 993 inline _LIBCPP_INLINE_VISIBILITY float       sqrt(float __x)       {return sqrtf(__x);}
 994 inline _LIBCPP_INLINE_VISIBILITY long double sqrt(long double __x) {return sqrtl(__x);}
 995 #endif
 996
 997 template <class _A1>
 998 inline _LIBCPP_INLINE_VISIBILITY
 999 typename enable_if<is_integral<_A1>::value, double>::type
1000 sqrt(_A1 __x) {return sqrt((double)__x);}
1001
1002 // tan
1003
1004 using ::tan;
1005 using ::tanf;
1006 #ifndef __sun__
1007
1008 #ifndef _MSC_VER
1009 inline _LIBCPP_INLINE_VISIBILITY float       tan(float __x)       {return tanf(__x);}
1010 inline _LIBCPP_INLINE_VISIBILITY long double tan(long double __x) {return tanl(__x);}
1011 #endif
1012
1013 template <class _A1>
1014 inline _LIBCPP_INLINE_VISIBILITY
1015 typename enable_if<is_integral<_A1>::value, double>::type
1016 tan(_A1 __x) {return tan((double)__x);}
1017
1018 // tanh
1019
1020 using ::tanh;
1021 using ::tanhf;
1022
1023 #ifndef _MSC_VER
1024 inline _LIBCPP_INLINE_VISIBILITY float       tanh(float __x)       {return tanhf(__x);}
1025 inline _LIBCPP_INLINE_VISIBILITY long double tanh(long double __x) {return tanhl(__x);}
1026 #endif
1027
1028 template <class _A1>
1029 inline _LIBCPP_INLINE_VISIBILITY
1030 typename enable_if<is_integral<_A1>::value, double>::type
1031 tanh(_A1 __x) {return tanh((double)__x);}
1032
1033 // acosh
1034
1035 #ifndef _MSC_VER
1036 using ::acosh;
1037 using ::acoshf;
1038
1039 inline _LIBCPP_INLINE_VISIBILITY float       acosh(float __x)       {return acoshf(__x);}
1040 inline _LIBCPP_INLINE_VISIBILITY long double acosh(long double __x) {return acoshl(__x);}
1041
1042 template <class _A1>
1043 inline _LIBCPP_INLINE_VISIBILITY
1044 typename enable_if<is_integral<_A1>::value, double>::type
1045 acosh(_A1 __x) {return acosh((double)__x);}
1046 #endif
1047
1048 // asinh
1049
1050 #ifndef _MSC_VER
1051 using ::asinh;
1052 using ::asinhf;
1053
1054 inline _LIBCPP_INLINE_VISIBILITY float       asinh(float __x)       {return asinhf(__x);}
1055 inline _LIBCPP_INLINE_VISIBILITY long double asinh(long double __x) {return asinhl(__x);}
1056
1057 template <class _A1>
1058 inline _LIBCPP_INLINE_VISIBILITY
1059 typename enable_if<is_integral<_A1>::value, double>::type
1060 asinh(_A1 __x) {return asinh((double)__x);}
1061 #endif
1062
1063 // atanh
1064
1065 #ifndef _MSC_VER
1066 using ::atanh;
1067 using ::atanhf;
1068
1069 inline _LIBCPP_INLINE_VISIBILITY float       atanh(float __x)       {return atanhf(__x);}
1070 inline _LIBCPP_INLINE_VISIBILITY long double atanh(long double __x) {return atanhl(__x);}
1071
1072 template <class _A1>
1073 inline _LIBCPP_INLINE_VISIBILITY
1074 typename enable_if<is_integral<_A1>::value, double>::type
1075 atanh(_A1 __x) {return atanh((double)__x);}
1076 #endif
1077
1078 // cbrt
1079
1080 #ifndef _MSC_VER
1081 using ::cbrt;
1082 using ::cbrtf;
1083
1084 inline _LIBCPP_INLINE_VISIBILITY float       cbrt(float __x)       {return cbrtf(__x);}
1085 inline _LIBCPP_INLINE_VISIBILITY long double cbrt(long double __x) {return cbrtl(__x);}
1086
1087 template <class _A1>
1088 inline _LIBCPP_INLINE_VISIBILITY
1089 typename enable_if<is_integral<_A1>::value, double>::type
1090 cbrt(_A1 __x) {return cbrt((double)__x);}
1091 #endif
1092
1093 // copysign
1094
1095 using ::copysign;
1096 using ::copysignf;
1097
1098 inline _LIBCPP_INLINE_VISIBILITY float       copysign(float __x, float __y)             {return copysignf(__x, __y);}
1099 inline _LIBCPP_INLINE_VISIBILITY long double copysign(long double __x, long double __y) {return copysignl(__x, __y);}
1100
1101 template <class _A1, class _A2>
1102 inline _LIBCPP_INLINE_VISIBILITY
1103 typename enable_if
1104 <
1105     is_arithmetic<_A1>::value &&
1106     is_arithmetic<_A2>::value,
1107     typename __promote<_A1, _A2>::type
1108 >::type
1109 copysign(_A1 __x, _A2 __y)
1110 {
1111     typedef typename __promote<_A1, _A2>::type __result_type;
1112     static_assert((!(is_same<_A1, __result_type>::value &&
1113                       is_same<_A2, __result_type>::value)), "");
1114     return copysign((__result_type)__x, (__result_type)__y);
1115 }
1116
1117 #ifndef _MSC_VER
1118
1119 // erf
1120
1121 using ::erf;
1122 using ::erff;
1123
1124 inline _LIBCPP_INLINE_VISIBILITY float       erf(float __x)       {return erff(__x);}
1125 inline _LIBCPP_INLINE_VISIBILITY long double erf(long double __x) {return erfl(__x);}
1126
1127 template <class _A1>
1128 inline _LIBCPP_INLINE_VISIBILITY
1129 typename enable_if<is_integral<_A1>::value, double>::type
1130 erf(_A1 __x) {return erf((double)__x);}
1131
1132 // erfc
1133
1134 using ::erfc;
1135 using ::erfcf;
1136
1137 inline _LIBCPP_INLINE_VISIBILITY float       erfc(float __x)       {return erfcf(__x);}
1138 inline _LIBCPP_INLINE_VISIBILITY long double erfc(long double __x) {return erfcl(__x);}
1139
1140 template <class _A1>
1141 inline _LIBCPP_INLINE_VISIBILITY
1142 typename enable_if<is_integral<_A1>::value, double>::type
1143 erfc(_A1 __x) {return erfc((double)__x);}
1144
1145 // exp2
1146
1147 using ::exp2;
1148 using ::exp2f;
1149
1150 inline _LIBCPP_INLINE_VISIBILITY float       exp2(float __x)       {return exp2f(__x);}
1151 inline _LIBCPP_INLINE_VISIBILITY long double exp2(long double __x) {return exp2l(__x);}
1152
1153 template <class _A1>
1154 inline _LIBCPP_INLINE_VISIBILITY
1155 typename enable_if<is_integral<_A1>::value, double>::type
1156 exp2(_A1 __x) {return exp2((double)__x);}
1157
1158 // expm1
1159
1160 using ::expm1;
1161 using ::expm1f;
1162
1163 inline _LIBCPP_INLINE_VISIBILITY float       expm1(float __x)       {return expm1f(__x);}
1164 inline _LIBCPP_INLINE_VISIBILITY long double expm1(long double __x) {return expm1l(__x);}
1165
1166 template <class _A1>
1167 inline _LIBCPP_INLINE_VISIBILITY
1168 typename enable_if<is_integral<_A1>::value, double>::type
1169 expm1(_A1 __x) {return expm1((double)__x);}
1170
1171 // fdim
1172
1173 using ::fdim;
1174 using ::fdimf;
1175
1176 inline _LIBCPP_INLINE_VISIBILITY float       fdim(float __x, float __y)             {return fdimf(__x, __y);}
1177 inline _LIBCPP_INLINE_VISIBILITY long double fdim(long double __x, long double __y) {return fdiml(__x, __y);}
1178
1179 template <class _A1, class _A2>
1180 inline _LIBCPP_INLINE_VISIBILITY
1181 typename enable_if
1182 <
1183     is_arithmetic<_A1>::value &&
1184     is_arithmetic<_A2>::value,
1185     typename __promote<_A1, _A2>::type
1186 >::type
1187 fdim(_A1 __x, _A2 __y)
1188 {
1189     typedef typename __promote<_A1, _A2>::type __result_type;
1190     static_assert((!(is_same<_A1, __result_type>::value &&
1191                       is_same<_A2, __result_type>::value)), "");
1192     return fdim((__result_type)__x, (__result_type)__y);
1193 }
1194
1195 // fma
1196
1197 inline _LIBCPP_INLINE_VISIBILITY float fmaf(float __x, float __y, float __z) {return (float)((double)__x*__y + __z);}
1198 #define FP_FAST_FMAF
1199
1200 using ::fma;
1201
1202 inline _LIBCPP_INLINE_VISIBILITY float       fma(float __x, float __y, float __z)                   {return fmaf(__x, __y, __z);}
1203 inline _LIBCPP_INLINE_VISIBILITY long double fma(long double __x, long double __y, long double __z) {return fmal(__x, __y, __z);}
1204
1205 template <class _A1, class _A2, class _A3>
1206 inline _LIBCPP_INLINE_VISIBILITY
1207 typename enable_if
1208 <
1209     is_arithmetic<_A1>::value &&
1210     is_arithmetic<_A2>::value &&
1211     is_arithmetic<_A3>::value,
1212     typename __promote<_A1, _A2, _A3>::type
1213 >::type
1214 fma(_A1 __x, _A2 __y, _A3 __z)
1215 {
1216     typedef typename __promote<_A1, _A2, _A3>::type __result_type;
1217     static_assert((!(is_same<_A1, __result_type>::value &&
1218                       is_same<_A2, __result_type>::value &&
1219                       is_same<_A3, __result_type>::value)), "");
1220     return fma((__result_type)__x, (__result_type)__y, (__result_type)__z);
1221 }
1222
1223 // fmax
1224
1225 using ::fmax;
1226 using ::fmaxf;
1227
1228 inline _LIBCPP_INLINE_VISIBILITY float       fmax(float __x, float __y)             {return fmaxf(__x, __y);}
1229 inline _LIBCPP_INLINE_VISIBILITY long double fmax(long double __x, long double __y) {return fmaxl(__x, __y);}
1230
1231 template <class _A1, class _A2>
1232 inline _LIBCPP_INLINE_VISIBILITY
1233 typename enable_if
1234 <
1235     is_arithmetic<_A1>::value &&
1236     is_arithmetic<_A2>::value,
1237     typename __promote<_A1, _A2>::type
1238 >::type
1239 fmax(_A1 __x, _A2 __y)
1240 {
1241     typedef typename __promote<_A1, _A2>::type __result_type;
1242     static_assert((!(is_same<_A1, __result_type>::value &&
1243                       is_same<_A2, __result_type>::value)), "");
1244     return fmax((__result_type)__x, (__result_type)__y);
1245 }
1246
1247 // fmin
1248
1249 using ::fmin;
1250 using ::fminf;
1251
1252 inline _LIBCPP_INLINE_VISIBILITY float       fmin(float __x, float __y)             {return fminf(__x, __y);}
1253 inline _LIBCPP_INLINE_VISIBILITY long double fmin(long double __x, long double __y) {return fminl(__x, __y);}
1254
1255 template <class _A1, class _A2>
1256 inline _LIBCPP_INLINE_VISIBILITY
1257 typename enable_if
1258 <
1259     is_arithmetic<_A1>::value &&
1260     is_arithmetic<_A2>::value,
1261     typename __promote<_A1, _A2>::type
1262 >::type
1263 fmin(_A1 __x, _A2 __y)
1264 {
1265     typedef typename __promote<_A1, _A2>::type __result_type;
1266     static_assert((!(is_same<_A1, __result_type>::value &&
1267                       is_same<_A2, __result_type>::value)), "");
1268     return fmin((__result_type)__x, (__result_type)__y);
1269 }
1270
1271 // hypot
1272
1273 using ::hypot;
1274 using ::hypotf;
1275
1276 inline _LIBCPP_INLINE_VISIBILITY float       hypot(float __x, float __y)             {return hypotf(__x, __y);}
1277 inline _LIBCPP_INLINE_VISIBILITY long double hypot(long double __x, long double __y) {return hypotl(__x, __y);}
1278
1279 template <class _A1, class _A2>
1280 inline _LIBCPP_INLINE_VISIBILITY
1281 typename enable_if
1282 <
1283     is_arithmetic<_A1>::value &&
1284     is_arithmetic<_A2>::value,
1285     typename __promote<_A1, _A2>::type
1286 >::type
1287 hypot(_A1 __x, _A2 __y)
1288 {
1289     typedef typename __promote<_A1, _A2>::type __result_type;
1290     static_assert((!(is_same<_A1, __result_type>::value &&
1291                       is_same<_A2, __result_type>::value)), "");
1292     return hypot((__result_type)__x, (__result_type)__y);
1293 }
1294
1295 // ilogb
1296
1297 using ::ilogb;
1298 using ::ilogbf;
1299
1300 inline _LIBCPP_INLINE_VISIBILITY int ilogb(float __x)       {return ilogbf(__x);}
1301 inline _LIBCPP_INLINE_VISIBILITY int ilogb(long double __x) {return ilogbl(__x);}
1302
1303 template <class _A1>
1304 inline _LIBCPP_INLINE_VISIBILITY
1305 typename enable_if<is_integral<_A1>::value, int>::type
1306 ilogb(_A1 __x) {return ilogb((double)__x);}
1307
1308 // lgamma
1309
1310 using ::lgamma;
1311 using ::lgammaf;
1312
1313 inline _LIBCPP_INLINE_VISIBILITY float       lgamma(float __x)       {return lgammaf(__x);}
1314 inline _LIBCPP_INLINE_VISIBILITY long double lgamma(long double __x) {return lgammal(__x);}
1315
1316
1317 template <class _A1>
1318 inline _LIBCPP_INLINE_VISIBILITY
1319 typename enable_if<is_integral<_A1>::value, double>::type
1320 lgamma(_A1 __x) {return lgamma((double)__x);}
1321
1322
1323 // llrint
1324
1325 using ::llrint;
1326 using ::llrintf;
1327
1328 inline _LIBCPP_INLINE_VISIBILITY long long llrint(float __x)       {return llrintf(__x);}
1329 inline _LIBCPP_INLINE_VISIBILITY long long llrint(long double __x) {return llrintl(__x);}
1330
1331 template <class _A1>
1332 inline _LIBCPP_INLINE_VISIBILITY
1333 typename enable_if<is_integral<_A1>::value, long long>::type
1334 llrint(_A1 __x) {return llrint((double)__x);}
1335
1336 // llround
1337
1338 using ::llround;
1339 using ::llroundf;
1340
1341 inline _LIBCPP_INLINE_VISIBILITY long long llround(float __x)       {return llroundf(__x);}
1342 inline _LIBCPP_INLINE_VISIBILITY long long llround(long double __x) {return llroundl(__x);}
1343
1344 template <class _A1>
1345 inline _LIBCPP_INLINE_VISIBILITY
1346 typename enable_if<is_integral<_A1>::value, long long>::type
1347 llround(_A1 __x) {return llround((double)__x);}
1348
1349 // log1p
1350
1351 using ::log1p;
1352 using ::log1pf;
1353
1354 inline _LIBCPP_INLINE_VISIBILITY float       log1p(float __x)       {return log1pf(__x);}
1355 inline _LIBCPP_INLINE_VISIBILITY long double log1p(long double __x) {return log1pl(__x);}
1356
1357 template <class _A1>
1358 inline _LIBCPP_INLINE_VISIBILITY
1359 typename enable_if<is_integral<_A1>::value, double>::type
1360 log1p(_A1 __x) {return log1p((double)__x);}
1361
1362 // log2
1363
1364 using ::log2;
1365 using ::log2f;
1366
1367 inline _LIBCPP_INLINE_VISIBILITY float       log2(float __x)       {return log2f(__x);}
1368 inline _LIBCPP_INLINE_VISIBILITY long double log2(long double __x) {return log2l(__x);}
1369
1370 template <class _A1>
1371 inline _LIBCPP_INLINE_VISIBILITY
1372 typename enable_if<is_integral<_A1>::value, double>::type
1373 log2(_A1 __x) {return log2((double)__x);}
1374
1375 // logb
1376
1377 using ::logb;
1378 using ::logbf;
1379
1380 inline _LIBCPP_INLINE_VISIBILITY float       logb(float __x)       {return logbf(__x);}
1381 inline _LIBCPP_INLINE_VISIBILITY long double logb(long double __x) {return logbl(__x);}
1382
1383 template <class _A1>
1384 inline _LIBCPP_INLINE_VISIBILITY
1385 typename enable_if<is_integral<_A1>::value, double>::type
1386 logb(_A1 __x) {return logb((double)__x);}
1387
1388 // lrint
1389
1390 using ::lrint;
1391 using ::lrintf;
1392
1393 inline _LIBCPP_INLINE_VISIBILITY long lrint(float __x)       {return lrintf(__x);}
1394 inline _LIBCPP_INLINE_VISIBILITY long lrint(long double __x) {return lrintl(__x);}
1395
1396 template <class _A1>
1397 inline _LIBCPP_INLINE_VISIBILITY
1398 typename enable_if<is_integral<_A1>::value, long>::type
1399 lrint(_A1 __x) {return lrint((double)__x);}
1400
1401 // lround
1402
1403 using ::lround;
1404 using ::lroundf;
1405
1406 inline _LIBCPP_INLINE_VISIBILITY long lround(float __x)       {return lroundf(__x);}
1407 inline _LIBCPP_INLINE_VISIBILITY long lround(long double __x) {return lroundl(__x);}
1408
1409 template <class _A1>
1410 inline _LIBCPP_INLINE_VISIBILITY
1411 typename enable_if<is_integral<_A1>::value, long>::type
1412 lround(_A1 __x) {return lround((double)__x);}
1413
1414 // nan
1415 #endif // _MSC_VER
1416 #endif // __sun__
1417 using ::nan;
1418 using ::nanf;
1419 #ifndef __sun__
1420 #ifndef _MSC_VER
1421
1422 // nearbyint
1423
1424 using ::nearbyint;
1425 using ::nearbyintf;
1426
1427 inline _LIBCPP_INLINE_VISIBILITY float       nearbyint(float __x)       {return nearbyintf(__x);}
1428 inline _LIBCPP_INLINE_VISIBILITY long double nearbyint(long double __x) {return nearbyintl(__x);}
1429
1430 template <class _A1>
1431 inline _LIBCPP_INLINE_VISIBILITY
1432 typename enable_if<is_integral<_A1>::value, double>::type
1433 nearbyint(_A1 __x) {return nearbyint((double)__x);}
1434
1435 // nextafter
1436
1437 using ::nextafter;
1438 using ::nextafterf;
1439
1440 inline _LIBCPP_INLINE_VISIBILITY float       nextafter(float __x, float __y)             {return nextafterf(__x, __y);}
1441 inline _LIBCPP_INLINE_VISIBILITY long double nextafter(long double __x, long double __y) {return nextafterl(__x, __y);}
1442
1443 template <class _A1, class _A2>
1444 inline _LIBCPP_INLINE_VISIBILITY
1445 typename enable_if
1446 <
1447     is_arithmetic<_A1>::value &&
1448     is_arithmetic<_A2>::value,
1449     typename __promote<_A1, _A2>::type
1450 >::type
1451 nextafter(_A1 __x, _A2 __y)
1452 {
1453     typedef typename __promote<_A1, _A2>::type __result_type;
1454     static_assert((!(is_same<_A1, __result_type>::value &&
1455                       is_same<_A2, __result_type>::value)), "");
1456     return nextafter((__result_type)__x, (__result_type)__y);
1457 }
1458
1459 // nexttoward
1460
1461 using ::nexttoward;
1462 using ::nexttowardf;
1463
1464 inline _LIBCPP_INLINE_VISIBILITY float       nexttoward(float __x, long double __y)       {return nexttowardf(__x, __y);}
1465 inline _LIBCPP_INLINE_VISIBILITY long double nexttoward(long double __x, long double __y) {return nexttowardl(__x, __y);}
1466
1467 template <class _A1>
1468 inline _LIBCPP_INLINE_VISIBILITY
1469 typename enable_if<is_integral<_A1>::value, double>::type
1470 nexttoward(_A1 __x, long double __y) {return nexttoward((double)__x, __y);}
1471
1472 // remainder
1473
1474 using ::remainder;
1475 using ::remainderf;
1476
1477 inline _LIBCPP_INLINE_VISIBILITY float       remainder(float __x, float __y)             {return remainderf(__x, __y);}
1478 inline _LIBCPP_INLINE_VISIBILITY long double remainder(long double __x, long double __y) {return remainderl(__x, __y);}
1479
1480 template <class _A1, class _A2>
1481 inline _LIBCPP_INLINE_VISIBILITY
1482 typename enable_if
1483 <
1484     is_arithmetic<_A1>::value &&
1485     is_arithmetic<_A2>::value,
1486     typename __promote<_A1, _A2>::type
1487 >::type
1488 remainder(_A1 __x, _A2 __y)
1489 {
1490     typedef typename __promote<_A1, _A2>::type __result_type;
1491     static_assert((!(is_same<_A1, __result_type>::value &&
1492                       is_same<_A2, __result_type>::value)), "");
1493     return remainder((__result_type)__x, (__result_type)__y);
1494 }
1495
1496 // remquo
1497
1498 using ::remquo;
1499 using ::remquof;
1500
1501 inline _LIBCPP_INLINE_VISIBILITY float       remquo(float __x, float __y, int* __z)             {return remquof(__x, __y, __z);}
1502 inline _LIBCPP_INLINE_VISIBILITY long double remquo(long double __x, long double __y, int* __z) {return remquol(__x, __y, __z);}
1503
1504 template <class _A1, class _A2>
1505 inline _LIBCPP_INLINE_VISIBILITY
1506 typename enable_if
1507 <
1508     is_arithmetic<_A1>::value &&
1509     is_arithmetic<_A2>::value,
1510     typename __promote<_A1, _A2>::type
1511 >::type
1512 remquo(_A1 __x, _A2 __y, int* __z)
1513 {
1514     typedef typename __promote<_A1, _A2>::type __result_type;
1515     static_assert((!(is_same<_A1, __result_type>::value &&
1516                       is_same<_A2, __result_type>::value)), "");
1517     return remquo((__result_type)__x, (__result_type)__y, __z);
1518 }
1519
1520 // rint
1521
1522 using ::rint;
1523 using ::rintf;
1524
1525 inline _LIBCPP_INLINE_VISIBILITY float       rint(float __x)       {return rintf(__x);}
1526 inline _LIBCPP_INLINE_VISIBILITY long double rint(long double __x) {return rintl(__x);}
1527
1528 template <class _A1>
1529 inline _LIBCPP_INLINE_VISIBILITY
1530 typename enable_if<is_integral<_A1>::value, double>::type
1531 rint(_A1 __x) {return rint((double)__x);}
1532
1533 // round
1534
1535 using ::round;
1536 using ::roundf;
1537
1538 inline _LIBCPP_INLINE_VISIBILITY float       round(float __x)       {return roundf(__x);}
1539 inline _LIBCPP_INLINE_VISIBILITY long double round(long double __x) {return roundl(__x);}
1540
1541 template <class _A1>
1542 inline _LIBCPP_INLINE_VISIBILITY
1543 typename enable_if<is_integral<_A1>::value, double>::type
1544 round(_A1 __x) {return round((double)__x);}
1545
1546 // scalbln
1547
1548 using ::scalbln;
1549 using ::scalblnf;
1550
1551 inline _LIBCPP_INLINE_VISIBILITY float       scalbln(float __x, long __y)       {return scalblnf(__x, __y);}
1552 inline _LIBCPP_INLINE_VISIBILITY long double scalbln(long double __x, long __y) {return scalblnl(__x, __y);}
1553
1554 template <class _A1>
1555 inline _LIBCPP_INLINE_VISIBILITY
1556 typename enable_if<is_integral<_A1>::value, double>::type
1557 scalbln(_A1 __x, long __y) {return scalbln((double)__x, __y);}
1558
1559 // scalbn
1560
1561 using ::scalbn;
1562 using ::scalbnf;
1563
1564 inline _LIBCPP_INLINE_VISIBILITY float       scalbn(float __x, int __y)       {return scalbnf(__x, __y);}
1565 inline _LIBCPP_INLINE_VISIBILITY long double scalbn(long double __x, int __y) {return scalbnl(__x, __y);}
1566
1567 template <class _A1>
1568 inline _LIBCPP_INLINE_VISIBILITY
1569 typename enable_if<is_integral<_A1>::value, double>::type
1570 scalbn(_A1 __x, int __y) {return scalbn((double)__x, __y);}
1571
1572 // tgamma
1573
1574 using ::tgamma;
1575 using ::tgammaf;
1576
1577 inline _LIBCPP_INLINE_VISIBILITY float       tgamma(float __x)       {return tgammaf(__x);}
1578 inline _LIBCPP_INLINE_VISIBILITY long double tgamma(long double __x) {return tgammal(__x);}
1579
1580 template <class _A1>
1581 inline _LIBCPP_INLINE_VISIBILITY
1582 typename enable_if<is_integral<_A1>::value, double>::type
1583 tgamma(_A1 __x) {return tgamma((double)__x);}
1584
1585 // trunc
1586
1587 using ::trunc;
1588 using ::truncf;
1589
1590 inline _LIBCPP_INLINE_VISIBILITY float       trunc(float __x)       {return truncf(__x);}
1591 inline _LIBCPP_INLINE_VISIBILITY long double trunc(long double __x) {return truncl(__x);}
1592
1593 template <class _A1>
1594 inline _LIBCPP_INLINE_VISIBILITY
1595 typename enable_if<is_integral<_A1>::value, double>::type
1596 trunc(_A1 __x) {return trunc((double)__x);}
1597
1598 #endif // !_MSC_VER
1599
1600 using ::acosl;
1601 using ::asinl;
1602 using ::atanl;
1603 using ::atan2l;
1604 using ::ceill;
1605 using ::cosl;
1606 using ::coshl;
1607 using ::expl;
1608 using ::fabsl;
1609 using ::floorl;
1610 using ::fmodl;
1611 using ::frexpl;
1612 using ::ldexpl;
1613 using ::logl;
1614 using ::log10l;
1615 using ::modfl;
1616 using ::powl;
1617 using ::sinl;
1618 using ::sinhl;
1619 using ::sqrtl;
1620 using ::tanl;
1621 #ifndef _MSC_VER
1622 using ::tanhl;
1623 using ::acoshl;
1624 using ::asinhl;
1625 using ::atanhl;
1626 using ::cbrtl;
1627 #endif  // !_MSC_VER
1628 using ::copysignl;
1629 #ifndef _MSC_VER
1630 using ::erfl;
1631 using ::erfcl;
1632 using ::exp2l;
1633 using ::expm1l;
1634 using ::fdiml;
1635 using ::fmal;
1636 using ::fmaxl;
1637 using ::fminl;
1638 using ::hypotl;
1639 using ::ilogbl;
1640 using ::lgammal;
1641 using ::llrintl;
1642 using ::llroundl;
1643 using ::log1pl;
1644 using ::log2l;
1645 using ::logbl;
1646 using ::lrintl;
1647 using ::lroundl;
1648 using ::nanl;
1649 using ::nearbyintl;
1650 using ::nextafterl;
1651 using ::nexttowardl;
1652 using ::remainderl;
1653 using ::remquol;
1654 using ::rintl;
1655 using ::roundl;
1656 using ::scalblnl;
1657 using ::scalbnl;
1658 using ::tgammal;
1659 using ::truncl;
1660 #endif // !_MSC_VER
1661
1662 #else
1663 using ::lgamma;
1664 using ::lgammaf;
1665 #endif // __sun__
1666 _LIBCPP_END_NAMESPACE_STD
1667
1668 #endif  // _LIBCPP_CMATH