code: def _random_pd_matrix(self, *shape):
      mat = rng.rand(*shape)
      chol = transform_diagonal.TransformDiagonal(softplus.Softplus())(
          mat
      )
      chol = tf.linalg.band_part(chol, -1, 0)
      return self.evaluate(tf.matmul(chol, chol, adjoint_b=True))
  
imports: import numpy
  from scipy import stats
  import tensorflow.compat.v2
  from tensorflow_probability.python.bijectors import softplus
  from tensorflow_probability.python.bijectors import transform_diagonal
  from tensorflow_probability.python.distributions import kullback_leibler
  from tensorflow_probability.python.distributions import mvn_full_covariance
  from tensorflow_probability.python.internal import tensorshape_util
  from tensorflow_probability.python.internal import test_util
label_desc: Use image node_id 1 for calling the MultivariateNormalFullCovarianceTest obj's underlying member method code with example usage: obj._random_pd_matrix() and returns: self