def METHOD_NAME(x1, x2):
    """
    Surrogate gradient of the squared exponential covariance function,
    with respect to `x1`.
    """
    # Corresponds to equation 11 in the paper.
    return cov_fun(x1, x2) * (x2 - x1) / loss_params["lengthscale"] ** 2