X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Facos.c;h=f9d28f07a2edf017f710028a567879085b91ecfb;hb=9b2923eb8edf5d0872c9fdd3561f00dbf8480840;hp=c38d7f2e7a991aa6395d106ec07c1bf0f5e6500f;hpb=0c7d465f8b04e9db61c2c9a32510799c19ec0503;p=libc-test diff --git a/src/math/acos.c b/src/math/acos.c index c38d7f2..f9d28f0 100644 --- a/src/math/acos.c +++ b/src/math/acos.c @@ -1,12 +1,38 @@ -#include "test.h" -#include +#include +#include +#include "util.h" -void bench_acos(int N) +static struct d_d t[] = { +#include "crlibm/acos.h" +#include "ucb/acos.h" +#include "sanity/acos.h" + +}; + +int main(void) { - int i; - volatile double y; + double y; + float d; + int e, i, err = 0; + struct d_d *p; - for (i = 0; i < N; i++) { - y = acos(0.3456); + for (i = 0; i < sizeof t/sizeof *t; i++) { + p = t + i; + setupfenv(p->r); + y = acos(p->x); + e = getexcept(); + if (!checkexcept(e, p->e, p->r)) { + printf("%s:%d: bad fp exception: %s acos(%a)=%a, want %s", + p->file, p->line, rstr(p->r), p->x, p->y, estr(p->e)); + printf(" got %s\n", estr(e)); + err++; + } + d = ulperr(y, p->y, p->dy); + if (!checkulp(d, p->r)) { + printf("%s:%d: %s acos(%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++; + } } + return !!err; }