diff options
Diffstat (limited to 'libc/sysdeps/ieee754/flt-32/e_powf.c')
| -rw-r--r-- | libc/sysdeps/ieee754/flt-32/e_powf.c | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/libc/sysdeps/ieee754/flt-32/e_powf.c b/libc/sysdeps/ieee754/flt-32/e_powf.c index 43069479a..12c408f93 100644 --- a/libc/sysdeps/ieee754/flt-32/e_powf.c +++ b/libc/sysdeps/ieee754/flt-32/e_powf.c @@ -141,6 +141,10 @@ __ieee754_powf(float x, float y) t2 = v-(t1-u); } else { float s2,s_h,s_l,t_h,t_l; + /* Avoid internal underflow for tiny y. The exact value + of y does not matter if |y| <= 2**-32. */ + if (iy < 0x2f800000) + SET_FLOAT_WORD (y, (hy & 0x80000000) | 0x2f800000); n = 0; /* take care subnormal number */ if(ix<0x00800000) |