bde 2008-02-28 16:22:36 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c
Log:
Fix and improve some magic numbers for the medium size case.
e_rem_pio2.c:
This case goes up to about 2**20pi/2, but the comment about it said that
bde 2008-02-25 11:43:20 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c k_rem_pio2.c
Log:
Fix some off-by-1 errors.
e_rem_pio2.c:
Float and double precision didn't work because init_jk[] was 1 too small.
It needs to be 2
bde 2008-02-25 18:28:58 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c
Log:
Use a temporary array instead of the arg array y[] for calling
__kernel_rem_pio2(). This simplifies analysis of aliasing and thus
results in better
bde 2008-02-23 12:53:21 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c
Log:
Optimize the 9pi/2 |x| = 2**19pi/2 case some more by avoiding an
fabs(), a conditional branch, and sign adjustments of 3 variables for
x 0 when the
bde 2008-02-22 15:55:15 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c
Log:
Optimize the 9pi/2 |x| = 2**19pi/2 case on amd64 and i386 by avoiding
the the double to int conversion operation which is very slow on these
arches.
bde 2008-02-22 18:43:23 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c
Log:
Avoid using FP-to-integer conversion for !(amd64 || i386) too. Use the
FP-to-FP method to round to an integer on all arches, and convert this
to an
das 2008-02-17 07:31:59 UTC
FreeBSD src repository
Modified files:
lib/msun/src e_rem_pio2.c e_rem_pio2f.c k_rem_pio2.c
math_private.h
Log:
Add more pi for long doubles. Also, avoid storing multiple copies
of the pi/2 array, as it is