Separable Concave Optimization Approximately Equals Piecewise-Linear Optimization

We study the problem of minimizing a nonnegative separable concave function over a compact feasible set. We approximate this problem to within a factor of 1+epsilon by a piecewise-linear minimization problem over the same feasible set. Our main result is that when the feasible set is a polyhedron, the number of resulting pieces is polynomial in the input size of the polyhedron and linear in 1/epsilon. For many practical concave cost problems, the resulting piecewise-linear cost problem can be formulated as a well-studied discrete optimization problem. As a result, a variety of polynomial-time exact algorithms, approximation algorithms, and polynomial-time heuristics for discrete optimization problems immediately yield fully polynomial-time approximation schemes, approximation algorithms, and polynomial-time heuristics for the corresponding concave cost problems. We illustrate our approach on two problems. For the concave cost multicommodity flow problem, we devise a new heuristic and study its performance using computational experiments. We are able to approximately solve significantly larger test instances than previously possible, and obtain solutions on average within 4.27% of optimality. For the concave cost facility location problem, we obtain a new 1.4991+epsilon approximation algorithm.Comment: Full pape

Separable Concave Optimization Approximately Equals Piecewise Linear Optimization

We show how to approximate a separable concave minimization problem over a general closed ground set by a single piecewise linear minimization problem. The approximation is to arbitrary 1 + ɛ precision in optimal cost. For polyhedral ground sets in R n + and nondecreasing cost functions, the number of pieces is polynomial in the input size and proportional to 1 / log(1 + ɛ). For general polyhedra, the number of pieces is polynomial in the input size and the size of the zeroes of the concave objective components. We illustrate our approach on the concave-cost uncapacitated multicommodity flow problem. By formulating the resulting piecewise linear approximation problem as a fixed charge, mixed integer model and using a dual ascent solution procedure, we solve randomly generated instances to within five to twenty percent of guaranteed optimality. The problem instances contain up to 50 nodes, 500 edges, 1,225 commodities, and 1,250,000 flow variables.

Global optimization with piecewise linear approximation

textGlobal optimization deals with the development of solution methodologies for nonlinear nonconvex optimization problems. These problems, which could arise in diverse situations ranging from optimizing hydro-power generation schedules to estimating coefficients of non-linear regression models, are difficult for traditional nonlinear solvers that iteratively search the neighborhood around a starting point. The Piecewise Linear Approximation (PLA) method that we study in this dissertation seeks to generate ‘good’ starting points, hopefully ones that lie in the basin of attraction of the globally optimal solution. In this approach, we approximate the non-linear functions in the optimization problem by piecewise linear functions defined over the vertices of a grid that partitions the domain of each nonlinear function into cells. Based on this approximation, we convert the original nonlinear program into a mixed integer program (MIP) and use the solution to this MIP as a starting point for a local nonlinear solver. In this dissertation, we validate the effectiveness of the PLA approach as a global optimization approach by applying it to a diverse set of continuous and discrete nonlinear optimization problems. Further, we develop various modeling and algorithmic strategies for enhancing the basic approach. Our computational results demonstrate that the PLA approach works well on non-convex problems and can, in some cases, provide better solutions than those provided by existing nonlinear solvers.Information, Risk, and Operations Management (IROM