Mandelbrot set: Difference between revisions
→Normalized Counting, Distance Estimation, Mercator Maps and Perturbation Theory: Small cuts and adjustments
(→{{header|Python}}: Duplicate code removed (see following subsection)) |
(→Normalized Counting, Distance Estimation, Mercator Maps and Perturbation Theory: Small cuts and adjustments) |
||
Line 6,521:
z = zero(c)
S = zeros(Complex{Float64}, n+1
for k in 1:n+1
S[k] = z
z = z ^ 2 + c
if abs2(z) > abs2(r)
error("
end
end
Line 6,536:
C = 5.0e-55 .* (A' .+ B .* im)
Z,
D = zeros(size(C))
Line 6,568:
z, dz = zero(c), zero(c)
S, dS = zeros(Complex{Float64}, n+1), zeros(Complex{Float64}, n+1
for k in 1:n+1
S[k], dS[k] = z, dz
z, dz = z ^ 2 + c, 2 * z * dz + 1
if abs2(z) > abs2(r)
error("
end
end
Line 6,583:
C = (.- 4.0) .* exp.((A' .+ B .* im) .* im)
Z,
D = zeros(size(C))
iteration(S,
for k in 1:n
M = abs2.(Z) .< abs2(r)
E[M], dE[M] = iteration(S[k],
Z[M], dZ[M] = S[k+1] .+ E[M], dS[k+1] .+ dE[M]
end
| |||