#ifndef _LIBC_MATH_H
#define _LIBC_MATH_H

#pragma once

#ifdef __cplusplus
extern "C" {
#endif

#define HUGE_VAL  __builtin_huge_val()
#define INFINITY  __builtin_inff()
#define NAN       __builtin_nanf("")

/*
 * X/Open math constants. Not in ISO C, but ported code assumes them freely
 * (NetSurf's about:chart uses M_PI and M_PI_2).
 */
#define M_E        2.7182818284590452354
#define M_LOG2E    1.4426950408889634074
#define M_LOG10E   0.43429448190325182765
#define M_LN2      0.69314718055994530942
#define M_LN10     2.30258509299404568402
#define M_PI       3.14159265358979323846
#define M_PI_2     1.57079632679489661923
#define M_PI_4     0.78539816339744830962
#define M_1_PI     0.31830988618379067154
#define M_2_PI     0.63661977236758134308
#define M_2_SQRTPI 1.12837916709551257390
#define M_SQRT2    1.41421356237309504880
#define M_SQRT1_2  0.70710678118654752440

/* C99 floating-point classification. */
#define FP_NAN       0
#define FP_INFINITE  1
#define FP_ZERO      2
#define FP_SUBNORMAL 3
#define FP_NORMAL    4

#define fpclassify(x) \
    __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL, FP_SUBNORMAL, \
                         FP_ZERO, x)
#define isnan(x)    __builtin_isnan(x)
#define isinf(x)    __builtin_isinf(x)
#define isfinite(x) __builtin_isfinite(x)
#define isnormal(x) __builtin_isnormal(x)
#define signbit(x)  __builtin_signbit(x)

double fabs(double x);
double frexp(double x, int *exp);

/* C89 modf and the C99 float variants assumed by hosted libstdc++
   (--with-newlib crossconfig). The float forms wrap the double
   implementations. */
double modf(double x, double *iptr);
float  modff(float x, float *iptr);
double hypot(double x, double y);
float  hypotf(float x, float y);
float  acosf(float x);
float  asinf(float x);
float  atanf(float x);
float  atan2f(float y, float x);
float  coshf(float x);
float  expf(float x);
float  fmodf(float x, float y);
float  frexpf(float x, int *exp);
float  ldexpf(float x, int exp);
float  logf(float x);
float  log10f(float x);
float  powf(float x, float y);
float  sinhf(float x);
float  tanf(float x);
float  tanhf(float x);
double ldexp(double x, int n);
long double ldexpl(long double x, int n);
double floor(double x);
double ceil(double x);
double sqrt(double x);
double sin(double x);
double cos(double x);
double tan(double x);
double asin(double x);
double acos(double x);
double atan(double x);
double atan2(double y, double x);
double pow(double base, double exp);
double log(double x);
double log2(double x);
double log10(double x);
double exp(double x);
double fmod(double x, double y);
double sinh(double x);
double cosh(double x);
double tanh(double x);
double round(double x);

float  floorf(float x);
float  ceilf(float x);
float  fabsf(float x);
float  sqrtf(float x);
float  sinf(float x);
float  cosf(float x);

#ifdef __cplusplus
}
#endif

#endif /* _LIBC_MATH_H */