diff options
Diffstat (limited to '')
| -rw-r--r-- | src/experimental/randomcmwc.cpp | 22 | 
1 files changed, 11 insertions, 11 deletions
| diff --git a/src/experimental/randomcmwc.cpp b/src/experimental/randomcmwc.cpp index 97f7262..a4e807e 100644 --- a/src/experimental/randomcmwc.cpp +++ b/src/experimental/randomcmwc.cpp | |||
| @@ -9,37 +9,37 @@ | |||
| 9 | 9 | ||
| 10 | #define PHI 0x9e3779b9 | 10 | #define PHI 0x9e3779b9 | 
| 11 | 11 | ||
| 12 | RandomCmwc::RandomCmwc() : | 12 | Bu::RandomCmwc::RandomCmwc() : | 
| 13 | q( NULL ), | 13 | q( 0 ), | 
| 14 | c( 362436 ) | 14 | c( 362436 ) | 
| 15 | { | 15 | { | 
| 16 | q = new uint32_t[4096]; | 16 | q = new uint32_t[4096]; | 
| 17 | } | 17 | } | 
| 18 | 18 | ||
| 19 | RandomCmwc::~RandomCmwc() | 19 | Bu::RandomCmwc::~RandomCmwc() | 
| 20 | { | 20 | { | 
| 21 | delete[] q; | 21 | delete[] q; | 
| 22 | } | 22 | } | 
| 23 | 23 | ||
| 24 | void RandomCmwc::seed( int32_t iSeed ) | 24 | void Bu::RandomCmwc::seed( int32_t iSeed ) | 
| 25 | { | 25 | { | 
| 26 | int i; | 26 | int i; | 
| 27 | 27 | ||
| 28 | Q[0] = iSeed; | 28 | q[0] = iSeed; | 
| 29 | Q[1] = iSeed + PHI; | 29 | q[1] = iSeed + PHI; | 
| 30 | Q[2] = iSeed + PHI + PHI; | 30 | q[2] = iSeed + PHI + PHI; | 
| 31 | 31 | ||
| 32 | for (i = 3; i < 4096; i++) | 32 | for (i = 3; i < 4096; i++) | 
| 33 | Q[i] = Q[i - 3] ^ Q[i - 2] ^ PHI ^ i; | 33 | q[i] = q[i - 3] ^ q[i - 2] ^ PHI ^ i; | 
| 34 | } | 34 | } | 
| 35 | 35 | ||
| 36 | int32_t RandomCmwc::rand() | 36 | int32_t Bu::RandomCmwc::rand() | 
| 37 | { | 37 | { | 
| 38 | uint64_t t, a = 18782LL; | 38 | uint64_t t, a = 18782LL; | 
| 39 | static uint32_t i = 4095; | 39 | static uint32_t i = 4095; | 
| 40 | uint32_t x, r = 0xfffffffe; | 40 | uint32_t x, r = 0xfffffffe; | 
| 41 | i = (i + 1) & 4095; | 41 | i = (i + 1) & 4095; | 
| 42 | t = a * Q[i] + c; | 42 | t = a * q[i] + c; | 
| 43 | c = (t >> 32); | 43 | c = (t >> 32); | 
| 44 | x = t + c; | 44 | x = t + c; | 
| 45 | if( x < c ) | 45 | if( x < c ) | 
| @@ -47,6 +47,6 @@ int32_t RandomCmwc::rand() | |||
| 47 | x++; | 47 | x++; | 
| 48 | c++; | 48 | c++; | 
| 49 | } | 49 | } | 
| 50 | return (Q[i] = r - x); | 50 | return (q[i] = r - x); | 
| 51 | } | 51 | } | 
| 52 | 52 | ||
