[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Reg :: My Gaussian Random Number Code

From: John D Lamb
Subject: Re: [Help-gsl] Reg :: My Gaussian Random Number Code
Date: Wed, 23 Jun 2010 07:33:56 +0100

On Tue, 2010-06-22 at 13:04 +0200, Sumit Adhikari wrote:

> P gauss_rand = sigma*y*sqrt(-2.0*log(rad)/rad) ;
>       return (gauss_rand);

The code looks correct. If you want a quick, simple test generate a few
hundred pairs of variates and produce a scatter plot.

The code is, however, not particularly efficient. The gsl_ran_gaussian
code should run faster and the gsl_ran_gaussian_ziggurat code faster
still. Even a simple Box–Müller generator should be quicker because you
can get two Gaussian variates from two random numbers.


reply via email to

[Prev in Thread] Current Thread [Next in Thread]