. This paper describes a software package, called LOQO, which implements a primaldual interior-point method for general nonlinear programming. We focus in this paper mainly on the algorithm as it applies to linear and quadratic programming with only brief mention of the extensions to convex and general nonlinear programming, since a detailed paper describing these extensions were published recently elsewhere. In particular, we emphasize the importance of establishing and maintaining symmetric quasidefiniteness of the reduced KKT system. We show that the industry standard MPS format can be nicely formulated in such a way to provide quasidefiniteness. Computational results are included for a variety of linear and quadratic programming problems. 1. INTRODUCTION LOQO is a software package for solving general (smooth) nonlinear optimization problems. It implements an infeasible-primal-dual path-following method. For linear programming, such methods were first proposed independently by Lust...

In the past 10 years the interior point methods (IPM) for linear programming have gained extraordinary interest as an alternative to the sparse simplex based methods. This has initiated a fruitful competition between the two types of algorithms which has lead to very efficient implementations on both sides. The significant difference between interior point and simplex based methods is reflected not only in the theoretical background but also in the practical implementation. In this paper we give an overview of the most important characteristics of advanced implementations of interior point methods. First, we present the infeasible-primal-dual algorithm which is widely considered the most efficient general purpose IPM. Our discussion includes various algorithmic enhancements of the basic algorithm. The only shortcoming of the "traditional" infeasible-primal-dual algorithm is to detect a possible primal or dual infeasibility of the linear program. We discuss how this problem can be solve...

1 Introduction Consider the linear programming (LP) problem in the standard form: (LP) minimize cT x

In this paper, we describe our implementation of a primal-dual infeasible-interior-point algorithm for large-scale linear programming under the MATLAB 1 environment. The resulting software is called LIPSOL -- Linear-programming Interior-Point SOLvers. LIPSOL is designed to take the advantages of MATLAB's sparse-matrix functions and external interface facilities, and of existing Fortran sparse Cholesky codes. Under the MATLAB environment, LIPSOL inherits a high degree of simplicity and versatility in comparison to its counterparts in Fortran or C language. More importantly, our extensive computational results demonstrate that LIPSOL also attains an impressive performance comparable with that of efficient Fortran or C codes in solving large-scale problems. In addition, we discuss in detail a technique for overcoming numerical instability in Cholesky factorization at the end-stage of iterations in interior-point algorithms. Keywords: Linear programming, Primal-Dual infeasible-interior-p...

This work concerns primal-dual interior-point methods for semidefinite programming (SDP) that use a search direction originally proposed by Helmberg-Rendl-Vanderbei-Wolkowicz [5] and Kojima-Shindoh-Hara [11], and recently rediscovered by Monteiro [15] in a more explicit form. In analyzing these methods, a number of basic equalities and inequalities were developed in [11] and also in [15] through different means and in different forms. In this paper, we give a concise derivation of the key equalities and inequalities for complexity analysis along the exact line used in linear programming (LP), producing basic relationships that have compact forms almost identical to their counterparts in LP. We also introduce a new formulation of the central path and variable-metric measures of centrality. These results provide convenient tools for deriving polynomiality results for primal-dual algorithms extended from LP to SDP using the aforementioned and related search directions. We present examples...

This work concerns primal-dual interior-point methods for semidefinite programming (SDP) that use a linearized complementarity equation originally proposed by Kojima, Shindoh and Hara [11], and recently rediscovered by Monteiro [15] in a more explicit form. In analyzing these methods, a number of basic equalities and inequalities were developed in [11] and also in [15] through different means and in different forms. In this paper, we give a very short derivation of the key equalities and inequalities along the exact line used in linear programming (LP), producing basic relationships that have highly compact forms almost identical to their counterparts in LP. We also introduce a new definition of the central path and variable-metric measures of centrality. These results provide convenient tools for extending existing polynomiality results for many, if not most, algorithms from LP to SDP with little complication. We present examples of such extensions, including the long-step infeasible-...

A new approach for the implementation of interior-point methods for solving linear programs is proposed. Its main feature is the iterative solution of the symmetric, but highly indefinite 2\Theta2-block systems of linear equations that arise within the interior-point algorithm. These linear systems are solved by a symmetric variant of the quasi-minimal residual (QMR) algorithm, which is an iterative solver for general linear systems. The symmetric QMR algorithm can be combined with indefinite preconditioners, which is crucial for the efficient solution of highly indefinite linear systems, yet it still fully exploits the symmetry of the linear systems to be solved. To support the use of the symmetric QMR iteration, a novel stable reduction of the original unsymmetric 3 \Theta 3-block systems to symmetric 2 \Theta 2-block systems is introduced, and a measure for a low relative accuracy for the solution of these linear systems within the interior-point algorithm is proposed. Some indefini...

. In this paper we introduce the truncated primal-infeasible dual-feasible interior point algorithm for linear programming and describe an implementation of this algorithm for solving the minimum cost network flow problem. In each iteration, the linear system that determines the search direction is computed inexactly, and the norm of the resulting residual vector is used in the stopping criteria of the iterative solver employed for the solution of the system. In the implementation, a preconditioned conjugate gradient method is used as the iterative solver. The details of the implementation are described and the code, pdnet, is tested on a large set of standard minimum cost network flow test problems. Computational results indicate that the implementation is competitive with state-of-the-art network flow codes. Key Words. Interior point method, linear programming, network flows, primal-infeasible dualfeasible, truncated Newton method, conjugate gradient, maximum flow, experimental test...

. In this paper, we propose primal-dual potential-reduction algorithms which can start from an infeasible interior point. We first describe two such algorithms and show that both are polynomial-time bounded. One of the algorithms decreases the Tanabe-Todd-Ye primal-dual potential function by a constant at each iteration under the condition that the duality gap decreases by at most the same ratio as the infeasibility. The other reduces a new potential function, which has one more term in the Tanabe-Todd-Ye potential function, by a fixed constant at each iteration without any other conditions on the step size. Finally, we describe modifications of these methods (incorporating centering steps) which dramatically decrease their computational complexity. Our algorithms also extend to the case of monotone linear complementarity problems. Key words. Polynomial Time, Linear Programming, Primal-Dual, Infeasible-Interior-Point Algorithm, Potential Function. AMS subject classifications. 90C05, ...

We present a generalization of a homogeneous self-dual linear programming (LP) algorithm to solving the monotone complementarity problem (MCP). The algorithm does not need to use any "big-M" parameter or two-phase method, and it generates either a solution converging towards feasibility and complementarity simultaneously or a certificate proving infeasibility. Moreover, if the MCP is polynomially solvable with an interior feasible starting point, then it can be polynomially solved without using or knowing such information at all. To our knowledge, this is the first interior-point and infeasible-starting algorithm for solving the MCP that possesses these desired features. Preliminary computational results are presented. Key words: Monotone complementarity problem, homogeneous and self-dual, infeasible-starting algorithm. Running head: A homogeneous algorithm for MCP. Department of Management, Odense University, Campusvej 55, DK-5230 Odense M, Denmark, email: eda@busieco.ou.dk. y De...