Fix a missing mask in a hi+lo decomposition. Thus bug made the extra

precision in software useless, so hypotf() had some errors in the 1-2
ulp range unless there is extra precision in hardware (as happens on
i386).

1 files changed, 1 insertions, 1 deletions

diff --git a/lib/msun/src/e_hypotf.c b/lib/msun/src/e_hypotf.c
index b100b75..b760cc6 100644
--- a/lib/msun/src/e_hypotf.c
+++ b/lib/msun/src/e_hypotf.c
@@ -72,7 +72,7 @@ __ieee754_hypotf(float x, float y)
 	    a  = a+a;
 	    SET_FLOAT_WORD(y1,hb&0xfffff000);
 	    y2 = b - y1;
-	    SET_FLOAT_WORD(t1,ha+0x00800000);
+	    SET_FLOAT_WORD(t1,(ha+0x00800000)&0xfffff000);
 	    t2 = a - t1;
 	    w  = __ieee754_sqrtf(t1*y1-(w*(-w)-(t1*y2+t2*b)));
 	}