math: logb needs to be correctly rounded

diff --git a/src/math/logb.c b/src/math/logb.c
index c9b90bc..b0ec283 100644 (file)
--- a/src/math/logb.c
+++ b/src/math/logb.c
@@ -33,7 +33,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 logb(%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/logbf.c b/src/math/logbf.c
index cb6d357..7d60218 100644 (file)
--- a/src/math/logbf.c
+++ b/src/math/logbf.c
@@ -33,7 +33,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 logbf(%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/logbl.c b/src/math/logbl.c
index ad0ccff..b314722 100644 (file)
--- a/src/math/logbl.c
+++ b/src/math/logbl.c
@@ -39,7 +39,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 logbl(%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++;