#ifndef __MATH__
#define __MATH__
#include "sys/cdefs.h"
#if (!defined(__WANT_LONG_DOUBLE_FORMAT__))
#if defined(__APPLE_CC__) && defined(__LONG_DOUBLE_128__)
#define __WANT_LONG_DOUBLE_FORMAT__ 128
#else
#define __WANT_LONG_DOUBLE_FORMAT__ 64
#endif
#endif
#if ( __WANT_LONG_DOUBLE_FORMAT__ - 0L == 128L )
#define __LIBMLDBL_COMPAT(sym) __asm("_" __STRING(sym) "$LDBL128")
#elif ( __WANT_LONG_DOUBLE_FORMAT__ - 0L == 64L )
#define __LIBMLDBL_COMPAT(sym)
#else
#define __LIBMLDBL_COMPAT(sym)
#endif
#ifdef __cplusplus
extern "C" {
#endif
#if __FLT_EVAL_METHOD__ == 0
typedef float float_t;
typedef double double_t;
#elif __FLT_EVAL_METHOD__ == 1
typedef double float_t;
typedef double double_t;
#elif __FLT_EVAL_METHOD__ == 2 || __FLT_EVAL_METHOD__ == -1
typedef long double float_t;
typedef long double double_t;
#else
#error "Unsupported value of __FLT_EVAL_METHOD__."
#endif
#if defined(__GNUC__)
#define HUGE_VAL __builtin_huge_val()
#define HUGE_VALF __builtin_huge_valf()
#define HUGE_VALL __builtin_huge_vall()
#define NAN __builtin_nanf("0x7fc00000")
#define __MATH_H_ALWAYS_INLINE__ __attribute__ ((always_inline))
#else
#define HUGE_VAL 1e500
#define HUGE_VALF 1e50f
#define HUGE_VALL 1e500L
#define NAN __nan( )
#define __MATH_H_ALWAYS_INLINE__
#endif
#define INFINITY HUGE_VALF
enum {
_FP_NAN = 1,
_FP_INFINITE = 2,
_FP_ZERO = 3,
_FP_NORMAL = 4,
_FP_SUBNORMAL = 5,
_FP_SUPERNORMAL = 6
};
#define FP_NAN _FP_NAN
#define FP_INFINITE _FP_INFINITE
#define FP_ZERO _FP_ZERO
#define FP_NORMAL _FP_NORMAL
#define FP_SUBNORMAL _FP_SUBNORMAL
#define FP_SUPERNORMAL _FP_SUPERNORMAL
#undef FP_FAST_FMA
#undef FP_FAST_FMAF
#undef FP_FAST_FMAL
#define FP_ILOGB0 (-2147483647)
#define FP_ILOGBNAN (2147483647)
#define MATH_ERRNO 1
#define MATH_ERREXCEPT 2
#define math_errhandling (__math_errhandling())
extern unsigned int __math_errhandling ( void );
#if (__WANT_LONG_DOUBLE_FORMAT__ - 0L == 128L)
#define fpclassify(x) \
( sizeof (x) == sizeof(float ) ? __fpclassifyf((float)x) \
: sizeof (x) == sizeof(double) ? __fpclassifyd((double)x) \
: __fpclassify ((long double)x))
extern int __fpclassifyf(float );
extern int __fpclassifyd(double );
extern int __fpclassify (long double);
#if defined( __GNUC__ )
#define isnormal(x) \
( sizeof (x) == sizeof(float ) ? __inline_isnormalf((float)(x)) \
: sizeof (x) == sizeof(double) ? __inline_isnormald((double)(x)) \
: __inline_isnormal ((long double)(x)))
#define isfinite(x) \
( sizeof (x) == sizeof(float ) ? __inline_isfinitef((float)(x)) \
: sizeof (x) == sizeof(double) ? __inline_isfinited((double)(x)) \
: __inline_isfinite ((long double)(x)))
#define isinf(x) \
( sizeof (x) == sizeof(float ) ? __inline_isinff((float)(x)) \
: sizeof (x) == sizeof(double) ? __inline_isinfd((double)(x)) \
: __inline_isinf ((long double)(x)))
#define isnan(x) \
( sizeof (x) == sizeof(float ) ? __inline_isnanf((float)(x)) \
: sizeof (x) == sizeof(double) ? __inline_isnand((double)(x)) \
: __inline_isnan ((long double)(x)))
#define signbit(x) \
( sizeof (x) == sizeof(float ) ? __inline_signbitf((float)(x)) \
: sizeof (x) == sizeof(double) ? __inline_signbitd((double)(x)) \
: __inline_signbit((long double)(x)))
static __inline__ int __inline_isfinitef (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isfinited (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isfinite (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinff (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinfd (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinf (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnanf (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnand (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnan (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnormalf (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnormald (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnormal (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_signbitf (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_signbitd (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_signbit (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinff( float __x ) { return __builtin_fabsf(__x) == __builtin_inff(); }
static __inline__ int __inline_isinfd( double __x ) { return __builtin_fabs(__x) == __builtin_inf(); }
static __inline__ int __inline_isinf( long double __x ) { return __builtin_fabsl(__x) == __builtin_infl(); }
static __inline__ int __inline_isfinitef( float __x ) { return __x == __x && __builtin_fabsf(__x) != __builtin_inff(); }
static __inline__ int __inline_isfinited( double __x ) { return __x == __x && __builtin_fabs(__x) != __builtin_inf(); }
static __inline__ int __inline_isfinite( long double __x ) { return __x == __x && __builtin_fabsl(__x) != __builtin_infl(); }
static __inline__ int __inline_isnanf( float __x ) { return __x != __x; }
static __inline__ int __inline_isnand( double __x ) { return __x != __x; }
static __inline__ int __inline_isnan( long double __x ) { return __x != __x; }
static __inline__ int __inline_signbitf( float __x ) { union{ float __f; unsigned int __u; }__u = {__x}; return (int)(__u.__u >> 31); }
static __inline__ int __inline_signbitd( double __x ) { union{ double __f; unsigned long long __u; }__u = {__x}; return (int)(__u.__u >> 63); }
static __inline__ int __inline_signbit( long double __x ){ union{ long double __ld; long long __p[2]; }__u = {__x}; return (int) ((unsigned long long) __u.__p[0] >> 63); }
static __inline__ int __inline_isnormalf( float __x ) { float fabsf = __builtin_fabsf(__x); if( __x != __x ) return 0; return fabsf < __builtin_inff() && fabsf >= __FLT_MIN__; }
static __inline__ int __inline_isnormald( double __x ) { double fabsf = __builtin_fabs(__x); if( __x != __x ) return 0; return fabsf < __builtin_inf() && fabsf >= __DBL_MIN__; }
static __inline__ int __inline_isnormal( long double __x ) { long double fabsf = __builtin_fabsl(__x); if( __x != __x ) return 0; return fabsf < __builtin_infl() && fabsf >= __LDBL_MIN__; }
#else
#define isnormal(x) \
( sizeof (x) == sizeof(float ) ? __isnormalf((float)(x)) \
: sizeof (x) == sizeof(double) ? __isnormald((double)(x)) \
: __isnormal ((long double)(x)))
#define isfinite(x) \
( sizeof (x) == sizeof(float ) ? __isfinitef((float)(x)) \
: sizeof (x) == sizeof(double) ? __isfinited((double)(x)) \
: __isfinite ((long double)(x)))
#define isinf(x) \
( sizeof (x) == sizeof(float ) ? __isinff((float)(x)) \
: sizeof (x) == sizeof(double) ? __isinfd((double)(x)) \
: __isinf ((long double)(x)))
#define isnan(x) \
( sizeof (x) == sizeof(float ) ? __isnanf((float)(x)) \
: sizeof (x) == sizeof(double) ? __isnand((double)(x)) \
: __isnan ((long double)(x)))
#define signbit(x) \
( sizeof (x) == sizeof(float ) ? __signbitf((float)(x)) \
: sizeof (x) == sizeof(double) ? __signbitd((double)(x)) \
: __signbitl((long double)(x)))
extern int __isnormalf (float );
extern int __isnormald (double );
extern int __isnormal (long double);
extern int __isfinitef (float );
extern int __isfinited (double );
extern int __isfinite (long double);
extern int __isinff (float );
extern int __isinfd (double );
extern int __isinf (long double);
extern int __isnanf (float );
extern int __isnand (double );
extern int __isnan (long double);
extern int __signbitf (float );
extern int __signbitd (double );
extern int __signbitl (long double);
#endif
#else
#define fpclassify(x) \
( sizeof (x) == sizeof(float ) ? __fpclassifyf((float)(x)) : __fpclassifyd((double)(x)) )
extern int __fpclassifyf(float );
extern int __fpclassifyd(double );
#if defined( __GNUC__ )
#define isnormal(x) \
( sizeof (x) == sizeof(float ) ? __inline_isnormalf((float)(x)) : __inline_isnormald((double)(x)) )
#define isfinite(x) \
( sizeof (x) == sizeof(float ) ? __inline_isfinitef((float)(x)) : __inline_isfinited((double)(x)) )
#define isinf(x) \
( sizeof (x) == sizeof(float ) ? __inline_isinff((float)(x)) : __inline_isinfd((double)(x)) )
#define isnan(x) \
( sizeof (x) == sizeof(float ) ? __inline_isnanf((float)(x)) : __inline_isnand((double)(x)) )
#define signbit(x) \
( sizeof (x) == sizeof(float ) ? __inline_signbitf((float)(x)) : __inline_signbitd((double)(x)) )
static __inline__ int __inline_isfinitef (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isfinited (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isfinite (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinff (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinfd (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinf (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnanf (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnand (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnan (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnormalf (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnormald (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isnormal (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_signbitf (float ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_signbitd (double ) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_signbit (long double) __MATH_H_ALWAYS_INLINE__;
static __inline__ int __inline_isinff( float __x ) { return __builtin_fabsf(__x) == __builtin_inff(); }
static __inline__ int __inline_isinfd( double __x ) { return __builtin_fabs(__x) == __builtin_inf(); }
static __inline__ int __inline_isinf( long double __x ) { return __builtin_fabsl(__x) == __builtin_infl(); }
static __inline__ int __inline_isfinitef( float __x ) { return __x == __x && __builtin_fabsf(__x) != __builtin_inff(); }
static __inline__ int __inline_isfinited( double __x ) { return __x == __x && __builtin_fabs(__x) != __builtin_inf(); }
static __inline__ int __inline_isfinite( long double __x ) { return __x == __x && __builtin_fabsl(__x) != __builtin_infl(); }
static __inline__ int __inline_isnanf( float __x ) { return __x != __x; }
static __inline__ int __inline_isnand( double __x ) { return __x != __x; }
static __inline__ int __inline_isnan( long double __x ) { return __x != __x; }
static __inline__ int __inline_signbitf( float __x ) { union{ float __f; unsigned int __u; }__u = {__x}; return (int)(__u.__u >> 31); }
static __inline__ int __inline_signbitd( double __x ) { union{ double __f; unsigned long long __u; }__u = {__x}; return (int)(__u.__u >> 63); }
static __inline__ int __inline_signbit( long double __x ){ union{ long double __ld; unsigned long long __p; }__u = {__x}; return (int) ((unsigned long long) __u.__p >> 63); }
static __inline__ int __inline_isnormalf( float __x ) { float fabsf = __builtin_fabsf(__x); if( __x != __x ) return 0; return fabsf < __builtin_inff() && fabsf >= __FLT_MIN__; }
static __inline__ int __inline_isnormald( double __x ) { double fabsf = __builtin_fabs(__x); if( __x != __x ) return 0; return fabsf < __builtin_inf() && fabsf >= __DBL_MIN__; }
static __inline__ int __inline_isnormal( long double __x ) { long double fabsf = __builtin_fabsl(__x); if( __x != __x ) return 0; return fabsf < __builtin_infl() && fabsf >= __LDBL_MIN__; }
#else
#define isnormal(x) \
( sizeof (x) == sizeof(float ) ? __isnormalf((float)(x)) : __isnormald((double)(x)) )
#define isfinite(x) \
( sizeof (x) == sizeof(float ) ? __isfinitef((float)(x)) : __isfinited((double)(x)) )
#define isinf(x) \
( sizeof (x) == sizeof(float ) ? __isinff((float)(x)) : __isinfd((double)(x)) )
#define isnan(x) \
( sizeof (x) == sizeof(float ) ? __isnanf((float)(x)) : __isnand((double)(x)) )
#define signbit(x) \
( sizeof (x) == sizeof(float ) ? __signbitf((float)(x)) : __signbitd((double)(x)) )
extern int __isnormalf (float );
extern int __isnormald (double );
extern int __isfinitef (float );
extern int __isfinited (double );
extern int __isinff (float );
extern int __isinfd (double );
extern int __isnanf (float );
extern int __isnand (double );
extern int __signbitf (float );
extern int __signbitd (double );
#endif
#endif
extern double acos( double );
extern float acosf( float );
extern double asin( double );
extern float asinf( float );
extern double atan( double );
extern float atanf( float );
extern double atan2( double, double );
extern float atan2f( float, float );
extern double cos( double );
extern float cosf( float );
extern double sin( double );
extern float sinf( float );
extern double tan( double );
extern float tanf( float );
extern double acosh( double );
extern float acoshf( float );
extern double asinh( double );
extern float asinhf( float );
extern double atanh( double );
extern float atanhf( float );
extern double cosh( double );
extern float coshf( float );
extern double sinh( double );
extern float sinhf( float );
extern double tanh( double );
extern float tanhf( float );
extern double exp( double );
extern float expf( float );
extern double exp2( double );
extern float exp2f( float );
extern double expm1( double );
extern float expm1f( float );
extern double log( double );
extern float logf( float );
extern double log10( double );
extern float log10f( float );
extern double log2( double );
extern float log2f( float );
extern double log1p( double );
extern float log1pf( float );
extern double logb( double );
extern float logbf( float );
extern double modf( double, double * );
extern float modff( float, float * );
extern double ldexp( double, int );
extern float ldexpf( float, int );
extern double frexp( double, int * );
extern float frexpf( float, int * );
extern int ilogb( double );
extern int ilogbf( float );
extern double scalbn( double, int );
extern float scalbnf( float, int );
extern double scalbln( double, long int );
extern float scalblnf( float, long int );
extern double fabs( double );
extern float fabsf( float );
extern double cbrt( double );
extern float cbrtf( float );
extern double hypot( double, double );
extern float hypotf( float, float );
extern double pow( double, double );
extern float powf( float, float );
extern double sqrt( double );
extern float sqrtf( float );
extern double erf( double );
extern float erff( float );
extern double erfc( double );
extern float erfcf( float );
extern double lgamma( double );
extern float lgammaf( float );
extern double tgamma( double );
extern float tgammaf( float );
extern double ceil( double );
extern float ceilf( float );
extern double floor( double );
extern float floorf( float );
extern double nearbyint( double );
extern float nearbyintf( float );
extern double rint( double );
extern float rintf( float );
extern long int lrint( double );
extern long int lrintf( float );
extern long long int llrint( double );
extern long long int llrintf( float );
extern double round( double );
extern float roundf( float );
extern long int lround( double );
extern long int lroundf( float );
extern long long int llround( double );
extern long long int llroundf( float );
extern double trunc( double );
extern float truncf( float );
extern double fmod( double, double );
extern float fmodf( float, float );
extern double remainder( double, double );
extern float remainderf( float, float );
extern double remquo( double, double, int * );
extern float remquof( float, float, int * );
extern double copysign( double, double );
extern float copysignf( float, float );
extern double nan( const char * );
extern float nanf( const char * );
extern double nextafter( double, double );
extern float nextafterf( float, float );
extern double fdim( double, double );
extern float fdimf( float, float );
extern double fmax( double, double );
extern float fmaxf( float, float );
extern double fmin( double, double );
extern float fminf( float, float );
extern double fma( double, double, double );
extern float fmaf( float, float, float );
#if ( __WANT_LONG_DOUBLE_FORMAT__ - 0L > 0L )
extern long double acosl( long double ) __LIBMLDBL_COMPAT(acosl);
extern long double asinl( long double ) __LIBMLDBL_COMPAT(asinl);
extern long double atanl( long double ) __LIBMLDBL_COMPAT(atanl);
extern long double atan2l( long double, long double ) __LIBMLDBL_COMPAT(atan2l);
extern long double cosl( long double ) __LIBMLDBL_COMPAT(cosl);
extern long double sinl( long double ) __LIBMLDBL_COMPAT(sinl);
extern long double tanl( long double ) __LIBMLDBL_COMPAT(tanl);
extern long double acoshl( long double ) __LIBMLDBL_COMPAT(acoshl);
extern long double asinhl( long double ) __LIBMLDBL_COMPAT(asinhl);
extern long double atanhl( long double ) __LIBMLDBL_COMPAT(atanhl);
extern long double coshl( long double ) __LIBMLDBL_COMPAT(coshl);
extern long double sinhl( long double ) __LIBMLDBL_COMPAT(sinhl);
extern long double tanhl( long double ) __LIBMLDBL_COMPAT(tanhl);
extern long double expl( long double ) __LIBMLDBL_COMPAT(expl);
extern long double exp2l( long double ) __LIBMLDBL_COMPAT(exp2l);
extern long double expm1l( long double ) __LIBMLDBL_COMPAT(expm1l);
extern long double logl( long double ) __LIBMLDBL_COMPAT(logl);
extern long double log10l( long double ) __LIBMLDBL_COMPAT(log10l);
extern long double log2l( long double ) __LIBMLDBL_COMPAT(log2l);
extern long double log1pl( long double ) __LIBMLDBL_COMPAT(log1pl);
extern long double logbl( long double ) __LIBMLDBL_COMPAT(logbl);
extern long double modfl( long double, long double * ) __LIBMLDBL_COMPAT(modfl);
extern long double ldexpl( long double, int ) __LIBMLDBL_COMPAT(ldexpl);
extern long double frexpl( long double, int * ) __LIBMLDBL_COMPAT(frexpl);
extern int ilogbl( long double ) __LIBMLDBL_COMPAT(ilogbl);
extern long double scalbnl( long double, int ) __LIBMLDBL_COMPAT(scalbnl);
extern long double scalblnl( long double, long int ) __LIBMLDBL_COMPAT(scalblnl);
extern long double fabsl( long double ) __LIBMLDBL_COMPAT(fabsl);
extern long double cbrtl( long double ) __LIBMLDBL_COMPAT(cbrtl);
extern long double hypotl( long double, long double ) __LIBMLDBL_COMPAT(hypotl);
extern long double powl( long double, long double ) __LIBMLDBL_COMPAT(powl);
extern long double sqrtl( long double ) __LIBMLDBL_COMPAT(sqrtl);
extern long double erfl( long double ) __LIBMLDBL_COMPAT(erfl);
extern long double erfcl( long double ) __LIBMLDBL_COMPAT(erfcl);
extern long double lgammal( long double ) __LIBMLDBL_COMPAT(lgammal);
extern long double tgammal( long double ) __LIBMLDBL_COMPAT(tgammal);
extern long double ceill( long double ) __LIBMLDBL_COMPAT(ceill);
extern long double floorl( long double ) __LIBMLDBL_COMPAT(floorl);
extern long double nearbyintl( long double ) __LIBMLDBL_COMPAT(nearbyintl);
extern long double rintl( long double ) __LIBMLDBL_COMPAT(rintl);
extern long int lrintl( long double ) __LIBMLDBL_COMPAT(lrintl);
extern long long int llrintl( long double ) __LIBMLDBL_COMPAT(llrintl);
extern long double roundl( long double ) __LIBMLDBL_COMPAT(roundl);
extern long int lroundl( long double ) __LIBMLDBL_COMPAT(lroundl);
extern long long int llroundl( long double ) __LIBMLDBL_COMPAT(llroundl);
extern long double truncl( long double ) __LIBMLDBL_COMPAT(truncl);
extern long double fmodl( long double, long double) __LIBMLDBL_COMPAT(fmodl);
extern long double remainderl( long double, long double ) __LIBMLDBL_COMPAT(remainderl);
extern long double remquol( long double, long double, int * ) __LIBMLDBL_COMPAT(remquol);
extern long double copysignl( long double, long double ) __LIBMLDBL_COMPAT(copysignl);
extern long double nanl( const char * ) __LIBMLDBL_COMPAT(nanl);
extern long double nextafterl( long double, long double ) __LIBMLDBL_COMPAT(nextafterl);
extern double nexttoward( double, long double ) __LIBMLDBL_COMPAT(nexttoward);
extern float nexttowardf( float, long double ) __LIBMLDBL_COMPAT(nexttowardf);
extern long double nexttowardl( long double, long double ) __LIBMLDBL_COMPAT(nexttowardl);
extern long double fdiml( long double, long double ) __LIBMLDBL_COMPAT(fdiml);
extern long double fmaxl( long double, long double ) __LIBMLDBL_COMPAT(fmaxl);
extern long double fminl( long double, long double ) __LIBMLDBL_COMPAT(fminl);
extern long double fmal( long double, long double, long double ) __LIBMLDBL_COMPAT(fmal);
#endif
#define isgreater(x, y) __builtin_isgreater ((x),(y))
#define isgreaterequal(x, y) __builtin_isgreaterequal ((x),(y))
#define isless(x, y) __builtin_isless ((x),(y))
#define islessequal(x, y) __builtin_islessequal ((x),(y))
#define islessgreater(x, y) __builtin_islessgreater ((x),(y))
#define isunordered(x, y) __builtin_isunordered ((x),(y))
extern double __inf( void );
extern float __inff( void );
extern float __nan( void );
#if !defined(_ANSI_SOURCE)
extern double j0 ( double );
extern double j1 ( double );
extern double jn ( int, double );
extern double y0 ( double );
extern double y1 ( double );
extern double yn ( int, double );
#if defined( __ppc__ )
#if defined( __GNUC__ )
extern double scalb ( double, double ) __asm("_scalb$UNIX2003" );
#else
extern double scalb ( double , int );
#endif
#else
extern double scalb ( double, double );
#endif
#define M_E 2.71828182845904523536028747135266250
#define M_LOG2E 1.44269504088896340735992468100189214
#define M_LOG10E 0.434294481903251827651128918916605082
#define M_LN2 0.693147180559945309417232121458176568
#define M_LN10 2.30258509299404568401799145468436421
#define M_PI 3.14159265358979323846264338327950288
#define M_PI_2 1.57079632679489661923132169163975144
#define M_PI_4 0.785398163397448309615660845819875721
#define M_1_PI 0.318309886183790671537767526745028724
#define M_2_PI 0.636619772367581343075535053490057448
#define M_2_SQRTPI 1.12837916709551257389615890312154517
#define M_SQRT2 1.41421356237309504880168872420969808
#define M_SQRT1_2 0.707106781186547524400844362104849039
#define MAXFLOAT ((float)3.40282346638528860e+38)
extern int signgam;
#endif
#if !defined(__NOEXTENSIONS__) && (!defined(_POSIX_C_SOURCE) || defined(_DARWIN_C_SOURCE))
#define __WANT_EXTENSIONS__
#endif
#ifdef __WANT_EXTENSIONS__
#define FP_SNAN FP_NAN
#define FP_QNAN FP_NAN
extern long int rinttol ( double );
extern long int roundtol ( double );
typedef struct __complex_s {
double Real;
double Imag;
} __complex_t;
#if !defined(_ANSI_SOURCE) && (!defined(_POSIX_C_SOURCE) || defined(_DARWIN_C_SOURCE))
#if (!defined(_XOPEN_SOURCE) || defined(_DARWIN_C_SOURCE))
enum fdversion {_fdlibm_ieee = -1, _fdlibm_svid, _fdlibm_xopen, _fdlibm_posix};
#define fdlibm_ieee _fdlibm_ieee
#define fdlibm_svid _fdlibm_svid
#define fdlibm_xopen _fdlibm_xopen
#define fdlibm_posix _fdlibm_posix
#define _LIB_VERSION_TYPE enum fdversion
#define _LIB_VERSION _fdlib_version
extern _LIB_VERSION_TYPE _LIB_VERSION;
#define _IEEE_ fdlibm_ieee
#define _SVID_ fdlibm_svid
#define _XOPEN_ fdlibm_xopen
#define _POSIX_ fdlibm_posix
#if !defined(__cplusplus)
struct exception {
int type;
char *name;
double arg1;
double arg2;
double retval;
};
#endif
#define HUGE MAXFLOAT
#define X_TLOSS 1.41484755040568800000e+16
#define DOMAIN 1
#define SING 2
#define OVERFLOW 3
#define UNDERFLOW 4
#define TLOSS 5
#define PLOSS 6
#endif
#endif
#if !defined(_ANSI_SOURCE) && (!defined(_POSIX_C_SOURCE) || defined(_DARWIN_C_SOURCE))
extern int finite ( double );
extern double gamma ( double );
#if (!defined(_XOPEN_SOURCE) || defined(_DARWIN_C_SOURCE))
#if !defined(__cplusplus)
extern int matherr ( struct exception * );
#endif
extern double significand ( double );
extern double drem ( double, double );
#ifdef _REENTRANT
extern double gamma_r ( double, int * );
extern double lgamma_r ( double, int * );
#endif
#endif
#endif
#endif
#if defined(__FAST_MATH__) && !defined(__cplusplus)
#define __FSELS(e,t,f) (((e) >= 0.0f) ? (t) : (f))
#define __FSEL(e,t,f) (((e) >= 0.0) ? (t) : (f))
static __inline__ float __fastmath_floorf( float f ) __attribute__((__always_inline__));
static __inline__ float __fastmath_floorf( float f )
{
float b, c, d, e, g, h, t;
c = __FSELS( f, -8388608.f, 8388608.f ); b = fabsf( f );
d = f - c; e = b - 8388608.f;
__asm__("" : "+f" (d));
d = d + c;
g = f - d;
h = __FSELS( g, 0.0f, 1.0f );
t = d - h;
return __FSELS( e, f, t );
}
static __inline__ float __fastmath_ceilf( float f ) __attribute__((__always_inline__));
static __inline__ float __fastmath_ceilf( float f )
{
float b, c, d, e, g, h, t;
c = __FSELS( f, -8388608.f, 8388608.f ); b = fabsf( f );
d = f - c; e = b - 8388608.f;
__asm__("" : "+f" (d));
d = d + c;
g = d - f;
h = __FSELS( g, 0.0f, 1.0f );
t = d + h;
return __FSELS( e, f, t );
}
static __inline__ double __fastmath_floor( double f ) __attribute__((__always_inline__));
static __inline__ double __fastmath_floor( double f )
{
double b, c, d, e, g, h, t;
c = __FSEL( f, -4503599627370496., 4503599627370496. ); b = fabs( f );
d = f - c; e = b - 4503599627370496.;
__asm__("" : "+f" (d));
d = d + c;
g = f - d;
h = __FSEL( g, 0.0, 1.0 );
t = d - h;
return __FSEL( e, f, t );
}
static __inline__ double __fastmath_ceil( double f ) __attribute__((__always_inline__));
static __inline__ double __fastmath_ceil( double f )
{
double b, c, d, e, g, h, t;
c = __FSEL( f, -4503599627370496., 4503599627370496. ); b = fabs( f );
d = f - c; e = b - 4503599627370496.;
__asm__("" : "+f" (d));
d = d + c;
g = d - f;
h = __FSEL( g, 0.0, 1.0 );
t = d + h;
return __FSEL( e, f, t );
}
#define floorf(x) __fastmath_floorf((x))
#define ceilf(x) __fastmath_ceilf((x))
#define floor(x) __fastmath_floor((x))
#define ceil(x) __fastmath_ceil((x))
#endif
#ifdef __cplusplus
}
#endif
#endif