code: def test_sample(self):
      num_samples = 2**20
  
      for bw in [1e-1, 5e-1, 1]:
          hp_kernel = hp_kernels.Gaussian(
              data=self.x_train, bandwidth=bw, fix_boundary=True
          )
  
          samples = hp_kernel.sample(num_samples=num_samples)
          phat1, x = np.histogram(samples, normed=True)
          phat2 = hp_kernel((x[1:] + x[:-1]) / 2).mean(axis=0)
  
          for p1, p2 in zip(phat1, phat2):
              self.assertAlmostEqual(p1, p2, delta=5e-2)
  
imports: import os
  import unittest
  import numpy
  from scipy.integrate import quad
  from statsmodels.nonparametric import kernels
  from hpbandster.optimizers.kde import kernels
  import ConfigSpace
  from pdb import set_trace
label_desc: Use image node_id 6 for calling the TestGaussian obj's underlying member method code with example usage: obj.test_sample() without return types