math: logb needs to be correctly rounded

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

diff --git a/src/math/lgammal_r.c b/src/math/lgammal_r.c
index 2ebe612..03e3901 100644 (file)
--- a/src/math/lgammal_r.c
+++ b/src/math/lgammal_r.c
@@ -4,10 +4,12 @@
 
 static struct l_li t[] = {
 #if LDBL_MANT_DIG == 53
-#include "sanity/lgammal_r.h"
+#include "sanity/lgamma_r.h"
+#include "special/lgamma_r.h"
 
 #elif LDBL_MANT_DIG == 64
 #include "sanity/lgammal_r.h"
+#include "special/lgammal_r.h"
 
 #endif
 };
@@ -38,7 +40,7 @@ int main(void)
                        err++;
                }
                d = ulperrl(y, p->y, p->dy);
-               if (!checkulp(d, p->r) || yi != p->i) {
+               if (!checkulp(d, p->r) || (!isnan(p->x) && p->x!=-inf && !(p->e&DIVBYZERO) && yi != p->i)) {
                        printf("%s:%d: %s lgammal_r(%La) want %La,%lld got %La,%d ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
                        err++;