math: correct-rounding check in remainder functions and 1ulp check in hypot

[libc-test] / src / math / y1.c

diff --git a/src/math/y1.c b/src/math/y1.c
index 95d971c..6bf5135 100644 (file)
--- a/src/math/y1.c
+++ b/src/math/y1.c
@@ -1,6 +1,8 @@
+#define _BSD_SOURCE 1
+#define _XOPEN_SOURCE 700
 #include <stdint.h>
 #include <stdio.h>
-#include "util.h"
+#include "mtest.h"
 
 static struct d_d t[] = {
 #include "sanity/y1.h"
@@ -33,7 +35,7 @@ int main(void)
                        err++;
                }
                d = ulperr(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if ((!(p->x < 0) && !checkulp(d, p->r)) || (p->x < 0 && !isnan(y) && y != -inf)) {
                        printf("%s:%d: %s y1(%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++;