code: def test_logpdf_overflow():
      # regression test for gh-12988; testing against linalg instability for
      # very high dimensionality kde
      np.random.seed(1)
      n_dimensions = 2500
      n_samples = 5000
      xn = np.array(
          [
              np.random.randn(n_samples) + (n)
              for n in range(0, n_dimensions)
          ]
      )
  
      # Default
      gkde = stats.gaussian_kde(xn)
  
      logpdf = gkde.logpdf(np.arange(0, n_dimensions))
      np.testing.assert_equal(np.isneginf(logpdf[0]), False)
      np.testing.assert_equal(np.isnan(logpdf[0]), False)
  
imports: from scipy import stats, linalg, integrate
  import numpy
  from numpy.testing import assert_almost_equal, assert_, assert_equal, assert_array_almost_equal, assert_array_almost_equal_nulp, assert_allclose
  import pytest
  from pytest import raises
label_desc: Use image node_id 18 for calling a global function with example usage: test_logpdf_overflow() without return types