Function run_sim

code/OO/collision_matrix.py:47–59

(N, n_iter)

    return sep < sum_r

def run_sim(N, n_iter):
  box_x, box_y = 10., 6.
  r_min, r_max = 0.1, 0.3
  print(f'n balls={N}')
  for i in range(n_iter):
    x = r_max + (box_x-2*r_max)*rand(N)
    y = r_max + (box_y-2*r_max)*rand(N)
    R = 0.1 + 0.8*rand(N)
    T_s = time.time()
    coll = collision_matrix(x, y, R)
    n_coll = (np.sum(coll) - N)//2
    print(f'n coll={n_coll} Hz = '
          f'{1/(time.time()-T_s):.2f}')