ジュリア集合をnumpyで描いてみる(そのl)

次のページを参考にする。 Pythonでジュリア集合を描画!条件変更で動画にしてみた

ほぼ書いているとおりですが、numpyがwarningを出すのでちょっと変更。 やっているのは「Zの初期化を@jit の関数の外に移す」とか「定数初期化をループの外に出す」とかそれぐらい。

• julia_cpu_2.py

jitで実行しているので0.67秒ですが、@jitを外すと15秒ぐらいなのでこれでも20倍速か。

GPGPU版は、時間切れなので宿題です。