code: def test_reference_implementation_consistency() -> None:
      q_mean = rng.rand(1, 1)
      q_cov = rng.rand(1, 1)
      p_mean = rng.rand(1, 1)
      p_cov = rng.rand(1, 1)
  
      multivariate_KL = multivariate_prior_KL(
          q_mean, p_mean, q_cov, p_cov
      )
      univariate_KL = univariate_prior_KL(
          q_mean.squeeze(),
          p_mean.squeeze(),
          q_cov.squeeze(),
          p_cov.squeeze(),
      )
  
      assert_allclose(univariate_KL, multivariate_KL.squeeze(), atol=4)
  
imports: from typing import Tuple, cast
  import numpy
  import pytest
  from check_shapes import check_shapes
  from numpy.testing import assert_allclose
  import gpflow
  from gpflow.base import AnyNDArray
  from gpflow.kernels import SquaredExponential
  from gpflow.likelihoods import Gaussian
  from tests.gpflow.kernels.reference import ref_rbf_kernel
label_desc: Use image node_id 5 for calling a global function with example usage: test_reference_implementation_consistency() without return types