def METHOD_NAME(
    call_result: Union[np.ndarray, np.floating],
    call_result_naive: Union[np.ndarray, np.floating],
):
    """Test whether the entries of the output of ``CovarianceFunction.__call__`` match
    the entries generated by the naive reference implementation."""
    np.testing.assert_allclose(
        call_result,
        call_result_naive,
        rtol=10**-12,
        atol=10**-12,
    )