double_t
static double_t roundtoint(double_t);
static int32_t converttoint(double_t);
double_t hz,z,r,w;
double_t r, w, z;
fn = (double_t)x*invpio2 + toint - toint;
double_t z,w,t,r,fn;
double_t fn;
fn = (double_t)x*invpio2 + toint - toint;
double_t z,r,v,w;
double_t r, s, w, z;
double_t z, r, v, w, s, a;
double_t z,r,w,s,t,u;
double_t p, q;
double_t p, q;
double_t w,s1,s2,z;
double_t y;
double_t r,s,t,w;
double_t r,T;
T = T*((double_t)x+x+r)/(x+r+r);
T = T*((double_t)x+x+r)/(x+r+r);
double_t dr, dt, dx;
double_t y;
static const double_t toint = 1/EPS;
double_t s,P,Q;
double_t s,R,S;
kd = (double_t)(int32_t)ki;
static inline double specialcase(double_t tmp, uint64_t sbits, uint64_t ki)
double_t scale, y;
double_t hi, lo;
double_t kd, z, r, r2, scale, tail, tmp;
static inline double specialcase(double_t tmp, uint64_t sbits, uint64_t ki)
double_t scale, y;
double_t hi, lo;
double_t kd, r, r2, scale, tail, tmp;
double_t kd, xd, z, r, r2, y, s;
xd = (double_t)x;
double_t kd, xd, z, r, r2, y, s;
xd = (double_t)x;
double_t y,hi,lo,c,t,e,hxs,hfx,r1,twopk;
double_t y;
static const double_t toint = 1/EPS;
double_t tiny = DBL_MIN/FLT_MIN * r;
static void sq(double_t *hi, double_t *lo, double x)
double_t xh, xl, xc;
xc = (double_t)x*SPLIT;
*hi = (double_t)x*x;
double_t hx, lx, hy, ly, z;
double_t z,r,s;
double_t s,r,z;
double_t z,r,s;
double_t s,r,z;
double_t t,y,z,nadj,p,p1,p2,p3,q,r,w;
double_t y;
double_t w, z, r, r2, r3, y, invc, logc, kd, hi, lo;
double_t rhi = r + w - w;
double_t rlo = r - rhi;
kd = (double_t)k;
double_t hfsq,f,s,z,R,w,t1,t2,dk,y,hi,lo,val_hi,val_lo;
double_t hfsq,f,c,s,z,R,w,t1,t2,dk;
double_t z, r, r2, r4, y, invc, logc, kd, hi, lo, t1, t2, t3, p;
double_t rhi, rlo;
kd = (double_t)k;
double_t rhi, rlo;
double_t z, r, r2, p, y, y0, invc, logc;
z = (double_t)asfloat(iz);
y0 = logc + (double_t)k;
double_t z, r, r2, y, y0, invc, logc;
z = (double_t)asfloat(iz);
y0 = logc + (double_t)k * Ln2;
double_t toint = asdouble(asuint64(1/EPS) | sign);
double_t y = x + toint - toint;
static inline double specialcase(double_t tmp, uint64_t sbits, uint64_t ki)
double_t scale, y;
double_t hi, lo, one = 1.0;
static inline double exp_inline(double_t x, double_t xtail, uint32_t sign_bias)
double_t kd, z, r, r2, scale, tail, tmp;
double_t one = WANT_ROUNDING ? 1.0 + x : 1.0;
kd = (double_t)(int32_t)ki;
double_t x2 = x * x;
double_t lo;
double_t hi = log_inline(ix, &lo);
double_t ehi, elo;
double_t yhi = asdouble(iy & -1ULL << 27);
double_t ylo = y - yhi;
double_t lhi = asdouble(asuint64(hi) & -1ULL << 27);
double_t llo = hi - lhi + lo;
static inline double_t log_inline(uint64_t ix, double_t *tail)
double_t z, r, y, invc, logc, logctail, kd, hi, t1, t2, lo, lo1, lo2, p;
kd = (double_t)k;
double_t zhi = asdouble((iz + (1ULL << 31)) & (-1ULL << 32));
double_t zlo = z - zhi;
double_t rhi = zhi * invc - 1.0;
double_t rlo = zlo * invc;
double_t ar, ar2, ar3, lo3, lo4;
double_t arhi = A[0] * rhi;
double_t arhi2 = rhi * arhi;
double_t logx = log2_inline(ix);
double_t ylogx = y * logx; /* cannot overflow, y is single prec. */
static inline double_t log2_inline(uint32_t ix)
double_t z, r, r2, r4, p, q, y, y0, invc, logc;
z = (double_t)asfloat(iz);
y0 = logc + (double_t)k;
static inline float exp2_inline(double_t xd, uint32_t sign_bias)
double_t kd, z, r, r2, y, s;
static const double_t toint = 1/EPS;
double_t y;
double_t y;
static const double_t toint = 1/EPS;
double_t y = x;
double_t t;
double_t dy, z, r;
double_t num = 0, den = 0;