1 // -*- C++ -*- C forwarding header.
3 // Copyright (C) 1997-2015 Free Software Foundation, Inc.
5 // This file is part of the GNU ISO C++ Library. This library is free
6 // software; you can redistribute it and/or modify it under the
7 // terms of the GNU General Public License as published by the
8 // Free Software Foundation; either version 3, or (at your option)
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 // GNU General Public License for more details.
16 // Under Section 7 of GPL version 3, you are granted additional
17 // permissions described in the GCC Runtime Library Exception, version
18 // 3.1, as published by the Free Software Foundation.
20 // You should have received a copy of the GNU General Public License and
21 // a copy of the GCC Runtime Library Exception along with this program;
22 // see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
23 // <http://www.gnu.org/licenses/>.
25 /** @file include/cmath
26 * This is a Standard C++ Library file. You should @c \#include this file
27 * in your programs, rather than any of the @a *.h implementation files.
29 * This is the C++ version of the Standard C Library header @c math.h,
30 * and its contents are (mostly) the same as that header, but are all
31 * contained in the namespace @c std (except for names which are defined
36 // ISO C++ 14882: 26.5 C library
39 #pragma GCC system_header
41 #include <bits/c++config.h>
42 #include <bits/cpp_type_traits.h>
43 #include <ext/type_traits.h>
46 #ifndef _GLIBCXX_CMATH
47 #define _GLIBCXX_CMATH 1
49 // Get rid of those macros defined in <math.h> in lieu of real functions.
75 namespace std _GLIBCXX_VISIBILITY(default)
77 _GLIBCXX_BEGIN_NAMESPACE_VERSION
79 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
80 inline _GLIBCXX_CONSTEXPR double
82 { return __builtin_fabs(__x); }
85 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
86 inline _GLIBCXX_CONSTEXPR float
88 { return __builtin_fabsf(__x); }
90 inline _GLIBCXX_CONSTEXPR long double
92 { return __builtin_fabsl(__x); }
95 template<typename _Tp>
96 inline _GLIBCXX_CONSTEXPR
97 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
100 { return __builtin_fabs(__x); }
104 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
105 inline _GLIBCXX_CONSTEXPR float
107 { return __builtin_acosf(__x); }
109 inline _GLIBCXX_CONSTEXPR long double
110 acos(long double __x)
111 { return __builtin_acosl(__x); }
114 template<typename _Tp>
115 inline _GLIBCXX_CONSTEXPR
116 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
119 { return __builtin_acos(__x); }
123 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
124 inline _GLIBCXX_CONSTEXPR float
126 { return __builtin_asinf(__x); }
128 inline _GLIBCXX_CONSTEXPR long double
129 asin(long double __x)
130 { return __builtin_asinl(__x); }
133 template<typename _Tp>
134 inline _GLIBCXX_CONSTEXPR
135 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
138 { return __builtin_asin(__x); }
142 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
143 inline _GLIBCXX_CONSTEXPR float
145 { return __builtin_atanf(__x); }
147 inline _GLIBCXX_CONSTEXPR long double
148 atan(long double __x)
149 { return __builtin_atanl(__x); }
152 template<typename _Tp>
153 inline _GLIBCXX_CONSTEXPR
154 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
157 { return __builtin_atan(__x); }
161 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
162 inline _GLIBCXX_CONSTEXPR float
163 atan2(float __y, float __x)
164 { return __builtin_atan2f(__y, __x); }
166 inline _GLIBCXX_CONSTEXPR long double
167 atan2(long double __y, long double __x)
168 { return __builtin_atan2l(__y, __x); }
171 template<typename _Tp, typename _Up>
172 inline _GLIBCXX_CONSTEXPR
173 typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
174 atan2(_Tp __y, _Up __x)
176 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
177 return atan2(__type(__y), __type(__x));
182 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
183 inline _GLIBCXX_CONSTEXPR float
185 { return __builtin_ceilf(__x); }
187 inline _GLIBCXX_CONSTEXPR long double
188 ceil(long double __x)
189 { return __builtin_ceill(__x); }
192 template<typename _Tp>
193 inline _GLIBCXX_CONSTEXPR
194 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
197 { return __builtin_ceil(__x); }
201 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
202 inline _GLIBCXX_CONSTEXPR float
204 { return __builtin_cosf(__x); }
206 inline _GLIBCXX_CONSTEXPR long double
208 { return __builtin_cosl(__x); }
211 template<typename _Tp>
212 inline _GLIBCXX_CONSTEXPR
213 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
216 { return __builtin_cos(__x); }
220 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
221 inline _GLIBCXX_CONSTEXPR float
223 { return __builtin_coshf(__x); }
225 inline _GLIBCXX_CONSTEXPR long double
226 cosh(long double __x)
227 { return __builtin_coshl(__x); }
230 template<typename _Tp>
231 inline _GLIBCXX_CONSTEXPR
232 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
235 { return __builtin_cosh(__x); }
239 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
240 inline _GLIBCXX_CONSTEXPR float
242 { return __builtin_expf(__x); }
244 inline _GLIBCXX_CONSTEXPR long double
246 { return __builtin_expl(__x); }
249 template<typename _Tp>
250 inline _GLIBCXX_CONSTEXPR
251 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
254 { return __builtin_exp(__x); }
258 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
259 inline _GLIBCXX_CONSTEXPR float
261 { return __builtin_fabsf(__x); }
263 inline _GLIBCXX_CONSTEXPR long double
264 fabs(long double __x)
265 { return __builtin_fabsl(__x); }
268 template<typename _Tp>
269 inline _GLIBCXX_CONSTEXPR
270 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
273 { return __builtin_fabs(__x); }
277 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
278 inline _GLIBCXX_CONSTEXPR float
280 { return __builtin_floorf(__x); }
282 inline _GLIBCXX_CONSTEXPR long double
283 floor(long double __x)
284 { return __builtin_floorl(__x); }
287 template<typename _Tp>
288 inline _GLIBCXX_CONSTEXPR
289 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
292 { return __builtin_floor(__x); }
296 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
297 inline _GLIBCXX_CONSTEXPR float
298 fmod(float __x, float __y)
299 { return __builtin_fmodf(__x, __y); }
301 inline _GLIBCXX_CONSTEXPR long double
302 fmod(long double __x, long double __y)
303 { return __builtin_fmodl(__x, __y); }
306 template<typename _Tp, typename _Up>
307 inline _GLIBCXX_CONSTEXPR
308 typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
309 fmod(_Tp __x, _Up __y)
311 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
312 return fmod(__type(__x), __type(__y));
317 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
319 frexp(float __x, int* __exp)
320 { return __builtin_frexpf(__x, __exp); }
323 frexp(long double __x, int* __exp)
324 { return __builtin_frexpl(__x, __exp); }
327 template<typename _Tp>
328 inline _GLIBCXX_CONSTEXPR
329 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
331 frexp(_Tp __x, int* __exp)
332 { return __builtin_frexp(__x, __exp); }
336 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
337 inline _GLIBCXX_CONSTEXPR float
338 ldexp(float __x, int __exp)
339 { return __builtin_ldexpf(__x, __exp); }
341 inline _GLIBCXX_CONSTEXPR long double
342 ldexp(long double __x, int __exp)
343 { return __builtin_ldexpl(__x, __exp); }
346 template<typename _Tp>
347 inline _GLIBCXX_CONSTEXPR
348 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
350 ldexp(_Tp __x, int __exp)
351 { return __builtin_ldexp(__x, __exp); }
355 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
356 inline _GLIBCXX_CONSTEXPR float
358 { return __builtin_logf(__x); }
360 inline _GLIBCXX_CONSTEXPR long double
362 { return __builtin_logl(__x); }
365 template<typename _Tp>
366 inline _GLIBCXX_CONSTEXPR
367 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
370 { return __builtin_log(__x); }
374 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
375 inline _GLIBCXX_CONSTEXPR float
377 { return __builtin_log10f(__x); }
379 inline _GLIBCXX_CONSTEXPR long double
380 log10(long double __x)
381 { return __builtin_log10l(__x); }
384 template<typename _Tp>
385 inline _GLIBCXX_CONSTEXPR
386 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
389 { return __builtin_log10(__x); }
393 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
395 modf(float __x, float* __iptr)
396 { return __builtin_modff(__x, __iptr); }
399 modf(long double __x, long double* __iptr)
400 { return __builtin_modfl(__x, __iptr); }
405 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
406 inline _GLIBCXX_CONSTEXPR float
407 pow(float __x, float __y)
408 { return __builtin_powf(__x, __y); }
410 inline _GLIBCXX_CONSTEXPR long double
411 pow(long double __x, long double __y)
412 { return __builtin_powl(__x, __y); }
414 #if __cplusplus < 201103L
415 // _GLIBCXX_RESOLVE_LIB_DEFECTS
416 // DR 550. What should the return type of pow(float,int) be?
418 pow(double __x, int __i)
419 { return __builtin_powi(__x, __i); }
422 pow(float __x, int __n)
423 { return __builtin_powif(__x, __n); }
426 pow(long double __x, int __n)
427 { return __builtin_powil(__x, __n); }
431 template<typename _Tp, typename _Up>
432 inline _GLIBCXX_CONSTEXPR
433 typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
434 pow(_Tp __x, _Up __y)
436 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
437 return pow(__type(__x), __type(__y));
442 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
443 inline _GLIBCXX_CONSTEXPR float
445 { return __builtin_sinf(__x); }
447 inline _GLIBCXX_CONSTEXPR long double
449 { return __builtin_sinl(__x); }
452 template<typename _Tp>
453 inline _GLIBCXX_CONSTEXPR
454 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
457 { return __builtin_sin(__x); }
461 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
462 inline _GLIBCXX_CONSTEXPR float
464 { return __builtin_sinhf(__x); }
466 inline _GLIBCXX_CONSTEXPR long double
467 sinh(long double __x)
468 { return __builtin_sinhl(__x); }
471 template<typename _Tp>
472 inline _GLIBCXX_CONSTEXPR
473 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
476 { return __builtin_sinh(__x); }
480 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
481 inline _GLIBCXX_CONSTEXPR float
483 { return __builtin_sqrtf(__x); }
485 inline _GLIBCXX_CONSTEXPR long double
486 sqrt(long double __x)
487 { return __builtin_sqrtl(__x); }
490 template<typename _Tp>
491 inline _GLIBCXX_CONSTEXPR
492 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
495 { return __builtin_sqrt(__x); }
499 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
500 inline _GLIBCXX_CONSTEXPR float
502 { return __builtin_tanf(__x); }
504 inline _GLIBCXX_CONSTEXPR long double
506 { return __builtin_tanl(__x); }
509 template<typename _Tp>
510 inline _GLIBCXX_CONSTEXPR
511 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
514 { return __builtin_tan(__x); }
518 #ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
519 inline _GLIBCXX_CONSTEXPR float
521 { return __builtin_tanhf(__x); }
523 inline _GLIBCXX_CONSTEXPR long double
524 tanh(long double __x)
525 { return __builtin_tanhl(__x); }
528 template<typename _Tp>
529 inline _GLIBCXX_CONSTEXPR
530 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
533 { return __builtin_tanh(__x); }
535 _GLIBCXX_END_NAMESPACE_VERSION
538 #if _GLIBCXX_USE_C99_MATH
539 #if !_GLIBCXX_USE_C99_FP_MACROS_DYNAMIC
541 // These are possible macros imported from C99-land.
549 #undef isgreaterequal
555 namespace std _GLIBCXX_VISIBILITY(default)
557 _GLIBCXX_BEGIN_NAMESPACE_VERSION
559 #if __cplusplus >= 201103L
561 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
563 fpclassify(float __x)
564 { return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
565 FP_SUBNORMAL, FP_ZERO, __x); }
568 fpclassify(double __x)
569 { return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
570 FP_SUBNORMAL, FP_ZERO, __x); }
573 fpclassify(long double __x)
574 { return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
575 FP_SUBNORMAL, FP_ZERO, __x); }
578 template<typename _Tp>
579 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
582 { return __x != 0 ? FP_NORMAL : FP_ZERO; }
584 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
587 { return __builtin_isfinite(__x); }
591 { return __builtin_isfinite(__x); }
594 isfinite(long double __x)
595 { return __builtin_isfinite(__x); }
598 template<typename _Tp>
599 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
604 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
607 { return __builtin_isinf(__x); }
611 { return __builtin_isinf(__x); }
614 isinf(long double __x)
615 { return __builtin_isinf(__x); }
618 template<typename _Tp>
619 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
624 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
627 { return __builtin_isnan(__x); }
631 { return __builtin_isnan(__x); }
634 isnan(long double __x)
635 { return __builtin_isnan(__x); }
638 template<typename _Tp>
639 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
644 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
647 { return __builtin_isnormal(__x); }
651 { return __builtin_isnormal(__x); }
654 isnormal(long double __x)
655 { return __builtin_isnormal(__x); }
658 template<typename _Tp>
659 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
662 { return __x != 0 ? true : false; }
664 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
665 // Note: middle-end/36757 is fixed, __builtin_signbit is type-generic.
668 { return __builtin_signbit(__x); }
672 { return __builtin_signbit(__x); }
675 signbit(long double __x)
676 { return __builtin_signbit(__x); }
679 template<typename _Tp>
680 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
683 { return __x < 0 ? true : false; }
685 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
687 isgreater(float __x, float __y)
688 { return __builtin_isgreater(__x, __y); }
691 isgreater(double __x, double __y)
692 { return __builtin_isgreater(__x, __y); }
695 isgreater(long double __x, long double __y)
696 { return __builtin_isgreater(__x, __y); }
699 template<typename _Tp, typename _Up>
701 __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
702 && __is_arithmetic<_Up>::__value), bool>::__type
703 isgreater(_Tp __x, _Up __y)
705 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
706 return __builtin_isgreater(__type(__x), __type(__y));
709 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
711 isgreaterequal(float __x, float __y)
712 { return __builtin_isgreaterequal(__x, __y); }
715 isgreaterequal(double __x, double __y)
716 { return __builtin_isgreaterequal(__x, __y); }
719 isgreaterequal(long double __x, long double __y)
720 { return __builtin_isgreaterequal(__x, __y); }
723 template<typename _Tp, typename _Up>
725 __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
726 && __is_arithmetic<_Up>::__value), bool>::__type
727 isgreaterequal(_Tp __x, _Up __y)
729 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
730 return __builtin_isgreaterequal(__type(__x), __type(__y));
733 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
735 isless(float __x, float __y)
736 { return __builtin_isless(__x, __y); }
739 isless(double __x, double __y)
740 { return __builtin_isless(__x, __y); }
743 isless(long double __x, long double __y)
744 { return __builtin_isless(__x, __y); }
747 template<typename _Tp, typename _Up>
749 __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
750 && __is_arithmetic<_Up>::__value), bool>::__type
751 isless(_Tp __x, _Up __y)
753 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
754 return __builtin_isless(__type(__x), __type(__y));
757 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
759 islessequal(float __x, float __y)
760 { return __builtin_islessequal(__x, __y); }
763 islessequal(double __x, double __y)
764 { return __builtin_islessequal(__x, __y); }
767 islessequal(long double __x, long double __y)
768 { return __builtin_islessequal(__x, __y); }
771 template<typename _Tp, typename _Up>
773 __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
774 && __is_arithmetic<_Up>::__value), bool>::__type
775 islessequal(_Tp __x, _Up __y)
777 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
778 return __builtin_islessequal(__type(__x), __type(__y));
781 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
783 islessgreater(float __x, float __y)
784 { return __builtin_islessgreater(__x, __y); }
787 islessgreater(double __x, double __y)
788 { return __builtin_islessgreater(__x, __y); }
791 islessgreater(long double __x, long double __y)
792 { return __builtin_islessgreater(__x, __y); }
795 template<typename _Tp, typename _Up>
797 __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
798 && __is_arithmetic<_Up>::__value), bool>::__type
799 islessgreater(_Tp __x, _Up __y)
801 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
802 return __builtin_islessgreater(__type(__x), __type(__y));
805 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
807 isunordered(float __x, float __y)
808 { return __builtin_isunordered(__x, __y); }
811 isunordered(double __x, double __y)
812 { return __builtin_isunordered(__x, __y); }
815 isunordered(long double __x, long double __y)
816 { return __builtin_isunordered(__x, __y); }
819 template<typename _Tp, typename _Up>
821 __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
822 && __is_arithmetic<_Up>::__value), bool>::__type
823 isunordered(_Tp __x, _Up __y)
825 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
826 return __builtin_isunordered(__type(__x), __type(__y));
831 template<typename _Tp>
832 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
836 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
837 return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
838 FP_SUBNORMAL, FP_ZERO, __type(__f));
841 template<typename _Tp>
842 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
846 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
847 return __builtin_isfinite(__type(__f));
850 template<typename _Tp>
851 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
855 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
856 return __builtin_isinf(__type(__f));
859 template<typename _Tp>
860 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
864 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
865 return __builtin_isnan(__type(__f));
868 template<typename _Tp>
869 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
873 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
874 return __builtin_isnormal(__type(__f));
877 template<typename _Tp>
878 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
882 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
883 return __builtin_signbit(__type(__f));
886 template<typename _Tp>
887 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
889 isgreater(_Tp __f1, _Tp __f2)
891 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
892 return __builtin_isgreater(__type(__f1), __type(__f2));
895 template<typename _Tp>
896 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
898 isgreaterequal(_Tp __f1, _Tp __f2)
900 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
901 return __builtin_isgreaterequal(__type(__f1), __type(__f2));
904 template<typename _Tp>
905 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
907 isless(_Tp __f1, _Tp __f2)
909 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
910 return __builtin_isless(__type(__f1), __type(__f2));
913 template<typename _Tp>
914 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
916 islessequal(_Tp __f1, _Tp __f2)
918 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
919 return __builtin_islessequal(__type(__f1), __type(__f2));
922 template<typename _Tp>
923 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
925 islessgreater(_Tp __f1, _Tp __f2)
927 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
928 return __builtin_islessgreater(__type(__f1), __type(__f2));
931 template<typename _Tp>
932 inline typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value,
934 isunordered(_Tp __f1, _Tp __f2)
936 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
937 return __builtin_isunordered(__type(__f1), __type(__f2));
942 _GLIBCXX_END_NAMESPACE_VERSION
945 #endif /* _GLIBCXX_USE_C99_FP_MACROS_DYNAMIC */
948 #if __cplusplus >= 201103L
950 #ifdef _GLIBCXX_USE_C99_MATH_TR1
1058 namespace std _GLIBCXX_VISIBILITY(default)
1060 _GLIBCXX_BEGIN_NAMESPACE_VERSION
1172 using ::nexttowardf;
1173 using ::nexttowardl;
1207 /// Additional overloads.
1208 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1211 { return __builtin_acoshf(__x); }
1213 constexpr long double
1214 acosh(long double __x)
1215 { return __builtin_acoshl(__x); }
1218 template<typename _Tp>
1219 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1222 { return __builtin_acosh(__x); }
1224 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1227 { return __builtin_asinhf(__x); }
1229 constexpr long double
1230 asinh(long double __x)
1231 { return __builtin_asinhl(__x); }
1234 template<typename _Tp>
1235 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1238 { return __builtin_asinh(__x); }
1240 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1243 { return __builtin_atanhf(__x); }
1245 constexpr long double
1246 atanh(long double __x)
1247 { return __builtin_atanhl(__x); }
1250 template<typename _Tp>
1251 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1254 { return __builtin_atanh(__x); }
1256 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1259 { return __builtin_cbrtf(__x); }
1261 constexpr long double
1262 cbrt(long double __x)
1263 { return __builtin_cbrtl(__x); }
1266 template<typename _Tp>
1267 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1270 { return __builtin_cbrt(__x); }
1272 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1274 copysign(float __x, float __y)
1275 { return __builtin_copysignf(__x, __y); }
1277 constexpr long double
1278 copysign(long double __x, long double __y)
1279 { return __builtin_copysignl(__x, __y); }
1282 template<typename _Tp, typename _Up>
1283 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1284 copysign(_Tp __x, _Up __y)
1286 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1287 return copysign(__type(__x), __type(__y));
1290 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1293 { return __builtin_erff(__x); }
1295 constexpr long double
1296 erf(long double __x)
1297 { return __builtin_erfl(__x); }
1300 template<typename _Tp>
1301 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1304 { return __builtin_erf(__x); }
1306 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1309 { return __builtin_erfcf(__x); }
1311 constexpr long double
1312 erfc(long double __x)
1313 { return __builtin_erfcl(__x); }
1316 template<typename _Tp>
1317 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1320 { return __builtin_erfc(__x); }
1322 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1325 { return __builtin_exp2f(__x); }
1327 constexpr long double
1328 exp2(long double __x)
1329 { return __builtin_exp2l(__x); }
1332 template<typename _Tp>
1333 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1336 { return __builtin_exp2(__x); }
1338 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1341 { return __builtin_expm1f(__x); }
1343 constexpr long double
1344 expm1(long double __x)
1345 { return __builtin_expm1l(__x); }
1348 template<typename _Tp>
1349 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1352 { return __builtin_expm1(__x); }
1354 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1356 fdim(float __x, float __y)
1357 { return __builtin_fdimf(__x, __y); }
1359 constexpr long double
1360 fdim(long double __x, long double __y)
1361 { return __builtin_fdiml(__x, __y); }
1364 template<typename _Tp, typename _Up>
1365 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1366 fdim(_Tp __x, _Up __y)
1368 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1369 return fdim(__type(__x), __type(__y));
1372 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1374 fma(float __x, float __y, float __z)
1375 { return __builtin_fmaf(__x, __y, __z); }
1377 constexpr long double
1378 fma(long double __x, long double __y, long double __z)
1379 { return __builtin_fmal(__x, __y, __z); }
1382 template<typename _Tp, typename _Up, typename _Vp>
1383 constexpr typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type
1384 fma(_Tp __x, _Up __y, _Vp __z)
1386 typedef typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type __type;
1387 return fma(__type(__x), __type(__y), __type(__z));
1390 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1392 fmax(float __x, float __y)
1393 { return __builtin_fmaxf(__x, __y); }
1395 constexpr long double
1396 fmax(long double __x, long double __y)
1397 { return __builtin_fmaxl(__x, __y); }
1400 template<typename _Tp, typename _Up>
1401 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1402 fmax(_Tp __x, _Up __y)
1404 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1405 return fmax(__type(__x), __type(__y));
1408 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1410 fmin(float __x, float __y)
1411 { return __builtin_fminf(__x, __y); }
1413 constexpr long double
1414 fmin(long double __x, long double __y)
1415 { return __builtin_fminl(__x, __y); }
1418 template<typename _Tp, typename _Up>
1419 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1420 fmin(_Tp __x, _Up __y)
1422 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1423 return fmin(__type(__x), __type(__y));
1426 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1428 hypot(float __x, float __y)
1429 { return __builtin_hypotf(__x, __y); }
1431 constexpr long double
1432 hypot(long double __x, long double __y)
1433 { return __builtin_hypotl(__x, __y); }
1436 template<typename _Tp, typename _Up>
1437 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1438 hypot(_Tp __x, _Up __y)
1440 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1441 return hypot(__type(__x), __type(__y));
1444 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1447 { return __builtin_ilogbf(__x); }
1450 ilogb(long double __x)
1451 { return __builtin_ilogbl(__x); }
1454 template<typename _Tp>
1456 typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1459 { return __builtin_ilogb(__x); }
1461 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1464 { return __builtin_lgammaf(__x); }
1466 constexpr long double
1467 lgamma(long double __x)
1468 { return __builtin_lgammal(__x); }
1471 template<typename _Tp>
1472 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1475 { return __builtin_lgamma(__x); }
1477 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1480 { return __builtin_llrintf(__x); }
1483 llrint(long double __x)
1484 { return __builtin_llrintl(__x); }
1487 template<typename _Tp>
1488 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1491 { return __builtin_llrint(__x); }
1493 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1496 { return __builtin_llroundf(__x); }
1499 llround(long double __x)
1500 { return __builtin_llroundl(__x); }
1503 template<typename _Tp>
1504 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1507 { return __builtin_llround(__x); }
1509 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1512 { return __builtin_log1pf(__x); }
1514 constexpr long double
1515 log1p(long double __x)
1516 { return __builtin_log1pl(__x); }
1519 template<typename _Tp>
1520 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1523 { return __builtin_log1p(__x); }
1525 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1529 { return __builtin_log2f(__x); }
1531 constexpr long double
1532 log2(long double __x)
1533 { return __builtin_log2l(__x); }
1536 template<typename _Tp>
1537 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1540 { return __builtin_log2(__x); }
1542 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1545 { return __builtin_logbf(__x); }
1547 constexpr long double
1548 logb(long double __x)
1549 { return __builtin_logbl(__x); }
1552 template<typename _Tp>
1553 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1556 { return __builtin_logb(__x); }
1558 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1561 { return __builtin_lrintf(__x); }
1564 lrint(long double __x)
1565 { return __builtin_lrintl(__x); }
1568 template<typename _Tp>
1569 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1572 { return __builtin_lrint(__x); }
1574 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1577 { return __builtin_lroundf(__x); }
1580 lround(long double __x)
1581 { return __builtin_lroundl(__x); }
1584 template<typename _Tp>
1585 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1588 { return __builtin_lround(__x); }
1590 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1592 nearbyint(float __x)
1593 { return __builtin_nearbyintf(__x); }
1595 constexpr long double
1596 nearbyint(long double __x)
1597 { return __builtin_nearbyintl(__x); }
1600 template<typename _Tp>
1601 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1604 { return __builtin_nearbyint(__x); }
1606 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1608 nextafter(float __x, float __y)
1609 { return __builtin_nextafterf(__x, __y); }
1611 constexpr long double
1612 nextafter(long double __x, long double __y)
1613 { return __builtin_nextafterl(__x, __y); }
1616 template<typename _Tp, typename _Up>
1617 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1618 nextafter(_Tp __x, _Up __y)
1620 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1621 return nextafter(__type(__x), __type(__y));
1624 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1626 nexttoward(float __x, long double __y)
1627 { return __builtin_nexttowardf(__x, __y); }
1629 constexpr long double
1630 nexttoward(long double __x, long double __y)
1631 { return __builtin_nexttowardl(__x, __y); }
1634 template<typename _Tp>
1635 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1637 nexttoward(_Tp __x, long double __y)
1638 { return __builtin_nexttoward(__x, __y); }
1640 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1642 remainder(float __x, float __y)
1643 { return __builtin_remainderf(__x, __y); }
1645 constexpr long double
1646 remainder(long double __x, long double __y)
1647 { return __builtin_remainderl(__x, __y); }
1650 template<typename _Tp, typename _Up>
1651 constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1652 remainder(_Tp __x, _Up __y)
1654 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1655 return remainder(__type(__x), __type(__y));
1658 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1660 remquo(float __x, float __y, int* __pquo)
1661 { return __builtin_remquof(__x, __y, __pquo); }
1664 remquo(long double __x, long double __y, int* __pquo)
1665 { return __builtin_remquol(__x, __y, __pquo); }
1668 template<typename _Tp, typename _Up>
1669 inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
1670 remquo(_Tp __x, _Up __y, int* __pquo)
1672 typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
1673 return remquo(__type(__x), __type(__y), __pquo);
1676 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1679 { return __builtin_rintf(__x); }
1681 constexpr long double
1682 rint(long double __x)
1683 { return __builtin_rintl(__x); }
1686 template<typename _Tp>
1687 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1690 { return __builtin_rint(__x); }
1692 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1695 { return __builtin_roundf(__x); }
1697 constexpr long double
1698 round(long double __x)
1699 { return __builtin_roundl(__x); }
1702 template<typename _Tp>
1703 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1706 { return __builtin_round(__x); }
1708 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1710 scalbln(float __x, long __ex)
1711 { return __builtin_scalblnf(__x, __ex); }
1713 constexpr long double
1714 scalbln(long double __x, long __ex)
1715 { return __builtin_scalblnl(__x, __ex); }
1718 template<typename _Tp>
1719 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1721 scalbln(_Tp __x, long __ex)
1722 { return __builtin_scalbln(__x, __ex); }
1724 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1726 scalbn(float __x, int __ex)
1727 { return __builtin_scalbnf(__x, __ex); }
1729 constexpr long double
1730 scalbn(long double __x, int __ex)
1731 { return __builtin_scalbnl(__x, __ex); }
1734 template<typename _Tp>
1735 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1737 scalbn(_Tp __x, int __ex)
1738 { return __builtin_scalbn(__x, __ex); }
1740 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1743 { return __builtin_tgammaf(__x); }
1745 constexpr long double
1746 tgamma(long double __x)
1747 { return __builtin_tgammal(__x); }
1750 template<typename _Tp>
1751 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1754 { return __builtin_tgamma(__x); }
1756 #ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO
1759 { return __builtin_truncf(__x); }
1761 constexpr long double
1762 trunc(long double __x)
1763 { return __builtin_truncl(__x); }
1766 template<typename _Tp>
1767 constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
1770 { return __builtin_trunc(__x); }
1772 _GLIBCXX_END_NAMESPACE_VERSION
1775 #endif // _GLIBCXX_USE_C99_MATH_TR1