code: def testSampleAndPdfConsistency4d(self):
      mean_dir = tf.math.l2_normalize(
          [[1.0, 2, 3, 4], [-2, -3, -1, 0]], axis=-1
      )
      concentration = [[0], [1e-4], [0.1], [1], [4]]
      pspherical = power_spherical.PowerSpherical(
          mean_direction=mean_dir,
          concentration=concentration,
          validate_args=True,
          allow_nan_stats=False,
      )
      self.VerifySampleAndPdfConsistency(pspherical)
      self.VerifyPdfWithNumpy(pspherical)
  
imports: import numpy
  from scipy import special
  import tensorflow.compat.v2
  from tensorflow_probability.python.distributions import kullback_leibler
  from tensorflow_probability.python.distributions import power_spherical
  from tensorflow_probability.python.distributions import spherical_uniform
  from tensorflow_probability.python.distributions import von_mises_fisher
  from tensorflow_probability.python.internal import dtype_util
  from tensorflow_probability.python.internal import tensorshape_util
  from tensorflow_probability.python.internal import test_util
  from tensorflow_probability.python.stats import sample_stats
label_desc: Use image node_id 10 for calling the _PowerSphericalTest obj's underlying member method code with example usage: obj.testSampleAndPdfConsistency4d() without return types