AuroraMiddleware/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2024-11-22 21:10:07 +00:00
Code Issues Packages Projects Releases Wiki Activity

* sysdeps/m68k/fpu/e_pow.c: Correct handling of some exceptional

Browse Source 
values. 
* sysdeps/m68k/fpu/e_scalb.c: Likewise.
2001-04-10  Andreas Schwab  <schwab@suse.de>

	* sysdeps/m68k/fpu/e_pow.c: Correct handling of some exceptional
	values.
	* sysdeps/m68k/fpu/e_scalb.c: Likewise.
This commit is contained in:
Andreas Schwab 2001-04-10 19:58:52 +00:00
parent 9deee4bb2a
commit d2784f806c
3 changed files with 16 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ChangeLog
View File

@ -1,3 +1,9 @@
2001-04-10 Andreas Schwab <schwab@suse.de>
* sysdeps/m68k/fpu/e_pow.c: Correct handling of some exceptional
values.
* sysdeps/m68k/fpu/e_scalb.c: Likewise.
2001-04-09 Ulrich Drepper <drepper@redhat.com>
* sysdeps/sparc/sparc32/sparcv9/Dist: Add hp-timing.c.

14
sysdeps/m68k/fpu/e_pow.c
View File

@ -41,23 +41,25 @@ s(__ieee754_pow) (float_type x, float_type y)
y_cond = __m81_test (y);
if (y_cond & __M81_COND_ZERO)
return 1.0;
if (y_cond & __M81_COND_NAN)
return x == 1.0 ? x : x + y;
x_cond = __m81_test (x);
if ((x_cond | y_cond) & __M81_COND_NAN)
if (x_cond & __M81_COND_NAN)
return x + y;
if (y_cond & __M81_COND_INF)
{
ax = s(fabs) (x);
if (ax == 1)
return y - y;
if (ax > 1)
if (ax == 1.0)
return ax;
if (ax > 1.0)
return y_cond & __M81_COND_NEG ? 0 : y;
else
return y_cond & __M81_COND_NEG ? -y : 0;
}
if (s(fabs) (y) == 1)
if (s(fabs) (y) == 1.0)
return y_cond & __M81_COND_NEG ? 1 / x : x;
if (y == 2)
@ -77,7 +79,7 @@ s(__ieee754_pow) (float_type x, float_type y)
}
ax = s(fabs) (x);
if (x_cond & (__M81_COND_INF | __M81_COND_ZERO) || ax == 1)
if (x_cond & (__M81_COND_INF | __M81_COND_ZERO) || ax == 1.0)
{
z = ax;
if (y_cond & __M81_COND_NEG)

5
sysdeps/m68k/fpu/e_scalb.c
View File

@ -18,6 +18,7 @@
Boston, MA 02111-1307, USA. */
#include <math.h>
#include "math_private.h"
#include "mathimpl.h"
#ifndef SUFF
@ -47,14 +48,12 @@ s(__ieee754_scalb) (float_type x, float_type fn)
return x * fn;
else if (x_cond & __M81_COND_ZERO)
return x;
else if (x_cond & __M81_COND_INF)
return 0.0/0.0;
else
return x / -fn;
}
if (m81(__rint) (fn) != fn)
return 0.0/0.0;
return (x - x) / (x - x);
__asm ("fscale%.x %1, %0" : "=f" (retval) : "f" (fn), "0" (x));
return retval;