AuroraMiddleware/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2025-01-10 19:30:10 +00:00
Code Issues Packages Projects Releases Wiki Activity

Fix uses of math_force_eval

Browse Source
This commit is contained in:
Andreas Schwab 2011-10-26 11:37:08 +02:00
parent f6ce9294b0
commit 804791474d
4 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
ChangeLog
View File

@ -1,5 +1,9 @@
2011-10-26 Andreas Schwab <schwab@redhat.com> 2011-10-26 Andreas Schwab <schwab@redhat.com>
* sysdeps/ieee754/flt-32/e_j0f.c: Fix use of math_force_eval.
* sysdeps/ieee754/dbl-64/s_round.c: Likewise.
* sysdeps/ieee754/flt-32/s_roundf.c: Likewise.
* math/math_private.h (math_force_eval): Allow non-addressable * math/math_private.h (math_force_eval): Allow non-addressable
arguments. arguments.
* sysdeps/i386/fpu/math_private.h (math_force_eval): Likewise. * sysdeps/i386/fpu/math_private.h (math_force_eval): Likewise.

6
sysdeps/ieee754/dbl-64/s_round.c
View File

@ -38,7 +38,7 @@ __round (double x)
{ {
if (j0 < 0) if (j0 < 0)
{ {
math_force_eval (huge + x > 0.0); math_force_eval (huge + x);
i0 &= 0x80000000; i0 &= 0x80000000;
if (j0 == -1) if (j0 == -1)
@ -51,7 +51,7 @@ __round (double x)
if (((i0 & i) | i1) == 0) if (((i0 & i) | i1) == 0)
/* X is integral. */ /* X is integral. */
return x; return x;
math_force_eval (huge + x > 0.0); math_force_eval (huge + x);
/* Raise inexact if x != 0. */ /* Raise inexact if x != 0. */
i0 += 0x00080000 >> j0; i0 += 0x00080000 >> j0;
@ -74,7 +74,7 @@ __round (double x)
/* X is integral. */ /* X is integral. */
return x; return x;
math_force_eval (huge + x > 0.0); math_force_eval (huge + x);
/* Raise inexact if x != 0. */ /* Raise inexact if x != 0. */
u_int32_t j = i1 + (1 << (51 - j0)); u_int32_t j = i1 + (1 << (51 - j0));

2
sysdeps/ieee754/flt-32/e_j0f.c
View File

@ -66,7 +66,7 @@ __ieee754_j0f(float x)
return z; return z;
} }
if(ix<0x39000000) { /* |x| < 2**-13 */ if(ix<0x39000000) { /* |x| < 2**-13 */
math_force_eval(huge+x>one); /* raise inexact if x != 0 */ math_force_eval(huge+x); /* raise inexact if x != 0 */
if(ix<0x32000000) return one; /* |x|<2**-27 */ if(ix<0x32000000) return one; /* |x|<2**-27 */
else return one - (float)0.25*x*x; else return one - (float)0.25*x*x;
} }

4
sysdeps/ieee754/flt-32/s_roundf.c
View File

@ -37,7 +37,7 @@ __roundf (float x)
{ {
if (j0 < 0) if (j0 < 0)
{ {
math_force_eval (huge + x > 0.0F); math_force_eval (huge + x);
i0 &= 0x80000000; i0 &= 0x80000000;
if (j0 == -1) if (j0 == -1)
@ -49,7 +49,7 @@ __roundf (float x)
if ((i0 & i) == 0) if ((i0 & i) == 0)
/* X is integral. */ /* X is integral. */
return x; return x;
math_force_eval (huge + x > 0.0F); math_force_eval (huge + x);
/* Raise inexact if x != 0. */ /* Raise inexact if x != 0. */
i0 += 0x00400000 >> j0; i0 += 0x00400000 >> j0;