]> git.ipfire.org Git - thirdparty/gcc.git/blobdiff - libstdc++-v3/include/c_global/cmath
Update copyright years.
[thirdparty/gcc.git] / libstdc++-v3 / include / c_global / cmath
index a333eb583c058dcd3905e1eaa022a868437581f9..b99aaf8df4041fbc9b7d87241ce100e332960bf4 100644 (file)
@@ -1,8 +1,6 @@
 // -*- C++ -*- C forwarding header.
 
-// Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005,
-// 2006, 2007, 2008, 2009, 2010
-// Free Software Foundation, Inc.
+// Copyright (C) 1997-2020 Free Software Foundation, Inc.
 //
 // This file is part of the GNU ISO C++ Library.  This library is free
 // software; you can redistribute it and/or modify it under the
 #include <bits/c++config.h>
 #include <bits/cpp_type_traits.h>
 #include <ext/type_traits.h>
-#include <math.h>
+#define _GLIBCXX_INCLUDE_NEXT_C_HEADERS
+#include_next <math.h>
+#undef _GLIBCXX_INCLUDE_NEXT_C_HEADERS
+#include <bits/std_abs.h>
 
 #ifndef _GLIBCXX_CMATH
 #define _GLIBCXX_CMATH 1
 
 // Get rid of those macros defined in <math.h> in lieu of real functions.
-#undef abs
 #undef div
 #undef acos
 #undef asin
 #undef tan
 #undef tanh
 
