所以我试图以非常大的放大倍数计算出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