RepoMirrors/musl
1
0
Fork 0
mirror of git://git.musl-libc.org/musl synced 2025-04-01 22:48:38 +00:00
Code Issues Packages Projects Releases Wiki Activity

Merge remote-tracking branch 'nsz/master'

Browse Source
This commit is contained in:
Rich Felker 2012-05-10 22:25:14 -04:00
commit 9cfa3065c5
3 changed files with 8 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/math/remquo.c
View File

@ -44,7 +44,7 @@ double remquo(double x, double y, int *quo)
goto fixup; goto fixup;
} }
if (lx == ly) { /* |x| = |y| return x*0 */ if (lx == ly) { /* |x| = |y| return x*0 */
*quo = 1; *quo = sxy ? -1 : 1;
return Zero[(uint32_t)sx>>31]; return Zero[(uint32_t)sx>>31];
} }
} }
@ -127,6 +127,7 @@ double remquo(double x, double y, int *quo)
/* convert back to floating value and restore the sign */ /* convert back to floating value and restore the sign */
if ((hx|lx) == 0) { /* return sign(x)*0 */ if ((hx|lx) == 0) { /* return sign(x)*0 */
q &= 0x7fffffff;
*quo = sxy ? -q : q; *quo = sxy ? -q : q;
return Zero[(uint32_t)sx>>31]; return Zero[(uint32_t)sx>>31];
} }
@ -144,10 +145,10 @@ double remquo(double x, double y, int *quo)
hx >>= n; hx >>= n;
} else if (n <= 31) { } else if (n <= 31) {
lx = (hx<<(32-n))|(lx>>n); lx = (hx<<(32-n))|(lx>>n);
hx = sx; hx = 0;
} else { } else {
lx = hx>>(n-32); lx = hx>>(n-32);
hx = sx; hx = 0;
} }
} }
fixup: fixup:

3
src/math/remquof.c
View File

@ -41,7 +41,7 @@ float remquof(float x, float y, int *quo)
q = 0; q = 0;
goto fixup; goto fixup;
} else if(hx==hy) { /* |x| = |y| return x*0*/ } else if(hx==hy) { /* |x| = |y| return x*0*/
*quo = 1; *quo = sxy ? -1 : 1;
return Zero[(uint32_t)sx>>31]; 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 */ /* convert back to floating value and restore the sign */
if (hx == 0) { /* return sign(x)*0 */ if (hx == 0) { /* return sign(x)*0 */
q &= 0x7fffffff;
*quo = sxy ? -q : q; *quo = sxy ? -q : q;
return Zero[(uint32_t)sx>>31]; return Zero[(uint32_t)sx>>31];
} }

3
src/math/remquol.c
View File

@ -94,7 +94,7 @@ long double remquol(long double x, long double y, int *quo)
goto fixup; /* |x|<|y| return x or x-y */ goto fixup; /* |x|<|y| return x or x-y */
} }
if (ux.bits.manh == uy.bits.manh && ux.bits.manl == uy.bits.manl) { if (ux.bits.manh == uy.bits.manh && ux.bits.manl == uy.bits.manl) {
*quo = 1; *quo = sxy ? -1 : 1;
return Zero[sx]; /* |x|=|y| return x*0*/ return Zero[sx]; /* |x|=|y| return x*0*/
} }
} }
@ -152,6 +152,7 @@ long double remquol(long double x, long double y, int *quo)
/* convert back to floating value and restore the sign */ /* convert back to floating value and restore the sign */
if ((hx|lx) == 0) { /* return sign(x)*0 */ if ((hx|lx) == 0) { /* return sign(x)*0 */
q &= 0x7fffffff;
*quo = sxy ? -q : q; *quo = sxy ? -q : q;
return Zero[sx]; return Zero[sx];
} }