> Bruce Carneal did some tests of robustness and speed for various normal
> generators. I don't know what his final tests showed for Box-Muller. IIRC,
> it had some failures but nothing spectacular. The tests were pretty
> stringent and based on using the erf to turn the normal distribution into a
On Wed, Dec 10, 2008 at 12:03 PM, Michael Gilbert <
[EMAIL PROTECTED]> wrote:
> Hello,
>
> I have been reading that there may be potential issues with the
> Box-Muller transform, which is used by the numpy.random.normal()
> function. Supposedly, since f*x1 and f*x2 are not independent variable
Wed, 10 Dec 2008 14:03:39 -0500, Michael Gilbert wrote:
> I have been reading that there may be potential issues with the
> Box-Muller transform, which is used by the numpy.random.normal()
> function. Supposedly, since f*x1 and f*x2 are not independent
> variables, then the individual elements
I think the use of a correct uniform generator will allow a good
normal distribution. Congruental generators are very basic generators,
everyone knows they should not be used. I think Numpy uses a Mersenne
Twisted generator, for which you can generate "independant" vectors
with several hundred valu
Hello,
I have been reading that there may be potential issues with the
Box-Muller transform, which is used by the numpy.random.normal()
function. Supposedly, since f*x1 and f*x2 are not independent variables, then
the individual elements (corresponding to f*x1 and f*x2 ) of the
distribution
