Lines Matching refs:dbl
231 ax.dbl = fabs (x);
237 z.dbl = ax.dbl; /* x is +-0,+-inf,+-1 */
240 z.dbl = one / z.dbl; /* z = (1 / |x|) */
246 z.dbl = NAN; /* (-1)**non-int is NaN */
250 z.dbl = -z.dbl; /* (x<0)**odd = -(|x|**odd) */
253 return z.dbl;
296 t.dbl = ax.dbl - one; /* t has 20 trailing zeros */
297 w = (t.dbl * t.dbl) * (0.5 - t.dbl * (0.3333333333333333333333 - t.dbl * 0.25));
298 u = ivln2_h * t.dbl; /* ivln2_h has 21 sig. bits */
299 v = t.dbl * ivln2_l - w * ivln2;
300 t1.dbl = u + v;
302 t2 = v - (t1.dbl - u);
313 ax.dbl *= two53;
338 u = ax.dbl - bp[k]; /* bp[0] = 1.0, bp[1] = 1.5 */
339 v = one / (ax.dbl + bp[k]);
341 s_h.dbl = ss;
344 t_h.dbl = zero;
346 t_l = ax.dbl - (t_h.dbl - bp[k]);
347 s_l = v * ((u - s_h.dbl * t_h.dbl) - s_h.dbl * t_l);
351 r += s_l * (s_h.dbl + ss);
352 s2 = s_h.dbl * s_h.dbl;
353 t_h.dbl = 3.0 + s2 + r;
355 t_l = r - ((t_h.dbl - 3.0) - s2);
357 u = s_h.dbl * t_h.dbl;
358 v = s_l * t_h.dbl + t_l * ss;
360 p_h.dbl = u + v;
362 p_l = v - (p_h.dbl - u);
363 z_h = cp_h * p_h.dbl; /* cp_h + cp_l = 2 / (3 * log2) */
364 z_l = cp_l * p_h.dbl + p_l * cp + dp_l[k];
366 t.dbl = (double) n;
367 t1.dbl = (((z_h + z_l) + dp_h[k]) + t.dbl);
369 t2 = z_l - (((t1.dbl - t.dbl) - dp_h[k]) - z_h);
373 y1.dbl = y;
375 p_l = (y - y1.dbl) * t1.dbl + y * t2;
376 p_h.dbl = y1.dbl * t1.dbl;
377 z.dbl = p_l + p_h.dbl;
388 if (p_l + ovt > z.dbl - p_h.dbl)
402 if (p_l <= z.dbl - p_h.dbl)
418 t.dbl = zero;
425 p_h.dbl -= t.dbl;
427 t.dbl = p_l + p_h.dbl;
429 u = t.dbl * lg2_h;
430 v = (p_l - (t.dbl - p_h.dbl)) * lg2 + t.dbl * lg2_l;
431 z.dbl = u + v;
432 w = v - (z.dbl - u);
433 t.dbl = z.dbl * z.dbl;
434 t1.dbl = z.dbl - t.dbl * (P1 + t.dbl * (P2 + t.dbl * (P3 + t.dbl * (P4 + t.dbl * P5))));
435 r = (z.dbl * t1.dbl) / (t1.dbl - two) - (w + z.dbl * w);
436 z.dbl = one - (r - z.dbl);
441 z.dbl = scalbn (z.dbl, n);
447 return s * z.dbl;