Computational performance of a projection and rescaling algorithm

This paper documents a computational implementation of a {\em projection and rescaling algorithm} for finding most interior solutions to the pair of feasibility problems \[ \text{find} \; x\in L\cap\mathbb{R}^n_{+} \;\;\;\; \text{ and } \; \;\;\;\; \text{find} \; \hat x\in L^\perp\cap\mathbb{R}^n_{+}, \] where $L$ denotes a linear subspace in $\mathbb{R}^n$ and $L^\perp$ denotes its orthogonal complement. The projection and rescaling algorithm is a recently developed method that combines a {\em basic procedure} involving only low-cost operations with a periodic {\em rescaling step.} We give a full description of a MATLAB implementation of this algorithm and present multiple sets of numerical experiments on synthetic problem instances with varied levels of conditioning. Our computational experiments provide promising evidence of the effectiveness of the projection and rescaling algorithm. Our MATLAB code is publicly available. Furthermore, the simplicity of the algorithm makes a computational implementation in other environments completely straightforward.

[1]  J. G. Wendel A Problem in Geometric Probability. , 1962 .

[2]  Michael J. Todd,et al.  Identifying the set of always-active constraints in a system of linear inequalities by a single linear program , 1985 .

[3]  Felipe Cucker,et al.  Unifying Condition Numbers for Linear Programming , 2003, Math. Oper. Res..

[4]  Santosh S. Vempala,et al.  An Efficient Re-Scaled Perceptron Algorithm for Conic Systems , 2006, Math. Oper. Res..

[5]  F. Cucker,et al.  Coverage processes on spheres and condition numbers for linear programming , 2007, 0712.2816.

[6]  Sergei Chubanov A strongly polynomial algorithm for linear systems having a binary solution , 2012, Math. Program..

[7]  Sergei Chubanov,et al.  A polynomial projection algorithm for linear feasibility problems , 2015, Math. Program..

[8]  T. Terlaky,et al.  A Polynomial Column-wise Rescaling von Neumann Algorithm , 2015 .

[9]  László A. Végh,et al.  Rescaled Coordinate Descent Methods for Linear Programming , 2016, IPCO.

[10]  Javier Peña,et al.  A deterministic rescaled perceptron algorithm , 2016, Math. Program..

[11]  Rebecca Hoberg,et al.  An Improved Deterministic Rescaling for Linear Programming Algorithms , 2016, IPCO.

[12]  Javier Peña,et al.  Solving Conic Systems via Projection and Rescaling , 2015, Math. Program..

[13]  Takashi Tsuchiya,et al.  An extension of Chubanov's polynomial-time linear programming algorithm to second-order cone programming , 2016, Optim. Methods Softw..

[14]  Kees Roos,et al.  An improved version of Chubanov's method for solving a homogeneous feasibility problem , 2018, Optim. Methods Softw..

[15]  Bruno F. Lourenço,et al.  An extension of Chubanov’s algorithm to symmetric cones , 2017, Mathematical Programming.

[16]  James Renegar,et al.  Accelerated first-order methods for hyperbolic programming , 2015, Mathematical Programming.

[17]  László A. Végh,et al.  Rescaling Algorithms for Linear Conic Feasibility , 2016, Math. Oper. Res..