about summary refs log tree commit diff
path: root/src/math/ldexp.c
diff options
context:
space:
mode:
authorSzabolcs Nagy <nsz@port70.net>2013-01-07 23:54:13 +0100
committerSzabolcs Nagy <nsz@port70.net>2013-01-07 23:54:13 +0100
commit121e3a38a1d0b324696fc1671029e4b815ea7fa4 (patch)
tree608923bcb097799f051cb509af138d8a49239cc5 /src/math/ldexp.c
parentd84923d89e4fe46877334cbaa5e4549f2c4cb8a6 (diff)
downloadmusl-121e3a38a1d0b324696fc1671029e4b815ea7fa4.tar.gz
musl-121e3a38a1d0b324696fc1671029e4b815ea7fa4.tar.xz
musl-121e3a38a1d0b324696fc1671029e4b815ea7fa4.zip
math: erf and erfc cleanup
common part of erf and erfc was put in a separate function which
saved some space and the new code is using unsigned arithmetics

erfcf had a bug: for some inputs in [7.95,8] the result had
more than 60ulp error: in expf(-z*z - 0.5625f) the argument
must be exact but not enough lowbits of z were zeroed,
-SET_FLOAT_WORD(z, ix&0xfffff000);
+SET_FLOAT_WORD(z, ix&0xffffe000);
fixed the issue
Diffstat (limited to 'src/math/ldexp.c')
0 files changed, 0 insertions, 0 deletions