$$\rho$$-regularization subproblems: strong duality and an eigensolver-based algorithm

Trust-region (TR) type method, based on a quadratic model such as the trust-region subproblem (TRS) and p-regularization subproblem (pRS), is arguably one of the most successful methods for unconstrained minimization. In this paper, we study a general regularized subproblem (named ρRS), which covers TRS and pRS as special cases. We derive a strong duality theorem for ρRS, and also its necessary and sufficient optimality condition under general assumptions on the regularization term. We then define the Rendl-Wolkowicz (RW) dual problem of ρRS, which is a maximization problem whose objective function is concave, and differentiable except possibly at two points. It is worth pointing out that our definition is based on an alternative derivation of the RW-dual problem for TRS. Then we propose an eigensolver-based algorithm for solving the RW-dual problem of ρRS. The algorithm is carried out by finding the smallest eigenvalue and its unit eigenvector of a certain matrix in each iteration. Finally, we present numerical results on randomly generated pRS’s, and on a new class of regularized problem that combines TRS and pRS, to illustrate our algorithm. Ting Kei Pong was supported in part by an internal funding, G-UAKK, of the Hong Kong Polytechnic University. Liaoyuan Zeng Department of Applied Mathematics, The Hong Kong Polytechnic University, Hong Kong, PRC E-mail: lyzeng@polyu.edu.hk Ting Kei Pong Department of Applied Mathematics, The Hong Kong Polytechnic University, Hong Kong, PRC E-mail: tk.pong@polyu.edu.hk 2 Liaoyuan Zeng, Ting Kei Pong

[1]  Nicholas I. M. Gould,et al.  Error estimates for iterative algorithms for minimizing regularized quadratic subproblems , 2020, Optim. Methods Softw..

[2]  Leon Hirsch,et al.  Fundamentals Of Convex Analysis , 2016 .

[3]  Adrian S. Lewis,et al.  Convex Analysis And Nonlinear Optimization , 2000 .

[4]  Akiko Takeda,et al.  Solving the Trust-Region Subproblem By a Generalized Eigenvalue Problem , 2017, SIAM J. Optim..

[5]  William W. Hager,et al.  Minimizing a Quadratic Over a Sphere , 2001, SIAM J. Optim..

[6]  Jorge J. Moré,et al.  Computing a Trust Region Step , 1983 .

[7]  Henry Wolkowicz,et al.  The trust region subproblem and semidefinite programming , 2004, Optim. Methods Softw..

[8]  Yurii Nesterov,et al.  Cubic regularization of Newton method and its global performance , 2006, Math. Program..

[9]  Henry Wolkowicz,et al.  Indefinite Trust Region Subproblems and Nonsymmetric Eigenvalue Perturbations , 1995, SIAM J. Optim..

[10]  Franz Rendl,et al.  A semidefinite framework for trust region subproblems with applications to large scale minimization , 1997, Math. Program..

[11]  Henry Wolkowicz,et al.  The generalized trust region subproblem , 2014, Comput. Optim. Appl..

[12]  J. Borwein,et al.  Convex Analysis And Nonlinear Optimization , 2000 .

[13]  Danny C. Sorensen,et al.  Algorithm 873: LSTRS: MATLAB software for large-scale trust-region subproblems and regularization , 2008, TOMS.

[14]  Ya-Xiang Yuan,et al.  Theory and application of p-regularized subproblems for p>2 , 2017, Optim. Methods Softw..

[15]  Felix Lieder Solving Large-Scale Cubic Regularization by a Generalized Eigenvalue Problem , 2020, SIAM J. Optim..

[16]  Nicholas I. M. Gould,et al.  Adaptive cubic regularisation methods for unconstrained optimization. Part I: motivation, convergence and numerical results , 2011, Math. Program..

[17]  Ya-Xiang Yuan,et al.  Theory and application of p-regularized subproblems for p>2 , 2017, Optim. Methods Softw..

[18]  Nicholas I. M. Gould,et al.  On solving trust-region and other regularised subproblems in optimization , 2010, Math. Program. Comput..

[19]  Yair Carmon,et al.  Analysis of Krylov Subspace Solutions of Regularized Non-Convex Quadratic Problems , 2018, NeurIPS.

[20]  Nicholas I. M. Gould,et al.  Solving the Trust-Region Subproblem using the Lanczos Method , 1999, SIAM J. Optim..

[21]  C. Zălinescu Convex analysis in general vector spaces , 2002 .

[22]  Nicholas I. M. Gould,et al.  Trust Region Methods , 2000, MOS-SIAM Series on Optimization.

[23]  Danny C. Sorensen,et al.  A New Matrix-Free Algorithm for the Large-Scale Trust-Region Subproblem , 2000, SIAM J. Optim..