X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flgammal.c;h=9050913cd9484b3d5bbea7a72fb4fbeff5a98589;hb=7da725cd0be7e880ff4971257720220d281c4ab3;hp=d282de2bcbc04c8f7f135277bc5a55c80e1736c4;hpb=d22f4f43182fe1cc70b343b6a52cb8e3a54d1f44;p=libc-test diff --git a/src/math/lgammal.c b/src/math/lgammal.c index d282de2..9050913 100644 --- a/src/math/lgammal.c +++ b/src/math/lgammal.c @@ -1,13 +1,17 @@ +#define _BSD_SOURCE 1 +#define _XOPEN_SOURCE 700 #include #include -#include "util.h" +#include "mtest.h" static struct l_li t[] = { #if LDBL_MANT_DIG == 53 #include "sanity/lgamma.h" +#include "special/lgamma.h" #elif LDBL_MANT_DIG == 64 #include "sanity/lgammal.h" +#include "special/lgammal.h" #endif }; @@ -39,7 +43,8 @@ int main(void) err++; } d = ulperrl(y, p->y, p->dy); - if (!checkulp(d, p->r) || 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++;