code: def perm(data, p):
      n = len(data)
      out = np.empty(n)
      # pragma omp simd
      for i in range(n):
          out[i] = data[p[i]] * 2
      return out
  
label_desc: Use image node_id 1 for calling a global function with example usage: perm(data, p) and returns: out