Using the Neumaier method, take a more accurate sum

Using the Neumaier method, take a more accurate sum

neumaierSum(numbers)

Arguments

numbers

A vector of numbers to sum.

Value

Sum of numbers, a bit more accurate than kahanSum

References

https://en.wikipedia.org/wiki/Kahan_summation_algorithm

Examples

sum(c(1,1e100,1,-1e100)) ## Should be 2, gives 0
#> [1] 0
neumaierSum(c(1,1e100,1,-1e100)) ## Gives 2
#> [1] 2