X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flgammal_r.c;h=848f48d464954303a24923ca6db72fbd3a968a24;hb=6abdab8bf1ccda7f28916222043bb498a1b9adf3;hp=03e3901e30b58d995f45cf54be47ffb403ab3c60;hpb=c886bb85233e642ca37731de72ad3e1f2f18a2f3;p=libc-test diff --git a/src/math/lgammal_r.c b/src/math/lgammal_r.c index 03e3901..848f48d 100644 --- a/src/math/lgammal_r.c +++ b/src/math/lgammal_r.c @@ -1,6 +1,7 @@ +#define _GNU_SOURCE 1 #include #include -#include "util.h" +#include "mtest.h" static struct l_li t[] = { #if LDBL_MANT_DIG == 53 @@ -40,7 +41,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_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++;