Julia Set - 自己的 Mathematica 实现的问题答案

【问题标题】：Julia Set - problem with own Mathematica implementationJulia Set - 自己的 Mathematica 实现的问题
【发布时间】：2021-01-09 21:51:13
【问题描述】：

我一直在研究我自己的 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

【问题讨论】：

请edit您的问题包含您拥有的完整源代码作为minimal reproducible example，其他人可以对其进行测试。不要发布源代码的屏幕截图。

标签： wolfram-mathematica fractals mandelbrot

【解决方案1】：

这个

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]

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

【讨论】：