More examples

examples/sieve.c

@@ -0,0 +1,38 @@
+#include <math.h>
+
+#define PFX bs
+#define SNAME bitset
+#include "cmc/bitset.h"
+
+// All the prime numbers up to 1 billion
+int main(void)
+{
+    const size_t n = 1000000000;
+
+    struct bitset *set = bs_new(n);
+
+    if (!set)
+        return 1;
+
+    bs_set_range(set, 0, n);
+
+    size_t up_to = sqrt((double)n) + 1;
+    for (size_t i = 2; i < up_to; ++i)
+    {
+        if (bs_get(set, i))
+        {
+            for (size_t j = 2 * i; j < n; j += i)
+                bs_clear(set, j);
+        }
+    }
+
+    FILE *file = fopen("primes.txt", "w");
+
+    for (size_t i = 2; i < n; i++)
+    {
+        if (bs_get(set, i))
+            fprintf(file, "%" PRIuMAX "\n", i);
+    }
+
+    fclose(file);
+}