| /src/lib/libm/src/ |
| e_pow.c | 102 double z,ax,z_h,z_l,p_h,p_l; local in function:__ieee754_pow
247 p_h = u+v;
248 SET_LOW_WORD(p_h,0);
249 p_l = v-(p_h-u);
250 z_h = cp_h*p_h; /* cp_h+cp_l = 2/(3*log2) */
251 z_l = cp_l*p_h+p_l*cp+dp_l[k];
263 p_h = yy1*t1;
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 *
[all...] |
| e_powf.c | 61 float z,ax,z_h,z_l,p_h,p_l; local in function:__ieee754_powf
191 p_h = u+v;
192 GET_FLOAT_WORD(is,p_h);
193 SET_FLOAT_WORD(p_h,is&0xfffff000);
194 p_l = v-(p_h-u);
195 z_h = cp_h*p_h; /* cp_h+cp_l = 2/(3*log2) */
196 z_l = cp_l*p_h+p_l*cp+dp_l[k];
209 p_h = yy1*t1;
210 z = p_l+p_h;
215 if(p_l+ovt>z-p_h) return sn*huge*huge; /* overflow *
[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
346 p_h = u + v;
347 o.value = p_h;
350 p_h = o.value;
351 p_l = v - (p_h - u);
352 z_h = cp_h * p_h; /* cp_h+cp_l = 2/(3*log2) */
353 z_l = cp_l * p_h + p_l * cp + dp_l[k];
375 p_h = yy1 * t1;
376 z = p_l + p_h;
387 if (p_l + ovt > z - p_h)
[all...] |
