#include <stdint.h>
#include <stdio.h>
#include <stdarg.h>
#include "util.h"
static int test_status;
#define error(...) print(__FILE__, __LINE__, __VA_ARGS__)
static void print(char *f, int l, char *fmt, ...)
{
test_status = 1;
va_list ap;
printf("%s:%d: ", f, l);
va_start(ap, fmt);
vprintf(fmt, ap);
va_end(ap);
}
#define F(n) {#n, n}
static struct {
char *name;
int i;
} te[] = {
#ifdef FE_DIVBYZERO
F(FE_DIVBYZERO),
#endif
#ifdef FE_INEXACT
F(FE_INEXACT),
#endif
#ifdef FE_INVALID
F(FE_INVALID),
#endif
#ifdef FE_OVERFLOW
F(FE_OVERFLOW),
#endif
#ifdef FE_UNDERFLOW
F(FE_UNDERFLOW),
#endif
};
static void test_except()
{
#pragma STDC FENV_ACCESS ON
int i,r;
for (i=0; i < sizeof te/sizeof*te; i++) {
feclearexcept(FE_ALL_EXCEPT);
r = feraiseexcept(te[i].i);
if (r)
error("feraiseexcept(%s) returned %d\n", te[i].name, r);
r = fetestexcept(FE_ALL_EXCEPT);
if (r != te[i].i)
error("feraiseexcept(%s) want %d got %d\n",
te[i].name, te[i].i, r);
}
}
static struct {
char *name;
int i;
} tr[] = {
F(FE_TONEAREST),
#ifdef FE_UPWARD
F(FE_UPWARD),
#endif
#ifdef FE_DOWNWARD
F(FE_DOWNWARD),
#endif
#ifdef FE_TOWARDZERO
F(FE_TOWARDZERO),
#endif
};
static void test_round()
{
#pragma STDC FENV_ACCESS ON
int i,r;
for (i=0; i < sizeof tr/sizeof*tr; i++) {
if (tr[i].i < 0)
error("%s (%d) < 0\n", tr[i].name, tr[i].i);
for (r=0; r < i; r++)
if (tr[r].i == tr[i].i)
error("%s (%d) == %s (%d)\n",
tr[r].name, tr[r].i, tr[i].name, tr[i].i);
}
for (i=0; i < sizeof tr/sizeof*tr; i++) {
r = fesetround(tr[i].i);
if (r != 0)
error("fesetround %d\n", r);
r = fegetround();
if (r != tr[i].i)
error("fegetround %x wanted %x\n", r, tr[i].i);
}
}
/* ieee double precision add operation */
static struct dd_d t[] = {
T(RN, 0x1p+0, 0x1p-52, 0x1.0000000000001p+0, 0x0p+0, 0)
T(RN, 0x1p+0, 0x1p-53, 0x1p+0, -0x1p-1, INEXACT)
T(RN, 0x1p+0, 0x1.01p-53, 0x1.0000000000001p+0, 0x1.fep-2, INEXACT)
T(RN, 0x1p+0, -0x1p-54, 0x1p+0, 0x1p-2, INEXACT)
T(RN, 0x1p+0, -0x1.01p-54, 0x1.fffffffffffffp-1, -0x1.fep-2, INEXACT)
T(RN, -0x1p+0, -0x1p-53, -0x1p+0, 0x1p-1, INEXACT)
T(RN, -0x1p+0, -0x1.01p-53, -0x1.0000000000001p+0, -0x1.fep-2, INEXACT)
T(RN, -0x1p+0, 0x1p-54, -0x1p+0, -0x1p-2, INEXACT)
T(RN, -0x1p+0, 0x1.01p-54, -0x1.fffffffffffffp-1, 0x1.fep-2, INEXACT)
T(RU, 0x1p+0, 0x1p-52, 0x1.0000000000001p+0, 0x0p+0, 0)
T(RU, 0x1p+0, 0x1p-53, 0x1.0000000000001p+0, 0x1p-1, INEXACT)
T(RU, 0x1p+0, 0x1.01p-53, 0x1.0000000000001p+0, 0x1.fep-2, INEXACT)
T(RU, 0x1p+0, -0x1p-54, 0x1p+0, 0x1p-2, INEXACT)
T(RU, 0x1p+0, -0x1.01p-54, 0x1p+0, 0x1.01p-2, INEXACT)
T(RU, -0x1p+0, -0x1p-53, -0x1p+0, 0x1p-1, INEXACT)
T(RU, -0x1p+0, -0x1.01p-53, -0x1p+0, 0x1.01p-1, INEXACT)
T(RU, -0x1p+0, 0x1p-54, -0x1.fffffffffffffp-1, 0x1p-1, INEXACT)
T(RU, -0x1p+0, 0x1.01p-54, -0x1.fffffffffffffp-1, 0x1.fep-2, INEXACT)
T(RD, 0x1p+0, 0x1p-52, 0x1.0000000000001p+0, 0x0p+0, 0)
T(RD, 0x1p+0, 0x1p-53, 0x1p+0, -0x1p-1, INEXACT)
T(RD, 0x1p+0, 0x1.01p-53, 0x1p+0, -0x1.01p-1, INEXACT)
T(RD, 0x1p+0, -0x1p-54, 0x1.fffffffffffffp-1, -0x1p-1, INEXACT)
T(RD, 0x1p+0, -0x1.01p-54, 0x1.fffffffffffffp-1, -0x1.fep-2, INEXACT)
T(RD, -0x1p+0, -0x1p-53, -0x1.0000000000001p+0, -0x1p-1, INEXACT)
T(RD, -0x1p+0, -0x1.01p-53, -0x1.0000000000001p+0, -0x1.fep-2, INEXACT)
T(RD, -0x1p+0, 0x1p-54, -0x1p+0, -0x1p-2, INEXACT)
T(RD, -0x1p+0, 0x1.01p-54, -0x1p+0, -0x1.01p-2, INEXACT)
T(RZ, 0x1p+0, 0x1p-52, 0x1.0000000000001p+0, 0x0p+0, 0)
T(RZ, 0x1p+0, 0x1p-53, 0x1p+0, -0x1p-1, INEXACT)
T(RZ, 0x1p+0, 0x1.01p-53, 0x1p+0, -0x1.01p-1, INEXACT)
T(RZ, 0x1p+0, -0x1p-54, 0x1.fffffffffffffp-1, -0x1p-1, INEXACT)
T(RZ, 0x1p+0, -0x1.01p-54, 0x1.fffffffffffffp-1, -0x1.fep-2, INEXACT)
T(RZ, -0x1p+0, -0x1p-53, -0x1p+0, 0x1p-1, INEXACT)
T(RZ, -0x1p+0, -0x1.01p-53, -0x1p+0, 0x1.01p-1, INEXACT)
T(RZ, -0x1p+0, 0x1p-54, -0x1.fffffffffffffp-1, 0x1p-1, INEXACT)
T(RZ, -0x1p+0, 0x1.01p-54, -0x1.fffffffffffffp-1, 0x1.fep-2, INEXACT)
};
static void test_round_add(void)
{
#pragma STDC FENV_ACCESS ON
double y;
float d;
int i;
struct dd_d *p;
for (i = 0; i < sizeof t/sizeof *t; i++) {
p = t + i;
if (p->r < 0)
continue;
fesetround(p->r);
y = p->x + p->x2;
d = ulperr(y, p->y, p->dy);
if (!checkcr(y, p->y, p->r)) {
printf("%s:%d: %s %a+%a want %a got %a ulperr %.3f = %a + %a\n",
p->file, p->line, rstr(p->r), p->x, p->x2, p->y, y, d, d-p->dy, p->dy);
test_status = 1;
}
}
}
int main(void)
{
test_except();
test_round();
test_round_add();
return test_status;
}