Abstract:We consider a quadratic program with a few negative eigenvalues (QP-r-NE) subject to linear and convex quadratic constraints that covers many applications and is known to be NP-hard even with one negative eigenvalue (QP1NE). In this paper, we first introduce a new global algorithm (ADMBB), which integrates several simple optimization techniques such as alternative direction method, and branch-and-bound, to find a globally optimal solution to the underlying QP within a pre-specified ϵ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\epsilon $$\end{document}-tolerance. We establish the convergence of the ADMBB algorithm and estimate its complexity. Second, we develop a global search algorithm (GSA) for QP1NE that can locate an optimal solution to QP1NE within ϵ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\epsilon $$\end{document}-tolerance and estimate the worst-case complexity bound of the GSA. Preliminary numerical results demonstrate that the ADMBB algorithm can effectively find a global optimal solution to large-scale QP-r-NE instances when r≤10\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r\le 10$$\end{document}, and the GSA outperforms the ADMBB for most of the tested QP1NE instances. The software reviewed as part of this submission was given the DOI (digital object identifier) 10.5281/zenodo.1344739.

Numerical methods for QBD matrix equations

Numerical Method for the Time Optimal Control Problem Governed by the Benjamin-Bona-Mahony Equation

On the Extreme Solutions of an Nonlinear Matrix Equation.

On the Iterative Methods for the Solution of Three Types of Nonlinear Matrix Equations

On Two Iteration Methods for the Quadratic Matrix Equations

Solutions and perturbation analysis of the matrix equation X - \sum_{i=1}^m A_i^* X^{-1} A_i = Q

The Extremal Solution Of The Matrix Equation Xs + A*X-Qa = I

QQMR: A Structure-Preserving Quaternion Quasi-Minimal Residual Method for Non-Hermitian Quaternion Linear Systems

Hermitian positive definite solution of nonlinear matrix equation X + A* X-2A = Q* 1

Towards positive definite solutions of a class of nonlinear matrix equations.

On Hermitian positive definite solution of nonlinear matrix equation X plus A* X-2 A = Q

Newton’s Iteration Method for Solving the Nonlinear Matrix Equation X+∑i=1mAi*X−1Ai=Q

On the Positive Definite Solutions of a Nonlinear Matrix Equation

ON THE PURE IMAGINARY QUATERNIONIC LEAST SQUARES SOLUTIONS OF MATRIX EQUATION

New global algorithms for quadratic programming with a few negative eigenvalues based on alternative direction method and convex relaxation

Numerical Solutions of Quaternionic Riccati Equations

A Modified Newton-based Matrix Splitting Iteration Method for Generalized Absolute Value Equations

Positive Solutions for a Class of Singular Boundary Value Problems with Fractional Q-Difference Equations

Special least squares solutions of the reduced biquaternion matrix equation with applications

LSQR Iterative Method for Generalized Coupled Sylvester Matrix Equations

Stable Solutions of Linear Systems Involving Long Chain of Matrix Multiplications