randomize() rather than system time seed?

The basics of how to develop, test, and use models.
Post Reply
abogaard
Posts: 8
Joined: Mon Feb 05, 2007 4:21 pm
Location: ann arbor, mi

randomize() rather than system time seed?

Post by abogaard »

Hello all,

I am attempting to get a true output of random numbers utilizing the c function randomize() in stdlib.h

I must be missing something critical about how neuron compiles everything because I cannot get it to work.

Any ideas?
ted
Site Admin
Posts: 6384
Joined: Wed May 18, 2005 4:50 pm
Location: Yale University School of Medicine
Contact:

Post by ted »

First, a word about terminology. Digital computers don't generate "real" random numbers.
They generate pseudorandom sequences that may pass various tests for statistical
independence, but they're not random in the same sense as, for example, the sequence
of times at which a radioactive material emits alpha particles.

Now for your question. The Random class makes it unnecessary to reach into the bowels
of a C library. If you need to generate different random sequences, use a different seed for
each sequence. To ensure reproducibility of resuts, keep a record of the seeds that were
used for each run. The highest quality generator is MCellRan4; to ensure statistical
independence of pseudorandom number streams, note that
each stream should be statistically independent as long as the highindex values differ by more than the eventual length of the stream.
--see
http://www.neuron.yale.edu/neuron/stati ... #MCellRan4
Post Reply