Optimization Online


Iterative Solution of Augmented Systems Arising in Interior Methods

Anders Forsgren (andersf***at***kth.se)
Philip Gill (pgill***at***ucsd.edu)
Joshua Griffin (jgriffi***at***sandia.gov)

Abstract: Iterative methods are proposed for certain augmented systems of linear equations that arise in interior methods for general nonlinear optimization. Interior methods define a sequence of KKT equations that represent the symmetrized (but indefinite) equations associated with Newton's method for a point satisfying the perturbed optimality conditions. These equations involve both the primal and dual variables and become increasingly ill-conditioned as the optimization proceeds. In this context, an iterative linear solver must not only handle the ill-conditioning but also detect the occurrence of KKT matrices with the wrong matrix inertia. A one-parameter family of equivalent linear equations is formulated that includes the KKT system as a special case. The discussion focuses on a particular system from this family, known as the ``doubly-augmented system'', that is positive definite with respect to both the primal and dual variables. This property means that a standard preconditioned conjugate-gradient method involving both primal and dual variables will either terminate successfully or detect if the KKT matrix has the wrong inertia. Constraint preconditioning is a well-known technique for preconditioning the conjugate-gradient method on augmented systems. A family of constraint preconditioners is proposed that provably eliminates the inherent ill-conditioning in the augmented system. A considerable benefit of combining constraint preconditioning with the doubly-augmented system is that the preconditioner need not be applied exactly. Two particular ``active-set'' constraint preconditioners are formulated that involve only a subset of the rows of the augmented system and thereby may be applied with considerably less work. Finally, some numerical experiments illustrate the numerical performance of the proposed preconditioners and highlight some theoretical properties of the preconditioned matrices.

Keywords: Large-scale nonlinear programming, nonconvex optimization, interior methods, augmented systems, KKT systems, iterative methods, conjugate-gradient method, preconditioning.

Category 1: Nonlinear Optimization (Constrained Nonlinear Optimization )

Citation: To appear in SIAM Journal on Optimization.


Entry Submitted: 08/31/2005
Entry Accepted: 08/31/2005
Entry Last Modified: 03/20/2007

Modify/Update this entry

  Visitors Authors More about us Links
  Subscribe, Unsubscribe
Digest Archive
Search, Browse the Repository


Coordinator's Board
Classification Scheme
Give us feedback
Optimization Journals, Sites, Societies
Mathematical Programming Society