scalb is buggy, do we need the *f and *l version?
scalbn:
needs a rewrite: signed int problems, wrong <-50000 check
-__ldexp_cexp:
- only complex in math/, move to cmath?
j0l, y0l,..:
declare long double bessel functions in math.h?
uniform inexact exception raising:
long double complex dc = sin(1.0+0*I);
long double complex fc = sin(1.0f+0*I);
- long double complex lc = sin(1.0l+0*I);
+// long double complex lc = sin(1.0l+0*I);
printf("i %2u %La\n", sizeof sin(1), i);
printf("d %2u %La\n", sizeof sin(1.0), d);
printf("dc %2u %La %La\n", sizeof sin(1.0+0*I), creal(dc), cimag(dc));
printf("fc %2u %La %La\n", sizeof sin(1.0f+0*I), creal(fc), cimag(fc));
- printf("lc %2u %La %La\n", sizeof sin(1.0l+0*I), creal(lc), cimag(lc));
+// printf("lc %2u %La %La\n", sizeof sin(1.0l+0*I), creal(lc), cimag(lc));
printf("sizeof pow(I,1) = %2u\n", sizeof pow(I,1));
printf("sizeof pow(I,1.0f) = %2u\n", sizeof pow(I,1.0f));