On the complexity of following the central path of linear programs by linear extrapolation II

A class of algorithms is proposed for solving linear programming problems (withm inequality constraints) by following the central path using linear extrapolation with a special adaptive choice of steplengths. The latter is based on explicit results concerning the convergence behaviour of Newton's method to compute points on the central pathx(r), r>0, and this allows to estimate the complexity, i.e. the total numberN = N(R, δ) of steps needed to go from an initial pointx(R) to a final pointx(δ), R>δ>0, by an integral of the local “weighted curvature” of the (primal—dual) path. Here, the central curve is parametrized with the logarithmic penalty parameterr↓0. It is shown that for large classes of problems the complexity integral, i.e. the number of stepsN, is not greater than constmα log(R/δ), whereα < 1/2 e.g.α = 1/4 orα = 3/8 (note thatα = 1/2 gives the complexity of zero order methods). We also provide a lower bound for the complexity showing that for some problems the above estimation can hold only forα ⩾ 1/3.As a byproduct, many analytical and structural properties of the primal—dual central path are obtained: there are, for instance, close relations between the weighted curvature and the logarithmic derivatives of the slack variables; the dependence of these quantities on the parameterr is described. Also, related results hold for a family of weighted trajectories, into which the central path can be embedded.

[1]  G. Sonnevend Sequential algorithms of optimal order global error for the uniform recovery of functions with monotone (r-1) derivatives , 1984 .

[2]  J. Stoer,et al.  An implementation of the method of analytic centers , 1988 .

[3]  M. Todd,et al.  Recent Developments and New Directions in Linear Programming , 1988 .

[4]  F. Jarre On the method of analytic centers for solving smooth convex programs , 1988 .

[5]  Aeneas Marxen Primal Barrier Methods for Linear Programming , 1989 .

[6]  C. C. Gonzaga,et al.  An Algorithm for Solving Linear Programming Problems in O(n 3 L) Operations , 1989 .

[7]  Paul D. Domich,et al.  Optimal 3-dimensional methods for linear programming , 1989 .

[8]  Paul D. Domich,et al.  Algorithmic Enhancements to the Method of Centers for Linear Programming Problems , 1989, INFORMS J. Comput..

[9]  Mauricio G. C. Resende,et al.  An implementation of Karmarkar's algorithm for linear programming , 1989, Math. Program..

[10]  J. Stoer,et al.  Global ellipsoidal approximations and homotopy methods for solving convex analytic programs , 1990 .

[11]  György Sonnevend Application of Analytic Centers to Feedback Design for Systems with Uncertainties , 1990 .

[12]  Masao Iri,et al.  Integrability of vector and multivector fields associated with interior point methods for linear programming , 1991, Math. Program..

[13]  I. Lustig,et al.  Computational experience with a primal-dual interior point method for linear programming , 1991 .

[14]  Shinji Mizuno,et al.  An $$O(\sqrt n L)$$ iteration potential reduction algorithm for linear complementarity problems , 1991, Math. Program..

[15]  Jean-Philippe Vial,et al.  A polynomial method of approximate centers for linear programming , 1992, Math. Program..

[16]  J. Goffin,et al.  Decomposition and nondifferentiable optimization with the projective algorithm , 1992 .