Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed

  • tests/PRNG.cfa

    re6d8d11 rae7a085c  
    88// Created On       : Wed Dec 29 09:38:12 2021
    99// Last Modified By : Peter A. Buhr
    10 // Last Modified On : Sun Apr 23 22:02:09 2023
    11 // Update Count     : 420
     10// Last Modified On : Wed Dec 21 20:39:59 2022
     11// Update Count     : 406
    1212//
    1313
    … …  
    1515#include <stdlib.hfa>                                                                   // PRNG
    1616#include <clock.hfa>
     17#include <thread.hfa>
    1718#include <limits.hfa>                                                                   // MAX
    1819#include <math.hfa>                                                                             // sqrt
    1920#include <malloc.h>                                                                             // malloc_stats
    2021#include <locale.h>                                                                             // setlocale
    21 #include <thread.hfa>
    2222#include <mutex_stmt.hfa>
    2323
    24 #define xstr(s) str(s)
    25 #define str(s) #s
    26 
    27 #if defined( __x86_64__ ) || defined( __aarch64__ )             // 64-bit architecture
     24#ifdef __x86_64__                                                                               // 64-bit architecture
    2825#define PRNG PRNG64
    2926#else                                                                                                   // 32-bit architecture
    3027#define PRNG PRNG32
    3128#endif // __x86_64__
    32 
    33 //#define TIME
    3429
    3530#ifdef TIME                                                                                             // use -O2 -nodebug
    … …  
    4338#endif // TIME
    4439
    45 static void avgstd( size_t trials, size_t buckets[] ) {
    46         size_t min = MAX, max = 0;
     40void avgstd( unsigned int buckets[] ) {
     41        unsigned int min = MAX, max = 0;
    4742        double sum = 0.0, diff;
    4843        for ( i; BUCKETS ) {
    … …  
    5954        } // for
    6055        double std = sqrt( sum / BUCKETS );
    61         mutex( sout ) sout | "trials"  | trials | "buckets" | BUCKETS
     56        mutex( sout ) sout | "trials"  | TRIALS | "buckets" | BUCKETS
    6257                | "min" | min | "max" | max
    6358                | "avg" | wd(0,1, avg) | "std" | wd(0,1, std) | "rstd" | wd(0,1, (avg == 0 ? 0.0 : std / avg * 100)) | "%";
    … …  
    6964thread T1 {};
    7065void main( T1 & ) {
    71         size_t * buckets = calloc( BUCKETS );                           // too big for task stack
    72         for ( TRIALS / 50 ) {
     66        unsigned int * buckets = calloc( BUCKETS );                     // too big for task stack
     67        for ( TRIALS / 100 ) {
    7368                buckets[rand() % BUCKETS] += 1;                                 // concurrent
    7469        } // for
    75         avgstd( TRIALS / 50, buckets );
     70        avgstd( buckets );
    7671        free( buckets );
    7772} // main
    … …  
    8176        PRNG prng;
    8277        if ( seed != 0 ) set_seed( prng, seed );
    83         size_t * buckets = calloc( BUCKETS );                           // too big for task stack
     78        unsigned int * buckets = calloc( BUCKETS );                     // too big for task stack
    8479        for ( TRIALS ) {
    8580                buckets[prng( prng ) % BUCKETS] += 1;                   // concurrent
    8681        } // for
    87         avgstd( TRIALS, buckets );
     82        avgstd( buckets );
    8883        free( buckets );
    8984} // main
    … …  
    9186thread T3 {};
    9287void main( T3 & th ) {
    93         size_t * buckets = calloc( BUCKETS );                           // too big for task stack
    94         for ( TRIALS / 5 ) {
     88        unsigned int * buckets = calloc( BUCKETS );                     // too big for task stack
     89        for ( TRIALS ) {
    9590                buckets[prng() % BUCKETS] += 1;                                 // concurrent
    9691        } // for
    97         avgstd( TRIALS / 5, buckets );
     92        avgstd( buckets );
    9893        free( buckets );
    9994} // main
    … …  
    10196thread T4 {};
    10297void main( T4 & th ) {
    103         size_t * buckets = calloc( BUCKETS );                           // too big for task stack
     98        unsigned int * buckets = calloc( BUCKETS );                     // too big for task stack
    10499        for ( TRIALS ) {
    105                 buckets[prng( th ) % BUCKETS] += 1;                             // concurrent
    106         } // for
    107         avgstd( TRIALS, buckets );
     100                buckets[prng( th ) % BUCKETS] += 1;     // concurrent
     101        } // for
     102        avgstd( buckets );
    108103        free( buckets );
    109104} // main
    … …  
    113108static void dummy( thread$ & th ) __attribute__(( noinline ));
    114109static void dummy( thread$ & th ) {
    115         size_t * buckets = (size_t *)calloc( BUCKETS, sizeof(size_t) ); // too big for task stack
    116         for ( size_t i = 0; i < TRIALS; i += 1 ) {
     110        unsigned int * buckets = (unsigned int *)calloc( BUCKETS, sizeof(unsigned int) ); // too big for task stack
     111        for ( unsigned int i = 0; i < TRIALS; i += 1 ) {
    117112                buckets[prng( th ) % BUCKETS] += 1;                             // sequential
    118113        } // for
    119         avgstd( TRIALS, buckets );
     114        avgstd( buckets );
    120115        free( buckets );
    121116} // dummy
    … …  
    123118
    124119int main() {
    125         // setlocale( LC_NUMERIC, getenv( "LANG" ) );           // causes leaked storage message
    126 
    127         // only works on the current pthread thread
     120        // causes leaked storage message
     121        // setlocale( LC_NUMERIC, getenv( "LANG" ) );                   // print digit separator
    128122        // locale_t loc = newlocale( LC_NUMERIC_MASK, getenv( "LANG" ), (locale_t)0p );
    129123        // if ( loc == (locale_t)0p ) abort( "newlocale" );
    … …  
    132126        enum { TASKS = 4 };
    133127        Time start;
    134 
    135128#ifdef TIME                                                                                             // too slow for test and generates non-repeatable results
    136129#if 1
    137         sout | "glib rand" | nl | nl;
    138 
    139         size_t rseed;
     130        unsigned int rseed;
    140131        if ( seed != 0 ) rseed = seed;
    141132        else rseed = rdtscl();
    … …  
    143134
    144135        sout | sepDisable;
    145         sout | nl | wd(26, "rand()" ) | wd(12, "rand(5)") | wd(12, "rand(0,5)" );
     136        sout | wd(26, "rand()" ) | wd(12, "rand(5)") | wd(12, "rand(0,5)" );
    146137        for ( 20 ) {
    147138                sout | wd(26, rand()) | nonl;
    … …  
    155146        STARTTIME;
    156147        {
    157                 size_t * buckets = calloc( BUCKETS );                   // too big for task stack
    158                 for ( i; TRIALS / 5 ) {
     148                unsigned int * buckets = calloc( BUCKETS );             // too big for task stack
     149                for ( i; TRIALS / 10 ) {
    159150                        buckets[rand() % BUCKETS] += 1;                         // sequential
    160151                } // for
    161                 avgstd( TRIALS / 5, buckets );
     152                avgstd( buckets );
    162153                free( buckets );
    163154        }
    164         ENDTIME( " x 5 " );
     155        ENDTIME( " x 10 " );
    165156
    166157        sout | nl | "Concurrent";
    … …  
    172163                } // wait for threads to complete
    173164        }
    174         ENDTIME( " x 50 " );
     165        ENDTIME( " x 100 " );
    175166#endif // 0
    176167#endif // TIME
    177 
    178         sout | nl | "CFA " xstr(PRNG_NAME);
    179 
    180168#if 1
    181169        PRNG prng;
    … …  
    196184        STARTTIME;
    197185        {
    198                 size_t * buckets = calloc( BUCKETS );                   // too big for task stack
     186                unsigned int * buckets = calloc( BUCKETS );             // too big for task stack
    199187                for ( TRIALS ) {
    200188                        buckets[prng( prng ) % BUCKETS] += 1;           // sequential
    201189                } // for
    202                 avgstd( TRIALS, buckets );
     190                avgstd( buckets );
    203191                free( buckets );
    204192        }
    … …  
    231219        STARTTIME;
    232220        {
    233                 size_t * buckets = calloc( BUCKETS );                   // too big for task stack
    234                 for ( TRIALS / 5 ) {
     221                unsigned int * buckets = calloc( BUCKETS );             // too big for task stack
     222                for ( TRIALS ) {
    235223                        buckets[prng() % BUCKETS] += 1;
    236224                } // for
    237                 avgstd( TRIALS / 5, buckets );
     225                avgstd( buckets );
    238226                free( buckets );
    239227        }
    240         ENDTIME( " x 5 " );
     228        ENDTIME();
    241229
    242230        sout | nl | "Concurrent";
    … …  
    248236                } // wait for threads to complete
    249237        }
    250         ENDTIME( " x 5 " );
     238        ENDTIME();
    251239#endif // 0
    252240#if 1
Note: See TracChangeset for help on using the changeset viewer.