def METHOD_NAME(multi,single): # range of sources or single source
  mean_single = single.mean()
  mean_multi = multi.mean()
  print("single source mean: %1.5g" % mean_single)
  print("multi source mean: %1.5g" % mean_multi)
  if np.allclose(mean_single, mean_multi): return True
  else:
    frac = mean_multi / mean_multi
    print("Means are off by a factor of %.6f" % frac)
    return False