Changeset aa8e24c3 for libcfa/src/stdlib.hfa

Timestamp:

Dec 29, 2021, 5:06:07 PM (2 years ago)

Author:

Peter A. Buhr <pabuhr@…>

Branches:

ADT, ast-experimental, enum, forall-pointer-decay, master, pthread-emulation, qualifiedEnum

Children:

aa0a1ad, db1ebed

Parents:

7a6ae53

git-author:

Peter A. Buhr <pabuhr@…> (12/29/21 17:00:15)

git-committer:

Peter A. Buhr <pabuhr@…> (12/29/21 17:06:07)

Message:

change range of integral random( u, l ) to [l,u], first attempt at specialized PRNG

File:

• libcfa/src/stdlib.hfa (modified) (3 diffs)

libcfa/src/stdlib.hfa

@@ -10 +10 @@
 // Created On       : Thu Jan 28 17:12:35 2016
 // Last Modified By : Peter A. Buhr
-// Last Modified On : Tue Apr 20 21:20:03 2021
-// Update Count     : 575
+// Last Modified On : Wed Dec 29 15:30:58 2021
+// Update Count     : 591
 //
 
@@ -363 +363 @@
 
 static inline {
-        long int random( long int l, long int u ) { if ( u < l ) [u, l] = [l, u]; return lrand48() % (u - l) + l; } // [l,u)
-        long int random( long int u ) { if ( u < 0 ) return random( u, 0 ); else return random( 0, u ); } // [0,u)
+        long int random( long int l, long int u ) { if ( u < l ) [u, l] = [l, u]; return lrand48() % (u - l + 1) + l; } // [l,u]
+        long int random( long int u ) { return random( 0, u - 1 ); } // [0,u)
         unsigned long int random( void ) { return lrand48(); }
         unsigned long int random( unsigned long int u ) { return lrand48() % u; } // [0,u)
-        unsigned long int random( unsigned long int l, unsigned long int u ) { if ( u < l ) [u, l] = [l, u]; return lrand48() % (u - l) + l; } // [l,u)
+        unsigned long int random( unsigned long int l, unsigned long int u ) { if ( u < l ) [u, l] = [l, u]; return lrand48() % (u - l + 1) + l; } // [l,u]
 
         char random( void ) { return (unsigned long int)random(); }
@@ -388 +388 @@
 //---------------------------------------
 
+struct PRNG {
+        uint32_t callcnt;                                                                       // call count
+        uint32_t seed;                                                                          // current seed
+        uint32_t state;                                                                         // random state
+}; // PRNG
+
+extern uint32_t prng( PRNG & prng ) __attribute__(( warn_unused_result )); // [0,UINT_MAX]
+static inline {
+        void set_seed( PRNG & prng, uint32_t seed_ ) with( prng ) { state = seed = seed_; } // set seed
+        void ?{}( PRNG & prng ) { set_seed( prng, rdtscl() ); } // random seed
+        void ?{}( PRNG & prng, uint32_t seed ) { set_seed( prng, seed ); } // fixed seed
+        uint32_t get_seed( PRNG & prng ) __attribute__(( warn_unused_result )) with( prng ) { return seed; } // get seed
+        uint32_t prng( PRNG & prng, uint32_t u ) __attribute__(( warn_unused_result )) { return prng( prng ) % u; } // [0,u)
+        uint32_t prng( PRNG & prng, uint32_t l, uint32_t u ) __attribute__(( warn_unused_result )) { return prng( prng, u - l + 1 ) + l; } // [l,u]
+        uint32_t calls( PRNG & prng ) __attribute__(( warn_unused_result )) with( prng ) { return callcnt; }
+} // distribution
+
+extern void set_seed( uint32_t seed );                                  // set per thread seed
+extern uint32_t get_seed();                                                             // get seed
+extern uint32_t prng( void ) __attribute__(( warn_unused_result )); // [0,UINT_MAX]
+static inline {
+        uint32_t prng( uint32_t u ) __attribute__(( warn_unused_result ));
+        uint32_t prng( uint32_t u ) { return prng() % u; }      // [0,u)
+        uint32_t prng( uint32_t l, uint32_t u ) __attribute__(( warn_unused_result ));
+        uint32_t prng( uint32_t l, uint32_t u ) { return prng( u - l + 1 ) + l; } // [l,u]
+} // distribution
+
+//---------------------------------------
+
 extern bool threading_enabled( void ) OPTIONAL_THREAD;