#include "test.h"

#include <math.h>

#include <stdint.h>

#include <fenv.h>

#include <stdio.h>

#include <float.h>

static struct {

	int flag;

	char *s;

} eflags[] = {

	{FE_INVALID, "FE_INVALID"},

	{FE_DIVBYZERO, "FE_DIVBYZERO"},

	{FE_OVERFLOW, "FE_OVERFLOW"},

	{FE_UNDERFLOW, "FE_UNDERFLOW"},

	{FE_INEXACT, "FE_INEXACT"},

};

static int ne = sizeof eflags / sizeof *eflags;

static struct {

	int flag;

	char *s;

} rflags[] = {

	{FE_TONEAREST,"FE_TONEAREST,"},

	{FE_DOWNWARD,"FE_DOWNWARD,"},

	{FE_UPWARD,"FE_UPWARD,"},

	{FE_TOWARDZERO,"FE_TOWARDZERO,"},

};

static int nr = sizeof rflags / sizeof *rflags;

void printexcept(int f) {

	int i, all=0;

	for (i = 0; i < ne; i++)

		if (f & eflags[i].flag) {

			printf("%s%s", all ? "|" : "", eflags[i].s);

			all |= eflags[i].flag;

		}

	if (all != f) {

		printf("%s%d", all ? "|" : "", f & ~all);

		all = f;

	}

	printf("%s,", all ? "" : "0");

}

void printround(int f) {

	int i;

	for (i = 0; i < nr; i++)

		if (f == rflags[i].flag) {

			printf("%s ", rflags[i].s);

			return;

		}

	printf("%d, ", f);

}

/*

struct {double x;} t[] = {

0.0,

0.25,

-0.25,

0.5,

-0.5,

0.75,

-0.75,

1.0,

-1.0,

1.25,

-1.25,

0x1p30,

-0x1p30,

0x1p31-1,

-0x1p31+1,

0x1p31,

-0x1p31,

0x1p31+1,

-0x1p31-1,

0x1p31-0.5,

-0x1p31+0.5,

0x1p31+0.5,

-0x1p31-0.5,

0x1p32,

-0x1p32,

0x1p32 - 0.5,

-0x1p32 + 0.5,

};

void test_gendata()

{

	int f, i, j;

	long n;

	for (i = 0; i < sizeof t/sizeof *t; i++) {

		for (j = 0; j < nr; j++) {

			fesetround(rflags[j].flag);

			feclearexcept(FE_ALL_EXCEPT);

			n = lrint(t[i].x);

			f = fetestexcept(FE_ALL_EXCEPT);

			if (f == FE_INVALID)

				n = 0;

			printround(rflags[j].flag);

			printf("%18a, %11ld, ", t[i].x, n);

			printexcept(f);

			printf("\n");

		}

	}

}

*/

static struct {

	int round;

	double x;

	long n;

