Preconditioning Saddle-Point Systems with Applications in Optimization

  • Authors:
  • H. Sue Dollar;Nicholas I. M. Gould;Martin Stoll;Andrew J. Wathen

  • Affiliations:
  • sue.dollar@stfc.ac.uk and nick.gould@stfc.ac.uk;-;martin.stoll@comlab.ox.ac.uk and andy.wathen@comlab.ox.ac.uk;-

  • Venue:
  • SIAM Journal on Scientific Computing
  • Year:
  • 2010

Quantified Score

Hi-index 0.00

Visualization

Abstract

Saddle-point systems arise in many applications areas, in fact in any situation where an extremum principle arises with constraints. The Stokes problem describing slow viscous flow of an incompressible fluid is a classic example coming from PDEs and in the area of optimization such problems are ubiquitous. In this paper we present a framework into which many well-known methods for solving saddle-point systems fit. Based on this description we show how new approaches for the solution of saddle-point systems arising in optimization can be derived from the Bramble-Pasciak conjugate gradient approach widely used in PDEs and more recent generalizations thereof. In particular we derive a class of new solution methods based on the use of preconditioned conjugate gradients in nonstandard inner products and demonstrate how these can be understood through more standard machinery. We show connections to constraint preconditioning and give the results of numerical computations on a number of standard optimization test examples.