ITK
6.0.0
Insight Toolkit
|
#include <itkMersenneTwisterRandomVariateGenerator.h>
MersenneTwisterRandom random variate generator.
It is recommended to create a separate object in each thread. By default, each instantiated class will have a different seed created by the GetNextSeed method. The creation of the initial seeds are initialized once from the time. For deterministic behavior, the individual instances' seeds should be manual set to separate values in each thread.
It is no longer recommended to use this class using a "Singleton-like" GetInstance method for the global instance of this class. This usage may result in unsafe concurrent access to the global instance.
This notice was included with the original implementation. The only changes made were to obfuscate the author's email addresses.
MersenneTwister.h Mersenne Twister random number generator – a C++ class MTRand Based on code by Makoto Matsumoto, Takuji Nishimura, and Shawn Cokus Richard J. Wagner v1.0 15 May 2003 rjwagner at writeme dot com
The Mersenne Twister is an algorithm for generating random numbers. It was designed with consideration of the flaws in various other generators. The period, 2^19937-1, and the order of equidistribution, 623 dimensions, are far greater. The generator is also fast; it avoids multiplication and division, and it benefits from caches and pipelines. For more information see the inventors' web page at http:*www.math.keio.ac.jp/~matumoto/emt.html
Reference M. Matsumoto and T. Nishimura, "Mersenne Twister: A 623-Dimensionally Equidistributed Uniform Pseudo-Random Number Generator", ACM Transactions on Modeling and Computer Simulation, Vol. 8, No. 1, January 1998, pp 3-30.
Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, Copyright (C) 2000 - 2003, Richard J. Wagner All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
The original code included the following notice:
When you use this, send an email to: matumoto at math dot keio dot ac dot jp with an appropriate reference to your work.
It would be nice to CC: rjwagner at writeme dot com and Cokus at math dot washington dot edu when you write.
Definition at line 127 of file itkMersenneTwisterRandomVariateGenerator.h.
Public Types | |
using | ConstPointer = SmartPointer< const Self > |
using | IntegerType = uint32_t |
using | Pointer = SmartPointer< Self > |
using | Self = MersenneTwisterRandomVariateGenerator |
using | Superclass = RandomVariateGeneratorBase |
Public Types inherited from itk::Statistics::RandomVariateGeneratorBase | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = RandomVariateGeneratorBase |
using | Superclass = Object |
Public Types inherited from itk::Object | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = Object |
using | Superclass = LightObject |
Public Types inherited from itk::LightObject | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = LightObject |
Public Member Functions | |
double | Get53BitVariate () |
IntegerType | GetIntegerVariate (const IntegerType &n) |
const char * | GetNameOfClass () const override |
IntegerType | GetSeed () const |
double | GetUniformVariate (const double a, const double b) |
double | GetVariate () override |
double | GetVariateWithClosedRange () |
double | GetVariateWithClosedRange (const double n) |
double | GetVariateWithOpenRange () |
double | GetVariateWithOpenRange (const double n) |
double | GetVariateWithOpenUpperRange () |
double | GetVariateWithOpenUpperRange (const double n) |
void | Initialize () |
void | Initialize (const IntegerType seed) |
double | operator() () |
IntegerType | GetIntegerVariate () |
double | GetNormalVariate (const double mean=0.0, const double variance=1.0) |
void | SetSeed (const IntegerType oneSeed) |
void | SetSeed () |
Public Member Functions inherited from itk::Statistics::RandomVariateGeneratorBase | |
const char * | GetNameOfClass () const override |
Public Member Functions inherited from itk::Object | |
unsigned long | AddObserver (const EventObject &event, Command *cmd) const |
unsigned long | AddObserver (const EventObject &event, std::function< void(const EventObject &)> function) const |
LightObject::Pointer | CreateAnother () const override |
virtual void | DebugOff () const |
virtual void | DebugOn () const |
Command * | GetCommand (unsigned long tag) |
bool | GetDebug () const |
MetaDataDictionary & | GetMetaDataDictionary () |
const MetaDataDictionary & | GetMetaDataDictionary () const |
virtual ModifiedTimeType | GetMTime () const |
virtual const TimeStamp & | GetTimeStamp () const |
bool | HasObserver (const EventObject &event) const |
void | InvokeEvent (const EventObject &) |
void | InvokeEvent (const EventObject &) const |
virtual void | Modified () const |
void | Register () const override |
void | RemoveAllObservers () |
void | RemoveObserver (unsigned long tag) const |
void | SetDebug (bool debugFlag) const |
void | SetReferenceCount (int) override |
void | UnRegister () const noexcept override |
void | SetMetaDataDictionary (const MetaDataDictionary &rhs) |
void | SetMetaDataDictionary (MetaDataDictionary &&rrhs) |
virtual void | SetObjectName (std::string _arg) |
virtual const std::string & | GetObjectName () const |
Public Member Functions inherited from itk::LightObject | |
Pointer | Clone () const |
virtual void | Delete () |
virtual int | GetReferenceCount () const |
void | Print (std::ostream &os, Indent indent=0) const |
Static Public Member Functions | |
static Pointer | GetInstance () |
static IntegerType | GetNextSeed () |
static Pointer | New () |
static void | ResetNextSeed () |
Static Public Member Functions inherited from itk::Object | |
static bool | GetGlobalWarningDisplay () |
static void | GlobalWarningDisplayOff () |
static void | GlobalWarningDisplayOn () |
static Pointer | New () |
static void | SetGlobalWarningDisplay (bool val) |
Static Public Member Functions inherited from itk::LightObject | |
static void | BreakOnError () |
static Pointer | New () |
Static Public Attributes | |
static constexpr IntegerType | StateVectorLength = 624 |
Protected Member Functions | |
IntegerType | hiBit (const IntegerType &u) const |
IntegerType | loBit (const IntegerType &u) const |
IntegerType | loBits (const IntegerType &u) const |
MersenneTwisterRandomVariateGenerator () | |
IntegerType | mixBits (const IntegerType &u, const IntegerType &v) const |
void | PrintSelf (std::ostream &os, Indent indent) const override |
void | reload () |
IntegerType | twist (const IntegerType &m, const IntegerType &s0, const IntegerType &s1) const |
~MersenneTwisterRandomVariateGenerator () override | |
Protected Member Functions inherited from itk::Statistics::RandomVariateGeneratorBase | |
RandomVariateGeneratorBase () | |
~RandomVariateGeneratorBase () override | |
Protected Member Functions inherited from itk::Object | |
Object () | |
bool | PrintObservers (std::ostream &os, Indent indent) const |
void | PrintSelf (std::ostream &os, Indent indent) const override |
virtual void | SetTimeStamp (const TimeStamp &timeStamp) |
~Object () override | |
Protected Member Functions inherited from itk::LightObject | |
virtual LightObject::Pointer | InternalClone () const |
LightObject () | |
virtual void | PrintHeader (std::ostream &os, Indent indent) const |
virtual void | PrintTrailer (std::ostream &os, Indent indent) const |
virtual | ~LightObject () |
Static Protected Member Functions | |
static IntegerType | hash (time_t t, clock_t c) |
Protected Attributes | |
int | m_Left {} |
IntegerType * | m_PNext {} |
std::atomic< IntegerType > | m_Seed {} |
IntegerType | state [StateVectorLength] |
Protected Attributes inherited from itk::LightObject | |
std::atomic< int > | m_ReferenceCount {} |
Static Protected Attributes | |
static constexpr unsigned int | M = 397 |
Private Member Functions | |
itkGetGlobalDeclarationMacro (MersenneTwisterGlobals, PimplGlobals) | |
Static Private Member Functions | |
static Pointer | CreateInstance () |
Private Attributes | |
std::mutex | m_InstanceMutex {} |
Static Private Attributes | |
static MersenneTwisterGlobals * | m_PimplGlobals |
using itk::Statistics::MersenneTwisterRandomVariateGenerator::ConstPointer = SmartPointer<const Self> |
Definition at line 134 of file itkMersenneTwisterRandomVariateGenerator.h.
using itk::Statistics::MersenneTwisterRandomVariateGenerator::IntegerType = uint32_t |
Definition at line 136 of file itkMersenneTwisterRandomVariateGenerator.h.
Definition at line 133 of file itkMersenneTwisterRandomVariateGenerator.h.
using itk::Statistics::MersenneTwisterRandomVariateGenerator::Self = MersenneTwisterRandomVariateGenerator |
Standard class type aliases.
Definition at line 131 of file itkMersenneTwisterRandomVariateGenerator.h.
using itk::Statistics::MersenneTwisterRandomVariateGenerator::Superclass = RandomVariateGeneratorBase |
Definition at line 132 of file itkMersenneTwisterRandomVariateGenerator.h.
|
protected |
|
overrideprotected |
|
staticprivate |
Internal method to actually create a new object.
|
inline |
Access to 53-bit random numbers (capacity of IEEE double precision) in the range [0,1)
Definition at line 503 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetIntegerVariate().
|
static |
Return the global Mersenne Twister instance.
This method returns a Singleton of the Mersenne Twister. The seed is initialized from the wall clock at first use, but can be globally set using the resulting instance's SetSeed().
|
inline |
Get an integer variate in [0, 2^32-1]
Definition at line 421 of file itkMersenneTwisterRandomVariateGenerator.h.
References m_Left, m_PNext, and reload().
Referenced by Get53BitVariate(), GetIntegerVariate(), GetVariateWithClosedRange(), GetVariateWithOpenRange(), and GetVariateWithOpenUpperRange().
|
inline |
Get an integer variate in [0, n] for n < 2^32
Definition at line 479 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetIntegerVariate().
|
overridevirtual |
Reimplemented from itk::LightObject.
|
static |
Return the next seed, derived as a sequence from the seed of the singleton instance.
|
inline |
Access to a normal random number distribution.
Definition at line 514 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariateWithOpenRange(), GetVariateWithOpenUpperRange(), and itk::Math::pi.
|
inline |
Return the current seed
Definition at line 414 of file itkMersenneTwisterRandomVariateGenerator.h.
References m_Seed.
|
inline |
Access to a uniform random number distribution
Definition at line 528 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariateWithOpenUpperRange().
|
inlineoverridevirtual |
Get a variate in the range [0, 1] Do NOT use for CRYPTOGRAPHY without securely hashing several returned values together, otherwise the generator state can be learned after reading 624 consecutive values.
Implements itk::Statistics::RandomVariateGeneratorBase.
Definition at line 536 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariateWithClosedRange().
Referenced by operator()().
|
inline |
Get a random variate in the range [0, 1]
Definition at line 438 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetIntegerVariate().
Referenced by GetVariate(), and GetVariateWithClosedRange().
|
inline |
Get a random variate in the range [0, n]
Definition at line 445 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariateWithClosedRange().
|
inline |
Get a range variate in the range (0, 1)
Definition at line 466 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetIntegerVariate().
Referenced by GetNormalVariate(), and GetVariateWithOpenRange().
|
inline |
Get a range variate in the range (0, n)
Definition at line 473 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariateWithOpenRange().
|
inline |
Get a range variate in the range [0, 1)
Definition at line 452 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetIntegerVariate().
Referenced by GetNormalVariate(), GetUniformVariate(), and GetVariateWithOpenUpperRange().
|
inline |
Get a range variate in the range [0, n)
Definition at line 459 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariateWithOpenUpperRange().
|
staticprotected |
Referenced by SetSeed().
|
inlineprotected |
Definition at line 283 of file itkMersenneTwisterRandomVariateGenerator.h.
|
inline |
Definition at line 393 of file itkMersenneTwisterRandomVariateGenerator.h.
References SetSeed().
Referenced by SetSeed().
|
inline |
Initialize with a simple IntegerType
Definition at line 344 of file itkMersenneTwisterRandomVariateGenerator.h.
References m_InstanceMutex, m_Seed, reload(), state, and StateVectorLength.
|
private |
Only used to synchronize the global variable across static libraries.
|
inlineprotected |
Definition at line 288 of file itkMersenneTwisterRandomVariateGenerator.h.
|
inlineprotected |
Definition at line 293 of file itkMersenneTwisterRandomVariateGenerator.h.
|
inlineprotected |
Definition at line 298 of file itkMersenneTwisterRandomVariateGenerator.h.
|
static |
Method for creation through the object factory.
This method allocates a new instance of a Mersenne Twister, and initializes it with the next seed from the global instance's seed.
Referenced by itk::RandomPermutation::RandomPermutation().
|
inline |
Same as GetVariate()
Definition at line 542 of file itkMersenneTwisterRandomVariateGenerator.h.
References GetVariate().
|
overrideprotectedvirtual |
Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes.
Reimplemented from itk::LightObject.
|
inlineprotected |
Reload array with N new values.
Definition at line 366 of file itkMersenneTwisterRandomVariateGenerator.h.
References M, m_Left, m_PNext, state, StateVectorLength, and twist().
Referenced by GetIntegerVariate(), and Initialize().
|
static |
Resets the internal data that is used to calculate the next seed. (Does not reset the initial seed.) Allows generating a reproducible sequence of pseudo-random numbers.
|
inline |
Re-seeding functions with same behavior as initializers
\note This method is thread-safe.
Definition at line 406 of file itkMersenneTwisterRandomVariateGenerator.h.
References hash().
Referenced by Initialize().
|
inline |
Re-seeding functions with same behavior as initializers
\note This method is thread-safe.
Definition at line 399 of file itkMersenneTwisterRandomVariateGenerator.h.
References Initialize().
|
inlineprotected |
Definition at line 304 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by reload().
|
staticconstexprprotected |
Period parameter
Definition at line 276 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by reload().
|
private |
Definition at line 333 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by Initialize().
|
protected |
Definition at line 319 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by GetIntegerVariate(), and reload().
|
staticprivate |
Definition at line 337 of file itkMersenneTwisterRandomVariateGenerator.h.
|
protected |
Definition at line 316 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by GetIntegerVariate(), and reload().
|
protected |
Definition at line 322 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by GetSeed(), and Initialize().
|
protected |
Definition at line 313 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by Initialize(), and reload().
|
staticconstexpr |
Length of state vector
Definition at line 169 of file itkMersenneTwisterRandomVariateGenerator.h.
Referenced by Initialize(), and reload().