code: def proximal_step(self, grad=None):
          
  
          if grad is None:
              grad = self.eval_grad()
  
          V = self.Y - (1. / self.L) * grad
  
          self.X = self.eval_proxop(V)
  
          return grad
docs: Compute proximal update (gradient descent + regularization).