#include <stdio.h>
#include <math.h>
#include "test.h"

static struct {
	double x;
	double yf;
	double yi;
} t[] = {
	0.7, 0x1.6666666666666p-1, 0.0,
	-0.7, -0x1.6666666666666p-1, -0.0,
	1.7, 0x1.6666666666666p-1, 1.0,
	-1.7, -0x1.6666666666666p-1, -1.0,
	2.0, 0.0, 2.0,
	-2.0, -0.0, -2.0,
	-0x1p99, -0.0, -0x1p99,
	-0.0, -0.0, -0.0,
	INFINITY, 0.0, INFINITY,
	-INFINITY, -0.0, -INFINITY,
	NAN, NAN, NAN,
};

static struct {
	float x;
	float yf;
	float yi;
} tf[] = {
	0.7, 0x1.666666p-1, 0.0,
	-0.7, -0x1.666666p-1, -0.0,
	1.7, 0x1.666668p-1, 1.0,
	-1.7, -0x1.666668p-1, -1.0,
	2.0, 0.0, 2.0,
	-2.0, -0.0, -2.0,
	-0x1p99, -0.0, -0x1p99,
	-0.0, -0.0, -0.0,
	INFINITY, 0.0, INFINITY,
	-INFINITY, -0.0, -INFINITY,
	NAN, NAN, NAN,
};

static struct {
	long double x;
	long double yf;
	long double yi;
} tl[] = {
	0.7L, 0x1.6666666666666666p-1L, 0.0,
	-0.7L, -0x1.6666666666666666p-1L, -0.0,
	1.7L, 0x1.6666666666666668p-1L, 1.0,
	-1.7L, -0x1.6666666666666668p-1L, -1.0,
	2.0, 0, 2.0,
	-2.0, -0.0, -2.0,
	-0x1p99, -0.0, -0x1p99,
	-0.0, -0.0, -0.0,
	INFINITY, 0.0, INFINITY,
	-INFINITY, -0.0, -INFINITY,
	NAN, NAN, NAN,
};

#define eq(a, b) ((isnan(a) && isnan(b)) || ((a) == (b) && signbit(a) == signbit(b)))

void test_modf()
{
	double yf, yi;
	int i;
	for (i = 0; i < sizeof t/sizeof *t; i++) {
		yf = modf(t[i].x, &yi);
		if (!eq(yf,t[i].yf) || !eq(yi,t[i].yi))
			error("modf(%a) want %a %a got %a %a\n", t[i].x, t[i].yf, t[i].yi, yf, yi);
	}
}
void test_modff()
{
	float yf, yi;
	int i;
	for (i = 0; i < sizeof tf/sizeof *tf; i++) {
		yf = modff(tf[i].x, &yi);
		if (!eq(yf,tf[i].yf) || !eq(yi,tf[i].yi))
			error("modff(%a) want %a %a got %a %a\n", tf[i].x, tf[i].yf, tf[i].yi, yf, yi);
	}
}
void test_modfl()
{
	long double yf, yi;
	int i;
	for (i = 0; i < sizeof tl/sizeof *tl; i++) {
		yf = modfl(tl[i].x, &yi);
		if (!eq(yf,tl[i].yf) || !eq(yi,tl[i].yi))
			error("modfl(%La) want %La %La got %La %La\n", tl[i].x, tl[i].yf, tl[i].yi, yf, yi);
	}
}

void bench_modf_small(int N)
{
	int i;
	double yi;
	volatile double yf;

	for (i = 0; i < N; i++)
		yf = modf(1234.5678, &yi);
}

void bench_modf_negint(int N)
{
	int i;
	double yi;
	volatile double yf;

	for (i = 0; i < N; i++)
		yf = modf(-1234.0, &yi);
}

void bench_modf_large(int N)
{
	int i;
	double yi;
	volatile double yf;

	for (i = 0; i < N; i++)
		yf = modf(1.2345678e300, &yi);
}


void bench_modff_small(int N)
{
	int i;
	float yi;
	volatile float yf;

	for (i = 0; i < N; i++)
		yf = modff(1234.5678, &yi);
}

void bench_modff_negint(int N)
{
	int i;
	float yi;
	volatile float yf;

	for (i = 0; i < N; i++)
		yf = modff(-1234.0, &yi);
}

void bench_modff_large(int N)
{
	int i;
	float yi;
	volatile float yf;

	for (i = 0; i < N; i++)
		yf = modff(1.2345678e100, &yi);
}