code: def test_predict_log_density() -> None:
      vgp0, vgp1, cvgp = _prepare_models()
  
      pred_ydensity0 = vgp0.predict_log_density(
          (Datum.Xtest, Datum.Ytest)
      )
      pred_ydensity_c0 = cvgp.predict_log_density(
          (Datum.Xtest_augmented0, Datum.Ytest_augmented0)
      )
      assert_allclose(pred_ydensity0, pred_ydensity_c0, atol=1e-2)
      pred_ydensity1 = vgp1.predict_log_density(
          (Datum.Xtest, Datum.Ytest)
      )
      pred_ydensity_c1 = cvgp.predict_log_density(
          (Datum.Xtest_augmented1, Datum.Ytest_augmented1)
      )
      assert_allclose(pred_ydensity1, pred_ydensity_c1, atol=1e-2)
  
file_path: GPflow/tests/gpflow/kernels/test_coregion.py
imports: from typing import Sequence, Tuple
  import numpy
  from check_shapes import ShapeChecker
  from numpy.testing import assert_allclose
  import gpflow
  from gpflow import set_trainable
  from gpflow.base import AnyNDArray
  from gpflow.mean_functions import Constant
  from gpflow.models import VGP
label_desc: Use image node_id 7 for calling a global function with example usage: test_predict_log_density() without return types