def julia_loops(cr, ci, N, bound=1.5, lim=1000., cutoff=1e6):
    julia = np.empty((N, N), dtype=np.uint32)
    grid_x = np.linspace(-bound, bound, N)
    for i, x in enumerate(grid_x):
        for j, y in enumerate(grid_x):
            julia[i,j] = kernel(x, y, cr, ci, lim, cutoff=cutoff)
    return julia
Pure Python calculation of the Julia set for a given `c`.  No NumPy
dataset/ETHPy150Open iskandr/parakeet/test/algorithms/test_julia.py/julia_loops