マンデルブロプロッターをもっと速くする方法は？

Question

私はPythonプログラムをより速くする方法を教えてください。このプログラムはマンデルブロセットを計算し、それをカメで描画します。私は問題がforループにあると思う。たぶん、ステップに時間がかかりすぎているかもしれません。

import numpy as np 
import turtle 
turtle.ht() 
turtle.pu() 
turtle.speed(0) 
turtle.delay(0) turtle.colormode(255) 

i= int(input("iteration = ")) 
g = int(input("accuracy = ")) 
xmin = float(input("X-min: ")) 
xmax = float(input("X-max: ")) 
ymin = float(input("Y-min: ")) 
ymax = float(input("Y-max: ")) 
cmode = int(255/i) 

input("PRESS TO START") 
for x in np.arange(xmin,xmax,1/g): 
    for y in np.arange(ymin,ymax,1/g): 
     c = x + y * 1j 
     z = 0 
     t = 1 
     for e in range(i): 
      z = z * z + c 
      if abs(z) > 3: 
       turtle.setx(g*c.real) 
       turtle.sety(g*c.imag) 
       turtle.dot(2,e*cmode,e*cmode,e*cmode) 
       t = 0 

     if t == 1: 
      turtle.setx(g*c.real) 
      turtle.sety(g*c.imag) 
      turtle.dot(2,"black") 

input("Calculated!") 
turtle.mainloop() 

Here is an example

Answer 1

次リワークは百倍速くあなたのオリジナルよりも次のようになります。

import numpy as np 
import turtle 

i = int(input("iteration = ")) 
g = int(input("accuracy = ")) 
xmin = float(input("X-min: ")) 
xmax = float(input("X-max: ")) 
ymin = float(input("Y-min: ")) 
ymax = float(input("Y-max: ")) 

cmode = int(255/i) 

input("PRESS TO START") 

turtle.hideturtle() 
turtle.penup() 
turtle.speed('fastest') 
turtle.colormode(255) 
turtle.setundobuffer(None) # turn off saving undo information 

turtle.tracer(0, 0) 

for x in np.arange(xmin, xmax, 1/g): 
    for y in np.arange(ymin, ymax, 1/g): 
     c = x + y * 1j 
     z = 0 
     t = True 

     for e in range(i): 
      z = z * z + c 

      if abs(z) > 3.0: 
       turtle.setposition(g * c.real, g * c.imag) 
       rgb = e * cmode 
       turtle.dot(2, rgb, rgb, rgb) 
       t = False 
       break 

     if t: 
      turtle.setposition(g * c.real, g * c.imag) 
      turtle.dot(2, "black") 

    turtle.update() 

print("Calculated!") 

turtle.mainloop() 

大きな変化が視覚的にすべてのプロットを回避するtracer()とupdate()の組み合わせの使用であります各垂直列が完成するたびに描画されます。