榕坚

你把问题文件传上来看看。还有一个可能就是sqrt函数。

柳烟

carr1407r.gsp (34.43 KB)

下载次数: 2199

2010-10-16 22:51

这是问题文件。sqrt是平方根函数。

榕坚

奇怪了，我今天做这个ahm.ufm系列中的分形。
PerpendicularMSet {
;This produces slice of the second order Mandelbulb directly perpendicular to the regular Mandelbrot set.
init:
Cx=Real(#pixel+@Perturbation)
Cy=Imag(#pixel+@Perturbation)
x=0
y=0
loop:
xTemp=x^2 - y^2 + Cx
y=-2*y*Abs(x) + Cy
x=xTemp
z=x+y*1i
bailout:
|x+y*(1i)|<=@Bailout
default:
center=(-0.5,0)
title="Perpendicular Mandelbrot"
float param Bailout
  default=4
endparam
switch:
type="PerpendicularJulia"
Start=#pixel
}
怎么做都是空白了一大块，逃逸区的部分全部数据溢出，莫名其妙的：

榕坚

还是弄不明白什么原因，今天下午再做一遍就全对了。柳老师的问题文件还没查出原因来，对数据溢出的原因还是没能完全明白：

柳烟

榕兄，这个M集与精典的M集的造法完全一样呀！就是与精典M集的迭代坐标变了点，你再造一遍，定能成功。

下载 (19.27 KB)

2010-10-17 15:50

柳烟

感谢你对问题文件的查找，不慌在一时，我暂时放开对那文件的关注，休息一下再来研究。
放大楼上的文件。

榕坚

柳老师，找到原因了：C点也要拉回来。你只拉回了Z点，不妨再试试。

柳烟

148# 榕坚
问题出在没将c拉回。我看了你的原文件，按你说的问题症结，重新作后，对了，谢谢。

榕坚

奇怪的UF代码：
arnewurzel {
; Arne Richter
init:
bool done = false
b = #pixel

loop:
a = @radikand / b
#z = (a+b) / 2
If abs (b-a) < @epsilon
done = true
   else
b = (a+b) / @quotient ; Parameter added to avoid all points to be
outside
endif

bailout:
done == false

default:
param @epsilon
default=(0.1,0)
endparam

param @quotient
default=(1,1)
endparam

param @radikand
default=(9,0)
endparam
}
其中红色部分的abs函数，按UF规定：abs(Z)还是复数，那么两个复数如何比较大小呢？

柳烟

帮助里没有，红色处确实费解。