def ystep(self):
"""Minimise Augmented Lagrangian with respect to :math:`\mathbf{y}`."""
self.Y[...,0:-1] = sl.shrink2(self.AX[...,0:-1] + self.U[...,0:-1],
(self.lmbda/self.rho)*self.Wtvna, axis=self.saxes)
self.Y[...,-1] = sl.shrink1(self.AX[...,-1] + self.U[...,-1] - self.S,
(1.0/self.rho)*self.Wdf)
def ystep(self):
r"""Minimise Augmented Lagrangian with respect to
:math:`\mathbf{y}`.
"""
self.Y = sl.shrink2(self.AX + self.U, (self.lmbda/self.rho)*self.Wtvna,
axis=self.saxes)
def ystep(self):
r"""Minimise Augmented Lagrangian with respect to
:math:`\mathbf{y}`."""
AXU = self.AX + self.U
self.Y[..., 0:-1] = sl.shrink2(AXU[..., 0:-1], self.mu / self.rho)
self.Y[..., -1] = sl.shrink1(AXU[..., -1],
(self.lmbda / self.rho) * self.Wl1)
def ystep(self):
"""Minimise Augmented Lagrangian with respect to y."""
self.Y[..., 0:-1] = sl.shrink2(self.AX[..., 0:-1] + self.U[..., 0:-1],
(self.lmbda / self.rho) * self.Wtvna)
self.Y[...,
-1] = sl.shrink1(self.AX[..., -1] + self.U[..., -1] - self.S,
(1.0 / self.rho))
def ystep(self):
r"""Minimise Augmented Lagrangian with respect to
:math:`\mathbf{y}`."""
AXU = self.AX + self.U
self.block_sep0(self.Y)[:] = sl.shrink1(
self.block_sep0(AXU), (self.lmbda/self.rho) * self.Wl1)
self.block_sep1(self.Y)[:] = sl.shrink2(
self.block_sep1(AXU), self.mu/self.rho, axis=(self.cri.axisC, -1))
def ystep(self):
"""Minimise Augmented Lagrangian with respect to y."""
self.Y = sl.shrink2(self.AX + self.U,
(self.lmbda / self.rho) * self.Wtvna)
def ystep(self):
"""Minimise Augmented Lagrangian with respect to y."""
self.Y = np.asarray(sl.shrink2(self.AX + self.U,
(self.lmbda / self.rho) * self.Wtvna),
dtype=self.dtype)
