#ifndef _MATH_H
#define _MATH_H

#include <sys/cdefs.h>

__BEGIN_DECLS

#define M_E		2.7182818284590452354	/* e */
#define M_LOG2E		1.4426950408889634074	/* log_2 e */
#define M_LOG10E	0.43429448190325182765	/* log_10 e */
#define M_LN2		0.69314718055994530942	/* log_e 2 */
#define M_LN10		2.30258509299404568402	/* log_e 10 */
#define M_PI		3.14159265358979323846	/* pi */
#define M_PI_2		1.57079632679489661923	/* pi/2 */
#define M_PI_4		0.78539816339744830962	/* pi/4 */
#define M_1_PI		0.31830988618379067154	/* 1/pi */
#define M_2_PI		0.63661977236758134308	/* 2/pi */
#define M_2_SQRTPI	1.12837916709551257390	/* 2/sqrt(pi) */
#define M_SQRT2		1.41421356237309504880	/* sqrt(2) */
#define M_SQRT1_2	0.70710678118654752440	/* 1/sqrt(2) */

#define M_El		2.7182818284590452353602874713526625L  /* e */
#define M_LOG2El	1.4426950408889634073599246810018922L  /* log_2 e */
#define M_LOG10El	0.4342944819032518276511289189166051L  /* log_10 e */
#define M_LN2l		0.6931471805599453094172321214581766L  /* log_e 2 */
#define M_LN10l		2.3025850929940456840179914546843642L  /* log_e 10 */
#define M_PIl		3.1415926535897932384626433832795029L  /* pi */
#define M_PI_2l		1.5707963267948966192313216916397514L  /* pi/2 */
#define M_PI_4l		0.7853981633974483096156608458198757L  /* pi/4 */
#define M_1_PIl		0.3183098861837906715377675267450287L  /* 1/pi */
#define M_2_PIl		0.6366197723675813430755350534900574L  /* 2/pi */
#define M_2_SQRTPIl	1.1283791670955125738961589031215452L  /* 2/sqrt(pi) */
#define M_SQRT2l	1.4142135623730950488016887242096981L  /* sqrt(2) */
#define M_SQRT1_2l	0.7071067811865475244008443621048490L  /* 1/sqrt(2) */

double sin(double d) __THROW __attribute__((__const__));
double cos(double d) __THROW __attribute__((__const__));
void sincos(double d, double *sin, double *cos) __THROW __attribute__((__const__));
double tan(double d) __THROW __attribute__((__const__));

double sinh(double d) __THROW __attribute__((__const__));
double cosh(double d) __THROW __attribute__((__const__));
double tanh(double d) __THROW __attribute__((__const__));

double asin(double d) __THROW __attribute__((__const__));
double acos(double d) __THROW __attribute__((__const__));
double atan(double d) __THROW __attribute__((__const__));
double atan2(double y, double x) __THROW __attribute__((__const__));

double asinh(double d) __THROW __attribute__((__const__));
double acosh(double d) __THROW __attribute__((__const__));
double atanh(double d) __THROW __attribute__((__const__));

float sinf(float d) __THROW __attribute__((__const__));
float cosf(float d) __THROW __attribute__((__const__));
void sincosf(float d, float *sin, float *cos) __THROW __attribute__((__const__));
float tanf(float d) __THROW __attribute__((__const__));

float sinhf(float d) __THROW __attribute__((__const__));
float coshf(float d) __THROW __attribute__((__const__));
float tanhf(float d) __THROW __attribute__((__const__));

float asinf(float d) __THROW __attribute__((__const__));
float acosf(float d) __THROW __attribute__((__const__));
float atanf(float d) __THROW __attribute__((__const__));
float atan2f(float y, float x) __THROW __attribute__((__const__));

float asinhf(float d) __THROW __attribute__((__const__));
float acoshf(float d) __THROW __attribute__((__const__));
float atanhf(float d) __THROW __attribute__((__const__));

