code: def run(n_bins):
      """
      Example of using openmp reduction to accelerate
      for loop
      """
  
      arr = np.zeros(n_bins)
  
      # omp parallel for reduction(+:arr)
      for i in range(10):
          for j in range(10):
              ij = i * 10 + j
              arr[int(ij % n_bins)] += 1
  
      return arr
  
label_desc: Use image node_id 2 for calling a global function with example usage: run(n_bins) and returns: arr