Skip to content

Instantly share code, notes, and snippets.

import collections
import mpmath as mp
import matplotlib.pyplot as plt
mp.dps = 100
alpha = beta = 2.2
sq2 = mp.sqrt(2)
gamma = (alpha * beta + mp.sqrt(alpha * alpha * beta * beta + 8 - 4 * (alpha * alpha+beta*beta)) )/ 2
@neozhaoliang
neozhaoliang / taichi_circle_packing.py
Last active October 6, 2023 08:01
Given an input image, convert it to a circle packing pattern
import cv2
import cairocffi as cairo
import numpy as np
import taichi as ti
ti.init(arch=ti.cpu)
scale = 5
@neozhaoliang
neozhaoliang / return_probability.py
Created November 29, 2022 10:28
Random walk return probability
# pip install taichi
import taichi as ti
ti.init(arch=ti.gpu)
d = 3
num_rounds = 100000
max_steps = 1000000
ivec = ti.types.vector(d, int)