X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flgammal.c;h=5a95266ea1aa81a1c6c4fcfc6b11dfc85f25c60d;hb=d05da15802af606664c6879a7463bae1c7d62223;hp=37047f968efb78ce45bca89f828f4657886a1cad;hpb=7b5ccd5b5297e4a0e413a1517699fbfb5bb575be;p=libc-test diff --git a/src/math/lgammal.c b/src/math/lgammal.c index 37047f9..5a95266 100644 --- a/src/math/lgammal.c +++ b/src/math/lgammal.c @@ -1,3 +1,4 @@ +#define _DEFAULT_SOURCE 1 #define _BSD_SOURCE 1 #define _XOPEN_SOURCE 700 #include @@ -43,7 +44,8 @@ int main(void) err++; } d = ulperrl(y, p->y, p->dy); - if (!checkulp(d, p->r) || (!isnan(p->x) && p->x!=-inf && !(p->e&DIVBYZERO) && yi != p->i)) { + // TODO: 2 ulp errors allowed + if ((p->r==RN && fabs(d) > 2) || (!isnan(p->x) && p->x!=-inf && !(p->e&DIVBYZERO) && yi != p->i)) { printf("%s:%d: %s lgammal(%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++;