diff options
| author | Adhemerval Zanella <azanella@linux.vnet.ibm.com> | 2012-04-11 16:30:13 -0300 |
|---|---|---|
| committer | Andreas Jaeger <aj@suse.de> | 2012-04-17 22:12:53 +0200 |
| commit | 76da7265320010c7a273ed99f53938c0f32d5fad (patch) | |
| tree | f094871a5f38652d2d320a11423fee602a8746e6 /sysdeps/generic | |
| parent | 0396e69dbefb3630651d93e0a7526c4326a071bb (diff) | |
| download | glibc-76da7265320010c7a273ed99f53938c0f32d5fad.zip glibc-76da7265320010c7a273ed99f53938c0f32d5fad.tar.gz glibc-76da7265320010c7a273ed99f53938c0f32d5fad.tar.bz2 | |
Fix ilogb exception and errno (bug 6794)
[BZ #6794]
Following Joseph comments about bug 6794, here is a proposed fix. It turned out
to be a large fix mainly because I had to move some file along to follow libm
files/names conventions.
Basically I have added wrappers (w_ilogb.c, w_ilogbf.c, w_ilogbl.c) that now calls
the symbol '__ieee754_ilogb'. The wrappers checks for '__ieee754_ilogb' output and
set the errno and raise exceptions as expected.
The '__ieee754_ilogb' is implemented in sysdeps. I have moved the 's_ilogb[f|l]' files
to e_ilogb[f|l] and renamed the '__ilogb[f|l]' to '__ieee754_ilogb[f|l]'.
I also found out a bug in i386 and x86-64 assembly coded ilogb implementation where
it raises a FE_DIVBYZERO when argument is '0.0'. I corrected this issue as well.
Finally I added the errno and FE_INVALID tests for 0.0, NaN and +-InF argument. Tested
on i386, x86-64, ppc32 and ppc64.
Diffstat (limited to 'sysdeps/generic')
| -rw-r--r-- | sysdeps/generic/math_private.h | 3 |
1 files changed, 3 insertions, 0 deletions
diff --git a/sysdeps/generic/math_private.h b/sysdeps/generic/math_private.h index e217224..5267ec3 100644 --- a/sysdeps/generic/math_private.h +++ b/sysdeps/generic/math_private.h @@ -213,6 +213,7 @@ extern double __ieee754_yn (int,double); extern double __ieee754_remainder (double,double); extern int32_t __ieee754_rem_pio2 (double,double*); extern double __ieee754_scalb (double,double); +extern int __ieee754_ilogb (double); /* fdlibm kernel function */ extern double __kernel_standard (double,double,int); @@ -260,6 +261,7 @@ extern float __ieee754_ynf (int,float); extern float __ieee754_remainderf (float,float); extern int32_t __ieee754_rem_pio2f (float,float*); extern float __ieee754_scalbf (float,float); +extern int __ieee754_ilogbf (float); /* float versions of fdlibm kernel functions */ @@ -305,6 +307,7 @@ extern long double __ieee754_ynl (int,long double); extern long double __ieee754_remainderl (long double,long double); extern int __ieee754_rem_pio2l (long double,long double*); extern long double __ieee754_scalbl (long double,long double); +extern int __ieee754_ilogbl (long double); /* long double versions of fdlibm kernel functions */ extern long double __kernel_sinl (long double,long double,int); |