code: def xstep(self):
          r
  
          self.X, self.ss = sp.prox_nuclear(self.S - self.Y - self.U,
                                            1/self.rho)
docs: r"""Minimise Augmented Lagrangian with respect to
          :math:`\mathbf{x}`.