summaryrefslogtreecommitdiffhomepage
path: root/examples/random.c
diff options
context:
space:
mode:
authorTyge Løvset <[email protected]>2020-11-26 15:08:07 +0100
committerTyge Løvset <[email protected]>2020-11-26 15:08:07 +0100
commitbdea38e432d109ba416d7f980f124da285d69a4f (patch)
treeeffeb89e993b01f374d5a44ce0b859f979ed7f27 /examples/random.c
parent65ad066312ed549472814e1576cec2f0dc0f51ea (diff)
downloadSTC-modified-bdea38e432d109ba416d7f980f124da285d69a4f.tar.gz
STC-modified-bdea38e432d109ba416d7f980f124da285d69a4f.zip
Removed stc/cfmt.h, as _Generic requires C11. I have made it a gist: https://gist.github.com/tylov/bcc956a4779f1d14204e66a14f17beb9
All examples are reverted to use printf() instead of c_print().
Diffstat (limited to 'examples/random.c')
-rw-r--r--examples/random.c28
1 files changed, 14 insertions, 14 deletions
diff --git a/examples/random.c b/examples/random.c
index 5649a930..7f1be3f5 100644
--- a/examples/random.c
+++ b/examples/random.c
@@ -1,13 +1,13 @@
#include <stdio.h>
#include <time.h>
#include <stc/crandom.h>
-#include <stc/cfmt.h>
+#include <stc/cstr.h>
int main()
{
enum {R = 30};
const size_t N = 1000000000;
- clock_t difference, before;
+ clock_t difference, before;
uint64_t sum = 0;
uint64_t seed = time(NULL);
@@ -15,36 +15,36 @@ int main()
uint32_t range = crand_i32(&pcg) & ((1u << 28) - 1);
crand_uniform_i32_t dist0 = crand_uniform_i32_init(0, range);
- c_print(1, "32 uniform: {}\n", dist0.range);
- sum = 0;
+ printf("32 uniform: %u\n", dist0.range);
+ double fsum = 0;
before = clock();
c_forrange (N) {
- sum += crand_uniform_i32(&pcg, &dist0);
+ fsum += (double) crand_uniform_i32(&pcg, &dist0) / dist0.range;
}
difference = clock() - before;
- c_print(1, "{} {}: {} secs\n", N, (double) sum / N / dist0.range, (float) difference / CLOCKS_PER_SEC);
+ printf("%zu %f: %f secs\n", N, fsum / N, (float) difference / CLOCKS_PER_SEC);
pcg = crand_rng32_init(seed);
dist0 = crand_uniform_i32_init(0, range);
puts("32 unbiased");
- sum = 0;
+ fsum = 0;
before = clock();
c_forrange (N) {
- sum += crand_unbiased_i32(&pcg, &dist0);
+ fsum += (double) crand_unbiased_i32(&pcg, &dist0) / dist0.range;
}
difference = clock() - before;
- c_print(1, "{} {}: {} secs\n", N, (double) sum / N / dist0.range, (float) difference / CLOCKS_PER_SEC);
+ printf("%zu %f: %f secs\n", N, fsum / N, (float) difference / CLOCKS_PER_SEC);
puts("64 uniform");
crand_rng64_t stc = crand_rng64_init(seed);
- crand_uniform_i64_t dist1 = crand_uniform_i64_init(0, dist0.range);
+ crand_uniform_i64_t dist1 = crand_uniform_i64_init(0, N);
sum = 0;
- before = clock();
+ before = clock();
c_forrange (N) {
sum += crand_uniform_i64(&stc, &dist1);
}
difference = clock() - before;
- c_print(1, "{} {}: {} secs\n", N, (double) sum / N / dist1.range, (float) difference / CLOCKS_PER_SEC);
+ printf("%zu %f: %f secs\n", N, (double) sum / N, (float) difference / CLOCKS_PER_SEC);
puts("normal distribution");
@@ -57,11 +57,11 @@ int main()
sum += n;
if (n >= 0 && n < R) ++hist[n];
}
-
+
cstr_t bar = cstr_init();
c_forrange (i, int, R) {
cstr_take(&bar, cstr_with_size(hist[i] * 25ull * R / N2, '*'));
- c_print(1, "{:3} {}\n", i, bar.str);
+ printf("%2d %s\n", i, bar.str);
}
cstr_del(&bar);
} \ No newline at end of file