From 9c6efc03eba4568ba9f1672a2bb0c50a047edd8e Mon Sep 17 00:00:00 2001
From: Szabolcs Nagy <nsz@port70.net>
Date: Thu, 1 Nov 2012 14:59:36 +0100
Subject: [PATCH] math: ldexp, modf, nearbyint, scalb are correctly rounded

---
 src/math/ldexp.c      | 2 +-
 src/math/ldexpf.c     | 2 +-
 src/math/ldexpl.c     | 2 +-
 src/math/modf.c       | 2 +-
 src/math/modff.c      | 2 +-
 src/math/modfl.c      | 2 +-
 src/math/nearbyint.c  | 2 +-
 src/math/nearbyintf.c | 2 +-
 src/math/nearbyintl.c | 2 +-
 src/math/scalb.c      | 2 +-
 src/math/scalbf.c     | 2 +-
 11 files changed, 11 insertions(+), 11 deletions(-)

diff --git a/src/math/ldexp.c b/src/math/ldexp.c
index 4ab011f..1d18160 100644
--- 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
--- 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
--- 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
--- 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
--- 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
--- 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
--- 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
--- 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
--- 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
--- 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
--- 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++;
-- 
2.20.1