Python中的快速,任意浮点

时间:2017-09-12 13:03:24

标签: python precision jit mandelbrot

所以我试图以非常大的放大倍数计算出mandelbrot的设置,并且遇到了一些砖墙。

使用pythons自己的浮动与jit配对给了我足够的速度。

但是在切换到使用十进制之后,我的函数刚开始为所有字段返回0。关闭Numba的jit解决了这个问题,但是我的软件的性能已经大大降低到无法使用(在64x64字段上进行100次迭代约30秒,浮点精度为6)。

我的代码:

def mandelbrot_point(creal, cimag, maxiter):
    real=creal
    imag=cimag
    for n in range(maxiter):
        real2=real*real
        imag2=imag*imag
        if real2+imag2>4.0:
            return n
        imag=2*real*imag+cimag
        real=real2-imag2+creal
    return 0

0 个答案:

没有答案