pairwiseSum {PreciseSums} | R Documentation |
Return an accurate floating point sum of values
Description
This was taken by NumPy and adapted for use here. It is more accurate than a standard sum, but still has numerical issues. Its main benefit is that it is about the same amount of time as a standard time with the added accuracy.
Usage
pairwiseSum(numbers)
Arguments
numbers |
A vector of numbers to sum. |
Value
A sum of numbers with a rounding error of O(lg n) instead of O(n).
Author(s)
Matthew Fidler (R implementation), Julian Taylor, Nathaniel J Smith, and NumPy team. #' @examples sum(c(1,1e100,1,-1e100)) ## Should be 2, gives 0 pairwiseSum(c(1,1e100,1,-1e100)) ## Should be 2, still 0
References
https://github.com/numpy/numpy/pull/3685
[Package PreciseSums version 0.6 Index]