An Efficient Parametric Linear Programming Solver and Application to Polyhedral Projection

Yu, Hang; Monniaux, David

doi:10.1007/978-3-030-32304-2_11

Cited by 6 publications

(4 citation statements)

References 16 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…It can be shown that the region of the ( 1 , … , k ) associated to a given optimum X * is a convex polyhedron (for C 0 , a convex polyhedral cone), and that these regions form a quasi-partition of the space of parameters (two regions may overlap at their boundary, but not in their interior). 14,19,27 The output of the parametric linear programming solver is this quasi-partition, 11 There exists another, dual, kind of parametric linear programming where the parameters are in the right-hand side B. We do not consider it here.…”

Section: Parametric Linear Programmingmentioning

confidence: 99%

“…We shall leave out here how polyhedral computations such as projection and convex hull can be reduced to parametric linear programming—this is covered in the literature 19,27 —and focus on solving the parametric linear programs.…”

Section: Sequential Algorithmsmentioning

confidence: 99%

“…Then, for given values of the

μ_{i}

, the problem is either unbounded, or there is one single optimal solution X ∗ . It can be shown that the region of the

(μ_{1}, \dots, μ_{k})

associated to a given optimum X ∗ is a convex polyhedron (for C 0 , a convex polyhedral cone), and that these regions form a quasi‐partition of the space of parameters (two regions may overlap at their boundary, but not in their interior) 14,19,27 . The output of the parametric linear programming solver is this quasi‐partition, and the associated optima—in our applications, the problem is always bounded in the optimization directions, so we do not deal with the unbounded case.…”

Section: Sequential Algorithmsmentioning

confidence: 99%

See 2 more Smart Citations

A task‐based approach to parallel parametric linear programming solving, and application to polyhedral computations

Coti

Monniaux

2020

Concurrency and Computation

Self Cite

View full text Add to dashboard Cite

Parametric linear programming is a central operation for polyhedral computations, as well as in certain control applications. Here, we propose a task-based scheme for parallelizing it, with quasi-linear speedup over large problems. This type of parallel applications is challenging, because several tasks might be computing the same region. In this article, we are presenting the algorithm itself with a parallel redundancy elimination algorithm, and conducting a thorough performance analysis. K E Y W O R D S parallel parametric linear programming, polyhedra 1 INTRODUCTION A convex polyhedron in dimension n is the solution set over Q n (or, equivalently, 1 R n) of a system of inequalities (with integer or rational coefficients). Since all the polyhedra we consider here are convex, we shall talk of polyhedra for short. There also exist integer polyhedra, defined over Z n , but they are very different in many respects; we shall not consider them here. 2 Computations over polyhedra arise in low dimension (n = 2 and n = 3) for modeling physical objects, but here we are interested in higher dimensions. Polyhedra in higher dimension are typically used to enclose the set of reachable states of systems whose state can be expressed, at least partially, as a vector of reals or rationals. For instance, one can study the flow of ordinary differential equations, or more generally the trajectories of hybrid systems, by enclosing the trajectories into a succession of convex polyhedra. If the internal state of a control system is expressed by a vector of n numeric variables, one may prove that this system never encounters a bad condition by exhibiting a polyhedron P such that the initial state belongs to P, no bad condition belongs to P, and all possible time steps of the control system leave P stable (i.e., it is not possible to execute a step starting in P and ending outside of P). One generally proves the correctness of programs by exhibiting inductive invariants-an inductive invariant for a loop is a set containing the precondition of the loop, stable by moving to the next loop iteration, and implying the desired postcondition. Since providing inductive invariants by hand is tedious, it is desirable to automate that process. One approach for doing so is abstract interpretation, where an ascending sequence of sets of states is computed until reaching a fixed point. One may, for instance, search such sets as products of intervals, an approach known as interval analysis, but the lack of relationships between the dimensions tends to severely limit the kind of properties that can be proved (e.g., one cannot have an invariant i < n where n is a parameter: this is neither an interval on i nor on n nor a combination thereof). Cousot and Halbwachs proposed searching for polyhedral inductive invariants. 1,2 The operations needed there are projection, more generally image by an affine transformation, convex hull, and inclusion (or equality) test, together with an extrapolation operator known as widening. 1 Whether emptiness tests, inclusio...

show abstract

Section: Parametric Linear Programmingmentioning

confidence: 99%

Section: Sequential Algorithmsmentioning

confidence: 99%

“…Then, for given values of the

μ_{i}

, the problem is either unbounded, or there is one single optimal solution X ∗ . It can be shown that the region of the

(μ_{1}, \dots, μ_{k})

Section: Sequential Algorithmsmentioning

confidence: 99%

See 1 more Smart Citation

A task‐based approach to parallel parametric linear programming solving, and application to polyhedral computations

Coti

Monniaux

2020

Concurrency and Computation

Self Cite

View full text Add to dashboard Cite

show abstract

“…The methods for doing so rely on lexicographic orderings or perturbations on the objective and/or constant term, or pivoting rules [20]. We have recently proposed a working solution to degeneracy [31] but there is still room for improvement.…”

Section: Parametric Linear Programmingmentioning

confidence: 99%