#ifndef _MATH_H
#define _MATH_H

#ifdef __cplusplus
extern "C" {
#endif

// NOTE(orca): not doing anything fancy for float_t and double_t
typedef float float_t;
typedef double double_t;

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

#define FP_NAN 0
#define FP_INFINITE 1
#define FP_ZERO 2
#define FP_SUBNORMAL 3
#define FP_NORMAL 4

int __fpclassify(double);
int __fpclassifyf(float);
int __fpclassifyl(long double);

static __inline unsigned __FLOAT_BITS(float __f)
{
    union
    {
        float __f;
        unsigned __i;
    } __u;

    __u.__f = __f;
    return __u.__i;
}

static __inline unsigned long long __DOUBLE_BITS(double __f)
{
    union
    {
        double __f;
        unsigned long long __i;
    } __u;

    __u.__f = __f;
    return __u.__i;
}

#define fpclassify(x) (                                                                           \
    sizeof(x) == sizeof(float) ? __fpclassifyf(x) : sizeof(x) == sizeof(double) ? __fpclassify(x) \
                                                                                : __fpclassifyl(x))

#define isinf(x) (                                                                                                                                              \
    sizeof(x) == sizeof(float) ? (__FLOAT_BITS(x) & 0x7fffffff) == 0x7f800000 : sizeof(x) == sizeof(double) ? (__DOUBLE_BITS(x) & -1ULL >> 1) == 0x7ffULL << 52 \
                                                                                                            : __fpclassifyl(x) == FP_INFINITE)

#define isnan(x) (                                                                                                                                            \
    sizeof(x) == sizeof(float) ? (__FLOAT_BITS(x) & 0x7fffffff) > 0x7f800000 : sizeof(x) == sizeof(double) ? (__DOUBLE_BITS(x) & -1ULL >> 1) > 0x7ffULL << 52 \
                                                                                                           : __fpclassifyl(x) == FP_NAN)

double fabs(double);
float fabsf(float);

double acos(double);
float acosf(float);

double cbrt(double);
float cbrtf(float);

double ceil(double);

double cos(double);
float cosf(float);

double floor(double);

double fmod(double, double);

double log(double);
float logf(float);
double log2(double);
float log2f(float);

double pow(double, double);
float powf(float, float);

double exp(double);
float expf(float);

double scalbn(double, int);

double sin(double);
float sinf(float);

double asin(double);
float asinf(float);

double tan(double);
float tanf(float);

double atan(double);
float atanf(float);
double atan2(double, double);
float atan2f(float, float);

double sqrt(double);
float sqrtf(float);

#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) */

//NOTE(orca) - implementation details
typedef unsigned uint32_t;
double __math_divzero(uint32_t sign);
float __math_divzerof(uint32_t sign);

#ifdef __cplusplus
}
#endif

#endif