diff options
Diffstat (limited to 'src')
| -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 | ||
