summaryrefslogtreecommitdiffhomepage
path: root/misc/examples/random.c
diff options
context:
space:
mode:
authorTyge Lovset <[email protected]>2022-12-20 23:31:51 +0100
committerTyge Lovset <[email protected]>2022-12-20 23:31:51 +0100
commit5f57d597cd27aef55adbcb3b452973b0c6e33667 (patch)
treedfd59c2fd0e36a6ef37912a9d0cc5a65970f1524 /misc/examples/random.c
parent1763be8c8cbbc0896477fcf924edd4180d1345a9 (diff)
downloadSTC-modified-5f57d597cd27aef55adbcb3b452973b0c6e33667.tar.gz
STC-modified-5f57d597cd27aef55adbcb3b452973b0c6e33667.zip
Restructured folders: examples, benchmarks, tests into misc folder.
Diffstat (limited to 'misc/examples/random.c')
-rw-r--r--misc/examples/random.c45
1 files changed, 45 insertions, 0 deletions
diff --git a/misc/examples/random.c b/misc/examples/random.c
new file mode 100644
index 00000000..fe64290d
--- /dev/null
+++ b/misc/examples/random.c
@@ -0,0 +1,45 @@
+#include <stdio.h>
+#include <time.h>
+#include <stc/crandom.h>
+
+int main()
+{
+ const size_t N = 1000000000;
+ const uint64_t seed = time(NULL), range = 1000000;
+ stc64_t rng = stc64_new(seed);
+
+ uint64_t sum;
+ clock_t diff, before;
+
+ printf("Compare speed of full and unbiased ranged random numbers...\n");
+ sum = 0;
+ before = clock();
+ c_forrange (N) {
+ sum += (uint32_t)stc64_rand(&rng);
+ }
+ diff = clock() - before;
+ printf("full range\t\t: %f secs, %" c_ZU ", avg: %f\n",
+ (float)diff / CLOCKS_PER_SEC, N, (double)sum / N);
+
+ stc64_uniform_t dist1 = stc64_uniform_new(0, range);
+ rng = stc64_new(seed);
+ sum = 0;
+ before = clock();
+ c_forrange (N) {
+ sum += stc64_uniform(&rng, &dist1); // unbiased
+ }
+ diff = clock() - before;
+ printf("unbiased 0-%" PRIu64 "\t: %f secs, %" c_ZU ", avg: %f\n",
+ range, (float)diff/CLOCKS_PER_SEC, N, (double)sum / N);
+
+ sum = 0;
+ rng = stc64_new(seed);
+ before = clock();
+ c_forrange (N) {
+ sum += stc64_rand(&rng) % (range + 1); // biased
+ }
+ diff = clock() - before;
+ printf("biased 0-%" PRIu64 " \t: %f secs, %" c_ZU ", avg: %f\n",
+ range, (float)diff / CLOCKS_PER_SEC, N, (double)sum / N);
+
+}