| /src/lib/libm/src/ |
| e_pow.c | 102 double z,ax,z_h,z_l,p_h,p_l; local in function:__ieee754_pow
249 p_l = v-(p_h-u);
251 z_l = cp_l*p_h+p_l*cp+dp_l[k];
262 p_l = (y-yy1)*t1+y*t2;
264 z = p_l+p_h;
270 if(p_l+ovt>z-p_h) return s*huge*huge; /* overflow */
276 if(p_l<=z-p_h) return s*tiny*tiny; /* underflow */
280 * compute 2**(p_h+p_l)
294 t = p_l+p_h;
297 v = (p_l-(t-p_h))*lg2+t*lg2_l
[all...] |
| e_powf.c | 61 float z,ax,z_h,z_l,p_h,p_l; local in function:__ieee754_powf
194 p_l = v-(p_h-u);
196 z_l = cp_l*p_h+p_l*cp+dp_l[k];
208 p_l = (y-yy1)*t1+y*t2;
210 z = p_l+p_h;
215 if(p_l+ovt>z-p_h) return sn*huge*huge; /* overflow */
220 if(p_l<=z-p_h) return sn*tiny*tiny; /* underflow */
223 * compute 2**(p_h+p_l)
236 t = p_l+p_h;
240 v = (p_l-(t-p_h))*lg2+t*lg2_l
[all...] |
| /src/lib/libm/ld128/ |
| e_powl.c | 148 long double z, ax, z_h, z_l, p_h, p_l; local in function:powl
351 p_l = v - (p_h - u);
353 z_l = cp_l * p_h + p_l * cp + dp_l[k];
374 p_l = (y - yy1) * t1 + y * t2;
376 z = p_l + p_h;
387 if (p_l + ovt > z - p_h)
400 if (p_l <= z - p_h)
404 /* compute 2**(p_h+p_l) */
414 t = p_l + p_h;
420 v = (p_l - (t - p_h)) * lg2 + t * lg2_l
[all...] |
