From 76da7265320010c7a273ed99f53938c0f32d5fad Mon Sep 17 00:00:00 2001 From: Adhemerval Zanella Date: Wed, 11 Apr 2012 16:30:13 -0300 Subject: 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. --- sysdeps/generic/math_private.h | 3 +++ 1 file changed, 3 insertions(+) (limited to 'sysdeps/generic') diff --git a/sysdeps/generic/math_private.h b/sysdeps/generic/math_private.h index e2172246f2..5267ec327b 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); -- cgit v1.2.3