Python中的快速任意浮点

Question

所以我试着在非常大的放大率下计算曼德布罗特集，遇到了一堵砖墙。你知道吗

使用pythons自己的float和jit给了我足够的速度。你知道吗

但是在切换到使用十进制之后，我的函数开始为所有字段返回0。关闭Numba的jit解决了这个问题，但是我的软件的性能急剧下降到无法使用（在一个浮点精度为6的64x64字段上，100次迭代大约30秒）。你知道吗

我的代码：

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