#ifndef _MATH
#define _MATH

#include <math.h>

#define PI 3.14159265358979323846

inline float sincf(const float x) {return (x==0)?1:sinf(x)/x;}

#define roundf(x) floorf((x)+0.5f)

#ifndef min
#define min(x,y) ((x<y)?x:y)
#endif

#ifndef max
#define max(x,y) ((x>y)?x:y)
#endif

inline float intAdjustf(const float x,const float diff=0.01f)
{
	register float f=roundf(x);
	return (fabsf(f-x)<diff)?f:x;
}

// code from Id Software
inline float rsqrtf(const float v)
{
	union {float vh; int i0;};
	union {float vr; int i1;};

	vh=v*0.5f;
	i1=0x5f3759df-(i0>>1);
	return vr*(1.5f-vh*vr*vr);
}

inline unsigned int getClosestPowerOfTwo(const unsigned int x)
{
	register unsigned int i=1;
	while (i<x) i+=i;
	if (4*x<3*i) i>>=1;
	return i;
}

inline unsigned int getUpperPowerOfTwo(const unsigned int x)
{
	register unsigned int i=1;
	while (i<x) i+=i;
	return i;
}

inline unsigned int getLowerPowerOfTwo(const unsigned int x)
{
	register unsigned int i=1;
	while (i<=x) i+=i;
	return i>>1;
}

inline int round(float x)
{
	return (int)(x>0?x+0.5:x-0.5);
}

#endif