double exp(double d) __THROW __attribute__((__const__));
double expm1(double d) __THROW __attribute__((__const__));
double exp2(double d) __THROW __attribute__((__const__));
double exp10(double d) __THROW __attribute__((__const__));
double pow10(double d) __THROW __attribute__((__const__));
double log(double d) __THROW __attribute__((__const__));
double log1p(double d) __THROW __attribute__((__const__));
double log2(double d) __THROW __attribute__((__const__));
double log10(double d) __THROW __attribute__((__const__));

float expf(float d) __THROW __attribute__((__const__));
float expm1f(float d) __THROW __attribute__((__const__));
float exp2f(float d) __THROW __attribute__((__const__));
float exp10f(float d) __THROW __attribute__((__const__));
float pow10f(float d) __THROW __attribute__((__const__));
float logf(float d) __THROW __attribute__((__const__));
float log1pf(float d) __THROW __attribute__((__const__));
float log2f(float d) __THROW __attribute__((__const__));
float log10f(float d) __THROW __attribute__((__const__));

double sqrt(double x) __THROW __attribute__((__const__));
double fabs(double x) __THROW __attribute__((__const__));
double fmod(double x, double y) __THROW __attribute__((__const__));
double pow(double x, double y) __THROW __attribute__((__const__));

float sqrtf(float x) __THROW __attribute__((__const__));
float fabsf(float x) __THROW __attribute__((__const__));
float fmodf(float x, float y) __THROW __attribute__((__const__));
float powf(float x, float y) __THROW __attribute__((__const__));

double floor(double x) __attribute__((__const__));
double ceil(double x) __attribute__((__const__));
double round(double x) __attribute__((__const__));
double trunc(double x) __attribute__((__const__));

float floorf(float x) __attribute__((__const__));
float ceilf(float x) __attribute__((__const__));
float roundf(float x) __attribute__((__const__));
float truncf(float x) __attribute__((__const__));

double remainder(double x, double y) __THROW __attribute__((__const__));
double drem(double x, double y) __THROW __attribute__((__const__));
double expm1(double x) __THROW __attribute__((__const__));
double hypot(double x, double y) __THROW __attribute__((__const__));
double atan2(double x, double y) __THROW __attribute__((__const__));
double copysign(double value, double sign) __attribute__((__const__));

float remainderf(float x, float y) __THROW __attribute__((__const__));
float dremf(float x, float y) __THROW __attribute__((__const__));
float expm1f(float x) __THROW __attribute__((__const__));
float hypotf(float x, float y) __THROW __attribute__((__const__));
float atan2f(float x, float y) __THROW __attribute__((__const__));
float copysignf(float value, float sign) __attribute__((__const__));

#  define HUGE_VAL \
  (__extension__ \
   ((union { unsigned long long __ll; double __d; }) \
    { __ll: 0x7ff0000000000000ULL }).__d)
#  define HUGE_VALF \
  (__extension__ \
   ((union { unsigned int __i; float __f; }) \
    { __i: 0x7f800000UL }).__f)
#define NAN (__builtin_nanf(""))

#ifdef _GNU_SOURCE
void sincos(double x, double* sinx, double* cosx);
double ipow (double mant, int expo);
#endif

int isnan(double d) __attribute__((__const__));
int isinf(double d) __attribute__((__const__));
int finite(double d) __attribute__((__const__));

double j0(double x);
double j1(double x);
double jn(int n, double x);
double y0(double x);
double y1(double x);
double yn(int n, double x);
double erf(double x);
double erfc(double x);
double lgamma(double x);

double rint(double x);
double fdim(double x, double y);
double fma(double x, double y, double z);
double modf(double x, double * iptr);

float  rintf(float x);
float  fdimf(float x, float y);
float  fmaf(float x, float y, float z);
float  modff(float x, float * iptr);

double frexp(double x, int *exp);
double ldexp(double x, int exp);
double logb(double x);
float  frexpf(float x, int *exp);
int    ilogb(double x);
float  ldexpf(float x, int exp);
float  logbf(float x);
int    ilogbf(float x);


__END_DECLS

#endif