Abstract
We present a primal interior point method for convex quadratic programming which is based upon a logarithmic barrier function approach. This approach generates a sequence of problems, each of which is approximately solved by taking a single Newton step. It is shown that the method requires\(O(\sqrt n L)\) iterations and O(n 3.5 L) arithmetic operations. By using modified Newton steps the number of arithmetic operations required by the algorithm can be reduced to O(n 3 L).
Similar content being viewed by others
References
E.R. Barnes, “A variation on Karmarkar's algorithm for solving linear programming problems,”Mathematical Programming 36 (1986) 174–182.
D.A. Bayer and J.C. Lagarias, “The nonlinear geometry of linear programming,” Preprints, AT&T Bell Laboratories (Murray Hill, NJ, 1986).
I.I. Dikin, “Iterative solution of problems of linear and quadratic programming,”Soviet Mathematics Doklady 8 (1967) 674–675.
A.V. Fiacco and G.P. McCormick,Nonlinear programming: Sequential unconstrained minimization techniques (John Wiley and Sons, New York, 1968).
K.R. Frisch, “The logarithmic potential method of convex programming,” Unpublished manuscript, University Institute of Economics (Oslo, Norway, 1955).
P.E. Gill, G.H. Golub, W. Murray and M.A. Saunders, “Methods for modifying matrix factorizations,”Mathematics of Computation 28 (1974) 505–535.
P.E. Gill, W. Murray, M.A. Saunders, J.A. Tomlin and M.H. Wright, “On projected Newton barrier methods for linear programming and an equivalence to Karmarkar's projective method,”Mathematical Programming 36 (1986) 183–209.
D. Goldfarb, “Matrix factorizations in optimization of nonlinear functions subject to linear constraints,”Mathematical Programming 10 (1975) 1–31.
D. Goldfarb and A. Idnani, “A numerically stable dual method for solving strictly convex quadratic programs,”Mathematical Programming 27 (1983) 1–33.
C.C. Gonzaga, “An algorithm for solving linear programming problems in O(n 3 L) operations,” in: N. Megiddo, ed.,Progress in Mathematical Programming (Springer, Berlin, 1988) pp. 1–28
A.S. Householder,The Theory of Matrices in Numerical Analysis (Bleisdel, Waltham, MA, 1964).
P. Huard, “Resolution of mathematical programming with nonlinear constraints by the method of centres,” in: J. Abadic, ed.,Nonlinear Programming (North-Holland, Amsterdam, 1967).
N. Karmarkar, “A new polynomial time algorithm for linear programming,”Combinatorica 4 (1984) 373–395.
S. Kapoor and P. Vaidya, “Fast algorithms for convex quadratic programming and multicommodity flows,”Proceedings of the 18th Annual ACM Symposium on the Theory of Computing (1986) 147–159.
M. Kojima, S. Mizuno and A. Yoshise, “A primal-dual interior point algorithm for linear programming,” in: N. Meggido, ed.,Progress in Mathematical Programming (Springer, Berlin) pp. 29–97.
M. Kojima, S. Mizuno and A. Yoshise, “A polynomial-time algorithm for a class of linear complementarity problems,”Mathematical Programming 44 (1989) 1–26.
M.K. Kozlov, S.P. Tarasov and L.G. Khachiyan, “Polynomial solvability of convex quadratic programming,”Doklady Akademiia Nauk SSSR 248 (1979). [Translated inSoviet Mathematics Doklady 20 (1979) 1108–1111.]
N. Meggido, “Pathways to the optimal set in linear programming,” in: N. Meggaddo, ed.,Progress in Mathematical Programming (Springer, Berlin, 1988) pp. 131–158.
S. Mehrotra and J. Sun, “An algorithm for convex quadratic programming that requires O(n 3.5 L) arithmetic operations,”Mathematics of Operations Research 15 (1990) 342–363.
R.D.C. Monteiro and I. Adler, “Interior path following primal—dual algorithms. Part I: Linear programming,”Mathematical Programming 44 (1989) 27–41.
R.D.C. Monteiro and I. Adler, “Interior path following primal—dual algorithms. Part II: Convex quadratic programming,”Mathematical Programming 44 (1989) 43–66.
J. Renegar, “A polynomial-time algorithm based on Newton's method for linear programming,”Mathematical Programming 40 (1988) 59–93.
N.Z. Shor, “Utilization of the operation of space dilatation in the minimization of convex functions,”Kibernetika 6, 6–12. [Translated inCybernetics 13 (1970) 94–96.]
P.M. Vaidya, “An algorithm for linear programming which requires O(((m+n)n 2+(m+n) 1.5 n)L) arithmetic operations,”Mathematical Programming 47 (1990) 175–201.
Y. Ye, “Interior algorithms for linear quadratic and linearly constrained convex programming,” Ph.D. Dissertation, Dept. of Engineering-Economic Systems, Stanford University (Stanford, CA, 1987).
Y. Ye, “Further development on the interior algorithm for convex quadratic programming,” Manuscript, Dept. of Engineering-Economic Systems, Stanford University (Stanford, CA, 1987).
Y. Ye and M. Kojima, “Recovering optimal dual solutions in Karmarkar's polynomial algorithm for linear programming,”Mathematical Programming 39 (1985) 305–318.
Y. Ye and E. Tse, “An extension of Karmarkar's projective algorithm for convex quadratic programming,”Mathematical Programming 44 (1989) 157–179.
D.B. Yudin and A.S. Nemirovski, “Informational complexity and effective methods of solution for convex extremal problems,”Ekonomika i Matematicheskie Metody 12 (1976) 357–369. [Translated inMatekon: Translations of Russian and East European Math. Economics, 13 (Spring 1977) 25–45.]
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Goldfarb, D., Liu, S. An O(n 3 L) primal interior point algorithm for convex quadratic programming. Mathematical Programming 49, 325–340 (1990). https://doi.org/10.1007/BF01588795
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01588795