code: def random_geometric(shape, rng):
      P = rng.integers(1, 1000, size=(shape[0], 2))
      Q = rng.integers(1, 1000, size=(shape[1], 2))
      return scipy.sparse.csr_matrix(cdist(P, Q, "sqeuclidean"))
  
file_path: scipy/benchmarks/benchmarks/sparse_csgraph_matching.py
imports: import numpy
  import scipy.sparse
  from scipy.spatial.distance import cdist
  from .common import Benchmark, safe_import
label_desc: Use image node_id 4 for calling a global function with example usage: random_geometric(shape, rng) and returns: scipy