Abstract
This paper proposes a novel approach for solving linear programs. We reformulate a primal-dual linear program as an unconstrained minimization of a convex and twice continuously differentiable merit function. When the optimal set of the primal-dual pair is nonempty, its optimal set is equal to the optimal set of the proposed merit function. Minimizing this merit function poses some challenges due to its Hessian being singular at some points in the domain, including the optimal solutions. We handle singular Hessians using the Newton method with Levenberg-Marquardt regularization. We show that the Newton method with Levenberg-Marquardt regularization yields global convergence to a solution of the primal-dual linear program in at most iterations requiring only the assumption that the optimal set of the primal-dual linear program is bounded. Testing on five NETLIB LP problems demonstrates promising results: for two instances, the method converges to very high accuracy, significantly faster than predicted by the derived worst-case bound.