code: def test_sample_pymc():
      distribs_pymc = [
          Geometric("G", 0.5),
          Poisson("P", 1),
          NegativeBinomial("N", 5, 0.4),
      ]
      size = 3
      pymc = import_module("pymc")
      if not pymc:
          skip("PyMC is not installed. Abort tests for _sample_pymc.")
      else:
          for X in distribs_pymc:
              samps = sample(X, size=size, library="pymc")
              for sam in samps:
                  assert sam in X.pspace.domain.set
          raises(
              NotImplementedError,
              lambda: sample(Skellam("S", 1, 1), library="pymc"),
          )
  
file_path: sympy/sympy/stats/sampling/tests/test_sample_discrete_rv.py
imports: from sympy.core.singleton import S
  from sympy.core.symbol import Symbol
  from sympy.external import import_module
  from sympy.stats import Geometric, Poisson, Zeta, sample, Skellam, DiscreteRV, Logarithmic, NegativeBinomial, YuleSimon
  from sympy.testing.pytest import skip, raises, slow
label_desc: Use image node_id 3 for calling a global function with example usage: test_sample_pymc() without return types