code: def ystep(self):
          r
  
          AXU = self.AX + self.U
          Y0 = sp.proj_l2(self.block_sep0(AXU) - self.S, self.epsilon,
                              axis=self.cri.axisN)
          Y1 = sp.prox_l1(self.block_sep1(AXU), self.wl1 / self.rho)
          self.Y = self.block_cat(Y0, Y1)
  
          super(ConvMinL1InL2Ball, self).ystep()
docs: r"""Minimise Augmented Lagrangian with respect to
          :math:`\mathbf{y}`.