code: def _check_multigammaln_array_result(a, d):
      # Test that the shape of the array returned by multigammaln
      # matches the input shape, and that all the values match
      # the value computed when multigammaln is called with a scalar.
      result = multigammaln(a, d)
      assert_array_equal(a.shape, result.shape)
      a1 = a.ravel()
      result1 = result.ravel()
      for i in range(a.size):
          assert_array_almost_equal_nulp(
              result1[i], multigammaln(a1[i], d)
          )
  
file_path: scipy/scipy/special/tests/test_spfun_stats.py
imports: import numpy
  from numpy.testing import assert_array_equal, assert_array_almost_equal_nulp, assert_almost_equal
  from pytest import raises
  from scipy.special import gammaln, multigammaln
label_desc: Use image node_id 1 for calling a global function with example usage: _check_multigammaln_array_result(a, d) without return types