#include <stdio.h>

#include <math.h>

#include "test.h"

static int check(double a, double b)

{

	double d = a - b;

	return fabs(d) <= 0x1p-52*fabs(a);

}

static int checkf(float a, float b)

{

	float d = a - b;

	return fabsf(d) <= 0x1p-23*fabsf(a);

}

static int checkl(long double a, long double b)

{

	long double d = a - b;

	return fabsl(d) <= 0x1p-63L*fabsl(a);

}

#define D(fx, y) do{ \

	volatile double yy = fx; \

	volatile long double yl = fx; \

	if (!check(yy, y)) \

		error("%s got %a = %.21e want %a = %.21e\n", #fx, yy, yy, y, y); \

	if (yy != yl) \

		error("float-store issue: (double)%s = %a (long double)%s = %La\n", #fx, yy, #fx, yl); \

}while(0)

#define F(fx, y) do{ \

	volatile float yy = fx; \

	volatile long double yl = fx; \

	if (!checkf(yy, y)) \

		error("%s got %a = %.21e want %a = %.21e\n", #fx, yy, yy, y, y); \

	if (yy != yl) \

		error("float-store issue: (float)%s = %a (long double)%s = %La\n", #fx, yy, #fx, yl); \

}while(0)

#define L(fx, y) do{ \

	long double yy = fx; \

	if (!checkl(yy, y)) \

		error("%s got %La = %.21Le want %La = %.21Le\n", #fx, yy, yy, y, y); \

}while(0)

#define I(fx, y) do{ \

	int yy = fx; \

	if (yy != y) \

		error("%s got %d want %d\n", #fx, yy, y); \

}while(0)

#define IL(fx, y) do{ \

	long yy = fx; \

	if (yy != y) \

		error("%s got %ld want %ld\n", #fx, yy, y); \

}while(0)

#define ILL(fx, y) do{ \

	long long yy = fx; \

	if (yy != y) \

		error("%s got %lld want %lld\n", #fx, yy, y); \

}while(0)

void test_math_sanity()

{

	int i;

	double q;

	float qf;

	long double ql;

D(acos(0.7), 0x1.973e83f5d5c9bp-1);

F(acosf(0.7f), 0x1.973e84p-1);

L(acosl(0.7L), 0x1.973e83f5d5c9aaf8p-1L);

D(acosh(1.7), 0x1.1f8c10d010fe6p+0);

F(acoshf(1.7f), 0x1.1f8c12p+0);

L(acoshl(1.7L), 0x1.1f8c10d010fe5d96p+0L);

D(asin(0.7), 0x1.8d00e692afd95p-1);

F(asinf(0.7f), 0x1.8d00e6p-1);

L(asinl(0.7L), 0x1.8d00e692afd95ddap-1L);

D(asinh(0.7), 0x1.4e2a4fe9085ddp-1);

F(asinhf(0.7f), 0x1.4e2a5p-1);

L(asinhl(0.7L), 0x1.4e2a4fe9085dd732p-1L);

D(atan(0.7), 0x1.38b112d7bd4adp-1);

F(atanf(0.7f), 0x1.38b112p-1);

L(atanl(0.7L), 0x1.38b112d7bd4ad786p-1L);

D(atan2(0.7, 1.0), 0x1.38b112d7bd4adp-1);

F(atan2f(0.7f, 1.0f), 0x1.38b112p-1);

L(atan2l(0.7L, 1.0L), 0x1.38b112d7bd4ad786p-1L);

D(atanh(0.7), 0x1.bc0ed0947fbe8p-1);

F(atanhf(0.7f), 0x1.bc0edp-1);

L(atanhl(0.7L), 0x1.bc0ed0947fbe9068p-1L);

D(cbrt(0.7), 0x1.c69b5a72f1a99p-1);

F(cbrtf(0.7f), 0x1.c69b5ap-1);

L(cbrtl(0.7L), 0x1.c69b5a72f1a99902p-1L);

D(ceil(0.7), 0x1p+0);

F(ceilf(0.7f), 0x1p+0);

L(ceill(0.7L), 0x1p+0L);

D(copysign(0.7, -1.0), -0x1.6666666666666p-1);

F(copysignf(0.7f, -1.0f), -0x1.666666p-1);

L(copysignl(0.7L, -1.0L), -0x1.6666666666666666p-1L);

D(cos(0.7), 0x1.87996529f9d93p-1);

F(cosf(0.7f), 0x1.879966p-1);

L(cosl(0.7L), 0x1.87996529f9d92618p-1L);

D(cosh(0.7), 0x1.4152c1862342fp+0);

F(coshf(0.7f), 0x1.4152c2p+0);

L(coshl(0.7L), 0x1.4152c1862342ef8ep+0L);

D(erf(0.7), 0x1.5b08c21171646p-1);

F(erff(0.7f), 0x1.5b08c2p-1);

L(erfl(0.7L), 0x1.5b08c21171646544p-1L);

D(erfc(0.7), 0x1.49ee7bdd1d374p-2);

F(erfcf(0.7f), 0x1.49ee7cp-2);

L(erfcl(0.7L), 0x1.49ee7bdd1d373576p-2L);

D(exp(0.7), 0x1.01c2a61268987p+1);

F(expf(0.7f), 0x1.01c2a6p+1);

L(expl(0.7L), 0x1.01c2a61268986bfep+1L);

D(exp2(0.7), 0x1.9fdf8bcce533dp+0);

F(exp2f(0.7f), 0x1.9fdf8cp+0);

L(exp2l(0.7L), 0x1.9fdf8bcce533d72p+0L);

D(expm1(0.7), 0x1.03854c24d130dp+0);

F(expm1f(0.7f), 0x1.03854cp+0);

L(expm1l(0.7L), 0x1.03854c24d130d7fep+0L);

D(fabs(-0.7), 0x1.6666666666666p-1);

F(fabsf(-0.7f), 0x1.666666p-1);

L(fabsl(-0.7L), 0x1.6666666666666666p-1L);

D(fdim(0.7, 0.5), 0x1.9999999999998p-3);

F(fdimf(0.7f, 0.5f), 0x1.999998p-3);

L(fdiml(0.7L, 0.5L), 0x1.9999999999999998p-3L);

D(floor(0.7), 0x0p+0);

F(floorf(0.7f), 0x0p+0);

L(floorl(0.7L), 0x0p+0L);

D(fma(0.7, 2.0, 0.1), 0x1.8p+0);

F(fmaf(0.7f, 2.0f, 0.1f), 0x1.8p+0);

L(fmal(0.7L, 2.0L, 0.1L), 0x1.8p+0L);

D(fmax(0.7, 0.5), 0x1.6666666666666p-1);

F(fmaxf(0.7f, 0.5), 0x1.666666p-1);

L(fmaxl(0.7L, 0.5), 0x1.6666666666666666p-1L);

D(fmin(0.7, 0.5), 0x1p-1);

F(fminf(0.7f, 0.5f), 0x1p-1);

L(fminl(0.7L, 0.5L), 0x1p-1L);

D(fmod(0.7, 0.5), 0x1.9999999999998p-3);

F(fmodf(0.7f, 0.5f), 0x1.999998p-3);

L(fmodl(0.7L, 0.5L), 0x1.9999999999999998p-3L);

D(frexp(0.7, &i), 0x1.6666666666666p-1);

F(frexpf(0.7f, &i), 0x1.666666p-1);

L(frexpl(0.7L, &i), 0x1.6666666666666666p-1L);

D(hypot(0.7, 1.0), 0x1.387ce204a35d2p+0);

F(hypotf(0.7f, 1.0f), 0x1.387ce2p+0);

L(hypotl(0.7L, 1.0L), 0x1.387ce204a35d1ff6p+0L);

I(ilogb(0.7), -1);

I(ilogbf(0.7f), -1);

I(ilogbl(0.7L), -1);

D(j0(0.7), 0x1.c32cc34b8cc59p-1);

F(j0f(0.7f), 0x1.c32cc4p-1);

D(j1(0.7), 0x1.50e44279c0546p-2);

F(j1f(0.7f), 0x1.50e442p-2);

D(jn(2, 0.7), 0x1.e195286f3b2fbp-5);

F(jnf(2, 0.7f), 0x1.e19528p-5);

D(ldexp(0.7, 3), 0x1.6666666666666p+2);

F(ldexpf(0.7f, 3), 0x1.666666p+2);

L(ldexpl(0.7L, 3), 0x1.6666666666666666p+2L);

D(lgamma(0.7), 0x1.0b20c891cde73p-2);

F(lgammaf(0.7f), 0x1.0b20cap-2);

L(lgammal(0.7L), 0x1.0b20c891cde72846p-2L);

D(lgamma_r(0.7, &i), 0x1.0b20c891cde73p-2);

F(lgammaf_r(0.7f, &i), 0x1.0b20cap-2);

L(lgammal_r(0.7L, &i), 0x1.0b20c891cde72846p-2L);

ILL(llrint(0.7), 1);

ILL(llrintf(0.7f), 1);

ILL(llrintl(0.7l), 1);

ILL(llround(0.7), 1);

ILL(llroundf(0.7f), 1);

ILL(llroundl(0.7L), 1);

D(log(0.7), -0x1.6d3c324e13f5p-2);

F(logf(0.7f), -0x1.6d3c34p-2);

L(logl(0.7L), -0x1.6d3c324e13f4ec54p-2L);

D(log10(0.7), -0x1.3d3d3d21ccf04p-3);

F(log10f(0.7f), -0x1.3d3d3ep-3);

L(log10l(0.7L), -0x1.3d3d3d21ccf035a6p-3L);

D(log1p(0.7), 0x1.0fae81914a991p-1);

F(log1pf(0.7f), 0x1.0fae82p-1);

L(log1pl(0.7L), 0x1.0fae81914a991308p-1L);

D(log2(0.7), -0x1.0776228967d13p-1);

F(log2f(0.7f), -0x1.077624p-1);

L(log2l(0.7L), -0x1.0776228967d1218cp-1L);

D(logb(0.7), -0x1p+0);

F(logbf(0.7f), -0x1p+0);

L(logbl(0.7L), -0x1p+0L);

IL(lrint(0.7), 1);

IL(lrintf(0.7f), 1);

IL(lrintl(0.7l), 1);

IL(lround(0.7), 1);

IL(lroundf(0.7f), 1);

IL(lroundl(0.7L), 1);

D(modf(0.7, &q), 0x1.6666666666666p-1);

F(modff(0.7f, &qf), 0x1.666666p-1);

L(modfl(0.7L, &ql), 0x1.6666666666666666p-1L);

D(nearbyint(0.7), 0x1p+0);

F(nearbyintf(0.7f), 0x1p+0);

L(nearbyintl(0.7L), 0x1p+0L);

D(nextafter(0.7, 1.0), 0x1.6666666666667p-1);

F(nextafterf(0.7f, 1.0f), 0x1.666668p-1);

L(nextafterl(0.7L, 1.0L), 0x1.6666666666666667p-1L);

D(nexttoward(0.7, 1.0L), 0x1.6666666666667p-1);

F(nexttowardf(0.7f, 1.0L), 0x1.666668p-1);

L(nexttowardl(0.7L, 1.0L), 0x1.6666666666666667p-1L);

D(pow(0.7, 1.5), 0x1.2bdbe460916ep-1);

F(powf(0.7f, 1.5f), 0x1.2bdbe4p-1);

L(powl(0.7L, 1.5L), 0x1.2bdbe460916e0b5p-1L);

D(remainder(0.7, 0.5), 0x1.9999999999998p-3);

F(remainderf(0.7f, 0.5f), 0x1.999998p-3);

L(remainderl(0.7L, 0.5L), 0x1.9999999999999998p-3L);

D(remquo(0.7, 0.5, &i), 0x1.9999999999998p-3);

F(remquof(0.7f, 0.5f, &i), 0x1.999998p-3);

L(remquol(0.7L, 0.5L, &i), 0x1.9999999999999998p-3L);

D(rint(0.7), 0x1p+0);

F(rintf(0.7f), 0x1p+0);

L(rintl(0.7L), 0x1p+0L);

D(round(0.7), 0x1p+0);

F(roundf(0.7f), 0x1p+0);

L(roundl(0.7L), 0x1p+0L);

D(scalb(0.7, 3), 0x1.6666666666666p+2);

F(scalbf(0.7f, 3), 0x1.666666p+2);

D(scalbln(0.7, 3), 0x1.6666666666666p+2);

F(scalblnf(0.7f, 3), 0x1.666666p+2);

L(scalblnl(0.7L, 3), 0x1.6666666666666666p+2L);

D(scalbn(0.7, 3), 0x1.6666666666666p+2);

F(scalbnf(0.7f, 3), 0x1.666666p+2);

L(scalbnl(0.7L, 3), 0x1.6666666666666666p+2L);

D(sin(0.7), 0x1.49d6e694619b8p-1);

F(sinf(0.7f), 0x1.49d6e6p-1);

L(sinl(0.7L), 0x1.49d6e694619b854ep-1L);

D(sinh(0.7), 0x1.8465153d5bdbdp-1);

F(sinhf(0.7f), 0x1.846514p-1);

L(sinhl(0.7L), 0x1.8465153d5bdbd0dep-1L);

D(sqrt(0.7), 0x1.ac5eb3f7ab2f8p-1);

F(sqrtf(0.7f), 0x1.ac5eb4p-1);

L(sqrtl(0.7L), 0x1.ac5eb3f7ab2f7de2p-1L);

D(tan(0.7), 0x1.af406c2fc78aep-1);

F(tanf(0.7f), 0x1.af406cp-1);

L(tanl(0.7L), 0x1.af406c2fc78ae54cp-1L);

D(tanh(0.7), 0x1.356fb17af2e91p-1);

F(tanhf(0.7f), 0x1.356fb2p-1);

L(tanhl(0.7L), 0x1.356fb17af2e9100ap-1L);

D(tgamma(0.7), 0x1.4c4d5ab21ea23p+0);

F(tgammaf(0.7f), 0x1.4c4d5cp+0);

L(tgammal(0.7L), 0x1.4c4d5ab21ea22798p+0L);

D(trunc(0.7), 0x0p+0);

F(truncf(0.7f), 0x0p+0);

L(truncl(0.7L), 0x0p+0L);

D(y0(0.7), -0x1.867b559ffc715p-3);

F(y0f(0.7f), -0x1.867b58p-3);

D(y1(0.7), -0x1.1a6e956728d35p+0);

F(y1f(0.7f), -0x1.1a6e96p+0);

D(yn(2, 0.7), -0x1.7b11b25df166ep+1);

F(ynf(2, 0.7f), -0x1.7b11b4p+1);

}