Table of Contents
Fetching ...

A Boosted-DCA with Power-Sum-DC Decomposition for Linearly Constrained Polynomial Programs

Hu Zhang, Yi-Shuai Niu

TL;DR

The Boosted DCA with Exact Line Search is introduced, and it is demonstrated that the exact line search equates to determining the roots of a univariate polynomial in an interval, with coefficients being computed explicitly based on the power-sum DC decompositions.

Abstract

This paper proposes a novel Difference-of-Convex (DC) decomposition for polynomials using a power-sum representation, achieved by solving a sparse linear system. We introduce the Boosted DCA with Exact Line Search (BDCAe) for addressing linearly constrained polynomial programs within the DC framework. Notably, we demonstrate that the exact line search equates to determining the roots of a univariate polynomial in an interval, with coefficients being computed explicitly based on the power-sum DC decompositions. The subsequential convergence of BDCAe to critical points is proven, and its convergence rate under the Kurdyka-Lojasiewicz property is established. To efficiently tackle the convex subproblems, we integrate the Fast Dual Proximal Gradient (FDPG) method by exploiting the separable block structure of the power-sum DC decompositions. We validate our approach through numerical experiments on the Mean-Variance-Skewness-Kurtosis (MVSK) portfolio optimization model and box-constrained polynomial optimization problems. Comparative analysis of BDCAe against DCA, BDCA with Armijo line search, UDCA, and UBDCA with projective DC decomposition, alongside standard nonlinear optimization solvers FMINCON and FILTERSD, substantiates the efficiency of our proposed approach.

A Boosted-DCA with Power-Sum-DC Decomposition for Linearly Constrained Polynomial Programs

TL;DR

The Boosted DCA with Exact Line Search is introduced, and it is demonstrated that the exact line search equates to determining the roots of a univariate polynomial in an interval, with coefficients being computed explicitly based on the power-sum DC decompositions.

Abstract

This paper proposes a novel Difference-of-Convex (DC) decomposition for polynomials using a power-sum representation, achieved by solving a sparse linear system. We introduce the Boosted DCA with Exact Line Search (BDCAe) for addressing linearly constrained polynomial programs within the DC framework. Notably, we demonstrate that the exact line search equates to determining the roots of a univariate polynomial in an interval, with coefficients being computed explicitly based on the power-sum DC decompositions. The subsequential convergence of BDCAe to critical points is proven, and its convergence rate under the Kurdyka-Lojasiewicz property is established. To efficiently tackle the convex subproblems, we integrate the Fast Dual Proximal Gradient (FDPG) method by exploiting the separable block structure of the power-sum DC decompositions. We validate our approach through numerical experiments on the Mean-Variance-Skewness-Kurtosis (MVSK) portfolio optimization model and box-constrained polynomial optimization problems. Comparative analysis of BDCAe against DCA, BDCA with Armijo line search, UDCA, and UBDCA with projective DC decomposition, alongside standard nonlinear optimization solvers FMINCON and FILTERSD, substantiates the efficiency of our proposed approach.
Paper Structure (22 sections, 10 theorems, 115 equations, 5 figures, 5 algorithms)

This paper contains 22 sections, 10 theorems, 115 equations, 5 figures, 5 algorithms.

Key Result

lemma 1

Let $\varphi(\boldsymbol{x})$ be a form in $\mathbb{H}_d[\boldsymbol{x}]$, expressed as $\varphi(\boldsymbol{x})=\sum_{\boldsymbol{\alpha}\in \mathcal{I}_{n,d}}c_{\boldsymbol{\alpha}}{\boldsymbol{x}}^{\boldsymbol{\alpha}},$ where $\boldsymbol{c}:=(c_{\boldsymbol{\alpha}^1},\ldots,c_{\boldsymbol{\alp where $\boldsymbol{\lambda}=(\lambda_{\boldsymbol{\alpha}^1},\ldots,\lambda_{\boldsymbol{\alpha}^{s

Figures (5)

  • Figure 1: The sparsity of the matrix $\widehat{\boldsymbol{V}}(n,d)$ changes with variable $n$ from $2$ to $16$ for some fixed degree $d\in\{3,4,5,6\}$. It is observed that, for $n>6$, the sparsity of $\widehat{\boldsymbol{V}}(n,d)$ for all $d\in\{3,4,5,6\}$ is more than $90\%$, and the sparsity goes to $1$ as $n\rightarrow\infty$ for any fixed $d$.
  • Figure 2: Log average CPU time along $n$ for results of DCA, BDCA and BDCA$_\text{e}$ for \ref{['eq:T-PSDC']} and \ref{['eq:HD-PSDC']} as well as FILTERSD and FMINCON reported in Tables \ref{['tab:tablet']}--\ref{['tab:tableo']}.
  • Figure 3: Comparison of DCA, BDCA and BDCA$_\text{e}$ applied to HD-PSDC decomposition with $\rho=1$, as well as FILTERSD and FMINCON for solving \ref{['eq:Box']} with $n\in\{10,20,30,40,50\},d=3, den\in\{0.1,0.4,0.7,1\}$.
  • Figure 4: Comparison of DCA, BDCA and BDCA$_\text{e}$ applied to HD-PSDC decomposition with $\rho=1$, as well as FILTERSD and FMINCON for solving \ref{['eq:Box']} with $n\in\{10,20,30,40,50\},d=4, den\in\{0.1,0.4,0.7,1\}$.
  • Figure 5: Log average CPU time along $n$ for results of DCA, BDCA and BDCA$_\text{e}$ for \ref{['eq:T-PSDC']} and \ref{['eq:HD-PSDC']} as well as FILTERSD and FMINCON reported in Table \ref{['tab:t_tablebox']} and Table \ref{['tab:hd_tablebox']}.

Theorems & Definitions (21)

  • lemma 1
  • proof
  • remark 1
  • proposition 1
  • proof
  • proposition 2
  • proof
  • lemma 2
  • proof
  • remark 2
  • ...and 11 more