diff options
| author | Tyge Lovset <[email protected]> | 2022-12-20 23:31:51 +0100 |
|---|---|---|
| committer | Tyge Lovset <[email protected]> | 2022-12-20 23:31:51 +0100 |
| commit | 5f57d597cd27aef55adbcb3b452973b0c6e33667 (patch) | |
| tree | dfd59c2fd0e36a6ef37912a9d0cc5a65970f1524 /misc/examples/random.c | |
| parent | 1763be8c8cbbc0896477fcf924edd4180d1345a9 (diff) | |
| download | STC-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.c | 45 |
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); + +} |
