hy &= 0x7fffffff; /* |y| */
/* purge off exception values */
- // FIXME: signed shift
if ((hy|ly) == 0 || hx >= 0x7ff00000 || /* y = 0, or x not finite */
(hy|((ly|-ly)>>31)) > 0x7ff00000) /* or y is NaN */
return (x*y)/(x*y);
goto fixup;
}
if (lx == ly) { /* |x| = |y| return x*0 */
- *quo = 1;
+ *quo = sxy ? -1 : 1;
return Zero[(uint32_t)sx>>31];
}
}
/* convert back to floating value and restore the sign */
if ((hx|lx) == 0) { /* return sign(x)*0 */
+ q &= 0x7fffffff;
*quo = sxy ? -q : q;
return Zero[(uint32_t)sx>>31];
}
hx >>= n;
} else if (n <= 31) {
lx = (hx<<(32-n))|(lx>>n);
- hx = sx;
+ hx = 0;
} else {
lx = hx>>(n-32);
- hx = sx;
+ hx = 0;
}
}
fixup: