--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct d_di t[] = {
+#include "sanity/frexp.h"
+
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct d_di *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = frexp(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s frexp(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperr(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s frexp(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct f_fi t[] = {
+#include "sanity/frexpf.h"
+
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct f_fi *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = frexpf(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s frexpf(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrf(y, p->y, p->dy);
+ if (!checkulp(d, p->r)) {
+ printf("%s:%d: %s frexpf(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct l_li t[] = {
+#if LDBL_MANT_DIG == 53
+#include "sanity/frexp.h"
+
+#elif LDBL_MANT_DIG == 64
+#include "sanity/frexpl.h"
+
+#endif
+};
+
+int main(void)
+{
+ int yi;
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct l_li *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = frexpl(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s frexpl(%La)=%La,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrl(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s frexpl(%La) want %La,%lld got %La,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
T(scalb, dd_d)
T(scalbf, ff_f)
T(j0f, f_f)
-T(j0l, l_l)
T(j1f, f_f)
-T(j1l, l_l)
T(y0f, f_f)
-T(y0l, l_l)
T(y1f, f_f)
-T(y1l, l_l)
T(exp10, d_d)
T(exp10f, f_f)
T(exp10l, l_l)
T(pow10, d_d)
T(pow10f, f_f)
T(pow10l, l_l)
-
T(frexp, d_di)
T(frexpf, f_fi)
T(frexpl, l_li)
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct d_di t[] = {
+HEADERS
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct d_di *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ___(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ___(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperr(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s ___(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct di_d t[] = {
+HEADERS
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct di_d *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ___(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ___(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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 ___(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct f_fi t[] = {
+HEADERS
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct f_fi *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ___(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ___(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrf(y, p->y, p->dy);
+ if (!checkulp(d, p->r)) {
+ printf("%s:%d: %s ___(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct fi_f t[] = {
+HEADERS
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct fi_f *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ___(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ___(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s ___(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct l_li t[] = {
+#if LDBL_MANT_DIG == 53
+DHEADERS
+#elif LDBL_MANT_DIG == 64
+HEADERS
+#endif
+};
+
+int main(void)
+{
+ int yi;
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct l_li *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ___(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ___(%La)=%La,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrl(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s ___(%La) want %La,%lld got %La,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct li_l t[] = {
+#if LDBL_MANT_DIG == 53
+DHEADERS
+#elif LDBL_MANT_DIG == 64
+HEADERS
+#endif
+};
+
+int main(void)
+{
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct li_l *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ___(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ___(%La, %lld)=%La, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s ___(%La, %lld) want %La got %La ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct di_d t[] = {
+#include "sanity/ldexp.h"
+
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct di_d *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ldexp(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ldexp(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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 ldexp(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct fi_f t[] = {
+#include "sanity/ldexpf.h"
+
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct fi_f *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ldexpf(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ldexpf(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s ldexpf(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct li_l t[] = {
+#if LDBL_MANT_DIG == 53
+#include "sanity/ldexp.h"
+
+#elif LDBL_MANT_DIG == 64
+#include "sanity/ldexpl.h"
+
+#endif
+};
+
+int main(void)
+{
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct li_l *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = ldexpl(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s ldexpl(%La, %lld)=%La, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s ldexpl(%La, %lld) want %La got %La ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct d_di t[] = {
+#include "sanity/lgamma.h"
+
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct d_di *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = lgamma(p->x);
+ yi = signgam;
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lgamma(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperr(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s lgamma(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct d_di t[] = {
+#include "sanity/lgamma_r.h"
+
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct d_di *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = lgamma_r(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lgamma_r(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperr(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s lgamma_r(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct f_fi t[] = {
+#include "sanity/lgammaf.h"
+
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct f_fi *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = lgammaf(p->x);
+ yi = signgam;
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lgammaf(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrf(y, p->y, p->dy);
+ if (!checkulp(d, p->r)) {
+ printf("%s:%d: %s lgammaf(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct f_fi t[] = {
+#include "sanity/lgammaf_r.h"
+
+};
+
+int main(void)
+{
+ int yi;
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct f_fi *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = lgammaf_r(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lgammaf_r(%a)=%a,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrf(y, p->y, p->dy);
+ if (!checkulp(d, p->r)) {
+ printf("%s:%d: %s lgammaf_r(%a) want %a,%lld got %a,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct l_li t[] = {
+#if LDBL_MANT_DIG == 53
+#include "sanity/lgamma.h"
+
+#elif LDBL_MANT_DIG == 64
+#include "sanity/lgammal.h"
+
+#endif
+};
+
+int main(void)
+{
+ int yi;
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct l_li *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = lgammal(p->x);
+ yi = signgam;
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lgammal(%La)=%La,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrl(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s lgammal(%La) want %La,%lld got %La,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct l_li t[] = {
+#if LDBL_MANT_DIG == 53
+#include "sanity/lgammal_r.h"
+
+#elif LDBL_MANT_DIG == 64
+#include "sanity/lgammal_r.h"
+
+#endif
+};
+
+int main(void)
+{
+ int yi;
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct l_li *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = lgammal_r(p->x, &yi);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s lgammal_r(%La)=%La,%lld, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
+ printf(" got %s\n", estr(e));
+ err++;
+ }
+ d = ulperrl(y, p->y, p->dy);
+ if (!checkulp(d, p->r) || yi != p->i) {
+ printf("%s:%d: %s lgammal_r(%La) want %La,%lld got %La,%d ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct di_d t[] = {
+#include "sanity/scalbln.h"
+
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct di_d *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = scalbln(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s scalbln(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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 scalbln(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct fi_f t[] = {
+#include "sanity/scalblnf.h"
+
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct fi_f *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = scalblnf(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s scalblnf(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s scalblnf(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct li_l t[] = {
+#if LDBL_MANT_DIG == 53
+#include "sanity/scalbln.h"
+
+#elif LDBL_MANT_DIG == 64
+#include "sanity/scalblnl.h"
+
+#endif
+};
+
+int main(void)
+{
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct li_l *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = scalblnl(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s scalblnl(%La, %lld)=%La, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s scalblnl(%La, %lld) want %La got %La ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct di_d t[] = {
+#include "sanity/scalbn.h"
+
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct di_d *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = scalbn(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s scalbn(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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 scalbn(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct fi_f t[] = {
+#include "sanity/scalbnf.h"
+
+};
+
+int main(void)
+{
+ double y;
+ float d;
+ int e, i, err = 0;
+ struct fi_f *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = scalbnf(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s scalbnf(%a, %lld)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s scalbnf(%a, %lld) want %a got %a, ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
--- /dev/null
+#include <stdint.h>
+#include <stdio.h>
+#include "util.h"
+
+static struct li_l t[] = {
+#if LDBL_MANT_DIG == 53
+#include "sanity/scalbn.h"
+
+#elif LDBL_MANT_DIG == 64
+#include "sanity/scalbnl.h"
+
+#endif
+};
+
+int main(void)
+{
+ long double y;
+ float d;
+ int e, i, err = 0;
+ struct li_l *p;
+
+ for (i = 0; i < sizeof t/sizeof *t; i++) {
+ p = t + i;
+ setupfenv(p->r);
+ y = scalbnl(p->x, p->i);
+ e = getexcept();
+ if (!checkexcept(e, p->e, p->r)) {
+ printf("%s:%d: bad fp exception: %s scalbnl(%La, %lld)=%La, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->i, 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:%d: %s scalbnl(%La, %lld) want %La got %La ulperr %.3f = %a + %a\n",
+ p->file, p->line, rstr(p->r), p->x, p->i, p->y, y, d, d-p->dy, p->dy);
+ err++;
+ }
+ }
+ return !!err;
+}
#undef RZ
#undef RD
#undef RU
+#ifdef FE_TONEAREST
#define RN FE_TONEAREST
+#else
+#define RN 0
+#endif
+#ifdef FE_TOWARDZERO
#define RZ FE_TOWARDZERO
+#else
+#define RZ -1
+#endif
+#ifdef FE_DOWNWARD
#define RD FE_DOWNWARD
+#else
+#define RD -1
+#endif
+#ifdef FE_UPWARD
#define RU FE_UPWARD
+#else
+#define RU -1
+#endif
#undef INEXACT
#undef INVALID
#undef DIVBYZERO
#undef UNDERFLOW
#undef OVERFLOW
+#ifdef FE_INEXACT
#define INEXACT FE_INEXACT
+#else
+#define INEXACT 0
+#endif
+#ifdef FE_INVALID
#define INVALID FE_INVALID
+#else
+#define INVALID 0
+#endif
+#ifdef FE_DIVBYZERO
#define DIVBYZERO FE_DIVBYZERO
+#else
+#define DIVBYZERO 0
+#endif
+#ifdef FE_UNDERFLOW
#define UNDERFLOW FE_UNDERFLOW
+#else
+#define UNDERFLOW 0
+#endif
+#ifdef FE_OVERFLOW
#define OVERFLOW FE_OVERFLOW
+#else
+#define OVERFLOW 0
+#endif
#undef inf
#undef nan
struct ff_f {POS int r; float x; float x2; float y; float dy; int e; };
struct dd_d {POS int r; double x; double x2; double y; float dy; int e; };
struct ll_l {POS int r; long double x; long double x2; long double y; float dy; int e; };
+struct d_di {POS int r; double x; double y; float dy; long long i; int e; };
+struct f_fi {POS int r; float x; float y; float dy; long long i; int e; };
+struct l_li {POS int r; long double x; long double y; float dy; long long i; int e; };
+struct di_d {POS int r; double x; long long i; double y; float dy; int e; };
+struct fi_f {POS int r; float x; long long i; float y; float dy; int e; };
+struct li_l {POS int r; long double x; long long i; long double y; float dy; int e; };
+#undef POS
char *estr(int);
char *rstr(int);