Crash Start of Interior Point Methods

Technical Report ERGO-2015-012

J. Gondzio

The starting point used by an interior point algorithm for linear and convex quadratic programming may significantly influence the behaviour of the method. A widely used heuristic to construct such a point consists of dropping variable nonnegativity constraints and computing a solution which minimizes the Euclidean norm of the primal (or dual) point while satisfying the appropriate primal (or dual) equality constraints, followed by shifting the variables so that all their components are positive and bounded away from zero. In this Short Communication a new approach for finding a starting point is proposed. It relies on a few inexact Newton steps performed at the start of the solution process. A formal justification of the new heuristic is given and computational results are presented to demonstrate its advantages in practice. Computational experience with a large collection of small- and medium-size test problems reveals that the new starting point is superior to the old one and saves 20-40\% of iterations needed by the primal-dual method. For larger and more difficult problems this translates into remarkable savings in the solution time.

Key words: Interior Point Methods, Initial Point, Crash Start, Linear Programming, Quadratic Programming.

PDF ERGO-15-012.pdf.

Written: October 5, 2015, revised January 6, 2016 and April 6, 2016.
European Journal of Operational Research (to appear).
Published online: June 1, 2016.

Related Software:
HOPDM Higher Order Primal Dual Method.