Minors on crandom and example.

1 files changed, 25 insertions, 17 deletions

diff --git a/examples/random.c b/examples/random.c
index 4834a974..0cf9da6f 100644
--- a/examples/random.c
+++ b/examples/random.c
@@ -1,43 +1,51 @@
 #include <stdio.h>

 #include <time.h>

+#include <math.h>

 #include <stc/crandom.h>

 #include <stc/cstr.h>

 

 int main()

 {

     enum {R = 30};

-    const size_t N = 100000000;

-    clock_t diff, before;    

-    uint64_t sum = 0;

+    const size_t N = 500000000;

     uint64_t seed = time(NULL);

-

     cstc64_t stc = cstc64_init(seed);

-    cstc64_uniform_t dist1 = cstc64_uniform_init(0, N);

-    sum = 0;

-    before = clock();    

-    c_forrange (N)  {

-        sum += cstc64_uniform(&stc, &dist1);

-    }

-    diff = clock() - before;

-    printf("uniform: %zu %f: %f secs\n", N, (double) sum / N, (float) diff / CLOCKS_PER_SEC);

+

+    uint64_t sum = 0;

 

     cstc64_normalf_t dist2 = cstc64_normalf_init(R / 2.0, R / 6.0);

     size_t N2 = 10000000;

     int hist[R] = {0};

     sum = 0;

-    before = clock();

     c_forrange (N2)  {

-        int n = (int) (cstc64_normalf(&stc, &dist2) + 0.5);

+        int n = round((cstc64_normalf(&stc, &dist2) + 0.5));

         sum += n;

         if (n >= 0 && n < R) ++hist[n];

     }

-    diff = clock() - before;

-    printf("normal : %zu %f: %f secs\n", N, (double) sum / N2, (float) diff / CLOCKS_PER_SEC);

-

     cstr_t bar = cstr_inits;

     c_forrange (i, int, R)  {

         cstr_resize(&bar, hist[i] * 25ull * R / N2, '*');

         printf("%3d %s\n", i, bar.str);

     }

+

+    clock_t diff, before;    

+

+    sum = 0;

+    before = clock();    

+    c_forrange (N)  {

+        sum += cstc64_rand(&stc);

+    }

+    diff = clock() - before;

+    printf("random : %f secs, %zu %f\n", (float) diff / CLOCKS_PER_SEC, N, (double) sum / N);

+

+    cstc64_uniform_t dist1 = cstc64_uniform_init(0, 1000);

+    sum = 0;

+    before = clock();    

+    c_forrange (N)  {

+        sum += cstc64_uniform(&stc, &dist1);

+    }

+    diff = clock() - before;

+    printf("uniform: %f secs, %zu %f\n", (float) diff / CLOCKS_PER_SEC, N, (double) sum / N);

+

     cstr_del(&bar);

 }
\ No newline at end of file