
Boost Users : 
Subject: [Boostusers] [accumulators] Kahan algorithmic variants of sum_impl, mean_impl?
From: Rhys Ulerich (rhys.ulerich_at_[hidden])
Date: 20110930 10:28:44
>From a look at boost::accumulators::impl::sum_impl [1] it appears a
straightforward, fast 'sum += value' is used to accumulate the running
sum. That can fall prey to numerical errors if the sum is large and
the value is small. A slower but more precise summation is possible
using the Kahan summation algorithm [2].
Has anyone implemented a Kahanbased variant for either sum_impl or
mean_impl [3]?
Assuming not, other than to mimic existing code, any implementation
suggestions? One gotcha for Kahan summation is that optimization
levels employing additive associativity can spoil the computation.
Are there other places in the Boost sources where this has become an
issue and where I might look for a portable way around it?
Thanks,
Rhys
[1] http://svn.boost.org/svn/boost/trunk/boost/accumulators/statistics/sum.hpp
[2] http://en.wikipedia.org/wiki/Kahan_summation_algorithm
[3] http://svn.boost.org/svn/boost/trunk/boost/accumulators/statistics/mean.hpp
Boostusers list run by williamkempf at hotmail.com, kalb at libertysoft.com, bjorn.karlsson at readsoft.com, gregod at cs.rpi.edu, wekempf at cox.net