prog8/examples/vm/bsieve.p8

41 lines
954 B
Plaintext
Raw Normal View History

2022-09-23 12:56:06 +00:00
%import textio
; The "Byte Sieve" test. https://en.wikipedia.org/wiki/Byte_Sieve
main {
sub start() {
const ubyte ITERS = 10
uword count
uword i
uword prime
uword k
const uword SIZEPL = 8191
uword @zp flags_ptr = memory("flags", SIZEPL, $100)
txt.print("calculating...\n")
repeat ITERS {
sys.memset(flags_ptr, SIZEPL, 1)
count = 0
for i in 0 to SIZEPL-1 {
if @(flags_ptr+i) {
prime = i + i + 3
k = i + prime
while k <= SIZEPL-1 {
@(flags_ptr + k) = false
k += prime
}
txt.print_uw(prime)
2022-09-24 14:00:25 +00:00
txt.nl()
2022-09-23 12:56:06 +00:00
count++
}
}
}
txt.nl()
2022-09-23 12:56:06 +00:00
txt.print_uw(count)
txt.print(" primes\n")
}
}