-#include "test.h"
-#include <math.h>
#include <stdint.h>
-#include <fenv.h>
#include <stdio.h>
-#include <float.h>
+#include "mtest.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");
-}
+static struct d_i t[] = {
+#include "sanity/lrint.h"
+#include "special/lrint.h"
-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 main(void)
{
- int f, i, j;
- long n;
+ #pragma STDC FENV_ACCESS ON
+ long long yi;
+ int e, i, err = 0;
+ struct d_i *p;
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");
- }
- }
-}
-*/
+ p = t + i;
-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, j;
- long n;
-
- for (i = 0; i < sizeof t/sizeof *t; i++) {
- fesetround(t[i].round);
+ if (p->r < 0)
+ continue;
+ fesetround(p->r);
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);
+ yi = lrint(p->x);
+ e = fetestexcept(INEXACT|INVALID|DIVBYZERO|UNDERFLOW|OVERFLOW);
+
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lrint(%a)=%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
}
- 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");
+ if (!(p->e&INVALID) && yi != p->i) {
+ printf("%s:%d: %s lrint(%a) want %lld got %lld\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, yi);
+ err++;
}
}
-}
-
-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);
- }
+ return !!err;
}