code: def pi_estimate(DARTS, rand, randsize):
      hits = 0
      "omp parallel for reduction(+:hits)"
      for i in range(0, DARTS):
          x = rand[i % randsize]
          y = rand[(randsize - i) % randsize]
          dist = sqrt(pow(x, 2) + pow(y, 2))
          if dist <= 1.0:
              hits += 1.0
      # hits / throws = 1/4 Pi
      pi = 4 * (hits / DARTS)
      return pi
  
label_desc: Use image node_id 1 for calling a global function with example usage: pi_estimate(DARTS, rand, randsize) and returns: pi