def METHOD_NAME(fun: Callable[..., tf.Tensor]) -> None:
    inducing_variable = mf.SharedIndependentInducingVariables(make_ip())
    kernel = mk.LinearCoregionalization(make_kernels(Datum.L), Datum.W)
    if fun is covariances.Kuu:
        t = tf.linalg.cholesky(fun(inducing_variable, kernel, jitter=1e-9))
        assert t.shape == (Datum.L, Datum.M, Datum.M)
    else:
        t = fun(inducing_variable, kernel, Datum.Xnew)
        assert t.shape == (Datum.L, Datum.M, Datum.N)