docs: def prox_lim(xy, step, boundary=None):
      """Proximal projection operator"""
      if boundary == "circle":
          return prox_circle(xy, step)
      if boundary == "line":
          return prox_line(xy, step)
      # default: do nothing
      return xy