[Bug-gsl] MISER subdivision bug

[Rudy Arthur]

Hi,

I found what I think is a bug in miser.c.
In the function estimate_corrmc, line 689 in the file
  sigma_l[i] = sqrt (fsum2_l[i] - fsum_l[i] * fsum_l[i] / hits_l[i]);
is supposed to be the variance in the left hand region with the division in
the i-axis ,
Shouldn't this be:
  sqrt( ( fsum2_l[i] / hits_l[i] - fsum_l[i] * fsum_l[i] ) / (hits_l[i] -
1) )
and correspondingly for the right hand variance? Or is there some subtlety
I am missing?

Rudy.