bug-gsl
[Top][All Lists]

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

 From: Brian Gough Subject: Re: [Bug-gsl] [PATCH] a one-pass run algorithm for sigma (histograms) Date: Tue, 01 Dec 2009 14:07:08 +0000 User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI)

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/