def prox_gradf_lim(xy, step, boundary=None):
    """Forward-backward step: gradient, followed by projection"""
    return prox_lim(prox_gradf(xy,step), step, boundary=boundary)
Forward-backward step: gradient, followed by projection