Tim Peters <[email protected]> added the comment:
Here's a "correct rounding" fail for the add_on approach:
xs = [16.000000000000004] * 9
decimal result = 48.00000000000001065814103642
which rounds to float 48.000000000000014
add_on result: 48.00000000000001
That's about 0.5000000000026 ulp too small - shocking ;-)
The fsum approach got this one right, but has similar failures on other inputs
of the same form; i.e.,
[i + ulp(i)] * 9
for various integers `i`.
----------
_______________________________________
Python tracker <[email protected]>
<https://bugs.python.org/issue41513>
_______________________________________
_______________________________________________
Python-bugs-list mailing list
Unsubscribe:
https://mail.python.org/mailman/options/python-bugs-list/archive%40mail-archive.com
