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
Update floating point tutorial
  • Loading branch information
@rhettinger
rhettinger committed Dec 22, 2022
commit af7613ebcf1d118291320f07e9fab8100266c13c
2 changes: 1 addition & 1 deletion Doc/tutorial/floatingpoint.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -192,7 +192,7 @@ added onto a running total. That can make a difference in overall accuracy
so that the errors do not accumulate to the point where they affect the
final total:

>>> sum([0.1] * 10) == 1.0
>>> 0.1 + 0.1 + 0.1 + 0.1 + 0.1 + 0.1 + 0.1 + 0.1 + 0.1 + 0.1 == 1.0
False
>>> math.fsum([0.1] * 10) == 1.0
True
Expand Down