14,539 research outputs found
A second derivative SQP method: theoretical issues
Sequential quadratic programming (SQP) methods form a class of highly efficient algorithms for solving nonlinearly constrained optimization problems. Although second derivative information may often be calculated, there is little practical theory that justifies exact-Hessian SQP methods. In particular, the resulting quadratic programming (QP) subproblems are often nonconvex, and thus finding their global solutions may be computationally nonviable. This paper presents a second-derivative SQP method based on quadratic subproblems that are either convex, and thus may be solved efficiently, or need not be solved globally. Additionally, an explicit descent-constraint is imposed on certain QP subproblems, which “guides” the iterates through areas in which nonconvexity is a concern. Global convergence of the resulting algorithm is established
A second derivative SQP method: local convergence
In [19], we gave global convergence results for a second-derivative SQP method for minimizing the exact â„“1-merit function for a fixed value of the penalty parameter. To establish this result, we used the properties of the so-called Cauchy step, which was itself computed from the so-called predictor step. In addition, we allowed for the computation of a variety of (optional) SQP steps that were intended to improve the efficiency of the algorithm. \ud
\ud
Although we established global convergence of the algorithm, we did not discuss certain aspects that are critical when developing software capable of solving general optimization problems. In particular, we must have strategies for updating the penalty parameter and better techniques for defining the positive-definite matrix Bk used in computing the predictor step. In this paper we address both of these issues. We consider two techniques for defining the positive-definite matrix Bk—a simple diagonal approximation and a more sophisticated limited-memory BFGS update. We also analyze a strategy for updating the penalty paramter based on approximately minimizing the ℓ1-penalty function over a sequence of increasing values of the penalty parameter.\ud
\ud
Algorithms based on exact penalty functions have certain desirable properties. To be practical, however, these algorithms must be guaranteed to avoid the so-called Maratos effect. We show that a nonmonotone varient of our algorithm avoids this phenomenon and, therefore, results in asymptotically superlinear local convergence; this is verified by preliminary numerical results on the Hock and Shittkowski test set
Inexact Convex Relaxations for AC Optimal Power Flow: Towards AC Feasibility
Convex relaxations of AC optimal power flow (AC-OPF) problems have attracted
significant interest as in several instances they provably yield the global
optimum to the original non-convex problem. If, however, the relaxation is
inexact, the obtained solution is not AC-feasible. The quality of the obtained
solution is essential for several practical applications of AC-OPF, but
detailed analyses are lacking in existing literature. This paper aims to cover
this gap. We provide an in-depth investigation of the solution characteristics
when convex relaxations are inexact, we assess the most promising AC
feasibility recovery methods for large-scale systems, and we propose two new
metrics that lead to a better understanding of the quality of the identified
solutions. We perform a comprehensive assessment on 96 different test cases,
ranging from 14 to 3120 buses, and we show the following: (i) Despite an
optimality gap of less than 1%, several test cases still exhibit substantial
distances to both AC feasibility and local optimality and the newly proposed
metrics characterize these deviations. (ii) Penalization methods fail to
recover an AC-feasible solution in 15 out of 45 cases, and using the proposed
metrics, we show that most failed test instances exhibit substantial distances
to both AC-feasibility and local optimality. For failed test instances with
small distances, we show how our proposed metrics inform a fine-tuning of
penalty weights to obtain AC-feasible solutions. (iii) The computational
benefits of warm-starting non-convex solvers have significant variation, but a
computational speedup exists in over 75% of the cases
- …