X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fsinhl.c;fp=src%2Fmath%2Fsinhl.c;h=d52adb311ad8851e829956d840066bc96a59022c;hb=f9d17902a35b6403b7c8354845e9f13f882c1c8e;hp=0000000000000000000000000000000000000000;hpb=1ee02029133b60651176619b0a6743a5d011c21b;p=libc-test diff --git a/src/math/sinhl.c b/src/math/sinhl.c new file mode 100644 index 0000000..d52adb3 --- /dev/null +++ b/src/math/sinhl.c @@ -0,0 +1,40 @@ +#include +#include +#include "util.h" + +static struct l_l t[] = { +#if LDBL_MANT_DIG == 53 +D +#include "sanity/sinhl.h" +#elif LDBL_MANT_DIG == 64 + +#include "sanity/sinhl.h" +#endif +}; + +int main(void) +{ + long double y; + float d; + int e, i, err = 0; + struct l_l *p; + + for (i = 0; i < sizeof t/sizeof *t; i++) { + p = t + i; + setupfenv(p->r); + y = sinhl(p->x); + e = getexcept(); + if (!checkexcept(e, p->e, p->r)) { + printf("%s sinhl(%La)==%La except: want %s", rstr(p->r), p->x, p->y, estr(p->e)); + printf(" got %s\n", estr(e)); + err++; + } + d = ulperrl(y, p->y, p->dy); + if (!checkulp(d, p->r)) { + printf("%s sinhl(%La) want %La got %La ulperr %.3f = %a + %a\n", + rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy); + err++; + } + } + return !!err; +}