English

Using incomplete indefinite $LDL^T$ preconditioning for inexact interior point methods for linear programming

Numerical Analysis 2017-08-17 v2

Abstract

Most linear algebra kernels in interior point methods for linear programming require the solution of linear systems of equation with the matrix N=ATD1AN = A^TD^{-1}A (or AD1ATAD^{-1}A^T), where AA denotes the constraint matrix of the linear program. This matrix NN arises from the reduced KKT system by block elimination. If the number of non-zeros in NN or in its Cholesky factorization N=LLTN= LL^T is very large, the computational cost and memory requirement to solve the linear systems of equations with NN may be prohibitively large. In this work we implement an interior point method described by R. Freund and F. Jarre. Forming the normal equation matrix NN is avoided altogether and we work with the reduced KKT system instead. We solve the linear systems for the Newton directions iteratively only to low accuracy using SQMR and an indefinite multilevel preconditioner. Preliminary numerical results are encouraging.

Keywords

Cite

@article{arxiv.1708.04298,
  title  = {Using incomplete indefinite $LDL^T$ preconditioning for inexact interior point methods for linear programming},
  author = {Robert Luce},
  journal= {arXiv preprint arXiv:1708.04298},
  year   = {2017}
}

Comments

APMOD 2012 extended abstract; posted for archivational purpose

R2 v1 2026-06-22T21:14:35.508Z