X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Facos.c;h=e8e7743ee81087ebaeacb774d72a988cdd7a627d;hb=4282b03a8c212b0e38af1876364be5f332673ce6;hp=c38d7f2e7a991aa6395d106ec07c1bf0f5e6500f;hpb=0c7d465f8b04e9db61c2c9a32510799c19ec0503;p=libc-test diff --git a/src/math/acos.c b/src/math/acos.c index c38d7f2..e8e7743 100644 --- a/src/math/acos.c +++ b/src/math/acos.c @@ -1,12 +1,37 @@ -#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 acos(%a)==%a except: want %s", 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 acos(%a) want %a got %a ulperr %.3f = %a + %a\n", + rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy); + err++; + } } + return !!err; }