X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flgamma.c;h=d76538f16550fa67a67a2c8bb0a98fce7af07963;hb=21dd45c44f745d3b85afab304362d01b3b45a5e4;hp=301a967628eaf2d34d9ab8e397563a732427c01e;hpb=75aab93ef75299e347be24f57223dc93f2af80a4;p=libc-test diff --git a/src/math/lgamma.c b/src/math/lgamma.c index 301a967..d76538f 100644 --- a/src/math/lgamma.c +++ b/src/math/lgamma.c @@ -9,6 +9,7 @@ static struct d_di t[] = { int main(void) { + #pragma STDC FENV_ACCESS ON int yi; double y; float d; @@ -17,10 +18,15 @@ int main(void) for (i = 0; i < sizeof t/sizeof *t; i++) { p = t + i; - setupfenv(p->r); + + if (p->r < 0) + continue; + fesetround(p->r); + feclearexcept(FE_ALL_EXCEPT); y = lgamma(p->x); yi = signgam; - e = getexcept(); + e = fetestexcept(INEXACT|INVALID|DIVBYZERO|UNDERFLOW|OVERFLOW); + if (!checkexcept(e, p->e, p->r)) { printf("%s:%d: bad fp exception: %s lgamma(%a)=%a,%lld, want %s", p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));