X-Git-Url: http://nsz.repo.hu/git/?p=libc-test;a=blobdiff_plain;f=src%2Fmath%2Fy0.c;h=b5b97b851a9908c19f63099d7c956e5f494bc4ff;hp=2ea45f74fb6e8126008028b0a3be914f6c64a3db;hb=5fc943786cdd4e93efe31f9f742f2097cec9f8a1;hpb=75543c20a003e3108b51cb63998752d27621f993 diff --git a/src/math/y0.c b/src/math/y0.c index 2ea45f7..b5b97b8 100644 --- a/src/math/y0.c +++ b/src/math/y0.c @@ -14,7 +14,7 @@ int main(void) #pragma STDC FENV_ACCESS ON double y; float d; - int e, i, err = 0; + int e, i, bad, err = 0; struct d_d *p; for (i = 0; i < sizeof t/sizeof *t; i++) { @@ -34,13 +34,14 @@ int main(void) err++; } d = ulperr(y, p->y, p->dy); - if ((!(p->x < 0) && !checkulp(d, p->r)) || (p->x < 0 && !isnan(y) && y != -inf)) { -// printf("%s:%d: %s y0(%a) want %a got %a ulperr %.3f = %a + %a\n", -// p->file, p->line, rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy); - err++; - // TODO: avoid spamming the output - printf(__FILE__ ": known to be broken near zeros\n"); - break; + bad = p->x < 0 && !isnan(y) && y != -inf; + if (bad || (!(p->x < 0) && !checkulp(d, p->r))) { + if (!bad && fabsf(d) < 0x1p52f) + printf("X "); + else + err++; + printf("%s:%d: %s y0(%a) want %a got %a ulperr %.3f = %a + %a\n", + p->file, p->line, rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy); } } return !!err;