Re: [Bug-gsl] [PATCH] a one-pass run algorithm for sigma (histograms)

[Brian Gough]

At Sat, 28 Nov 2009 07:57:49 +1300,
Johannes Buchner wrote:
> I looked into the algorithm in the mentioned reference.
> However, the formula 1.6b, for n values at x_k, boils down to 
>  \sum_(1+k)^(n+k){1/l}. This can be calculated if there is a discrete
> weight (n values at x_k). I was unable to find a shorthand formula for
> this sum (so it can be generalized from n to a continuous value).
> Thus I conclude that it is not possible to use the mentioned algorithm
> for continuous histogram weights.

Ok, thanks for looking into it.  Since the current algorithm is
working ok I will leave it alone, as the speed is not a major issue.

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/