diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2013-01-07 23:54:13 +0100 |
|---|---|---|
| committer | Szabolcs Nagy <nsz@port70.net> | 2013-01-07 23:54:13 +0100 |
| commit | 121e3a38a1d0b324696fc1671029e4b815ea7fa4 (patch) | |
| tree | 608923bcb097799f051cb509af138d8a49239cc5 /src/math/modff.c | |
| parent | d84923d89e4fe46877334cbaa5e4549f2c4cb8a6 (diff) | |
| download | musl-121e3a38a1d0b324696fc1671029e4b815ea7fa4.tar.gz musl-121e3a38a1d0b324696fc1671029e4b815ea7fa4.tar.bz2 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/modff.c')
0 files changed, 0 insertions, 0 deletions