fsum: Return an accurate floating point sum of values
Description
This method avoids loss of precision by tracking multiple
intermediate partial sums. Based on python's math.fsum
Usage
fsum(numbers)
Value
Sum of numbers without loss of precision
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 builds, the underlying C library uses extended
precision addition and may occasionally double-round an
intermediate sum causing it to be off in its least significant bit.
Arguments
numbers
A vector of numbers to sum.
Author
Matthew Fidler (R implementation), Raymond Hettinger,
Jonathan Shewchuk, Python Team