code: def test_identities2():
      # test the identity loggamma(z + 1) = log(z) + loggamma(z)
      x = np.array([-99.5, -9.5, -0.5, 0.5, 9.5, 99.5])
      y = x.copy()
      x, y = np.meshgrid(x, y)
      z = (x + 1j * y).flatten()
      dataset = np.vstack((z, np.log(z) + loggamma(z))).T
  
      def f(z):
          return loggamma(z + 1)
  
      FuncData(f, dataset, 0, 1, rtol=1e-14, atol=1e-14).check()
  
file_path: scipy/scipy/special/tests/test_loggamma.py
imports: import numpy
  from numpy.testing import assert_allclose, assert_
  from scipy.special._testutils import FuncData
  from scipy.special import gamma, gammaln, loggamma
label_desc: Use image node_id 2 for calling a global function with example usage: test_identities2() and returns: loggamma