math: fix remquo.c when x==-y and a subnormal remainder bug as well

backported fix from freebsd:
http://svnweb.FreeBSD.org/base?view=revision&revision=233973

1 files changed, 2 insertions, 1 deletions

diff --git a/src/math/remquof.c b/src/math/remquof.c
index 11569ce8..536a050a 100644
--- a/src/math/remquof.c
+++ b/src/math/remquof.c
@@ -41,7 +41,7 @@ float remquof(float x, float y, int *quo)
 		q = 0;
 		goto fixup;
 	} else if(hx==hy) {  /* |x| = |y| return x*0*/
-		*quo = 1;
+		*quo = sxy ? -1 : 1;
 		return Zero[(uint32_t)sx>>31];
 	}
 
@@ -92,6 +92,7 @@ float remquof(float x, float y, int *quo)
 
 	/* convert back to floating value and restore the sign */
 	if (hx == 0) {                             /* return sign(x)*0 */
+		q &= 0x7fffffff;
 		*quo = sxy ? -q : q;
 		return Zero[(uint32_t)sx>>31];
 	}