A brief history of linear and mixed-integer programming computation
For many of us, modern-day linear programming (LP) started with the work of George Dantzig in 1947. However, it must be said that many other scientists have also made seminal contributions to the subject, and some would argue that the origins of LP predate Dantzig’s contribution. It is matter open to debate [36]. However, what is not open to debate is Dantzig’s key contribution to LP computation. In contrast to the economists of his time, Dantzig viewed LP not just as a qualitative tool in the analysis of economic phenomena, but as a method that could be used to compute actual answers to specific real-world problems. Consistent with that view, he proposed an algorithm for solving LPs, the simplex algorithm [12]. To this day the simplex algorithm remains a primary computational tool in linear and mixed-integer programming (MIP). In [11] it is reported that the first application of Dantzig’s simplex algorithm to the solution of a non-trivial LP was Laderman’s solution of a 21 constraint, 77 variable instance of the classical Stigler Diet Problem [41]. It is reported that the total computation time was 120 man-days! The first computer implementation of an at-least modestly general version of the simplex algorithm is reported to have been on the SEAC computer at the then National Bureau of Standards [25]. (There were apparently some slightly earlier implementations for dealing with models that were “triangular”, that is, where all the linear systems could be solved by simple addition and subtraction.) Orchard-Hays [35] reports that several small instances having as many as 10 constraints and 20 variables were solved with this implementation. The first systematic development of computer codes for the simplex algorithm began very shortly thereafter at the RAND Corporation in Santa Monica, California. Dantzig’s initial LP work occurred at the Air Force following
- Conference Article
5
- 10.2118/17777-ms
- Jun 27, 1988
A common type of mathematical optimization is Linear Programming (LP). An LP solution of aquifer influence functions has recently been reported by Gadjica, etal.1 (1987) and Targac, etal.2 Their LP matrices were large and sparse (only 3% of the elements were non-zero) and were solved on main frame computers. Another recent application of LP is equation-of-state matching of laboratory PVT data3. This problem leads to a smaller, denser LP matrix. Three methods of LP solution were investigated on microcomputers: (1) the simplex method, (2) the revised simplex method, and (3) the symmetric method. These methods were run on several LP problems ranging from a small dense matrix to large sparse matrices. The different methods have different characteristics which affect the speed, storage requirements and simplicity of coding. The simplex method is straightforward, but usually is slower and requires more storage than the other methods. The results of this study are tabulated with running times and storage requirements for the various LP methods and microcomputers. The computers range from the IBM XT to the Compaq 386. This information serves as a documentation of the LP codes and should be useful for an engineer interested in using LP codes on a microcomputer.
- Single Book
6
- 10.1201/9781003256052
- Nov 29, 2021
This new volume provides the information needed to understand the simplex method, the revised simplex method, dual simplex method, and more for solving linear programming problems.Following a logical order, the book first gives a mathematical model of the linear problem programming and describes the usual assumptions under which the problem is solved. It gives a brief description of classic algorithms for solving linear programming problems as well as some theoretical results. It goes on to explain the definitions and solutions of linear programming problems, outlining the simplest geometric methods and showing how they can be implemented. Practical examples are included along the way. The book concludes with a discussion of multi-criteria decision-making methods.Advances in Optimization and Linear Programming is a highly useful guide to linear programming for professors and students in optimization and linear programming.
- Single Book
3524
- 10.1002/9781119454816
- Oct 22, 2019
Engineering Optimization Theory and Practice
- Research Article
36
- 10.6100/ir735328
- Nov 18, 2015
- Data Archiving and Networked Services (DANS)
Secure multiparty computation is a basic concept of growing interest in modern cryptography. It allows a set of mutually distrusting parties to perform a computation on their private information in such a way that as little as possible is revealed about each private input. The early results of multiparty computation have only theoretical signi??cance since they are not able to solve computationally complex functions in a reasonable amount of time. Nowadays, e??ciency of secure multiparty computation is an important topic of cryptographic research. As a case study we apply multiparty computation to solve the problem of secure linear programming. The results enable, for example in the context of the EU-FP7 project SecureSCM, collaborative supply chain management. Collaborative supply chain management is about the optimization of the supply and demand con??guration of a supply chain. In order to optimize the total bene??t of the entire chain, parties should collaborate by pooling their sensitive data. With the focus on e??ciency we design protocols that securely solve any linear program using the simplex algorithm. The simplex algorithm is well studied and there are many variants of the simplex algorithm providing a simple and e??cient solution to solving linear programs in practice. However, the cryptographic layer on top of any variant of the simplex algorithm imposes restrictions and new complexity measures. For example, hiding the number of iterations of the simplex algorithm has the consequence that the secure implementations have a worst case number of iterations. Then, since the simplex algorithm has exponentially many iterations in the worst case, the secure implementations have exponentially many iterations in all cases. To give a basis for understanding the restrictions, we review the basic theory behind the simplex algorithm and we provide a set of cryptographic building blocks used to implement secure protocols evaluating basic variants of the simplex algorithm. We show how to balance between privacy and e??ciency; some protocols reveal data about the internal state of the simplex algorithm, such as the number of iterations, in order to improve the expected running times. For the sake of simplicity and e??ciency, the protocols are based on Shamir's secret sharing scheme. We combine and use the results from the literature on secure random number generation, secure circuit evaluation, secure comparison, and secret indexing to construct e??cient building blocks for secure simplex. The solutions for secure linear programming in this thesis can be split into two categories. On the one hand, some protocols evaluate the classical variants of the simplex algorithm in which numbers are truncated, while the other protocols evaluate the variants of the simplex algorithms in which truncation is avoided. On the other hand, the protocols can be separated by the size of the tableaus. Theoretically there is no clear winner that has both the best security properties and the best performance.
- Research Article
12
- 10.1016/j.cor.2010.07.001
- Jul 10, 2010
- Computers & Operations Research
Hybrid-LP: Finding advanced starting points for simplex, and pivoting LP methods
- Research Article
5728
- 10.2307/253980
- Jul 1, 2000
- The Journal of the Operational Research Society
Introduction and Preliminaries. Problems, Algorithms, and Complexity. LINEAR ALGEBRA. Linear Algebra and Complexity. LATTICES AND LINEAR DIOPHANTINE EQUATIONS. Theory of Lattices and Linear Diophantine Equations. Algorithms for Linear Diophantine Equations. Diophantine Approximation and Basis Reduction. POLYHEDRA, LINEAR INEQUALITIES, AND LINEAR PROGRAMMING. Fundamental Concepts and Results on Polyhedra, Linear Inequalities, and Linear Programming. The Structure of Polyhedra. Polarity, and Blocking and Anti--Blocking Polyhedra. Sizes and the Theoretical Complexity of Linear Inequalities and Linear Programming. The Simplex Method. Primal--Dual, Elimination, and Relaxation Methods. Khachiyana s Method for Linear Programming. The Ellipsoid Method for Polyhedra More Generally. Further Polynomiality Results in Linear Programming. INTEGER LINEAR PROGRAMMING. Introduction to Integer Linear Programming. Estimates in Integer Linear Programming. The Complexity of Integer Linear Programming. Totally Unimodular Matrices: Fundamental Properties and Examples. Recognizing Total Unimodularity. Further Theory Related to Total Unimodularity. Integral Polyhedra and Total Dual Integrality. Cutting Planes. Further Methods in Integer Linear Programming. References. Indexes.
- Research Article
3
- 10.1287/mnsc.2021.3955
- Mar 19, 2021
- Management Science
Optimization has been one of the most fundamental and extensive contributions of management science/operations research, with an enormous number of contributions and subfields developed by many researchers and practitioners. When the journal Management Science launched in 1954, little was known about optimization, including some results in nonlinear optimization and the simplex method and duality developed for linear programming. However, linear programming computations were limited to problems with at most 101 linear constraints. Then some early contributions by seminal researchers began to develop foundations for the field. I will review a few of these early contributions, focusing on the traveling sales problem and integer programming, decomposition, and column generation. I will summarize some research and applied contributions since then, including the enormous development of computations. I will focus on linear and integer programs with some material on combinatorial optimization. This paper was accepted by David Simchi-Levi, Special Section of Management Science: 65th Anniversary.
- Research Article
42
- 10.1016/j.amc.2005.11.134
- Jan 24, 2006
- Applied Mathematics and Computation
Linear bilevel programming with upper level constraints depending on the lower level solution
- Book Chapter
1
- 10.1201/9781003089636-3
- Jan 19, 2022
In Operations Research, linear-fractional programming is considered as the generalization of linear programming problem. While in a linear programming the objective function is a linear function, and in a linear-fractional programming the objective function is the ratio of two linear or non-linear functions. The majority of the algorithm used for solving the linear fractional programming problem relies upon the classical simplex method. In this paper, we have proposed a new algorithm for solving a linear fractional programming problem in which the objective function is a combination of linear fractional function, while constraint functions are in the form of linear inequalities. Our proposed algorithm is based on the extension of the method, which is used to solve linear programming problems with linear constraints. The primary intent behind developing this method is that we did not need to transform the linear fractional programming problem into linear programming problem, and also it helps in finding out the feasible region via a sequence of points in the direction that improves the feasibility of the fractional objective function. Numerical examples are given to illustrate the use of these proposed methods. Lastly, to demonstrate the efficacy of the proposed algorithm, we have compared the findings obtained with other approaches to display our algorithm's efficacy.
- Research Article
8
- 10.1016/s0305-0548(99)00075-1
- Jun 2, 2000
- Computers and Operations Research
Piecewise linear programming via interior points
- Research Article
16
- 10.1002/cjce.5450700222
- Apr 1, 1992
- The Canadian Journal of Chemical Engineering
Multi‐period planning problems in the oil and refinery industry are typically large, sparse, staircase/band diagonal structured and nonlinear optimization problems. Successive linear programming (SLP) type methods have been widely used for solving these planning problems. But, it has long been recognized that the simplex method used in solving linear programs requires a large number of iterations for staircase/band diagonal structured problems. In this paper, we report results of an application of a recently developed interior point method that promises to be many times faster than the simplex method for multi‐period planning problems. However, to facilitate the use of interior point method in the current SLP algorithms a hybrid method combining the interior point method and the simplex method is developed. Therefore, the results determined with this hybrid method are qualitatively equivalent to that obtained with the simplex method alone. The CPU times corresponding to the hybrid method are compared with the CPU times of simplex and dual affine methods. The new hybrid method generates a basic feasible solution of the linear programming problem and is approximately 7 times faster than the simplex method on the tested planning problems. Moreover, the interior point and hybrid methods become faster as the problem size increases.
- Conference Article
- 10.1145/2491845.2491856
- Sep 19, 2013
Linear Programming (LP) is a significant research area in the field of operations research. The simplex algorithm is the most widely used and well-studied method for solving Linear Programming problems (LPs). Many algorithms have been proposed for the solution of LPs. The vast majority of these algorithms belong to three main categories: (i) Simplex-type or pivoting algorithms, (ii) interior-point methods (IPMs) and (iii) exterior point simplex type algorithms (EPSA). The aim of this paper is to present an implementation of a hybrid simplex algorithm that begins to solve the LP using an IPM and after a number of iterations continues with a primal-dual EPSA algorithm. This hybrid approach aims to take advantage of: (i) IPM strengths, which is the fast convergence in the first iterations, and (ii) EPSA strengths, i.e. the fast convergence when making steps in directions that are linear combinations of attractive directions. The idea of combining different types of linear programming algorithms is not new; to the best of our knowledge, this is the first time that interior point methods and exterior point algorithms are combined. The interior point that is calculated by IPM after a number of iterations can lead to such attractive directions. In order to gain an insight into the practical behavior of the proposed algorithm, we have performed some computational experiments over sparse randomly generated optimal LPs. Finally, in the computational study that we have conducted, we investigate the adequate number of iterations that IPM should run in order to decrease the CPU time and the iterations of the proposed algorithm.
- Research Article
- 10.29977/jciie.200605.0005
- May 1, 2006
- Journal of The Chinese Institute of Industrial Engineers
This paper establishes a framework for solving some optimization problems with linear constraints using simplex-type methods. The problems include those found in linear programs, linear fractional programs, and generalized linear fractional programs. In this study, these problems refer to a standard form of minimizing a single parameter subject to parameterized linear equations. Based on the analysis of parameterized basis-based solutions, a unified simplex-type approach is proposed. The adaptability of the parameterized model and that of the solution procedure are discussed. In particular, the proposed algorithm can prevent cycling when compared with the conventional simplex method used for solving linear programs.
- Research Article
29
- 10.1016/j.ejor.2006.02.021
- Apr 1, 2007
- European Journal of Operational Research
A simplex algorithm for piecewise-linear fractional programming problems
- Research Article
194
- 10.1137/0115113
- Sep 1, 1967
- SIAM Journal on Applied Mathematics
: So far the study of stochastic programs with recourse has been limited to the case (called by G. Dantzig programming under uncertainty) when only the right-hand sides or resources of the problem are random. In this paper the authors extend the theory to the general case when essentially all the parameters involved are random. This generalization immediately raises the problem of attributing a precise meaning to the stochastic constraints. They examine a probability formulation (satisfying the constraints almost surely) and a possibility formulation (satisfying the constraints for all values of the random parameters in the support of their joint distribution) and show them equivalent under a rather weak but curious W-condition. Finally, they prove that without restriction the equivalent deterministic form of a stochastic program with recourse is a convex program for which we obtain some additional properties when some of the parameters of the original problem are constant. The applications of the theoretical results of this paper to certain classes of stochastic programs which have arisen from practical problems will be presented in a separate paper: 'Stochastic Programs with Recourse: Special Forms.' (Author)