Difference between revisions of "Proposals:RandomNumbers"
|Line 44:||Line 44:|
Latest revision as of 17:15, 28 March 2010
Resolve platform dependancies with current random number generation. Many registration algorithms depend on random sampling, exactly where I came across the problem.
It would be great to standardize the random number generator to be sufficiently random, fast and generate the same sequences (from a given seed) across platforms. The discussion culminated with a post: http://www.itk.org/mailman/private/insight-developers/2005-January/006220.html I'm not sure any action happened on the suggestions.
More information on mailing list under "random number generator".
It may be worth following the thread http://www.itk.org/mailman/private/insight-developers/2005-January/006189.html A stable portable may also be obtained from http://www-personal.engin.umich.edu/~wagnerr/MersenneTwister.html
The Boost site has a good library for random generators:
Among others, it contains an implementation of the Mersenne Twister. So it may be a helpful source for writing an itkMersenneTwister class, or something like that.
In a previous email to the ITK users mailing list I described some tests to check the actual 'randomness' of two random generators, vnl_sample (which is currently used in ITK) and vnl_random (which is in a more recent version of VXL):
It appeared that vnl_random is much better. What I didn't mention in this email, is that I also tested the MersenneTwister from the Boost library. It gave similar results as vnl_random.
So, I would say that the MersenneTwister indeed is a good option.
MersenneTwister has been added to ITK. Its faster and from the tests that I have, much more random. I've committed some tests to ITK. See for instance http://www.itk.org/Testing/Sites/bathurst.bwh.harvard.edu/Linux-icpc-i686-release/20050610-0100-Nightly/Results/__Testing_Code_Numerics_Statistics_itkMersenneTwisterRandomVariateGeneratorTest_cxx.html In general, its anywhere from 1 to 30 times faster than the vnl counterparts depending on the platform.
I suggest future uses of random number generator requiring algorithms try to use MersenneTwister. It provides methods to get both Normal and uniform distributed random variables. ALso it will not suffer from re-seeding issues on cygwin etc as in http://public.kitware.com/pipermail/insight-users/2005-May/012962.html