X-Git-Url: http://nsz.repo.hu/git/?p=libc-test;a=blobdiff_plain;f=src%2Fmath%2Fscalbf.c;fp=src%2Fmath%2Fscalbf.c;h=249032582aa2690d8a7409667d1331301d308218;hp=0000000000000000000000000000000000000000;hb=f9d17902a35b6403b7c8354845e9f13f882c1c8e;hpb=1ee02029133b60651176619b0a6743a5d011c21b diff --git a/src/math/scalbf.c b/src/math/scalbf.c new file mode 100644 index 0000000..2490325 --- /dev/null +++ b/src/math/scalbf.c @@ -0,0 +1,34 @@ +#include +#include +#include "util.h" + +static struct ff_f t[] = { + +}; + +int main(void) +{ + float y; + float d; + int e, i, err = 0; + struct ff_f *p; + + for (i = 0; i < sizeof t/sizeof *t; i++) { + p = t + i; + setupfenv(p->r); + y = scalbf(p->x, p->x2); + e = getexcept(); + if (!checkexcept(e, p->e, p->r)) { + printf("%s scalbf(%a,%a)==%a except: want %s", rstr(p->r), p->x, p->x2, p->y, estr(p->e)); + printf(" got %s\n", estr(e)); + err++; + } + d = ulperrf(y, p->y, p->dy); + if (!checkulp(d, p->r)) { + printf("%s scalbf(%a,%a) want %a got %a ulperr %.3f = %a + %a\n", + rstr(p->r), p->x, p->x2, p->y, y, d, d-p->dy, p->dy); + err++; + } + } + return !!err; +}