+2012-11-17 David S. Miller <davem@davemloft.net>
+
+ [BZ #14811]
+ * sysdeps/ieee754/ldbl-128/e_powl.c (__ieee754_powl): Saturate
+ nonzero exponents with absolute value below 0x1p-128 to +/-
+ 0x1p-128.
+
2012-11-17 Joseph Myers <joseph@codesourcery.com>
* sysdeps/unix/sysv/linux/fxstatat.c: Include <string.h>.
{
long double z, ax, z_h, z_l, p_h, p_l;
long double y1, t1, t2, r, s, t, u, v, w;
- long double s2, s_h, s_l, t_h, t_l;
+ long double s2, s_h, s_l, t_h, t_l, ay;
int32_t i, j, k, yisint, n;
u_int32_t ix, iy;
int32_t hx, hy;
return (hy > 0) ? huge * huge : tiny * tiny;
}
+ ay = y > 0 ? y : -y;
+ if (ay < 0x1p-128)
+ y = y < 0 ? -0x1p-128 : 0x1p-128;
+
n = 0;
/* take care subnormal number */
if (ix < 0x00010000)