Connecting Tech Pros Worldwide Forums | Help | Site Map

random number generator

asdf
Guest
 
Posts: n/a
#1: Oct 21 '06
I want a random number generator, the random number should be subject a
uniform distribution in [0,1]. Could you please give me some hints?
Thanks.


Daniel T.
Guest
 
Posts: n/a
#2: Oct 21 '06

re: random number generator


"asdf" <likemursili@gmail.comwrote:
Quote:
I want a random number generator, the random number should be subject a
uniform distribution in [0,1]. Could you please give me some hints?
Thanks.
http://members.cox.net/srice1/random/crandom.html

--
There are two things that simply cannot be doubted, logic and perception.
Doubt those, and you no longer*have anyone to discuss your doubts with,
nor any ability to discuss them.
Peter Jansson
Guest
 
Posts: n/a
#3: Oct 21 '06

re: random number generator


asdf wrote:
Quote:
I want a random number generator, the random number should be subject a
uniform distribution in [0,1]. Could you please give me some hints?
Thanks.
>
You might be interested in the Mersenne Twister:
http://www.math.sci.hiroshima-u.ac.j...at/MT/emt.html

Which have a link to:
http://www.math.sci.hiroshima-u.ac.j...ES/mt19937ar.c


Sincerely,

Peter Jansson
http://www.p-jansson.com/
http://www.jansson.net/
Pete Becker
Guest
 
Posts: n/a
#4: Oct 22 '06

re: random number generator


asdf wrote:
Quote:
I want a random number generator, the random number should be subject a
uniform distribution in [0,1]. Could you please give me some hints?
Thanks.
>
If you've got a copy of TR1, there are random number generators galore
in it. Dinkumware implements all of the TR. Boost has, among other
pieces, the random number generators (although I haven't checked to see
if they're up to date with the changes made for TR1).

--

-- Pete

Author of "The Standard C++ Library Extensions: a Tutorial and
Reference." For more information about this book, see
www.petebecker.com/tr1book.
noone
Guest
 
Posts: n/a
#5: Oct 23 '06

re: random number generator


On Sat, 21 Oct 2006 07:53:59 -0700, asdf wrote:
Quote:
I want a random number generator, the random number should be subject a
uniform distribution in [0,1]. Could you please give me some hints?
Thanks.
Are you sure you want (1) included in your range?

Others have given you pointers to RNG literature. One thing I've noticed
that seems to be an almost universal misuse of random number functions
such as the srand() and rand() functions is that programmers who
use them to generate numbers [0..1) virtually NEVER check to see if the
returned value of rand()==RAND_MAX. While the probability of this
happening is 1/RAND_MAX, you should check because rand()/RAND_MAX could
end up outside of your intended range.


For simple non-cryptographic uniform distributions the rand() function is
often adequate but for more stringent requirements the BOOST library
provides several RNGs. Of course as with many open source projects, the
documentation is lacking...at least it was when I last used them.


Steve Pope
Guest
 
Posts: n/a
#6: Oct 23 '06

re: random number generator


noone <noone@all.comwrote:
Quote:
One thing I've noticed that seems to be an almost universal
misuse of random number functions such as the srand() and rand()
functions is that programmers who use them to generate numbers
[0..1) virtually NEVER check to see if the returned value of
rand()==RAND_MAX. While the probability of this happening is
1/RAND_MAX, you should check because rand()/RAND_MAX could end
up outside of your intended range.
I usually do something like

int x = rand() & 0x3fffffff;
double y = (double) x / (double) 0x40000000;

y is now uniform in [0,1).

A test for RAND_MAX >= 0x3fffffff would be nice as you say.
On my system RAND_MAX is pow(2,31)-1 which is what you'd expect,
so it works.

Steve
Pete Becker
Guest
 
Posts: n/a
#7: Oct 23 '06

re: random number generator


noone wrote:
Quote:
On Sat, 21 Oct 2006 07:53:59 -0700, asdf wrote:
>
Quote:
>I want a random number generator, the random number should be subject a
>uniform distribution in [0,1]. Could you please give me some hints?
>Thanks.
>
Are you sure you want (1) included in your range?
>
Others have given you pointers to RNG literature. One thing I've noticed
that seems to be an almost universal misuse of random number functions
such as the srand() and rand() functions is that programmers who
use them to generate numbers [0..1) virtually NEVER check to see if the
returned value of rand()==RAND_MAX. While the probability of this
happening is 1/RAND_MAX, you should check because rand()/RAND_MAX could
end up outside of your intended range.
>
Better yet, don't check, but map the range so that this isn't a problem:

rand() / ((double)RAND_MAX + 1)
Quote:
>
For simple non-cryptographic uniform distributions the rand() function is
often adequate but for more stringent requirements the BOOST library
provides several RNGs. Of course as with many open source projects, the
documentation is lacking...at least it was when I last used them.
>
>
There's documentation for the TR1 random number generators in chapter 13
of my book, "The Standard C++ Library Extensions." The Boost generators
are close to what's in TR1.

--

-- Pete

Author of "The Standard C++ Library Extensions: a Tutorial and
Reference." For more information about this book, see
www.petebecker.com/tr1book.
Martin Steen
Guest
 
Posts: n/a
#8: Oct 23 '06

re: random number generator


asdf wrote:
Quote:
I want a random number generator, the random number should be subject a
uniform distribution in [0,1]. Could you please give me some hints?
Thanks.
>
double r = (double) rand() / RAND_MAX;

http://www.cplusplus.com/ref/cstdlib/rand.html

-Martin




Pete Becker
Guest
 
Posts: n/a
#9: Oct 23 '06

re: random number generator


Pete Becker wrote:
Quote:
noone wrote:
Quote:
>On Sat, 21 Oct 2006 07:53:59 -0700, asdf wrote:
>>
Quote:
>>I want a random number generator, the random number should be subject a
>>uniform distribution in [0,1]. Could you please give me some hints?
>>Thanks.
>>
>Are you sure you want (1) included in your range?
>>
>Others have given you pointers to RNG literature. One thing I've noticed
>that seems to be an almost universal misuse of random number functions
>such as the srand() and rand() functions is that programmers who
>use them to generate numbers [0..1) virtually NEVER check to see if the
>returned value of rand()==RAND_MAX. While the probability of this
>happening is 1/RAND_MAX, you should check because rand()/RAND_MAX could
>end up outside of your intended range.
>>
>
Better yet, don't check, but map the range so that this isn't a problem:
>
rand() / ((double)RAND_MAX + 1)
>
Just a clarification (since I confused myself): this gives a uniform
distribution over the half-open range [0.0, 1.0), as the message from
"noone" suggests. It does not answer the original question of how to get
a uniform distribution over the closed range [0.0, 1.0].

--

-- Pete

Author of "The Standard C++ Library Extensions: a Tutorial and
Reference." For more information about this book, see
www.petebecker.com/tr1book.
Closed Thread