-
-
Couldn't load subscription status.
- Fork 5.7k
Closed
Description
This is my fault in #4039, I think: although cumsum implements a "pairwise" algorithm, unlike the sum function it doesn't actually get any accuracy benefit from this, and is in fact equivalent to naive summation:
x = rand(10^7)
foldl(+, 0.0, x) == cumsum(x)[end]produces true. Is there an easy way to rearrange cumsum_pairwise so that it is as accurate as sum?
Metadata
Metadata
Assignees
Labels
No labels