musl/src/math, branch v0.9.0 musl - an implementation of the standard library for Linux-based systems add FORCE_EVAL macro to evaluate float expr for their side effect 2012-05-06T19:24:28+00:00 nsz nsz@port70.net 2012-05-06T19:24:28+00:00 6ab8136b4477fd75381c06fa0e7fa93c89c712a1 updated nextafter* to use FORCE_EVAL, it can be used in many other places in the math code to improve readability. 
updated nextafter* to use FORCE_EVAL, it can be used in many other
places in the math code to improve readability.
fix unused variable warnings in new nextafter/nexttoward code 2012-05-06T18:48:20+00:00 Rich Felker dalias@aerifal.cx 2012-05-06T18:48:20+00:00 4e597feef0595caefa39ab43c813734a1244fa84 apparently initializing a variable is not "using" it but assigning to it is "using" it. i don't really like this fix, but it's better than trying to make a bigger cleanup just before a release, and it should work fine (tested against nsz's math tests). 
apparently initializing a variable is not "using" it but assigning to
it is "using" it. i don't really like this fix, but it's better than
trying to make a bigger cleanup just before a release, and it should
work fine (tested against nsz's math tests).
math: nextafter and nexttoward cleanup 2012-05-06T11:08:59+00:00 nsz nsz@port70.net 2012-05-06T11:08:59+00:00 6cf865dba69bab6346dc268d9173609af36b984e make nexttoward, nexttowardf independent of long double representation. fix nextafterl: it did not raise underflow flag when the result was 0. 
make nexttoward, nexttowardf independent of long double representation.
fix nextafterl: it did not raise underflow flag when the result was 0.
math: change the formula used for acos.s 2012-05-04T23:11:56+00:00 nsz nsz@port70.net 2012-05-04T23:11:56+00:00 f697d66b81912af59128ac1b96bc0e2a4514b758 old: 2*atan2(sqrt(1-x),sqrt(1+x)) new: atan2(fabs(sqrt((1-x)*(1+x))),x) improvements: * all edge cases are fixed (sign of zero in downward rounding) * a bit faster (here a single call is about 131ns vs 162ns) * a bit more precise (at most 1ulp error on 1M uniform random samples in [0,1), the old formula gave some 2ulp errors as well) 
old: 2*atan2(sqrt(1-x),sqrt(1+x))
new: atan2(fabs(sqrt((1-x)*(1+x))),x)
improvements:
* all edge cases are fixed (sign of zero in downward rounding)
* a bit faster (here a single call is about 131ns vs 162ns)
* a bit more precise (at most 1ulp error on 1M uniform random
samples in [0,1), the old formula gave some 2ulp errors as well)
support alternate glibc name pow10 for exp10 2012-05-01T04:07:37+00:00 Rich Felker dalias@aerifal.cx 2012-05-01T04:07:37+00:00 a917c03706d94564082d5ccff032a3e2c3dad537 






add C stub for sqrtl (already implemented in asm on i386 and x86_64)
2012-05-01T01:32:19+00:00

Rich Felker
dalias@aerifal.cx

2012-05-01T01:32:19+00:00

da5d89d42fb47c648fb83645e5e4a55bae907032












first try at writing an efficient and "correct" exp10
2012-04-30T07:26:53+00:00

Rich Felker
dalias@aerifal.cx

2012-04-30T07:26:53+00:00

f6819755779a084bf2f82cb90175a4d9a018de73

this is a nonstandard function so it's not clear what conditions it
should satisfy. my intent is that it be fast and exact for positive
integral exponents when the result fits in the destination type, and
fast and correctly rounded for small negative integral exponents.
otherwise we aim for at most 1ulp error; it seems to differ from pow
by at most 1ulp and it's often 2-5 times faster than pow.





this is a nonstandard function so it's not clear what conditions it
should satisfy. my intent is that it be fast and exact for positive
integral exponents when the result fits in the destination type, and
fast and correctly rounded for small negative integral exponents.
otherwise we aim for at most 1ulp error; it seems to differ from pow
by at most 1ulp and it's often 2-5 times faster than pow.







fix typo in the x86_64 rounding asm
2012-04-30T00:36:32+00:00

Rich Felker
dalias@aerifal.cx

2012-04-30T00:36:32+00:00

e5a9b50e973aadff2e6303d01be8807b1fb9c325












new math asm (abs/rounding) for x86_64
2012-04-30T00:31:46+00:00

Rich Felker
dalias@aerifal.cx

2012-04-30T00:31:46+00:00

0e4a995213620630005a14f32c089879f7c0315e

untested





untested







math: fix x86 asin accuracy
2012-04-04T15:34:28+00:00

nsz
nsz@port70.net

2012-04-04T15:34:28+00:00

37eaec3ad396ea4502f4aea2395b3b5f5f69a1c6

use (1-x)*(1+x) instead of (1-x*x) in asin.s
the later can be inaccurate with upward rounding when x is close to 1





use (1-x)*(1+x) instead of (1-x*x) in asin.s
the later can be inaccurate with upward rounding when x is close to 1