Skip to content

GH-100425: Improve accuracy of builtin sum() for float inputs #100426

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
rhettinger merged 13 commits into from
Dec 23, 2022
Merged

GH-100425: Improve accuracy of builtin sum() for float inputs #100426

Changes from 1 commit
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
916ab7e
Neumaier compensated summation
rhettinger Dec 22, 2022
787f2e5
Kahan summation
rhettinger Dec 22, 2022
bc308aa
Revert "Kahan summation"
rhettinger Dec 22, 2022
0c2ba04
Restore the bogus sign tests
rhettinger Dec 22, 2022
89e46c8
Add tests and comments
rhettinger Dec 22, 2022
a9cb7d2
Add tests and docs
rhettinger Dec 22, 2022
839cc05
Add blurb
rhettinger Dec 22, 2022
66db937
Remove out of date example for fsum()
rhettinger Dec 22, 2022
af7613e
Update floating point tutorial
rhettinger Dec 22, 2022
4f4b731
Correct handling of overflowed or infinite sums
rhettinger Dec 23, 2022
4da0b8f
Move local declarations inside the block where they are used.
rhettinger Dec 23, 2022
6c7894c
For consistent results, apply compensation on the slow path as well.
rhettinger Dec 23, 2022
3cf9536
Mark as an implementation detail.
rhettinger Dec 23, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Remove out of date example for fsum()
  • Loading branch information
@rhettinger
rhettinger committed Dec 22, 2022
commit 66db9377f5c0317149b9c3c1bd31b68d449a92ae
7 changes: 1 addition & 6 deletions Doc/library/math.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -108,12 +108,7 @@ Number-theoretic and representation functions
.. function:: fsum(iterable)

Return an accurate floating point sum of values in the iterable. Avoids
loss of precision by tracking multiple intermediate partial sums:

>>> sum([.1, .1, .1, .1, .1, .1, .1, .1, .1, .1])
0.9999999999999999
>>> fsum([.1, .1, .1, .1, .1, .1, .1, .1, .1, .1])
1.0
loss of precision by tracking multiple intermediate partial sums.

The algorithm's accuracy depends on IEEE-754 arithmetic guarantees and the
typical case where the rounding mode is half-even. On some non-Windows
Expand Down