An adaptive stochastic sequential quadratic programming with differentiable exact augmented lagrangians

• Published in: Mathematical programming Vol. 199; no. 1-2; pp. 721 - 791
• Main Authors: Na, Sen, Anitescu, Mihai, Kolar, Mladen
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.05.2023; Springer
• Subjects: Algorithms; Calculus of Variations and Optimal Control; Optimization; Combinatorics; Full Length Paper; Mathematical and Computational Physics; Mathematical Methods in Physics; Mathematics; Mathematics and Statistics; Mathematics of Computing; Numerical Analysis; Theoretical; Exact penalty; 90C30; Sequential quadratic programming; 90C55; Augmented Lagrangian; 90-08; 90C90; 90C15; 90C26; Stochastic optimization
• ISSN: 0025-5610; 1436-4646; 1436-4646
• DOI: 10.1007/s10107-022-01846-z

We consider solving nonlinear optimization problems with a stochastic objective and deterministic equality constraints. We assume for the objective that its evaluation, gradient, and Hessian are inaccessible, while one can compute their stochastic estimates by, for example, subsampling. We propose a stochastic algorithm based on sequential quadratic programming (SQP) that uses a differentiable exact augmented Lagrangian as the merit function. To motivate our algorithm design, we first revisit and simplify an old SQP method Lucidi (J. Optim. Theory Appl. 67(2): 227–245, 1990) developed for solving deterministic problems, which serves as the skeleton of our stochastic algorithm. Based on the simplified deterministic algorithm, we then propose a non-adaptive SQP for dealing with stochastic objective, where the gradient and Hessian are replaced by stochastic estimates but the stepsizes are deterministic and prespecified. Finally, we incorporate a recent stochastic line search procedure Paquette and Scheinberg (SIAM J. Optim. 30(1): 349–376 2020) into the non-adaptive stochastic SQP to adaptively select the random stepsizes, which leads to an adaptive stochastic SQP. The global “almost sure” convergence for both non-adaptive and adaptive SQP methods is established. Numerical experiments on nonlinear problems in CUTEst test set demonstrate the superiority of the adaptive algorithm.