2

仅用 100 次迭代来渲染这个 mandelbrot 集需要一个多小时，而 10,000 次迭代需要 10 个小时。有没有办法让它更快：

from graphics import *

width = 700
height = 700
win = GraphWin("Mandelbrot",width,height)
spacing = 1
zoom = 0.1
xOffset = -0.171

yOffset = 0.61
win.setBackground('black')
for x in range(0,width,spacing):
    for y in range(1,height,spacing):
        a = ((x / width) * zoom) - xOffset
        b = ((y / height) * zoom) - yOffset

        pt = Point(x,y)


        n = 0
        ca = a
        cb = b
        while(n<10000):
            aa = a * a - b * b
            bb = 2 * a * b
            a = aa + ca
            b = bb + cb
            n+=1
            if(abs(a+b) > 2000):
                break
            if(n < 2000):
               pt.setFill('black')
            if(n>5000):
                pt.setFill('grey')
            if(n>1000):
                pt.setFill('white')

        pt.draw(win)
4

2 回答 2

4

numpy可能是最快的方法。有关此方法的详细信息，请参阅“如何在 Python 中快速计算 Mandelbrot 集” 。

对于纯 Python，使用本机复数来加速循环。还可以使用abs()函数快速计算复数的大小：

>>> def mandle(c, boundary=2.0, maxloops=10000):
        # https://en.wikipedia.org/wiki/Mandelbrot_set
        z = 0.0j
        for i in range(maxloops):
            z = z * z + c
            if abs(z) > boundary:
                break
        return i

>>> mandle(0.04 + 0.65j)
21
>>> mandle(0.04 + 0.66j)
16
>>> mandle(0.04 + 0.67j)
12

渲染本身不太可能是程序的慢部分（10,000 个循环可能会使绘制点的时间相形见绌）。也就是说，如果您想加快渲染速度，通常唯一的选择是在每次调用图形库时绘制多个点。

最后，考虑您是否真的希望最大迭代次数为 10,000。最多 200 次迭代可以获得良好的结果。

于 2017-05-14T21:14:24.487 回答
0

除了@RaymondHettinger 建议的复数之外，我们还可以在 Zelle 图形方面做几件事来加快速度。首先是不使用Point()，它的开销太大。它自己的win实例有一个plot()位图操作的方法，它没有 的开销Point()，即不能取消绘制它，不能移动它。

第二个是关闭自动刷新并在每列上执行我们自己的屏幕刷新。最后，简单地避免做任何你不需要的计算——例如ca可以在外循环中计算，而不是在内循环中。颜色可以在最内层循环之外计算，等等。

这是我上面的返工——对 70 x 70 图像计时，它比原始代码快 7 倍：

from graphics import *

spacing = 1
zoom = 0.1

xOffset, yOffset = -0.171, 0.61

width, height = 700, 700

win = GraphWin('Mandelbrot', width, height, autoflush=False)

win.setBackground('black')

xzoom, yzoom = zoom / width, zoom / height

for real in range(0, width, spacing):

    ca = real * xzoom - xOffset

    for imaginary in range(0, height, spacing):

        c, z = complex(ca, imaginary * yzoom - yOffset), 0j

        n = 0

        while n < 10000:

            if abs(z) > 2000:
                break

            z = z * z + c

            n += 1

        color = 'black'

        if n > 5000:
            color = 'grey'
        elif n > 1000:
            color = 'white'

        if color != 'black':
            win.plot(real, imaginary, color=color)

    win.flush()

不是我们可能希望的完整数量级，但减少七个小时的周转时间仍然很重要！

最后，您的代码中存在一个错误，它使像素永远不会变灰——我已经修复了这个问题。

于 2017-05-15T00:10:37.327 回答