random number generator 
October 21st, 2006, 03:35 PM
| | | |
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. | 
October 21st, 2006, 04:25 PM
| | | | 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. | 
October 22nd, 2006, 03:55 PM
| | | | 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. | 
October 23rd, 2006, 03:35 AM
| | | | 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. | 
October 23rd, 2006, 07:15 AM
| | | | 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 | 
October 23rd, 2006, 11:35 AM
| | | | 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. | 
October 23rd, 2006, 09:05 PM
| | | | 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 | 
October 23rd, 2006, 10:15 PM
| | | | 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. |  | | | | /bytes/about
We are a network of experts and professionals in IT and software development that help one another with answers to tough questions and share insights.
Get the best answers to your questions from over 225,662 network members.
|