Skip to content

Instantly share code, notes, and snippets.

@rexwangcc
Last active August 30, 2021 01:09
Show Gist options
  • Save rexwangcc/80eca823500f932be64e4ec37a0879dd to your computer and use it in GitHub Desktop.
Save rexwangcc/80eca823500f932be64e4ec37a0879dd to your computer and use it in GitHub Desktop.
A Julia Set fractal implementation in Taichi language. Originally created by @yuanming-hu.
import taichi as ti
ti.init(arch=ti.gpu)
n = 320
pixels = ti.field(dtype=float, shape=(n * 2, n))
@ti.func
def complex_sqr(z):
return ti.Vector([z[0]**2 - z[1]**2, z[1] * z[0] * 2])
@ti.kernel
def paint(t: float):
for i, j in pixels: # Parallized over all pixels
c = ti.Vector([-0.8, ti.cos(t) * 0.2])
z = ti.Vector([i / n - 1, j / n - 0.5]) * 2
iterations = 0
while z.norm() < 20 and iterations < 50:
z = complex_sqr(z) + c
iterations += 1
pixels[i, j] = 1 - iterations * 0.02
gui = ti.GUI("Julia Set", res=(n * 2, n))
for i in range(1000000):
paint(i * 0.03)
gui.set_image(pixels)
gui.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment