Julia Set - 自己的 Mathematica 实现问题

问题描述

我一直在研究我自己的 Julia Set Plot 实现。我不想使用 JuliaSetPlot,(不过我很想使用 JuliaSetIterationPoints 和 JuliaSetCount,我只是不知道如何使用)。

我想出了这样的方法,但我有一个问题,我不知道哪里出了问题,为什么它不起作用。

有人可以帮忙吗?

'''数学

firstFun= Function[ {Typed[pixel0,"ComplexReal64"]},Module[{i = 1,maksi=100,pixel = pixel0},While[i < maksi && (Abs[pixel])^2 < 2,temp = (Re[pixel])^2 - (Im[pixel])^2
    Re[pixel] = 2 * Re[pixel] * Im[pixel] - 0.8\[Iota] * Im[pixel0]
    Im[pixel] = temp - 0.8\[Iota]* Re[pixel0];
     i++ ];
    i]];  

''' my code

解决方法

这个

firstFun=Function[{Typed[pixel0,"ComplexReal64"]},Module[{i=1,maksi=100,pixel=pixel0},While[i<maksi&&Abs[pixel]^2<2,pixel=2*Re[pixel]*Im[pixel]-0.8*I*Im[pixel0]+
       I*(Re[pixel]^2-Im[pixel]^2-0.8*I*Re[pixel0]);
      i++];
    i]];
compFun[c_]=FunctionCompile[firstFun]

编译时没有任何编译时错误消息。

如果我没有弄错那么我认为你的像素计算可以简化为

pixel=I*Conjugate[pixel]^2+0.8*Conjugate[pixel0]

请非常仔细地测试所有这些以确保它是正确的。