aboutsummaryrefslogtreecommitdiff
path: root/docs/history/15-o3/code/examples/primes.zc
diff options
context:
space:
mode:
Diffstat (limited to 'docs/history/15-o3/code/examples/primes.zc')
-rw-r--r--docs/history/15-o3/code/examples/primes.zc38
1 files changed, 0 insertions, 38 deletions
diff --git a/docs/history/15-o3/code/examples/primes.zc b/docs/history/15-o3/code/examples/primes.zc
deleted file mode 100644
index 26df0b5..0000000
--- a/docs/history/15-o3/code/examples/primes.zc
+++ /dev/null
@@ -1,38 +0,0 @@
-// primes below 60000: an o3 sieve over a fixed array, the primes
-// gathered into a growable one, then a few questions asked of them.
-// with ZYTRACE=1 the o3 lines show each region being lifted and compiled.
-
-// marks the composites; f[i] ends up 1 for every prime i
-o3 int sieve(fixed []int f)
-{
- int n = len(f), c = 0;
- for (int i = 2; i < n; i++) f[i] = 1;
- for (int i = 2; i * i < n; i++)
- if (f[i]) for (int j = i * i; j < n; j += i) f[j] = 0;
- for (int i = 2; i < n; i++) c += f[i];
- return c;
-}
-// the largest gap between neighbours, and where it starts
-o3 int gap(fixed []int f, *int at)
-{
- int n = len(f), last = 2, best = 0, where = 0;
- for (int i = 3; i < n; i++)
- if (f[i]) { if (i - last > best) { best = i - last; where = last; } last = i; }
- *at = where;
- return best;
-}
-fn main()
-{
- fixed []int f = new(fixed []int, 60000);
- []int ps = new([]int);
- int c = sieve(f), at = 0, twins = 0;
- for (int i = 2; i < len(f); i++) if (f[i]) push(ps, i);
- o3 { for (int k = 1; k < len(ps); k++) if (ps[k] - ps[k - 1] == 2) twins++; }
- puts("primes below 60000: "); putn(c); putc(10);
- puts("the last five:"); for (int k = len(ps) - 5; k < len(ps); k++) { putc(' '); putn(ps[k]); } putc(10);
- puts("twin pairs: "); putn(twins); putc(10);
- int g = gap(f, &at);
- puts("widest gap: "); putn(g); puts(", after "); putn(at); putc(10);
- free(d*f); free(d*ps);
- return 0;
-}