Displaying 2 results from an estimated 2 matches for "neumaier".
Did you mean:
neubauer
2019 Feb 21
1
code for sum function
Specifically: https://svn.r-project.org/R/trunk/src/main/summary.c
And if you don't want to deal with Subversion, you can look at the
read-only github mirror:
https://github.com/wch/r-source/blob/e5b21d0397c607883ff25cca379687b86933d730/src/main/summary.c#L115-L131
On Thu, Feb 21, 2019 at 11:57 AM David Winsemius <dwinsemius at comcast.net> wrote:
>
>
> On 2/20/19 2:55 PM,
2019 Feb 19
4
code for sum function
The algorithm does make a differece. You can use Kahan's summation
algorithm (https://en.wikipedia.org/wiki/Kahan_summation_algorithm) to
reduce the error compared to the naive summation algorithm. E.g., in R
code:
naiveSum <-
function(x) {
s <- 0.0
for(xi in x) s <- s + xi
s
}
kahanSum <- function (x)
{
s <- 0.0
c <- 0.0 # running compensation for lost