C1
C1 = 4.16666666666666019037e-02, /* 0x3FA55555, 0x5555554C */
r = z*(C1+z*(C2+z*C3)) + w*w*(C4+z*(C5+z*C6));
C1 = 0x155553e1053a42.0p-57, /* 0.0416666233237390631894 */
return ((1.0+z*C0) + w*C1) + (w*z)*r;
C1 = 0.0416666666666666666136L; /* 0xaaaaaaaaaaaaaa9b.0p-68 */
#define POLY(z) (z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*C7)))))))
C1 = 0.04166666666666666666666666666666658424671L,
#define POLY(z) (z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*(C7+ \
tmp = tail + r * C1 + r2 * (C2 + r * C3) + r2 * r2 * (C4 + r * C5);
x -= px * C1;
C1 = 6.93145751953125E-1L,
xx = C1 + C2;
z = z + e * C1;
z = z + e * C1;
static const long double C1 = 6.9314575195312500000000E-1L;
static const long double C1 = 6.9314575195312500000000E-1L;
z = z + e * C1;
z = z + e * C1; /* This sum has an error of 1/2 lsb. */