code: def compute_julia(cr, ci, N, bound=1.5, lim=1000.0, cutoff=1e6):
      """Pure Python calculation of the Julia set for a given `c`.  No NumPy
      array operations are used.
      """
      julia = np.empty((N, N), np.uint32)
      grid_x = np.linspace(-bound, bound, N)
      t0 = time()
      "omp parallel for default(none) shared(grid_x, cr, ci, lim, cutoff, julia)"
      for i, x in enumerate(grid_x):
          for j, y in enumerate(grid_x):
              julia[i, j] = kernel(x, y, cr, ci, lim, cutoff)
      return julia, time() - t0
  
imports: import numpy
  from time import time
label_desc: Use image node_id 2 for calling a global function with example usage: compute_julia(cr, ci, N, bound, lim, cutoff) and returns: julia,  unknown