+extern "C++"
+{
 namespace std _GLIBCXX_VISIBILITY(default)
 {
 _GLIBCXX_BEGIN_NAMESPACE_VERSION
 
-  inline double
-  abs(double __x)
-  { return __builtin_fabs(__x); }
-
-  inline float
-  abs(float __x)
-  { return __builtin_fabsf(__x); }
-
-  inline long double
-  abs(long double __x)
-  { return __builtin_fabsl(__x); }
-
-  template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
-    abs(_Tp __x)
-    { return __builtin_fabs(__x); }
-
   using ::acos;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   acos(float __x)
   { return __builtin_acosf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   acos(long double __x)
   { return __builtin_acosl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     acos(_Tp __x)
     { return __builtin_acos(__x); }
 
   using ::asin;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   asin(float __x)
   { return __builtin_asinf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   asin(long double __x)
   { return __builtin_asinl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     asin(_Tp __x)
     { return __builtin_asin(__x); }
 
   using ::atan;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   atan(float __x)
   { return __builtin_atanf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   atan(long double __x)
   { return __builtin_atanl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     atan(_Tp __x)
     { return __builtin_atan(__x); }
 
   using ::atan2;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   atan2(float __y, float __x)
   { return __builtin_atan2f(__y, __x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   atan2(long double __y, long double __x)
   { return __builtin_atan2l(__y, __x); }
+#endif
 
   template<typename _Tp, typename _Up>
-    inline
-    typename __gnu_cxx::__promote_2<
-    typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value
-                                   && __is_arithmetic<_Up>::__value,
-                                   _Tp>::__type, _Up>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     atan2(_Tp __y, _Up __x)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
@@ -168,112 +160,142 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
 
   using ::ceil;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   ceil(float __x)
   { return __builtin_ceilf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   ceil(long double __x)
   { return __builtin_ceill(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     ceil(_Tp __x)
     { return __builtin_ceil(__x); }
 
   using ::cos;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   cos(float __x)
   { return __builtin_cosf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   cos(long double __x)
   { return __builtin_cosl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     cos(_Tp __x)
     { return __builtin_cos(__x); }
 
   using ::cosh;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   cosh(float __x)
   { return __builtin_coshf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   cosh(long double __x)
   { return __builtin_coshl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     cosh(_Tp __x)
     { return __builtin_cosh(__x); }
 
   using ::exp;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   exp(float __x)
   { return __builtin_expf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   exp(long double __x)
   { return __builtin_expl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     exp(_Tp __x)
     { return __builtin_exp(__x); }
 
   using ::fabs;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   fabs(float __x)
   { return __builtin_fabsf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   fabs(long double __x)
   { return __builtin_fabsl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     fabs(_Tp __x)
     { return __builtin_fabs(__x); }
 
   using ::floor;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   floor(float __x)
   { return __builtin_floorf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   floor(long double __x)
   { return __builtin_floorl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     floor(_Tp __x)
     { return __builtin_floor(__x); }
 
   using ::fmod;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   fmod(float __x, float __y)
   { return __builtin_fmodf(__x, __y); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   fmod(long double __x, long double __y)
   { return __builtin_fmodl(__x, __y); }
+#endif
+
+  template<typename _Tp, typename _Up>
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    fmod(_Tp __x, _Up __y)
+    {
+      typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
+      return fmod(__type(__x), __type(__y));
+    }
 
   using ::frexp;
 
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
   inline float
   frexp(float __x, int* __exp)
   { return __builtin_frexpf(__x, __exp); }
@@ -281,63 +303,75 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   inline long double
   frexp(long double __x, int* __exp)
   { return __builtin_frexpl(__x, __exp); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     frexp(_Tp __x, int* __exp)
     { return __builtin_frexp(__x, __exp); }
 
   using ::ldexp;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   ldexp(float __x, int __exp)
   { return __builtin_ldexpf(__x, __exp); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   ldexp(long double __x, int __exp)
   { return __builtin_ldexpl(__x, __exp); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
-  ldexp(_Tp __x, int __exp)
-  { return __builtin_ldexp(__x, __exp); }
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
+    ldexp(_Tp __x, int __exp)
+    { return __builtin_ldexp(__x, __exp); }
 
   using ::log;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   log(float __x)
   { return __builtin_logf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   log(long double __x)
   { return __builtin_logl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     log(_Tp __x)
     { return __builtin_log(__x); }
 
   using ::log10;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   log10(float __x)
   { return __builtin_log10f(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   log10(long double __x)
   { return __builtin_log10l(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                    double>::__type
     log10(_Tp __x)
     { return __builtin_log10(__x); }
 
   using ::modf;
 
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
   inline float
   modf(float __x, float* __iptr)
   { return __builtin_modff(__x, __iptr); }
@@ -345,18 +379,20 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   inline long double
   modf(long double __x, long double* __iptr)
   { return __builtin_modfl(__x, __iptr); }
+#endif
 
   using ::pow;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   pow(float __x, float __y)
   { return __builtin_powf(__x, __y); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   pow(long double __x, long double __y)
   { return __builtin_powl(__x, __y); }
 
-#ifndef __GXX_EXPERIMENTAL_CXX0X__
+#if __cplusplus < 201103L
   // _GLIBCXX_RESOLVE_LIB_DEFECTS
   // DR 550. What should the return type of pow(float,int) be?
   inline double
@@ -370,14 +406,12 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   inline long double
   pow(long double __x, int __n)
   { return __builtin_powil(__x, __n); }
+#endif
 #endif
 
   template<typename _Tp, typename _Up>
-    inline
-    typename __gnu_cxx::__promote_2<
-    typename __gnu_cxx::__enable_if<__is_arithmetic<_Tp>::__value
-                                   && __is_arithmetic<_Up>::__value,
-                                   _Tp>::__type, _Up>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     pow(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
@@ -386,87 +420,99 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
 
   using ::sin;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   sin(float __x)
   { return __builtin_sinf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   sin(long double __x)
   { return __builtin_sinl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     sin(_Tp __x)
     { return __builtin_sin(__x); }
 
   using ::sinh;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   sinh(float __x)
   { return __builtin_sinhf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   sinh(long double __x)
   { return __builtin_sinhl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     sinh(_Tp __x)
     { return __builtin_sinh(__x); }
 
   using ::sqrt;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   sqrt(float __x)
   { return __builtin_sqrtf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   sqrt(long double __x)
   { return __builtin_sqrtl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     sqrt(_Tp __x)
     { return __builtin_sqrt(__x); }
 
   using ::tan;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   tan(float __x)
   { return __builtin_tanf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   tan(long double __x)
   { return __builtin_tanl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     tan(_Tp __x)
     { return __builtin_tan(__x); }
 
   using ::tanh;
 
-  inline float
+#ifndef __CORRECT_ISO_CPP_MATH_H_PROTO
+  inline _GLIBCXX_CONSTEXPR float
   tanh(float __x)
   { return __builtin_tanhf(__x); }
 
-  inline long double
+  inline _GLIBCXX_CONSTEXPR long double
   tanh(long double __x)
   { return __builtin_tanhl(__x); }
+#endif
 
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
-                                          double>::__type
+    inline _GLIBCXX_CONSTEXPR
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    double>::__type
     tanh(_Tp __x)
     { return __builtin_tanh(__x); }
 
-_GLIBCXX_END_NAMESPACE_VERSION
-} // namespace
-
 #if _GLIBCXX_USE_C99_MATH
 #if !_GLIBCXX_USE_C99_FP_MACROS_DYNAMIC
 
@@ -484,136 +530,171 @@ _GLIBCXX_END_NAMESPACE_VERSION
 #undef islessgreater
 #undef isunordered
 
-namespace std _GLIBCXX_VISIBILITY(default)
-{
-_GLIBCXX_BEGIN_NAMESPACE_VERSION
+#if __cplusplus >= 201103L
 
-#ifdef __GXX_EXPERIMENTAL_CXX0X__
-  inline int
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr int
   fpclassify(float __x)
   { return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
                                FP_SUBNORMAL, FP_ZERO, __x); }
 
-  inline int
+  constexpr int
   fpclassify(double __x)
   { return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
                                FP_SUBNORMAL, FP_ZERO, __x); }
 
-  inline int
+  constexpr int
   fpclassify(long double __x)
   { return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL,
                                FP_SUBNORMAL, FP_ZERO, __x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          int>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                              int>::__type
     fpclassify(_Tp __x)
     { return __x != 0 ? FP_NORMAL : FP_ZERO; }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isfinite(float __x)
   { return __builtin_isfinite(__x); }
 
-  inline bool
+  constexpr bool
   isfinite(double __x)
   { return __builtin_isfinite(__x); }
 
-  inline bool
+  constexpr bool
   isfinite(long double __x)
   { return __builtin_isfinite(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          bool>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                              bool>::__type
     isfinite(_Tp __x)
     { return true; }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isinf(float __x)
   { return __builtin_isinf(__x); }
 
-  inline bool
+#if _GLIBCXX_HAVE_OBSOLETE_ISINF \
+  && !_GLIBCXX_NO_OBSOLETE_ISINF_ISNAN_DYNAMIC
+  using ::isinf;
+#else
+  constexpr bool
   isinf(double __x)
   { return __builtin_isinf(__x); }
+#endif
 
-  inline bool
+  constexpr bool
   isinf(long double __x)
   { return __builtin_isinf(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          bool>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                              bool>::__type
     isinf(_Tp __x)
     { return false; }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isnan(float __x)
   { return __builtin_isnan(__x); }
 
-  inline bool
+#if _GLIBCXX_HAVE_OBSOLETE_ISNAN \
+  && !_GLIBCXX_NO_OBSOLETE_ISINF_ISNAN_DYNAMIC
+  using ::isnan;
+#else
+  constexpr bool
   isnan(double __x)
   { return __builtin_isnan(__x); }
+#endif
 
-  inline bool
+  constexpr bool
   isnan(long double __x)
   { return __builtin_isnan(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          bool>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                              bool>::__type
     isnan(_Tp __x)
     { return false; }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isnormal(float __x)
   { return __builtin_isnormal(__x); }
 
-  inline bool
+  constexpr bool
   isnormal(double __x)
   { return __builtin_isnormal(__x); }
 
-  inline bool
+  constexpr bool
   isnormal(long double __x)
   { return __builtin_isnormal(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          bool>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                              bool>::__type
     isnormal(_Tp __x)
     { return __x != 0 ? true : false; }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  // Note: middle-end/36757 is fixed, __builtin_signbit is type-generic.
+  constexpr bool
   signbit(float __x)
   { return __builtin_signbit(__x); }
 
-  inline bool
+  constexpr bool
   signbit(double __x)
   { return __builtin_signbit(__x); }
 
-  inline bool
+  constexpr bool
   signbit(long double __x)
   { return __builtin_signbit(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
-                                          bool>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value,
+                                              bool>::__type
     signbit(_Tp __x)
     { return __x < 0 ? true : false; }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isgreater(float __x, float __y)
   { return __builtin_isgreater(__x, __y); }
 
-  inline bool
+  constexpr bool
   isgreater(double __x, double __y)
   { return __builtin_isgreater(__x, __y); }
 
-  inline bool
+  constexpr bool
   isgreater(long double __x, long double __y)
   { return __builtin_isgreater(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename
+    constexpr typename
     __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
                            && __is_arithmetic<_Up>::__value), bool>::__type
     isgreater(_Tp __x, _Up __y)
@@ -621,21 +702,25 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return __builtin_isgreater(__type(__x), __type(__y));
     }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isgreaterequal(float __x, float __y)
   { return __builtin_isgreaterequal(__x, __y); }
 
-  inline bool
+  constexpr bool
   isgreaterequal(double __x, double __y)
   { return __builtin_isgreaterequal(__x, __y); }
 
-  inline bool
+  constexpr bool
   isgreaterequal(long double __x, long double __y)
   { return __builtin_isgreaterequal(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename
+    constexpr typename
     __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
                            && __is_arithmetic<_Up>::__value), bool>::__type
     isgreaterequal(_Tp __x, _Up __y)
@@ -643,21 +728,25 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return __builtin_isgreaterequal(__type(__x), __type(__y));
     }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isless(float __x, float __y)
   { return __builtin_isless(__x, __y); }
 
-  inline bool
+  constexpr bool
   isless(double __x, double __y)
   { return __builtin_isless(__x, __y); }
 
-  inline bool
+  constexpr bool
   isless(long double __x, long double __y)
   { return __builtin_isless(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename
+    constexpr typename
     __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
                            && __is_arithmetic<_Up>::__value), bool>::__type
     isless(_Tp __x, _Up __y)
@@ -665,21 +754,25 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return __builtin_isless(__type(__x), __type(__y));
     }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   islessequal(float __x, float __y)
   { return __builtin_islessequal(__x, __y); }
 
-  inline bool
+  constexpr bool
   islessequal(double __x, double __y)
   { return __builtin_islessequal(__x, __y); }
 
-  inline bool
+  constexpr bool
   islessequal(long double __x, long double __y)
   { return __builtin_islessequal(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename
+    constexpr typename
     __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
                            && __is_arithmetic<_Up>::__value), bool>::__type
     islessequal(_Tp __x, _Up __y)
@@ -687,21 +780,25 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return __builtin_islessequal(__type(__x), __type(__y));
     }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   islessgreater(float __x, float __y)
   { return __builtin_islessgreater(__x, __y); }
 
-  inline bool
+  constexpr bool
   islessgreater(double __x, double __y)
   { return __builtin_islessgreater(__x, __y); }
 
-  inline bool
+  constexpr bool
   islessgreater(long double __x, long double __y)
   { return __builtin_islessgreater(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename
+    constexpr typename
     __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
                            && __is_arithmetic<_Up>::__value), bool>::__type
     islessgreater(_Tp __x, _Up __y)
@@ -709,21 +806,25 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return __builtin_islessgreater(__type(__x), __type(__y));
     }
+#endif
 
-  inline bool
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr bool
   isunordered(float __x, float __y)
   { return __builtin_isunordered(__x, __y); }
 
-  inline bool
+  constexpr bool
   isunordered(double __x, double __y)
   { return __builtin_isunordered(__x, __y); }
 
-  inline bool
+  constexpr bool
   isunordered(long double __x, long double __y)
   { return __builtin_isunordered(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename
+    constexpr typename
     __gnu_cxx::__enable_if<(__is_arithmetic<_Tp>::__value
                            && __is_arithmetic<_Up>::__value), bool>::__type
     isunordered(_Tp __x, _Up __y)
@@ -731,6 +832,7 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return __builtin_isunordered(__type(__x), __type(__y));
     }
+#endif
 
 #else
 
@@ -843,15 +945,11 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       return __builtin_isunordered(__type(__f1), __type(__f2));
     }
 
-#endif
-
-_GLIBCXX_END_NAMESPACE_VERSION
-} // namespace
-
+#endif // C++11
 #endif /* _GLIBCXX_USE_C99_FP_MACROS_DYNAMIC */
-#endif
+#endif /* _GLIBCXX_USE_C99_MATH */
 
-#ifdef __GXX_EXPERIMENTAL_CXX0X__
+#if __cplusplus >= 201103L
 
 #ifdef _GLIBCXX_USE_C99_MATH_TR1
 
@@ -903,12 +1001,14 @@ _GLIBCXX_END_NAMESPACE_VERSION
 #undef lgamma
 #undef lgammaf
 #undef lgammal
+#ifndef _GLIBCXX_NO_C99_ROUNDING_FUNCS
 #undef llrint
 #undef llrintf
 #undef llrintl
 #undef llround
 #undef llroundf
 #undef llroundl
+#endif
 #undef log1p
 #undef log1pf
 #undef log1pl
@@ -961,10 +1061,6 @@ _GLIBCXX_END_NAMESPACE_VERSION
 #undef truncf
 #undef truncl
 
-namespace std _GLIBCXX_VISIBILITY(default)
-{
-_GLIBCXX_BEGIN_NAMESPACE_VERSION
-
   // types
   using ::double_t;
   using ::float_t;
@@ -1034,6 +1130,7 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   using ::lgammaf;
   using ::lgammal;
 
+#ifndef _GLIBCXX_NO_C99_ROUNDING_FUNCS
   using ::llrint;
   using ::llrintf;
   using ::llrintl;
@@ -1041,6 +1138,7 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   using ::llround;
   using ::llroundf;
   using ::llroundl;
+#endif
 
   using ::log1p;
   using ::log1pf;
@@ -1111,439 +1209,511 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   using ::truncl;
 
   /// Additional overloads.
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   acosh(float __x)
   { return __builtin_acoshf(__x); }
 
-  inline long double
+  constexpr long double
   acosh(long double __x)
   { return __builtin_acoshl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     acosh(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return acosh(__type(__x));
-    }
+    { return __builtin_acosh(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   asinh(float __x)
   { return __builtin_asinhf(__x); }
 
-  inline long double
+  constexpr long double
   asinh(long double __x)
   { return __builtin_asinhl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     asinh(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return asinh(__type(__x));
-    }
+    { return __builtin_asinh(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   atanh(float __x)
   { return __builtin_atanhf(__x); }
 
-  inline long double
+  constexpr long double
   atanh(long double __x)
   { return __builtin_atanhl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     atanh(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return atanh(__type(__x));
-    }
+    { return __builtin_atanh(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   cbrt(float __x)
   { return __builtin_cbrtf(__x); }
 
-  inline long double
+  constexpr long double
   cbrt(long double __x)
   { return __builtin_cbrtl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     cbrt(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return cbrt(__type(__x));
-    }
+    { return __builtin_cbrt(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   copysign(float __x, float __y)
   { return __builtin_copysignf(__x, __y); }
 
-  inline long double
+  constexpr long double
   copysign(long double __x, long double __y)
   { return __builtin_copysignl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     copysign(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return copysign(__type(__x), __type(__y));
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   erf(float __x)
   { return __builtin_erff(__x); }
 
-  inline long double
+  constexpr long double
   erf(long double __x)
   { return __builtin_erfl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     erf(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return erf(__type(__x));
-    }
+    { return __builtin_erf(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   erfc(float __x)
   { return __builtin_erfcf(__x); }
 
-  inline long double
+  constexpr long double
   erfc(long double __x)
   { return __builtin_erfcl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     erfc(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return erfc(__type(__x));
-    }
+    { return __builtin_erfc(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   exp2(float __x)
   { return __builtin_exp2f(__x); }
 
-  inline long double
+  constexpr long double
   exp2(long double __x)
   { return __builtin_exp2l(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     exp2(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return exp2(__type(__x));
-    }
+    { return __builtin_exp2(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   expm1(float __x)
   { return __builtin_expm1f(__x); }
 
-  inline long double
+  constexpr long double
   expm1(long double __x)
   { return __builtin_expm1l(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     expm1(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return expm1(__type(__x));
-    }
+    { return __builtin_expm1(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   fdim(float __x, float __y)
   { return __builtin_fdimf(__x, __y); }
 
-  inline long double
+  constexpr long double
   fdim(long double __x, long double __y)
   { return __builtin_fdiml(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     fdim(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return fdim(__type(__x), __type(__y));
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   fma(float __x, float __y, float __z)
   { return __builtin_fmaf(__x, __y, __z); }
 
-  inline long double
+  constexpr long double
   fma(long double __x, long double __y, long double __z)
   { return __builtin_fmal(__x, __y, __z); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up, typename _Vp>
-    inline typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type
+    constexpr typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type
     fma(_Tp __x, _Up __y, _Vp __z)
     {
       typedef typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type __type;
       return fma(__type(__x), __type(__y), __type(__z));
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   fmax(float __x, float __y)
   { return __builtin_fmaxf(__x, __y); }
 
-  inline long double
+  constexpr long double
   fmax(long double __x, long double __y)
   { return __builtin_fmaxl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     fmax(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return fmax(__type(__x), __type(__y));
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   fmin(float __x, float __y)
   { return __builtin_fminf(__x, __y); }
 
-  inline long double
+  constexpr long double
   fmin(long double __x, long double __y)
   { return __builtin_fminl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     fmin(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return fmin(__type(__x), __type(__y));
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   hypot(float __x, float __y)
   { return __builtin_hypotf(__x, __y); }
 
-  inline long double
+  constexpr long double
   hypot(long double __x, long double __y)
   { return __builtin_hypotl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     hypot(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return hypot(__type(__x), __type(__y));
     }
+#endif
 
-  inline int
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr int
   ilogb(float __x)
   { return __builtin_ilogbf(__x); }
 
-  inline int
+  constexpr int
   ilogb(long double __x)
   { return __builtin_ilogbl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline int
+    constexpr
+    typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                    int>::__type
     ilogb(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return ilogb(__type(__x));
-    }
+    { return __builtin_ilogb(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   lgamma(float __x)
   { return __builtin_lgammaf(__x); }
 
-  inline long double
+  constexpr long double
   lgamma(long double __x)
   { return __builtin_lgammal(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     lgamma(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return lgamma(__type(__x));
-    }
+    { return __builtin_lgamma(__x); }
+#endif
 
-  inline long long
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr long long
   llrint(float __x)
   { return __builtin_llrintf(__x); }
 
-  inline long long
+  constexpr long long
   llrint(long double __x)
   { return __builtin_llrintl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline long long
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              long long>::__type
     llrint(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return llrint(__type(__x));
-    }
+    { return __builtin_llrint(__x); }
+#endif
 
-  inline long long
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr long long
   llround(float __x)
   { return __builtin_llroundf(__x); }
 
-  inline long long
+  constexpr long long
   llround(long double __x)
   { return __builtin_llroundl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline long long
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              long long>::__type
     llround(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return llround(__type(__x));
-    }
+    { return __builtin_llround(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   log1p(float __x)
   { return __builtin_log1pf(__x); }
 
-  inline long double
+  constexpr long double
   log1p(long double __x)
   { return __builtin_log1pl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     log1p(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return log1p(__type(__x));
-    }
+    { return __builtin_log1p(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
   // DR 568.
-  inline float
+  constexpr float
   log2(float __x)
   { return __builtin_log2f(__x); }
 
-  inline long double
+  constexpr long double
   log2(long double __x)
   { return __builtin_log2l(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     log2(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return log2(__type(__x));
-    }
+    { return __builtin_log2(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   logb(float __x)
   { return __builtin_logbf(__x); }
 
-  inline long double
+  constexpr long double
   logb(long double __x)
   { return __builtin_logbl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     logb(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return logb(__type(__x));
-    }
+    { return __builtin_logb(__x); }
+#endif
 
-  inline long
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr long
   lrint(float __x)
   { return __builtin_lrintf(__x); }
 
-  inline long
+  constexpr long
   lrint(long double __x)
   { return __builtin_lrintl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline long
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              long>::__type
     lrint(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return lrint(__type(__x));
-    }
+    { return __builtin_lrint(__x); }
+#endif
 
-  inline long
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr long
   lround(float __x)
   { return __builtin_lroundf(__x); }
 
-  inline long
+  constexpr long
   lround(long double __x)
   { return __builtin_lroundl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline long
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              long>::__type
     lround(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return lround(__type(__x));
-    }
+    { return __builtin_lround(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   nearbyint(float __x)
   { return __builtin_nearbyintf(__x); }
 
-  inline long double
+  constexpr long double
   nearbyint(long double __x)
   { return __builtin_nearbyintl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     nearbyint(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return nearbyint(__type(__x));
-    }
+    { return __builtin_nearbyint(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   nextafter(float __x, float __y)
   { return __builtin_nextafterf(__x, __y); }
 
-  inline long double
+  constexpr long double
   nextafter(long double __x, long double __y)
   { return __builtin_nextafterl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     nextafter(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return nextafter(__type(__x), __type(__y));
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   nexttoward(float __x, long double __y)
   { return __builtin_nexttowardf(__x, __y); }
 
-  inline long double
+  constexpr long double
   nexttoward(long double __x, long double __y)
   { return __builtin_nexttowardl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     nexttoward(_Tp __x, long double __y)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return nexttoward(__type(__x), __y);
-    }
+    { return __builtin_nexttoward(__x, __y); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   remainder(float __x, float __y)
   { return __builtin_remainderf(__x, __y); }
 
-  inline long double
+  constexpr long double
   remainder(long double __x, long double __y)
   { return __builtin_remainderl(__x, __y); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
-    inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
+    constexpr typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     remainder(_Tp __x, _Up __y)
     {
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return remainder(__type(__x), __type(__y));
     }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
   inline float
   remquo(float __x, float __y, int* __pquo)
   { return __builtin_remquof(__x, __y, __pquo); }
@@ -1551,7 +1721,9 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
   inline long double
   remquo(long double __x, long double __y, int* __pquo)
   { return __builtin_remquol(__x, __y, __pquo); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp, typename _Up>
     inline typename __gnu_cxx::__promote_2<_Tp, _Up>::__type
     remquo(_Tp __x, _Up __y, int* __pquo)
@@ -1559,108 +1731,202 @@ _GLIBCXX_BEGIN_NAMESPACE_VERSION
       typedef typename __gnu_cxx::__promote_2<_Tp, _Up>::__type __type;
       return remquo(__type(__x), __type(__y), __pquo);
     }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   rint(float __x)
   { return __builtin_rintf(__x); }
 
-  inline long double
+  constexpr long double
   rint(long double __x)
   { return __builtin_rintl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     rint(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return rint(__type(__x));
-    }
+    { return __builtin_rint(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   round(float __x)
   { return __builtin_roundf(__x); }
 
-  inline long double
+  constexpr long double
   round(long double __x)
   { return __builtin_roundl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     round(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return round(__type(__x));
-    }
+    { return __builtin_round(__x); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   scalbln(float __x, long __ex)
   { return __builtin_scalblnf(__x, __ex); }
 
-  inline long double
+  constexpr long double
   scalbln(long double __x, long __ex)
   { return __builtin_scalblnl(__x, __ex); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     scalbln(_Tp __x, long __ex)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return scalbln(__type(__x), __ex);
-    }
+    { return __builtin_scalbln(__x, __ex); }
+#endif
  
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   scalbn(float __x, int __ex)
   { return __builtin_scalbnf(__x, __ex); }
 
-  inline long double
+  constexpr long double
   scalbn(long double __x, int __ex)
   { return __builtin_scalbnl(__x, __ex); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     scalbn(_Tp __x, int __ex)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return scalbn(__type(__x), __ex);
-    }
+    { return __builtin_scalbn(__x, __ex); }
+#endif
 
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   tgamma(float __x)
   { return __builtin_tgammaf(__x); }
 
-  inline long double
+  constexpr long double
   tgamma(long double __x)
   { return __builtin_tgammal(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     tgamma(_Tp __x)
-    {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return tgamma(__type(__x));
-    }
+    { return __builtin_tgamma(__x); }
+#endif
  
-  inline float
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_FP
+  constexpr float
   trunc(float __x)
   { return __builtin_truncf(__x); }
 
-  inline long double
+  constexpr long double
   trunc(long double __x)
   { return __builtin_truncl(__x); }
+#endif
 
+#ifndef __CORRECT_ISO_CPP11_MATH_H_PROTO_INT
   template<typename _Tp>
-    inline typename __gnu_cxx::__promote<_Tp>::__type 
+    constexpr typename __gnu_cxx::__enable_if<__is_integer<_Tp>::__value, 
+                                              double>::__type
     trunc(_Tp __x)
+    { return __builtin_trunc(__x); }
+#endif
+
+#endif // _GLIBCXX_USE_C99_MATH_TR1
+#endif // C++11
+
+#if __cplusplus > 201402L
+
+  // [c.math.hypot3], three-dimensional hypotenuse
+#define __cpp_lib_hypot 201603
+
+  template<typename _Tp>
+    inline _Tp
+    __hypot3(_Tp __x, _Tp __y, _Tp __z)
     {
-      typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
-      return trunc(__type(__x));
+      __x = std::abs(__x);
+      __y = std::abs(__y);
+      __z = std::abs(__z);
+      if (_Tp __a = __x < __y ? __y < __z ? __z : __y : __x < __z ? __z : __x)
+       return __a * std::sqrt((__x / __a) * (__x / __a)
+                              + (__y / __a) * (__y / __a)
+                              + (__z / __a) * (__z / __a));
+      else
+       return {};
+    }
+
+  inline float
+  hypot(float __x, float __y, float __z)
+  { return std::__hypot3<float>(__x, __y, __z); }
+
+  inline double
+  hypot(double __x, double __y, double __z)
+  { return std::__hypot3<double>(__x, __y, __z); }
+
+  inline long double
+  hypot(long double __x, long double __y, long double __z)
+  { return std::__hypot3<long double>(__x, __y, __z); }
+
+  template<typename _Tp, typename _Up, typename _Vp>
+    typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type
+    hypot(_Tp __x, _Up __y, _Vp __z)
+    {
+      using __type = typename __gnu_cxx::__promote_3<_Tp, _Up, _Vp>::__type;
+      return std::__hypot3<__type>(__x, __y, __z);
+    }
+#endif // C++17
+
+#if __cplusplus > 201703L
+  // linear interpolation
+# define __cpp_lib_interpolate 201902L
+
+  template<typename _Fp>
+    constexpr _Fp
+    __lerp(_Fp __a, _Fp __b, _Fp __t) noexcept
+    {
+      if (__a <= 0 && __b >= 0 || __a >= 0 && __b <= 0)
+       return __t * __b + (1 - __t) * __a;
+
+      if (__t == 1)
+       return __b;                        // exact
+
+      // Exact at __t=0, monotonic except near __t=1,
+      // bounded, determinate, and consistent:
+      const _Fp __x = __a + __t * (__b - __a);
+      return __t > 1 == __b > __a
+       ? (__b < __x ? __x : __b)
+       : (__b > __x ? __x : __b);  // monotonic near __t=1
     }
 
+  constexpr float
+  lerp(float __a, float __b, float __t) noexcept
+  { return std::__lerp(__a, __b, __t); }
+
+  constexpr double
+  lerp(double __a, double __b, double __t) noexcept
+  { return std::__lerp(__a, __b, __t); }
+
+  constexpr long double
+  lerp(long double __a, long double __b, long double __t) noexcept
+  { return std::__lerp(__a, __b, __t); }
+#endif // C++20
+
 _GLIBCXX_END_NAMESPACE_VERSION
 } // namespace
 
-#endif // _GLIBCXX_USE_C99_MATH_TR1
+#if _GLIBCXX_USE_STD_SPEC_FUNCS
+#  include <bits/specfun.h>
+#endif
 
-#endif // __GXX_EXPERIMENTAL_CXX0X__
+} // extern "C++"
 
 #endif