Mandelbrot set: Difference between revisions

m
Line 8,977:
 
E, Z, dZ = np.zeros_like(C), np.zeros_like(C), np.zeros_like(C)
D, I, J = np.zeros(C.shape), np.zeros(C.shape, dtype=np.int64), np.zeros(C.shape, dtype=np.int64)
 
for k in range(n): # you can also try range(2*n) due to rebasing
R = Z.real ** 2 + Z.imag ** 2 < E.real ** 2 + E.imag ** 2
E[R], I[R] = Z[R], J[R]0 # rebase if z is closer to zero
M = Z.real ** 2 + Z.imag ** 2 < r ** 2
E[M], I[M] = (2 * S[I[M]] + E[M]) * E[M] + C[M], I[M] + 1
305

edits