A Partial Exact Penalty Function Approach for Constrained Optimization

In this paper, we focus on a class of constrained nonlinear optimization problems (NLP), where some of its equality constraints define a closed embedded submanifold $\mathcal{M}$ in $\mathbb{R}^n$. Although NLP can be solved directly by various existing approaches for constrained optimization in Euclidean space, these approaches usually fail to recognize the manifold structure of $\mathcal{M}$. To achieve better efficiency by utilizing the manifold structure of $\mathcal{M}$ in directly applying these existing optimization approaches, we propose a partial penalty function approach for NLP. In our proposed penalty function approach, we transform NLP into the corresponding constraint dissolving problem (CDP) in the Euclidean space, where the constraints that define $\mathcal{M}$ are eliminated through exact penalization. We establish the relationships on the constraint qualifications between NLP and CDP, and prove that NLP and CDP have the same stationary points and KKT points in a neighborhood of the feasible region under mild conditions. Therefore, various existing optimization approaches developed for constrained optimization in the Euclidean space can be directly applied to solve NLP through CDP. Preliminary numerical experiments demonstrate that by dissolving the constraints that define $\mathcal{M}$, CDP gains superior computational efficiency when compared to directly applying existing optimization approaches to solve NLP, especially in high dimensional scenarios.

[1]  Nicolas Boumal An Introduction to Optimization on Smooth Manifolds , 2023 .

[2]  K. Toh,et al.  Solving graph equipartition SDPs on an algebraic variety , 2021, Math. Program..

[3]  K. Toh,et al.  An Improved Unconstrained Approach for Bilevel Optimization , 2022, SIAM Journal on Optimization.

[4]  K. Toh,et al.  A Constraint Dissolving Approach for Nonsmooth Optimization over the Stiefel Manifold , 2022, 2205.10500.

[5]  Tatjana Stykel,et al.  Computing Symplectic Eigenpairs of Symmetric Positive-Definite Matrices via Trace Minimization and Riemannian Optimization , 2021, SIAM J. Matrix Anal. Appl..

[6]  Ying Cui,et al.  Clustering by Orthogonal NMF Model and Non-Convex Penalty Optimization , 2019, IEEE Transactions on Signal Processing.

[7]  A. Schiela,et al.  An SQP method for equality constrained optimization on manifolds , 2020, 2005.06844.

[8]  Suvrit Sra,et al.  Nonconvex stochastic optimization on manifolds via Riemannian Frank-Wolfe methods , 2019, ArXiv.

[9]  Nicolas Boumal,et al.  Simple Algorithms for Optimization on Riemannian Manifolds with Constraints , 2019, Applied Mathematics & Optimization.

[10]  Roland Herzog,et al.  Intrinsic Formulation of KKT Conditions and Constraint Qualifications on Smooth Manifolds , 2018, SIAM J. Optim..

[11]  Song Mei,et al.  Analysis of Sequential Quadratic Programming through the Lens of Riemannian Optimization , 2018 .

[12]  Tim Mitchell,et al.  A BFGS-SQP method for nonsmooth, nonconvex, constrained optimization and its evaluation using relative minimization profiles , 2017, Optim. Methods Softw..

[13]  Xiaojun Chen,et al.  Complexity analysis of interior point algorithms for non-Lipschitz and nonconvex minimization , 2015, Math. Program..

[14]  S. Fang,et al.  On constraint qualifications: motivation, design and inter-relations , 2013 .

[15]  Ruben E. Perez,et al.  pyOpt: a Python-based object-oriented framework for nonlinear constrained optimization , 2011, Structural and Multidisciplinary Optimization.

[16]  B. Afsari Riemannian Lp center of mass: existence, uniqueness, and convexity , 2011 .

[17]  Levent Tunçel,et al.  Optimization algorithms on matrix manifolds , 2009, Math. Comput..

[18]  José Mario Martínez,et al.  Augmented Lagrangian methods under the constant positive linear dependence constraint qualification , 2007, Math. Program..

[19]  José Mario Martínez,et al.  On Augmented Lagrangian Methods with General Lower-Level Constraints , 2007, SIAM J. Optim..

[20]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[21]  R. Osserman,et al.  A Panoramic View of Riemannian Geometry. , 2005 .

[22]  Maher Moakher,et al.  To appear in: SIAM J. MATRIX ANAL. APPL. MEANS AND AVERAGING IN THE GROUP OF ROTATIONS∗ , 2002 .

[23]  Sven Leyffer,et al.  Nonlinear programming without a penalty function , 2002, Math. Program..

[24]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[25]  Jorge Nocedal,et al.  An Interior Point Algorithm for Large-Scale Nonlinear Programming , 1999, SIAM J. Optim..

[26]  Franz Rendl,et al.  QAPLIB – A Quadratic Assignment Problem Library , 1997, J. Glob. Optim..

[27]  Olvi L. Mangasarian,et al.  Computable numerical bounds for lagrange multipliers of stationary points of non-convex differentiable non-linear programs , 1985 .

[28]  F. Clarke Optimization And Nonsmooth Analysis , 1983 .

[29]  H. Karcher Riemannian center of mass and mollifier smoothing , 1977 .

[30]  M. Hestenes Multiplier and gradient methods , 1969 .

[31]  M. Powell A method for nonlinear constraints in minimization problems , 1969 .