2019-06-26 05:26:20 +01:00
|
|
|
# Compute the Mandelbrot set, to test complex numbers
|
|
|
|
|
2020-03-23 02:26:08 +00:00
|
|
|
|
2019-06-26 05:26:20 +01:00
|
|
|
def mandelbrot(w, h):
|
|
|
|
def in_set(c):
|
|
|
|
z = 0
|
|
|
|
for i in range(32):
|
|
|
|
z = z * z + c
|
|
|
|
if abs(z) > 100:
|
|
|
|
return i
|
|
|
|
return 0
|
|
|
|
|
|
|
|
img = bytearray(w * h)
|
|
|
|
|
2020-03-23 02:26:08 +00:00
|
|
|
xscale = (w - 1) / 2.4
|
|
|
|
yscale = (h - 1) / 3.2
|
2019-06-26 05:26:20 +01:00
|
|
|
for v in range(h):
|
2020-03-23 02:26:08 +00:00
|
|
|
line = memoryview(img)[v * w : v * w + w]
|
2019-06-26 05:26:20 +01:00
|
|
|
for u in range(w):
|
|
|
|
c = in_set(complex(v / yscale - 2.3, u / xscale - 1.2))
|
|
|
|
line[u] = c
|
|
|
|
|
|
|
|
return img
|
|
|
|
|
2020-03-23 02:26:08 +00:00
|
|
|
|
2019-06-26 05:26:20 +01:00
|
|
|
bm_params = {
|
|
|
|
(100, 100): (20, 20),
|
|
|
|
(1000, 1000): (80, 80),
|
|
|
|
(5000, 1000): (150, 150),
|
|
|
|
}
|
|
|
|
|
2020-03-23 02:26:08 +00:00
|
|
|
|
2019-06-26 05:26:20 +01:00
|
|
|
def bm_setup(ps):
|
|
|
|
return lambda: mandelbrot(ps[0], ps[1]), lambda: (ps[0] * ps[1], None)
|