2020-09-09 20:49:32 +00:00
|
|
|
%target cx16
|
2020-09-20 21:49:36 +00:00
|
|
|
%import textio
|
|
|
|
%import floats
|
2020-08-29 21:55:26 +00:00
|
|
|
%zeropage basicsafe
|
|
|
|
|
|
|
|
main {
|
|
|
|
const uword width = 256
|
|
|
|
const uword height = 200
|
2020-08-31 21:00:53 +00:00
|
|
|
const ubyte max_iter = 16 ; 32 actually looks pretty nice but takes longer
|
2020-08-29 21:55:26 +00:00
|
|
|
|
|
|
|
sub start() {
|
2020-08-30 17:31:20 +00:00
|
|
|
initialize()
|
|
|
|
mandel()
|
|
|
|
repeat {
|
|
|
|
; do nothing
|
|
|
|
}
|
|
|
|
}
|
2020-08-29 21:55:26 +00:00
|
|
|
|
2020-08-30 17:31:20 +00:00
|
|
|
sub mandel() {
|
2020-08-30 21:50:31 +00:00
|
|
|
const float XL=-2.000
|
|
|
|
const float XU=0.500
|
|
|
|
const float YL=-1.100
|
|
|
|
const float YU=1.100
|
|
|
|
float dx = (XU-XL)/width
|
|
|
|
float dy = (YU-YL)/height
|
2020-08-29 21:55:26 +00:00
|
|
|
ubyte pixelx
|
|
|
|
ubyte pixely
|
|
|
|
|
|
|
|
for pixely in 0 to height-1 {
|
2020-08-30 22:36:40 +00:00
|
|
|
float yy = YL+dy*(pixely as float)
|
2020-08-29 21:55:26 +00:00
|
|
|
|
|
|
|
cx16.r0 = 0
|
|
|
|
cx16.r1 = pixely
|
|
|
|
cx16.FB_cursor_position()
|
|
|
|
|
|
|
|
for pixelx in 0 to width-1 {
|
2020-08-30 22:36:40 +00:00
|
|
|
float xx = XL+dx*(pixelx as float)
|
2020-08-29 21:55:26 +00:00
|
|
|
|
|
|
|
float xsquared = 0.0
|
|
|
|
float ysquared = 0.0
|
|
|
|
float x = 0.0
|
|
|
|
float y = 0.0
|
|
|
|
ubyte iter = 0
|
|
|
|
|
|
|
|
while iter<max_iter and xsquared+ysquared<4.0 {
|
|
|
|
y = x*y*2.0 + yy
|
|
|
|
x = xsquared - ysquared + xx
|
|
|
|
xsquared = x*x
|
|
|
|
ysquared = y*y
|
|
|
|
iter++
|
|
|
|
}
|
|
|
|
cx16.FB_set_pixel(max_iter-iter)
|
|
|
|
}
|
2020-08-30 17:31:20 +00:00
|
|
|
|
|
|
|
print_time()
|
2020-08-29 21:55:26 +00:00
|
|
|
}
|
|
|
|
}
|
2020-08-30 17:31:20 +00:00
|
|
|
|
|
|
|
sub initialize() {
|
|
|
|
void cx16.screen_set_mode($80)
|
|
|
|
|
|
|
|
txt.plot(32, 5)
|
|
|
|
txt.print("256*200")
|
|
|
|
txt.plot(32, 6)
|
|
|
|
txt.print("mandel-")
|
|
|
|
txt.plot(33, 7)
|
|
|
|
txt.print("brot")
|
|
|
|
txt.plot(32, 9)
|
|
|
|
txt.print("floats")
|
|
|
|
txt.plot(32, 10)
|
2020-08-30 21:50:31 +00:00
|
|
|
txt.print_b(max_iter)
|
|
|
|
txt.print(" iter")
|
2020-08-30 17:31:20 +00:00
|
|
|
|
|
|
|
cx16.r0 = 0
|
|
|
|
cx16.r1 = 0
|
|
|
|
cx16.r2 = 0
|
|
|
|
cx16.r3 = 0
|
|
|
|
cx16.clock_set_date_time()
|
|
|
|
|
|
|
|
cx16.r0=0
|
|
|
|
cx16.FB_init()
|
|
|
|
}
|
|
|
|
|
|
|
|
sub print_time() {
|
|
|
|
cx16.clock_get_date_time()
|
|
|
|
txt.plot(33, 12)
|
|
|
|
if lsb(cx16.r2) < 10
|
|
|
|
c64.CHROUT('0')
|
|
|
|
txt.print_ub(lsb(cx16.r2))
|
|
|
|
c64.CHROUT(':')
|
|
|
|
if msb(cx16.r2) < 10
|
|
|
|
c64.CHROUT('0')
|
|
|
|
txt.print_ub(msb(cx16.r2))
|
|
|
|
}
|
2020-08-29 21:55:26 +00:00
|
|
|
}
|