The nonlinear geometry of linear programming. I. Affine and projective scaling trajectories

Authors:
D. A. Bayer and J. C. Lagarias

Journal:
Trans. Amer. Math. Soc. **314** (1989), 499-526

MSC:
Primary 90C05; Secondary 34A99, 52A40

DOI:
https://doi.org/10.1090/S0002-9947-1989-1005525-6

MathSciNet review:
1005525

Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: This series of papers studies a geometric structure underlying Karmarkar's projective scaling algorithm for solving linear programming problems. A basic feature of the projective scaling algorithm is a vector field depending on the objective function which is defined on the interior of the polytope of feasible solutions of the linear program. The geometric structure studied is the set of trajectories obtained by integrating this vector field, which we call -trajectories. We also study a related vector field, the affine scaling vector field, and its associated trajectories, called -trajectories. The affine scaling vector field is associated to another linear programming algorithm, the affine scaling algorithm. Affine and projective scaling vector fields are each defined for linear programs of a special form, called strict standard form and canonical form, respectively.

This paper derives basic properties of -trajectories and -trajectories. It reviews the projective and affine scaling algorithms, defines the projective and affine scaling vector fields, and gives differential equations for -trajectories and -trajectories. It shows that projective transformations map -trajectories into -trajectories. It presents Karmarkar's interpretation of -trajectories as steepest descent paths of the objective function with respect to the Riemannian geometry restricted to the relative interior of the polytope of feasible solutions. -trajectories of a canonical form linear program are radial projections of -trajectories of an associated standard form linear program. As a consequence there is a polynomial time linear programming algorithm using the affine scaling vector field of this associated linear program: This algorithm is essentially Karmarkar's algorithm.

These trajectories are studied in subsequent papers by two nonlinear changes of variables called Legendre transform coordinates and projective Legendre transform coordinates, respectively. It will be shown that -trajectories have an algebraic and a geometric interpretation. They are algebraic curves, and they are geodesics (actually distinguished chords) of a geometry isometric to a Hubert geometry on a polytope combinatorially dual to the polytope of feasible solutions. The -trajectories of strict standard form linear programs have similar interpretations: They are algebraic curves, and are geodesics of a geometry isometric to Euclidean geometry.

