Changes in libcfa/src/stdlib.cfa [1959528:15b05c8]
- File:
-
- 1 edited
-
libcfa/src/stdlib.cfa (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
libcfa/src/stdlib.cfa
r1959528 r15b05c8 10 10 // Created On : Thu Jan 28 17:10:29 2016 11 11 // Last Modified By : Peter A. Buhr 12 // Last Modified On : Wed Jan 12 16:25:36202213 // Update Count : 5 7812 // Last Modified On : Mon Jan 3 09:36:27 2022 13 // Update Count : 519 14 14 // 15 15 16 16 #include "stdlib.hfa" 17 //#include "concurrency/kernel/fwd.hfa"18 #include "concurrency/invoke.h" // random_state19 17 20 18 //--------------------------------------- … … 223 221 //--------------------------------------- 224 222 225 // Pipelined to allow OoO overlap with reduced dependencies. Critically, return the current value, and compute and store 226 // the next value. 223 static uint32_t seed = 0; // current seed 224 static thread_local uint32_t state; // random state 225 226 void set_seed( uint32_t seed_ ) { state = seed = seed_; } 227 uint32_t get_seed() { return seed; } 227 228 228 229 #define GENERATOR LCG 229 230 230 static inline uint32_t MarsagliaXor( uint32_t & state ) { 231 uint32_t ret = state; 231 inline uint32_t MarsagliaXor( uint32_t & state ) { 232 if ( unlikely( seed == 0 ) ) set_seed( rdtscl() ); 233 else if ( unlikely( state == 0 ) ) state = seed; 232 234 state ^= state << 6; 233 235 state ^= state >> 21; 234 236 state ^= state << 7; 235 return ret;237 return state; 236 238 } // MarsagliaXor 237 239 238 static inline uint32_t LCG( uint32_t & state ) {// linear congruential generator239 uint32_t ret = state;240 state = 36969 * (state & 65535) + (state >> 16); // 36969 is NOT prime! No not change it!241 return ret;240 inline uint32_t LCG( uint32_t & state ) { // linear congruential generator 241 if ( unlikely( seed == 0 ) ) set_seed( rdtscl() ); 242 else if ( unlikely( state == 0 ) ) state = seed; 243 return state = 36969 * (state & 65535) + (state >> 16); // 36969 is NOT prime! 242 244 } // LCG 243 245 244 uint32_t __global_random_seed;245 246 void set_seed( PRNG & prng, uint32_t seed_ ) with( prng ) { state = seed = seed_; GENERATOR( state ); } // set seed247 246 uint32_t prng( PRNG & prng ) with( prng ) { callcnt += 1; return GENERATOR( state ); } 248 247 249 void set_seed( uint32_t seed ) { 250 active_thread()->random_state = __global_random_seed = seed; 251 GENERATOR( active_thread()->random_state ); 252 } // set_seed 253 uint32_t get_seed() { return __global_random_seed; } 254 uint32_t prng( void ) { return GENERATOR( active_thread()->random_state ); } // [0,UINT_MAX] 248 uint32_t prng( void ) { return GENERATOR( state ); } 255 249 256 250 //---------------------------------------
Note:
See TracChangeset
for help on using the changeset viewer.