"sum([1.0, 1e100, 1.0, -1e100])" fails

[johnmyleswhite]

Relevant?

Julia equivalent of Python's "fsum" for floating point summation

I have been looking for a Julia equivalent of Python’s fsum and have thus far been unable to find one. Is there an implementation already available? If not I will write my own and put it on Github for others to use. From the Python docs: math.fsum ( iterable ) Return an accurate floating point sum of values in the iterable. Avoids 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…