Convex quartic problems: homogenized gradient method and preconditioning
Radu Dragomir,Yurii Nesterov +1 more
- 30 Jun 2023
TL;DR: In this article , the authors consider a convex minimization problem for which the objective is the sum of a homogeneous polynomial of degree four and a linear term, and they design a first-order method called Homogenized Gradient, along with an accelerated version, which enjoy fast convergence rates.
read more
Abstract: We consider a convex minimization problem for which the objective is the sum of a homogeneous polynomial of degree four and a linear term. Such task arises as a subproblem in algorithms for quadratic inverse problems with a difference-of-convex structure. We design a first-order method called Homogenized Gradient, along with an accelerated version, which enjoy fast convergence rates of respectively $\mathcal{O}(\kappa^2/K^2)$ and $\mathcal{O}(\kappa^2/K^4)$ in relative accuracy, where $K$ is the iteration counter. The constant $\kappa$ is the quartic condition number of the problem. Then, we show that for a certain class of problems, it is possible to compute a preconditioner for which this condition number is $\sqrt{n}$, where $n$ is the problem dimension. To establish this, we study the more general problem of finding the best quadratic approximation of an $\ell_p$ norm composed with a quadratic map. Our construction involves a generalization of the so-called Lewis weights.
read more
Chat with Paper
AI Agents for this Paper
Find similar papers on Google Scholar, PubMed and Arxiv
Write a critical review of this paper
Analyze citations of this paper to find unaddressed research gaps
Figures

Figure 2: Performance of the homogenized gradient method (Algorithm 1) on Problem (6.1) with three different preconditioners: HomGDuses the standard Euclidean norm (B = In), HomGD + optimal preconditioneruses the norm induced by B∗ as computed by Algorithm 3, while HomGD + initial preconditioneruses the norm induced by B(0), the initial iterate of Algorithm 3. As predicted by Section 5.3, the impact of using B∗ over B(0) is significant only for highly coherent matrices. 
Figure 1: Comparison of different methods on the synthetic quartic problem (6.1) for different condition numbers of A. We plot the objective gap f(xk) − f∗, where f∗ is estimated as the minimal objective value computed by all methods. We observe that the restart scheme used in Accel-HomGD (Algorithm 2) is beneficial.
References
•Book
Problem complexity and method efficiency in optimization
John Darzentas
- 01 Jan 1983
TL;DR: In this article, problem complexity and method efficiency in optimisation are discussed in terms of problem complexity, method efficiency, and method complexity in the context of OO optimization, respectively.
Proximal alternating linearized minimization for nonconvex and nonsmooth problems
TL;DR: A self-contained convergence analysis framework is derived and it is established that each bounded sequence generated by PALM globally converges to a critical point.
Phase Retrieval via Wirtinger Flow: Theory and Algorithms
TL;DR: This paper develops a nonconvex formulation of the phase retrieval problem as well as a concrete solution algorithm that is shown to rigorously allow the exact retrieval of phase information from a nearly minimal number of random measurements.
1.4K