code: def ystep(self):
          r
  
          self.Y = np.asarray(sp.prox_l1(self.S - self.AX - self.U,
                                         self.lmbda/self.rho), dtype=self.dtype)
docs: r"""Minimise Augmented Lagrangian with respect to
          :math:`\mathbf{y}`.