code: def simulator_moment(op, rv, *inputs):
      sim_inputs = inputs[3:]
      # Take the mean of 10 draws
      multiple_sim = rv.owner.op(
          *sim_inputs, size=pt.concatenate([[10], rv.shape])
      )
      return pt.mean(multiple_sim, axis=0)
  
imports: import logging
  import numpy
  import pytensor
  import pytensor.tensor
  from pytensor.graph.op import Apply, Op
  from pytensor.tensor.random.op import RandomVariable
  from pytensor.tensor.variable import TensorVariable
  from scipy.spatial import cKDTree
  from pymc.distributions.distribution import Distribution, _moment
  from pymc.logprob.abstract import _logprob
  from pymc.pytensorf import floatX
label_desc: Use image node_id 1 for calling a global function with example usage: simulator_moment(op, rv) and returns: pt