	int except;

} t[] = {

FE_TONEAREST,             0x0p+0,           0, 0,

FE_DOWNWARD,              0x0p+0,           0, 0,

FE_UPWARD,                0x0p+0,           0, 0,

FE_TOWARDZERO,            0x0p+0,           0, 0,

FE_TONEAREST,             0x1p-2,           0, FE_INEXACT,

FE_DOWNWARD,              0x1p-2,           0, FE_INEXACT,

FE_UPWARD,                0x1p-2,           1, FE_INEXACT,

FE_TOWARDZERO,            0x1p-2,           0, FE_INEXACT,

FE_TONEAREST,            -0x1p-2,           0, FE_INEXACT,

FE_DOWNWARD,             -0x1p-2,          -1, FE_INEXACT,

FE_UPWARD,               -0x1p-2,           0, FE_INEXACT,

FE_TOWARDZERO,           -0x1p-2,           0, FE_INEXACT,

FE_TONEAREST,             0x1p-1,           0, FE_INEXACT,

FE_DOWNWARD,              0x1p-1,           0, FE_INEXACT,

FE_UPWARD,                0x1p-1,           1, FE_INEXACT,

FE_TOWARDZERO,            0x1p-1,           0, FE_INEXACT,

FE_TONEAREST,            -0x1p-1,           0, FE_INEXACT,

FE_DOWNWARD,             -0x1p-1,          -1, FE_INEXACT,

FE_UPWARD,               -0x1p-1,           0, FE_INEXACT,

FE_TOWARDZERO,           -0x1p-1,           0, FE_INEXACT,

FE_TONEAREST,           0x1.8p-1,           1, FE_INEXACT,

FE_DOWNWARD,            0x1.8p-1,           0, FE_INEXACT,

FE_UPWARD,              0x1.8p-1,           1, FE_INEXACT,

FE_TOWARDZERO,          0x1.8p-1,           0, FE_INEXACT,

FE_TONEAREST,          -0x1.8p-1,          -1, FE_INEXACT,

FE_DOWNWARD,           -0x1.8p-1,          -1, FE_INEXACT,

FE_UPWARD,             -0x1.8p-1,           0, FE_INEXACT,

FE_TOWARDZERO,         -0x1.8p-1,           0, FE_INEXACT,

FE_TONEAREST,             0x1p+0,           1, 0,

FE_DOWNWARD,              0x1p+0,           1, 0,

FE_UPWARD,                0x1p+0,           1, 0,

FE_TOWARDZERO,            0x1p+0,           1, 0,

FE_TONEAREST,            -0x1p+0,          -1, 0,

FE_DOWNWARD,             -0x1p+0,          -1, 0,

FE_UPWARD,               -0x1p+0,          -1, 0,

FE_TOWARDZERO,           -0x1p+0,          -1, 0,

FE_TONEAREST,           0x1.4p+0,           1, FE_INEXACT,

FE_DOWNWARD,            0x1.4p+0,           1, FE_INEXACT,

FE_UPWARD,              0x1.4p+0,           2, FE_INEXACT,

FE_TOWARDZERO,          0x1.4p+0,           1, FE_INEXACT,

FE_TONEAREST,          -0x1.4p+0,          -1, FE_INEXACT,

FE_DOWNWARD,           -0x1.4p+0,          -2, FE_INEXACT,

FE_UPWARD,             -0x1.4p+0,          -1, FE_INEXACT,

FE_TOWARDZERO,         -0x1.4p+0,          -1, FE_INEXACT,

FE_TONEAREST,            0x1p+30,  1073741824, 0,

FE_DOWNWARD,             0x1p+30,  1073741824, 0,

FE_UPWARD,               0x1p+30,  1073741824, 0,

FE_TOWARDZERO,           0x1p+30,  1073741824, 0,

FE_TONEAREST,           -0x1p+30, -1073741824, 0,

FE_DOWNWARD,            -0x1p+30, -1073741824, 0,

FE_UPWARD,              -0x1p+30, -1073741824, 0,

FE_TOWARDZERO,          -0x1p+30, -1073741824, 0,

FE_TONEAREST,   0x1.fffffffcp+30,  2147483647, 0,

FE_DOWNWARD,    0x1.fffffffcp+30,  2147483647, 0,

FE_UPWARD,      0x1.fffffffcp+30,  2147483647, 0,

FE_TOWARDZERO,  0x1.fffffffcp+30,  2147483647, 0,

FE_TONEAREST,  -0x1.fffffffcp+30, -2147483647, 0,

FE_DOWNWARD,   -0x1.fffffffcp+30, -2147483647, 0,

FE_UPWARD,     -0x1.fffffffcp+30, -2147483647, 0,

FE_TOWARDZERO, -0x1.fffffffcp+30, -2147483647, 0,

FE_TONEAREST,            0x1p+31,           0, FE_INVALID,

FE_DOWNWARD,             0x1p+31,           0, FE_INVALID,

FE_UPWARD,               0x1p+31,           0, FE_INVALID,

FE_TOWARDZERO,           0x1p+31,           0, FE_INVALID,

FE_TONEAREST,           -0x1p+31, -2147483648, 0,

FE_DOWNWARD,            -0x1p+31, -2147483648, 0,

FE_UPWARD,              -0x1p+31, -2147483648, 0,

FE_TOWARDZERO,          -0x1p+31, -2147483648, 0,

FE_TONEAREST,   0x1.00000002p+31,           0, FE_INVALID,

FE_DOWNWARD,    0x1.00000002p+31,           0, FE_INVALID,

FE_UPWARD,      0x1.00000002p+31,           0, FE_INVALID,

FE_TOWARDZERO,  0x1.00000002p+31,           0, FE_INVALID,

FE_TONEAREST,  -0x1.00000002p+31,           0, FE_INVALID,

FE_DOWNWARD,   -0x1.00000002p+31,           0, FE_INVALID,

FE_UPWARD,     -0x1.00000002p+31,           0, FE_INVALID,

FE_TOWARDZERO, -0x1.00000002p+31,           0, FE_INVALID,

FE_TONEAREST,   0x1.fffffffep+30,           0, FE_INVALID,

FE_DOWNWARD,    0x1.fffffffep+30,  2147483647, FE_INEXACT,

FE_UPWARD,      0x1.fffffffep+30,           0, FE_INVALID,

FE_TOWARDZERO,  0x1.fffffffep+30,  2147483647, FE_INEXACT,

FE_TONEAREST,  -0x1.fffffffep+30, -2147483648, FE_INEXACT,

FE_DOWNWARD,   -0x1.fffffffep+30, -2147483648, FE_INEXACT,

FE_UPWARD,     -0x1.fffffffep+30, -2147483647, FE_INEXACT,

FE_TOWARDZERO, -0x1.fffffffep+30, -2147483647, FE_INEXACT,

FE_TONEAREST,   0x1.00000001p+31,           0, FE_INVALID,

FE_DOWNWARD,    0x1.00000001p+31,           0, FE_INVALID,

FE_UPWARD,      0x1.00000001p+31,           0, FE_INVALID,

FE_TOWARDZERO,  0x1.00000001p+31,           0, FE_INVALID,

FE_TONEAREST,  -0x1.00000001p+31, -2147483648, FE_INEXACT,

FE_DOWNWARD,   -0x1.00000001p+31,           0, FE_INVALID,

FE_UPWARD,     -0x1.00000001p+31, -2147483648, FE_INEXACT,

FE_TOWARDZERO, -0x1.00000001p+31, -2147483648, FE_INEXACT,

FE_TONEAREST,            0x1p+32,           0, FE_INVALID,

FE_DOWNWARD,             0x1p+32,           0, FE_INVALID,

FE_UPWARD,               0x1p+32,           0, FE_INVALID,

FE_TOWARDZERO,           0x1p+32,           0, FE_INVALID,

FE_TONEAREST,           -0x1p+32,           0, FE_INVALID,

FE_DOWNWARD,            -0x1p+32,           0, FE_INVALID,

FE_UPWARD,              -0x1p+32,           0, FE_INVALID,

FE_TOWARDZERO,          -0x1p+32,           0, FE_INVALID,

FE_TONEAREST,   0x1.ffffffffp+31,           0, FE_INVALID,

FE_DOWNWARD,    0x1.ffffffffp+31,           0, FE_INVALID,

FE_UPWARD,      0x1.ffffffffp+31,           0, FE_INVALID,

FE_TOWARDZERO,  0x1.ffffffffp+31,           0, FE_INVALID,

FE_TONEAREST,  -0x1.ffffffffp+31,           0, FE_INVALID,

FE_DOWNWARD,   -0x1.ffffffffp+31,           0, FE_INVALID,

FE_UPWARD,     -0x1.ffffffffp+31,           0, FE_INVALID,

FE_TOWARDZERO, -0x1.ffffffffp+31,           0, FE_INVALID,

};

