| Age | Commit message (Collapse) | Author |
|---|
|
|
|
[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.
|
|
|
|
We can reduce the number of STMXCSR, and often we can avoid the
call to __feraiseexcept.
|
|
|
|
|
|
For code generation to stay identical on x86_64, this requires that
we define the fp word manipulation macros before including the
generic header.
|
|
so that double arithmetic in s_sin is done in 53 bit (without extend i386 double precision)
|
|
This reverts commit 60d6f5a6f50d838bcb4240fcc0223cac445c6c83.
|
