diff options
| author | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2018-07-03 15:05:31 +0300 |
|---|---|---|
| committer | Corinna Vinschen <corinna@vinschen.de> | 2018-07-06 11:29:01 +0300 |
| commit | 2805b07fa1f54ef001c2ce78618e8efc87155232 (patch) | |
| tree | 1d54b22c61463070bd3ee28fa9014c19962685e4 /newlib/libm/common/log_data.c | |
| parent | 6a85e1a4e5a01dec3b423eb172b3ad9155487e75 (diff) | |
Fix large ulp error in pow without fma very near 1.0
The !HAVE_FAST_FMA code path split r = z/c - 1 into r = rhi + rlo such
that when z = 1-tiny and c = 1 then rlo and rhi could have much larger
magnitude than r which later caused large rounding errors.
So do a nearest rounding instead of truncation at the split.
In newlib with default settings this was observable on some arm targets
that enable the new math code but has no fma.
Diffstat (limited to 'newlib/libm/common/log_data.c')
0 files changed, 0 insertions, 0 deletions