[PATCH 0/3] Fix powf inaccuracies up to ~169 ULP reported by Paul Zimmermann
Fabian Schriever <[email protected]>
| Newsgroups | gmane.comp.lib.newlib |
|---|---|
| Message-ID | <[email protected]> |
This patch series fixes the powf accuracy errors reported by Paul Zimmermann on the 6th of August.
All three patches are taken from FreeBSD and authored by Bruce Evans, not me. The original commits are referenced in the individual patches. The patches also contain the original messages.
Having run our test setup the error is reduced from 169 ULP to just shy of 1 ULP (0.9998...). Note however that, while quite exhaustive, this is not a full range test, and there still may be unobserved errors above 1 ULP.
Fabian Schriever (3):
powf: Fix the hi+lo decomposition for 2/(3ln2) (FreeBSD)
powf: Fixed 2 bugs in the computation /* t_h=ax+bp[k] High */.
(FreeBSD)
powf: Fixed another precision bug in powf() (FreeBSD)
newlib/libm/math/ef_pow.c | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
--
2.33.0.windows.1