code: def test_gamma(x: float, shape: float, scale: float) -> None:
      gpf = logdensities.gamma(x, shape, scale).numpy()
      sps = scipy.stats.gamma(a=shape, loc=0.0, scale=scale).logpdf(x)
      np.testing.assert_allclose(gpf, sps)
  
file_path: GPflow/tests/gpflow/test_logdensities.py
imports: import numpy
  import pytest
  import scipy.stats
  from numpy.random import randn
  from numpy.testing import assert_allclose
  from scipy.stats import multivariate_normal
  from gpflow import logdensities
  from gpflow.base import AnyNDArray
  from gpflow.utilities import to_default_float
label_desc: Use image node_id 6 for calling a global function with example usage: test_gamma(x, shape, scale) without return types