dmj-NovaIJulia {
;
; This is the Nova fractal (Julia form), a
; modified Newtonian-style fractal. The formula
; was first shown to me by Paul Derbyshire (who
; named it Nova). It has also appeared elsewhere
; under other names. If you leave the Julia
; seed at the default (0,0), you can use this as
; a general Newton-style fractal as in FractInt.
;
; This variant uses Kerry Mitchell's inverted
; computation so that it works with coloring
; methods expecting divergent z.
;
init:
complex zsquared = (0,0)
complex zcubed = (0,0)
complex zcurrent = #pixel
z = (0,0)
loop:
IF (@power == (3,0)); special optimized routine for power 3
zsquared = sqr(zcurrent)
zcubed = zsquared * zcurrent
z = @relax * (zcubed-1) / (3*zsquared) - @seed
zcurrent = zcurrent - z
z = 1/z
ELSE
z = @relax * (zcurrent^@power-1) / (@power * zcurrent^(@power-1)) - @seed
zcurrent = zcurrent - z
z = 1/z
ENDIF
IF (@fudge == true &&\
|z| > 4); fudging angle
z = z * zcurrent/cabs(zcurrent)
ENDIF
bailout:
|z| < @bailout
default:
title = "Nova-I (Julia)"
helpfile = "dmj-pub\dmj-pub-uf-ni.htm"
maxiter = 1000
periodicity = 0
center = (0,0)
magn = 1.5
param seed
caption = "Julia Seed"
default = (0,0)
hint = "This is the Julia seed, a constant parameter which \
defines the shape of the fractal."
endparam
param power
caption = "Exponent"
default = (3,0)
hint = "Overall exponent for the equation. (3,0) gives \
the classic NovaM type."
endparam
param bailout
caption = "Bailout"
default = 10000.0
hint = "Bailout value; larger values will cause more \
iterations to be done for each point."
endparam
param relax
caption = "Relaxation"
default = (1,0)
hint = "This can be used to slow down the convergence of \
the formula."
endparam
param fudge
caption = "Fudge z Angle"
default = false
hint = "Modifies angle of z based on starting point. \
Turning this on will make decomposition more \
consistent between iterations, regardless of \
the root converged on."
endparam
switch:
type = "dmj-NovaIMandel"
power = @power
bailout = @bailout
relax = @relax
fudge = @fudge
}
位于dmj.ufm中,代码看似不长,算起来长。数据我已弄好,可就是扫出的没多少东西,不对劲。不知是不是fudge这个变量我没用对,实不知该如何用。
Nova-I (Julia).gsp (65.23 KB)
|