prog8/examples/vm/bsieve.p8
2024-12-15 22:54:06 +01:00

44 lines
968 B
Lua

; The "Byte Sieve" test. https://en.wikipedia.org/wiki/Byte_Sieve
%import textio
%zeropage basicsafe
%option no_sysinit
main {
sub start() {
uword count
uword i
uword prime
uword k
const uword SIZEPL = 8191
uword @zp flags_ptr = memory("flags", SIZEPL, $100)
txt.print("calculating...\n")
sys.memset(flags_ptr, SIZEPL, 1)
count = 1
for i in 0 to SIZEPL-1 {
if flags_ptr[i]!=0 {
prime = i + i + 3
k = i + prime
while k <= SIZEPL-1 {
flags_ptr[k] = 0 ; false
k += prime
}
; txt.print_uw(prime)
; txt.nl()
count++
}
}
txt.nl()
txt.print("last prime: ")
txt.print_uw(prime)
txt.print("\nnumber of primes: ")
txt.print_uw(count)
txt.nl()
}
}