Ignore:
Timestamp:
Jan 12, 2022, 6:34:58 PM (4 years ago)
Author:
Peter A. Buhr <pabuhr@…>
Branches:
ADT, ast-experimental, enum, forall-pointer-decay, master, pthread-emulation, qualifiedEnum
Children:
00f5fde
Parents:
a10f6b4
Message:

third attempt at specialized PRNG

Location:
libcfa/src/concurrency
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • libcfa/src/concurrency/kernel/startup.cfa

    ra10f6b4 r1959528  
    101101extern void __wake_proc(processor *);
    102102extern int cfa_main_returned;                                                   // from interpose.cfa
     103extern uint32_t __global_random_seed;
    103104
    104105//-----------------------------------------------------------------------------
     
    489490        preferred = ready_queue_new_preferred();
    490491        last_proc = 0p;
     492        random_state = __global_random_seed;
    491493        #if defined( __CFA_WITH_VERIFY__ )
    492494                canary = 0x0D15EA5E0D15EA5Ep;
  • libcfa/src/concurrency/thread.cfa

    ra10f6b4 r1959528  
    1010// Created On       : Tue Jan 17 12:27:26 2017
    1111// Last Modified By : Peter A. Buhr
    12 // Last Modified On : Mon Jan 10 17:05:35 2022
    13 // Update Count     : 28
     12// Last Modified On : Wed Jan 12 18:28:18 2022
     13// Update Count     : 35
    1414//
    1515
     
    2727uint64_t thread_rand();
    2828
    29 extern uint32_t __thread_seed;                                                  // global thread seed
     29extern uint32_t __global_random_seed;
    3030
    3131//-----------------------------------------------------------------------------
     
    4141        self_mon.owner = &this;
    4242        self_mon.recursion = 1;
    43         random_state = __thread_seed;
    4443        self_mon_p = &self_mon;
    4544        curr_cluster = &cl;
     
    4847        preferred = ready_queue_new_preferred();
    4948        last_proc = 0p;
     49        random_state = __global_random_seed;
    5050        #if defined( __CFA_WITH_VERIFY__ )
    5151                canary = 0x0D15EA5E0D15EA5Ep;
     
    180180        return ret;
    181181}
    182 
     182 
    183183#define GENERATOR LCG
    184184
    185 inline uint32_t MarsagliaXor( uint32_t & state ) {
     185static inline uint32_t MarsagliaXor( uint32_t & state ) {
     186        uint32_t ret = state;
    186187        state ^= state << 6;
    187188        state ^= state >> 21;
    188189        state ^= state << 7;
    189         return state;
     190        return ret;
    190191} // MarsagliaXor
    191192
    192 inline uint32_t LCG( uint32_t & state ) {                               // linear congruential generator
    193         return state = 36969 * (state & 65535) + (state >> 16); // 36969 is NOT prime!
     193static inline uint32_t LCG( uint32_t & state ) {                // linear congruential generator
     194        uint32_t ret = state;
     195        state = 36969 * (state & 65535) + (state >> 16);        // 36969 is NOT prime! No not change it!
     196        return ret;
    194197} // LCG
    195198
    196 void set_seed( uint32_t seed ) { active_thread()->random_state = seed; __thread_seed = seed; }
     199void set_seed( uint32_t seed ) { active_thread()->random_state = seed; __global_random_seed = seed; }
    197200uint32_t prng( void ) { return GENERATOR( active_thread()->random_state ); } // [0,UINT_MAX]
    198201
Note: See TracChangeset for help on using the changeset viewer.