math: ldexp, modf, nearbyint, scalb are correctly rounded

diff --git a/src/math/ldexp.c b/src/math/ldexp.c
index 4ab011f..1d18160 100644 (file)
--- a/src/math/ldexp.c
+++ b/src/math/ldexp.c
@@ -32,7 +32,7 @@ int main(void)
                        err++;
                }
                d = ulperr(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s ldexp(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
                        err++;

diff --git a/src/math/ldexpf.c b/src/math/ldexpf.c
index 01b037c..37a7c03 100644 (file)
--- a/src/math/ldexpf.c
+++ b/src/math/ldexpf.c
@@ -32,7 +32,7 @@ int main(void)
                        err++;
                }
                d = ulperrf(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s ldexpf(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
                        err++;

diff --git a/src/math/ldexpl.c b/src/math/ldexpl.c
index f1c6439..c45e9c5 100644 (file)
--- a/src/math/ldexpl.c
+++ b/src/math/ldexpl.c
@@ -37,7 +37,7 @@ int main(void)
                        err++;
                }
                d = ulperrl(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s ldexpl(%La, %lld) want %La got %La ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
                        err++;

diff --git a/src/math/modf.c b/src/math/modf.c
index 86be04f..f30582d 100644 (file)
--- a/src/math/modf.c
+++ b/src/math/modf.c
@@ -32,7 +32,7 @@ int main(void)
                }
                d = ulperr(y, p->y, p->dy);
                di = ulperr(yi, p->y2, p->dy2);
-               if (!checkulp(d, p->r) || !checkulp(di, p->r)) {
+               if (!checkcr(y, p->y, p->r) || !checkcr(yi, p->y2, p->r)) {
                        printf("%s:%d: %s modf(%a) want %a,%a got %a,%a, ulperr %.3f = %a + %a, %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, p->y2, y, yi, d, d-p->dy, p->dy, di, di-p->dy2, p->dy2);
                        err++;

diff --git a/src/math/modff.c b/src/math/modff.c
index e6ed290..0791dc0 100644 (file)
--- a/src/math/modff.c
+++ b/src/math/modff.c
@@ -32,7 +32,7 @@ int main(void)
                }
                d = ulperr(y, p->y, p->dy);
                di = ulperr(yi, p->y2, p->dy2);
-               if (!checkulp(d, p->r) || !checkulp(di, p->r)) {
+               if (!checkcr(y, p->y, p->r) || !checkcr(yi, p->y2, p->r)) {
                        printf("%s:%d: %s modff(%a) want %a,%a got %a,%a, ulperr %.3f = %a + %a, %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, p->y2, y, yi, d, d-p->dy, p->dy, di, di-p->dy2, p->dy2);
                        err++;

diff --git a/src/math/modfl.c b/src/math/modfl.c
index 72d3b1d..c555cd4 100644 (file)
--- a/src/math/modfl.c
+++ b/src/math/modfl.c
@@ -36,7 +36,7 @@ int main(void)
                }
                d = ulperr(y, p->y, p->dy);
                di = ulperr(yi, p->y2, p->dy2);
-               if (!checkulp(d, p->r) || !checkulp(di, p->r)) {
+               if (!checkcr(y, p->y, p->r) || !checkcr(yi, p->y2, p->r)) {
                        printf("%s:%d: %s modf(%La) want %La,%La got %La,%La, ulperr %.3f = %a + %a, %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, p->y2, y, yi, d, d-p->dy, p->dy, di, di-p->dy2, p->dy2);
                        err++;

diff --git a/src/math/nearbyint.c b/src/math/nearbyint.c
index 927f129..db5a58c 100644 (file)
--- a/src/math/nearbyint.c
+++ b/src/math/nearbyint.c
@@ -32,7 +32,7 @@ int main(void)
                        err++;
                }
                d = ulperr(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s nearbyint(%a) want %a got %a ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy);
                        err++;

diff --git a/src/math/nearbyintf.c b/src/math/nearbyintf.c
index 13de8ed..4e8be0a 100644 (file)
--- a/src/math/nearbyintf.c
+++ b/src/math/nearbyintf.c
@@ -32,7 +32,7 @@ int main(void)
                        err++;
                }
                d = ulperrf(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s nearbyintf(%a) want %a got %a ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy);
                        err++;

diff --git a/src/math/nearbyintl.c b/src/math/nearbyintl.c
index 9bbc94e..a86e83b 100644 (file)
--- a/src/math/nearbyintl.c
+++ b/src/math/nearbyintl.c
@@ -37,7 +37,7 @@ int main(void)
                        err++;
                }
                d = ulperrl(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s nearbyintl(%La) want %La got %La ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy);
                        err++;

diff --git a/src/math/scalb.c b/src/math/scalb.c
index 35ef2ce..95a5533 100644 (file)
--- a/src/math/scalb.c
+++ b/src/math/scalb.c
@@ -32,7 +32,7 @@ int main(void)
                        err++;
                }
                d = ulperr(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s scalb(%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);
                        err++;

diff --git a/src/math/scalbf.c b/src/math/scalbf.c
index 233c14d..8c190cc 100644 (file)
--- a/src/math/scalbf.c
+++ b/src/math/scalbf.c
@@ -32,7 +32,7 @@ int main(void)
                        err++;
                }
                d = ulperrf(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s scalbf(%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);
                        err++;