void test_lrint()

{

	int f, i;

	long n;

	for (i = 0; i < sizeof t/sizeof *t; i++) {

		fesetround(t[i].round);

		feclearexcept(FE_ALL_EXCEPT);

		n = lrint(t[i].x);

		f = fetestexcept(FE_ALL_EXCEPT);

		if (t[i].except != FE_INVALID && n != t[i].n) {

			error("round=");

			printround(t[i].round);

			printf("lrint(%a) want %ld got %ld\n", t[i].x, t[i].n, n);

		}

		if (f != t[i].except) {

			error("round=");

			printround(t[i].round);

			printf("lrint(%a)==%ld want except=", t[i].x, t[i].n);

			printexcept(t[i].except);

			printf(" got except=");

			printexcept(f);

			printf("\n");

		}

	}

}

void bench_lrint_simple(int N)

{

	int i;

	volatile int n;

	for (i = 0; i < N; i++) {

		n = lrint(1.25);

	}

}

void bench_lrint_hard(int N)

{

	int i;

	volatile int n;

	for (i = 0; i < N; i++) {

//		feclearexcept(FE_ALL_EXCEPT);

//		n = lrint(1.5);

//		n = lrint(0x1p32);

//		n = lrint(-0x1p31);

		n = lrint(0x1p31+0.5);

	}

}