**[AKRV]**I. Adler, N. Karmarkar, G. Resende and S. Veiga,*An implementation of Karmarkar's algorithm for linear programming*, preprint, Univ. of California, Berkeley, 1986.**[A]**Kurt M. Anstreicher,*A monotonic projective algorithm for fractional linear programming*, Algorithmica**1**(1986), no. 4, 483–498. MR**880734**, https://doi.org/10.1007/BF01840458**[Ar]**V. I. Arnol′d,*Mathematical methods of classical mechanics*, Springer-Verlag, New York-Heidelberg, 1978. Translated from the Russian by K. Vogtmann and A. Weinstein; Graduate Texts in Mathematics, 60. MR**0690288****[B]**Earl R. Barnes,*A variation on Karmarkar’s algorithm for solving linear programming problems*, Math. Programming**36**(1986), no. 2, 174–182. MR**866987**, https://doi.org/10.1007/BF02592024**[BL2]**D. A. Bayer and J. C. Lagarias,*The nonlinear geometry of linear programming. II. Legendre transform coordinates and central trajectories*, Trans. Amer. Math. Soc.**314**(1989), no. 2, 527–581. MR**1005526**, https://doi.org/10.1090/S0002-9947-1989-1005526-8**[BP]**H. Busemann and B. B. Phadke,*A general version of Beltrami’s theorem in the large*, Pacific J. Math.**115**(1984), no. 2, 299–315. MR**765189****[Bu]**Herbert Busemann,*The geometry of geodesics*, Academic Press Inc., New York, N. Y., 1955. MR**0075623****[Bu2]**Herbert Busemann,*Spaces with distinguished shortest joins*, A spectrum of mathematics (Essays presented to H. G. Forder), Auckland Univ. Press, Auckland, 1971, pp. 108–120. MR**0431075****[CH]**R. Courant and D. Hubert,*Methods of mathematical physics*, Vols. I, II, Wiley, New York, 1962.**[D1]**I. I. Dikin,*Iterative solution of problems of linear and quadratic programming*, Dokl. Akad. Nauk SSSR**174**(1967), 747–748 (Russian). MR**0221850****[D2]**-,*About the convergence of an iterative process*, Controllable Systems IM IK SO AN SSR 1974, No. 12, 54-60. (Russian)**[F]**W. Fenchel,*On conjugate convex functions*, Canadian J. Math.**1**(1949), 73–77. MR**0028365****[FM]**Anthony V. Fiacco and Garth P. McCormick,*Nonlinear programming: Sequential unconstrained minimization techniques*, John Wiley and Sons, Inc., New York-London-Sydney, 1968. MR**0243831****[Fl]**Wendell H. Fleming,*Functions of several variables*, Addison-Wesley Publishing Co., Inc., Reading, MAss.-London, 1965. MR**0174675****[GZ]**C. B. Garcia and W. I. Zangwill,*Pathways to solutions, fixed points and equilibria*, Prentice-Hall, Englewood Cliffs, N. J., 1981.**[Go]**Clovis C. Gonzaga,*An algorithm for solving linear programming problems in 𝑂(𝑛³𝐿) operations*, Progress in mathematical programming (Pacific Grove, CA, 1987) Springer, New York, 1989, pp. 1–28. MR**982713****[H]**D. Hilbert,*Grundlagen der Geometrie*, 7th ed., Leipzig, 1930. (English transl.,*Foundations of geometry*.)**[Ho]**J. Hooker,*The projective linear programming algorithm*, Interfaces**16**(1986), 75-90.**[II]**Masao Iri and Hiroshi Imai,*A multiplicative barrier function method for linear programming*, Algorithmica**1**(1986), no. 4, 455–482. MR**880733**, https://doi.org/10.1007/BF01840457**[K]**N. Karmarkar,*A new polynomial-time algorithm for linear programming*, Combinatorica**4**(1984), no. 4, 373–395. MR**779900**, https://doi.org/10.1007/BF02579150**[KLSW]**Narendra K. Karmarkar, Jeffrey C. Lagarias, Lev Slutsman, and Pyng Wang,*Power series variants of Karmarkar-type algorithms*, AT&T Tech. J.**68**(1989), no. 3, 20–36. MR**1021069****[KV]**S. Kapoor and P. M. Vaidya,*Fast algorithms for convex quadratic programming and multicommodity flows*, Proc. 18th ACM Sympos. on Theory of Computing, 1986, pp. 147-159.**[KMY]**Masakazu Kojima, Shinji Mizuno, and Akiko Yoshise,*A primal-dual interior point algorithm for linear programming*, Progress in mathematical programming (Pacific Grove, CA, 1987) Springer, New York, 1989, pp. 29–47. MR**982714****[L3]**J. C. Lagarias,*The nonlinear geometry of linear programming. III. Projective Legendre transform coordinates and Hilbert geometry*, Trans. Amer. Math. Soc.**320**(1990), no. 1, 193–225. MR**1058199**, https://doi.org/10.1090/S0002-9947-1990-1058199-0**[Ln]**Cornelius Lanczos,*The Variational Principles of Mechanics*, Mathematical Expositions, no. 4, University of Toronto Press, Toronto, Ont., 1949. MR**0034139****[M1]**Nimrod Megiddo,*On the complexity of linear programming*, Advances in economic theory (Cambridge, MA, 1985) Econom. Soc. Monogr., vol. 12, Cambridge Univ. Press, Cambridge, 1989, pp. 225–268. MR**1117042****[M2]**Nimrod Megiddo,*Pathways to the optimal set in linear programming*, Progress in mathematical programming (Pacific Grove, CA, 1987) Springer, New York, 1989, pp. 131–158. MR**982720****[MS]**N. Megiddo and M. Shub,*Boundary behavior of interior point algorithms in linear programming*, IBM Research Report RJ-5319, Sept. 1986.**[N]**J. L. Nazareth,*Homotopy techniques in linear programming*, Algorithmica**1**(1986), no. 4, 529–535. MR**880737**, https://doi.org/10.1007/BF01840461**[Re]**James Renegar,*A polynomial-time algorithm, based on Newton’s method, for linear programming*, Math. Programming**40**(1988), no. 1, (Ser. A), 59–93. MR**923697**, https://doi.org/10.1007/BF01580724**[R1]**R. T. Rockafellar,*Conjugates and Legendre transforms of convex functions*, Canad. J. Math.**19**(1967), 200–205. MR**0213496**, https://doi.org/10.4153/CJM-1967-012-4**[R2]**-,*Convex analysis*, Princeton Univ. Press, Princeton, N. J., 1970.**[Sh]**Michael Shub,*On the asymptotic behavior of the projective rescaling algorithm for linear programming*, J. Complexity**3**(1987), no. 3, 258–269. MR**919676**, https://doi.org/10.1016/0885-064X(87)90015-X**[So1]**Gy. Sonnevend,*An “analytical centre” for polyhedrons and new classes of global algorithms for linear (smooth, convex) programming*, System modelling and optimization (Budapest, 1985) Lect. Notes Control Inf. Sci., vol. 84, Springer, Berlin, 1986, pp. 866–875. MR**903521**, https://doi.org/10.1007/BFb0043914**[So2]**-,*A new method for solving a set of linear (convex) inequalities and its application for identification and optimization*, Proc. Sympos. on Dynamic Modelling, IFAC-IFORS, Budapest, June 1986.**[SW]**Josef Stoer and Christoph Witzgall,*Convexity and optimization in finite dimensions. I*, Die Grundlehren der mathematischen Wissenschaften, Band 163, Springer-Verlag, New York-Berlin, 1970. MR**0286498****[Va]**P. Vaidya,*An algorithm for linear programming which requires**arithmetic operations*, Proc. 19th ACM Sympos. on Theory of Computing, 1987, pp. 29-38.**[VMF]**Robert J. Vanderbei, Marc S. Meketon, and Barry A. Freedman,*A modification of Karmarkar’s linear programming algorithm*, Algorithmica**1**(1986), no. 4, 395–407. MR**880730**, https://doi.org/10.1007/BF01840454

Retrieve articles in *Transactions of the American Mathematical Society*
with MSC:
90C05,
34A99,
52A40

Retrieve articles in all journals with MSC: 90C05, 34A99, 52A40

Additional Information

DOI:
https://doi.org/10.1090/S0002-9947-1989-1005525-6

Article copyright:
© Copyright 1989
American Mathematical Society