Integer Programming Reformulations: Dantzig-Wolfe & Benders Decomposition the Coluna Software Platform
|
|
- Ursula Jackson
- 5 years ago
- Views:
Transcription
1 Integer Programming Reformulations: Dantzig-Wolfe & Benders Decomposition the Coluna Software Platform François Vanderbeck B. Detienne, F. Clautiaux, R. Griset, T. Leite, G. Marques, V. Nesello, A. Pessoa, E. Uchoa, R. Sadykov, H. Sen, I. Tahiri, L. Wolsey ROADEF, February 2019 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 1 / 54
2 Formulation Combinatorial Optimization Problem (CO) min{c(s) : s S} where S is the discrete set of feasible solutions. (7) (6) (5) S (4) François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 2 / 54
3 Formulation Combinatorial Optimization Problem (CO) min{c(s) : s S} where S is the discrete set of feasible solutions. Integer Programming Formulation A polyhedron P = {x n : Ax a} is a formulation for (CO) iff min{c(s) : s S} min{cx : x P I = P n }. x2 x2 (7) 5 5 (6) (5) 3 (7) 3 (7) (6) (5) (6) (5) S (4) 2 (4) 2 (4) x x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 2 / 54
4 Alternative formulations A formulation is typically not unique P and P can be alternative formulations for (CO) if (CO) min{cx : x P n } min{c x : x P n } x2 5 3 (7) (6) (5) 2 1 (4) x1 More generally: can expressed in different variable-spaces. François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 3 / 54
5 Quality of Formulations Stronger formulation (in the same space) Formulation P n is a stronger than P n if P P. Then, min{cx : x P } min{cx : x P} x2 5 3 (7) (6) (5) 2 1 (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 4 / 54
6 Extended Formulation Given an initial compact formulation: x2 5 3 (7) (6) (5) 2 1 (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 5 / 54
7 Extended Formulation w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 6 / 54
8 Extended Formulation w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 7 / 54
9 Extended Formulation w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 8 / 54
10 Projection The Projection of Q = {(x, w) n+e : Gx + Hw d} on the x-space is: proj x (Q) := {x n : w e such that(x, w) Q}. w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 9 / 54
11 Take away messages Integer Programming Reformulations A powerful way to exploit the combinatorial structure. Benders & Dantzig Wolfe decomposition are generic schemes to derive & handle reformulations Size can be coped with using dynamic generation: a small % of variables and constraints are needed; hence it scales up to real-life applications. With efficiency enhancement features, these approaches can be highly competitive Implementation can be generic, with tools such as Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 10 / 54
12 Outline 1 Sub-gradient linearization : Dantzig & Benders 2 Solution Methods 3 Efficient Implementation 4 Integer Solutions 5 Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 11 / 54
13 Subgradient linearization of a convex function f (x) max f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 12 / 54
14 Subgradient linearization of convex functions η max η η f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 13 / 54
15 Subgradient linearization of convex functions η max η η f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 13 / 54
16 Subgradient linearization of convex functions η max η η f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 13 / 54
17 Subgradient linearization of convex functions η max η η f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 13 / 54
18 Subgradient linearization of convex functions η max η η f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 13 / 54
19 Subgradient linearization of convex functions η max η η f(x) g(x) 0 x 0 where f (.) and/or g(.) are defined by a bundle of linear pieces that are generated dynamically during the optimization process. Separation: Evaluate the true fonction f at the current point x yields a proof of optimality or a violated constraint. x François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 13 / 54
20 Outline 1 Sub-gradient linearization : Dantzig & Benders 2 Solution Methods 3 Efficient Implementation 4 Integer Solutions 5 Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 14 / 54
21 Dantzig-Wolfe Decomposition Assume a bounded integer integer problem with structure: [F] min c x : A x a x X = { B x b x n } François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 15 / 54
22 Dantzig-Wolfe Decomposition Assume a bounded integer integer problem with structure: [F] min c x : A x a x X = { B x b where B is typically block diagonal x n } [F B ] min c 1 x 1 + c 2 x c K x K A 1 x 1 + A 2 x A K x K a B 1 x 1 b 1 B 2 x 2 b B K x K x 1 p 1 x 2 p 2..., x K p K b K François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 15 / 54
23 Dantzig-Wolfe Decomposition Assume a bounded integer integer problem with structure: [F] min c x : A x a x X = { B x b x n } Assume that subproblem [SP] min{c x : x X} (1) is relatively easy to solve compared to problem [F]. Then, X = {x q } q Q François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 15 / 54
24 Dantzig-Wolfe Decomposition Assume a bounded integer integer problem with structure: [F] min c x : A x a x X = { B x b x n } Assume that subproblem [SP] min{c x : x X} (1) is relatively easy to solve compared to problem [F]. Then, X = {x q } q Q conv(x) = {x n + : x = x q λ q, λ q = 1, λ q 0 q Q} q Q q Q François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 15 / 54
25 Dantzig-Wolfe Decomposition Assume a bounded integer integer problem with structure: [F] min c x : A x a x X = { B x b x n } Assume that subproblem [SP] min{c x : x X} (1) is relatively easy to solve compared to problem [F]. Then, X = {x q } q Q conv(x) = {x n + : x = x q λ q, λ q = 1, λ q 0 q Q} q Q q Q [F] min{ cx q λ q : Ax q λ q a, λ q = 1, x q λ q n } q Q q Q q Q q François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 15 / 54
26 The Cutting Stock Problem (CSP) 7 x waste : 0 min x waste : x waste : x waste : 0 8 x waste : 0 4 x waste : 48 5 x waste : 54 4 x waste : x waste : 50 6 x waste : x waste : 0 s.t. K k=1 y k K x ik d i i k=1 w i x ik W y k k i x ik i, k y k {0, 1} k 8 x waste : x waste : x waste : x waste : x waste : 44 2 x waste : 368 Number of distinct cutting patterns : 17 Number of rolls : 88 Total waste : 5090 min q λ q x q i λ q d i i = 1,, n q λ q 0 q François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 16 / 54
27 The Vehicle Routing Problem (VRP) min s.t. c e x r e e,r e δ(i),r x r e 2 i x r X r r 6, 4, 3 min e,r,q Q r c e x q e λ q e δ(i),r,q Q r x q e λ q 2 i q Q r λ q U r r λ q 0 r, q Q r François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 17 / 54
28 Lagrangian Duality [F] = min {c x : A x a; B x b; x n } L(π) = min q Q {c xq + π (a Ax q )} [LD] = max L(π) π m + η η L(π) π o π π François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 18 / 54
29 The Lagrangian Dual as an LP η [LD] max min {π a + (c π m + q Q πa)xq }; max{η, η cx q + π(a Ax q ) q Q, π m +, η 1 }; η L(π) π L(π) π François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 19 / 54
30 The Lagrangian Dual as an LP η [LD] max min {π a + (c π m + q Q πa)xq }; max{η, η cx q + π(a Ax q ) q Q, π m +, η 1 }; min{ (cx q )λ q, q Q (Ax q )λ q a, q Q λ q = 1, q Q λ q 0 q Q}; η L(π) L(π) François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 19 / 54 π π
31 Cut generation in the dual L(u) x2 x1 L(u) x2 x1 x3 x3 x4 x4 x5 x5 u u L(u) x2 x1 L(u) x2 x1 x3 x3 x4 x4 x5 x5 u u François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 20 / 54
32 Column generation in the primal Xlp Xlp Xlp Xlp François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 21 / 54
33 Properties Restricted Master: DUAL DM t max L t () where L t (π) = min q Q t{πa + (c πa)x q } PRIMAL PM t min{cx : Ax a, x conv({x q } q Q t)}. Valid upper bound: UB η L(π) (π t, η t ) π François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 22 / 54
34 Properties Restricted Master: DUAL DM t max L t () where L t (π) = min q Q t{πa + (c πa)x q } PRIMAL PM t min{cx : Ax a, x conv({x q } q Q t)}. Valid upper bound: UB η L(π) (π t, η t ) SubProb: Computing L(π t ) yields a 1 most violated constr. for DM t 2 most neg. red. cost col. for PM t 3 a correction on L t (.) at point π t 4 a sub-gradient g t = (a A x t ) of L(.) π Valid lower bound: LB François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 22 / 54
35 The benefits of the approach Better Dual Bound: Lagrangian duality Mastering the size of the formulation: projection w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 23 / 54
36 The benefits of the approach Better Dual Bound: Lagrangian duality Mastering the size of the formulation: projection Divide and conquer & Allow to Exploit problem structure complexity size François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 23 / 54
37 The benefits of the approach Better Dual Bound: Lagrangian duality Mastering the size of the formulation: projection Divide and conquer & Allow to Exploit problem structure Avoid symmetry between identical blocks & Allow a compact representation of a solution min cx 1 + c x c x K A x 1 + A x A x K a B x 1 b 1 B x 2 b B x K x 1 p 1 x 2 p 2..., x K p K If identical subproblems, solution x q is taken λ q times b K François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 23 / 54
38 Dual Decompositions On a subset of constraints (Dantzig-Wolfe): multi-resources min c 1 x 1 + c 2 x c K x K A x 1 + A x A x K a B 1 x 1 b 1 B 2 x 2 b B K x K b K On a subset of variables (Benders): multi-decision-levels min a x + b 1 y 1 + b 2 y b K y K A x + B 1 y 1 c 1 A x + B 2 y 2 c 2 A x A x + B K y K c K François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 24 / 54
39 Capacitated Facility Location P fx + Pi i i s.t. y Pj ij i yij yij xi min François Vanderbeck (U. de Bordeaux & Atoptima) P i,j cij yij si xi i dj j 0 i, j {0, 1}, i Integer Programming Reformulations 25 / 54
40 Multi Commodity Network Design min i,j f ijx ij + i,j,k c ijky ijk s.t. k y ijk u ij x ij i, j j y ijk j y jik = b ik i, k y ijk 0 i, j, k x ij {0, 1} i, j François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 26 / 54
41 Benders Decomposition max a x + b 1 y 1 + b 2 y b K y K A x + B 1 y 1 c 1 A x + B 2 y 2 c 2 A x A x + B K y K c K max { a x + b y : A x + B y c, x n, y p + }, François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 27 / 54
42 Benders Decomposition max { a x + b y : A x + B y c, φ(x) := max {b y : B y c A x, x n, y p + }, y p + } First stage: max x,η {η : η a x + φ(x), x n, η }, François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 28 / 54
43 Benders Decomposition max { a x + b y : A x + B y c, φ(x) := max {b y : B y c A x, x n, y p + }, y p + } First stage: max x,η {η : η a x + φ(x), x n, η }, Second stage: min ν 0? b y + ν (η a x) B y + 1 ν (A x c) (y, ν) p+1 + François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 28 / 54
44 Benders Decomposition max { a x + b y : A x + B y c, φ(x) := max {b y : B y c A x, x n, y p + }, y p + } First stage: max x,η {η : η a x + φ(x), x n, η }, u q 0 (η a x) + uq (A x c) 0 q Second stage: min ν 0? b y + ν (η a x) B y + 1 ν (A x c) (y, ν) p+1 + max u 0 (η a x) + u (A x c) 0? u 0 b u B 0 u 0 + 1x 1 (u 0, u) 1+m 2 + {u 1,..., u T } are its extreme dual solutions. François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 28 / 54
45 Benders Decomposition max { a x + b y : A x + B y c, φ(x) := max {b y : B y c A x, x n, y p + }, y p + } First stage: Relaxed Master UB max x,η {η : η a x + φ(x), x n, η }, Second stage: min ν 0? b y + ν (η a x) B y + 1 ν (A x c) (y, ν) p+1 + u q 0 (η a x) + uq (A x c) 0 q Sub-Problem Lagrangian LB max u 0 (η a x) + u (A x c) 0? u 0 b u B 0 u 0 + 1x 1 (u 0, u) 1+m 2 + {u 1,..., u T } are its extreme dual solutions. François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 28 / 54
46 Benders Relaxed Master PRIMAL: [M t ] UB max x,η {η : x q 0 η + (xq A x q 0 a) x cxq q Q t x n +, η }, (x t, η t ) DUAL: [DM t ] min{ (c x q ) λ q : q Q t (x q A x q 0 a) λ q 0, q Q t x q 0 λ q = 1, q λ q 0 q Q t }. (π t, σ t ) = q (uq, u q 0 )λ q SP: L(x, η) = η + min q {c u q (u q A u q 0 a) x uq 0 η} u, LB François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 29 / 54
47 Outline 1 Sub-gradient linearization : Dantzig & Benders 2 Solution Methods 3 Efficient Implementation 4 Integer Solutions 5 Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 30 / 54
48 Method s Classification 1 Ascent: max π {L(π)} Subgradient Volume Conjugate SubGradient Subgradient: π t+1 = (π t + s t η L(π) g t g t )+ + simple unconstraint optimization + low memory use - myoptic : use local info. - slow convergence tail π François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 31 / 54
49 Method s Classification 1 Ascent: Subgradient Volume Conjugate SubGradient f (x) 2 Polyhedral - Active Set: Kelley Bundle x + use global info + restricted formulation - combinatorics of active set - oscillation bt extreme sol. François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 31 / 54
50 Method s Classification 1 Ascent: Subgradient Volume Conjugate SubGradient 2 Polyhedral - Active Set: Kelley Bundle 3 Polyhedral - Interior Point: Analytic Center Cutting Plane Primal Dual Col Gen + nb iterations less sensitive to size - full size solution François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 31 / 54
51 Comparing Methods Best Methods? Ascent methods are computationally light but only allow a rough optimization (good for: warm start, huge problems [Nesterov-14]) Active set methods are best when combinatoric is low Extreme solutions are good to yield integer solutions Otherwise, Interior Point approaches (ACCPM, PDCGM) take over Quality of the implementation François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 32 / 54
52 Outline 1 Sub-gradient linearization: Dantzig & Benders 2 Solution Methods 3 Efficient Implementation 4 Integer Solutions 5 Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 33 / 54
53 Convergence Issue {π t } t π Π Oscillations: π t (bang-bang), {L(π t )} (yo-yo) and π t π > π t 1 π Tailing-off effect: cuts co-linear with objective (inherent to cutting plane proc.) step size very small. Multiple optima François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 34 / 54
54 Stabilization - Acceleration Techniques 1 Proximality of a stability center (to avoid oscillations): Volume: π t = ( ˆπ + s t dt d t )+. Bundle: π t = argmax π {πa + σ 1 2θ π ˆπ 2 } ( ˆπ, ˆL) L(π) François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 35 / 54
55 Stabilization - Acceleration Techniques 1 Proximality of a stability center (to avoid oscillations): 2 Central Point Cutting (to avoid tailing and degeneracy): Optimal face center Chebyshev center Analytic center (of trust polyhedra) Interior Point LP-Solver η = η t ( ˆπ, ˆL) (π t, η t ) η ˆL η ˆL ( ˆπ, ˆL) François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 35 / 54
56 Stabilization - Acceleration Techniques 1 Proximality of a stability center (to avoid oscillations): 2 Central Point Cutting (to avoid tailing and degeneracy): 3 Smoothing (to avoid oscillations): Primal Smoothing Volume: d t = α d t 1 + (1 α) g t Dual Smoothing (also contributes to centralize) [Neame-99] π t = α π t 1 + (1 α) π t [Wentges-97] π t = α ˆπ + (1 α) π t G x y0 G x(a)=x a G OUT SEP d y2 y1 IN f(x)=c1 f(x)=c3<c2 f(x)=c2<c1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 35 / 54
57 Outline 1 Sub-gradient linearization: Dantzig & Benders 2 Solution Methods 3 Efficient Implementation 4 Integer Solutions 5 Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 36 / 54
58 Integrality enforcement in Branch-and-Price w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 37 / 54
59 Integrality enforcement in Branch-and-Price w x2 (5) (4) x1 François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 38 / 54
60 Integrality enforcement in Branch-and-Price Assume a bounded integer integer problem with structure: Assume that subproblem [F] min c x : A x a x Z = { B x b x n } [SP] min{c x : x Z} (2) with conv(z) = {x n + : x = x q λ q, λ q = 1, λ q 0 q Q} q Q q Q The reformulation of [F] is : min{ cx q λ q : Ax q λ q a, λ q = 1, q Q q Q q Q q xq λ q n } François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 39 / 54
61 Integrality enforcement in Branch-and-Price Assume a bounded integer integer problem with structure: Assume that subproblem [F] min c x : A x a x Z = { B x b x n } [SP] min{c x : x Z} (2) with conv(z) = {x n + : x = x q λ q, λ q = 1, λ q 0 q Q} q Q q Q The reformulation of [F] is : min{ cx q λ q : Ax q λ q a, λ q = 1, q Q But in the case of identical sub-problems... q Q q Q q xq λ q n } François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 39 / 54
62 Integrality enforcement in B&P: identical SP A Mapping that partitions the aggregate SP solutions: λ 1 = 1 2 λ 2 = 1 2 λ 3 = 1 2 λ 4 = x 1 = x 2 = François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 40 / 54
63 Integrality enforcement in B&P: identical SP A Mapping that partitions the aggregate SP solutions: λ q x x x x X 1 = {x X : x 1 1} and X 12 = {x X : x 1 1, x 2 1} Branching: [V, MP 2009] q X 12 λ q 2 or q (X 1 \X 12 ) λ q 2 or q (X\X 1 ) λ q 3 dynamically increases the # of SP, up to K non-empty nested SP François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 41 / 54
64 Integer Benders Algorithm (text-book) Oracle Multi-cuts: separation for each integer solution of the pool François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 42 / 54
65 Benders using Lazy Constraints Oracle Oracle Oracle Oracle Oracle François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 43 / 54
66 Benders using Cut-and-Branch Oracle Oracle Oracle Oracle Oracle Oracle François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 44 / 54
67 Benders using Branch-and-Cut Oracle Oracle Oracle Oracle Oracle Oracle Oracle Oracle Oracle François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 45 / 54
68 Outline 1 Sub-gradient linearization: Dantzig & Benders 2 Solution Methods 3 Efficient Implementation 4 Integer Solutions 5 Coluna.jl François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 46 / 54
69 Software What: Uses: Coluna.jl a Column-and-Row generation code Open-Source: MPLicence in Julia 1.0 Extending JuMP / JuliaOpt and BlockDecomposition.jl Dantzig-Wolfe & Benders decomposition Robust & Stochastic Optimization Machine Learning First release: Jan. 19 Support: Math Optimization Society (MOS) & JuliaOpt URL François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 47 / 54
70 Coluna Environment François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 48 / 54
71 Cutting Stock Problem: identical subproblems François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 49 / 54
72 Vehicle Routing Problem: Q-route subproblem oracle François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 50 / 54
73 Vehicle Routing Problem: Cut & Pricing oracles François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 51 / 54
74 Production Problem: combining Dantzig & Benders François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 52 / 54
75 Coluna s Future Present Modeling JuMP Extension Automated Reformul. Algorithms Column Generation Pricing Call-backs Branch-and-Price Preprocessing Primal Heuristics Applications RCSP Oracle VRP Application (MIT L) François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 53 / 54
76 Coluna s Future Present Coming Soon Modeling JuMP Extension Automated Reformul. Algorithms Column Generation Pricing Call-backs Branch-and-Price Preprocessing Primal Heuristics Applications RCSP Oracle VRP Application (MIT L) Efficiency Stabilization Strong Branching Further Primal Heuristics Pricing Strategies Parallelisation Cutomization Cutting Callback Branching Callback Functionality Ad-hoc Branching Benders Decomposition Nested Decomposition François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 53 / 54
77 Coluna & the community Welcome β-users Integration Algorithms: Subgradient, Interior Point, Primal Heuristics Oracles: Pricing, Separation Code review Co-development François Vanderbeck (U. de Bordeaux & Atoptima) Integer Programming Reformulations 54 / 54
Stabilization in Column Generation: numerical study
1 / 26 Stabilization in Column Generation: numerical study Artur Pessoa 3 Ruslan Sadykov 1,2 Eduardo Uchoa 3 François Vanderbeck 2,1 1 INRIA Bordeaux, France 2 Univ. Bordeaux I, France 3 Universidade Federal
More informationReformulation and Decomposition of Integer Programs
Reformulation and Decomposition of Integer Programs François Vanderbeck 1 and Laurence A. Wolsey 2 (Reference: CORE DP 2009/16) (1) Université Bordeaux 1 & INRIA-Bordeaux (2) Université de Louvain, CORE.
More informationColumn Generation for Extended Formulations
1 / 28 Column Generation for Extended Formulations Ruslan Sadykov 1 François Vanderbeck 2,1 1 INRIA Bordeaux Sud-Ouest, France 2 University Bordeaux I, France ISMP 2012 Berlin, August 23 2 / 28 Contents
More informationExtended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications
Extended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications François Vanderbeck University of Bordeaux INRIA Bordeaux-Sud-Ouest part : Defining Extended Formulations
More informationPedro Munari - COA 2017, February 10th, University of Edinburgh, Scotland, UK 2
Pedro Munari [munari@dep.ufscar.br] - COA 2017, February 10th, University of Edinburgh, Scotland, UK 2 Outline Vehicle routing problem; How interior point methods can help; Interior point branch-price-and-cut:
More informationLecture 9: Dantzig-Wolfe Decomposition
Lecture 9: Dantzig-Wolfe Decomposition (3 units) Outline Dantzig-Wolfe decomposition Column generation algorithm Relation to Lagrangian dual Branch-and-price method Generated assignment problem and multi-commodity
More informationFeasibility Pump Heuristics for Column Generation Approaches
1 / 29 Feasibility Pump Heuristics for Column Generation Approaches Ruslan Sadykov 2 Pierre Pesneau 1,2 Francois Vanderbeck 1,2 1 University Bordeaux I 2 INRIA Bordeaux Sud-Ouest SEA 2012 Bordeaux, France,
More informationOutline. Relaxation. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Lagrangian Relaxation. Lecture 12 Single Machine Models, Column Generation
Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING 1. Lagrangian Relaxation Lecture 12 Single Machine Models, Column Generation 2. Dantzig-Wolfe Decomposition Dantzig-Wolfe Decomposition Delayed Column
More informationAcceleration and Stabilization Techniques for Column Generation
Acceleration and Stabilization Techniques for Column Generation Zhouchun Huang Qipeng Phil Zheng Department of Industrial Engineering & Management Systems University of Central Florida Sep 26, 2014 Outline
More informationLagrangian Relaxation in MIP
Lagrangian Relaxation in MIP Bernard Gendron May 28, 2016 Master Class on Decomposition, CPAIOR2016, Banff, Canada CIRRELT and Département d informatique et de recherche opérationnelle, Université de Montréal,
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective Vehicle routing problem 1 / 33 Cutting stock
More informationBranching in Branch-and-Price: a Generic Scheme
Branching in Branch-and-Price: a Generic Scheme François Vanderbeck To cite this version: François Vanderbeck. Branching in Branch-and-Price: a Generic Scheme. Mathematical Programming, Series A, Springer,
More informationColumn generation for extended formulations
EURO J Comput Optim (2013) 1:81 115 DOI 10.1007/s13675-013-0009-9 ORIGINAL PAPER Column generation for extended formulations Ruslan Sadykov François Vanderbeck Received: 5 May 2012 / Accepted: 21 December
More informationA generic view of Dantzig Wolfe decomposition in mixed integer programming
Operations Research Letters 34 (2006) 296 306 Operations Research Letters www.elsevier.com/locate/orl A generic view of Dantzig Wolfe decomposition in mixed integer programming François Vanderbeck a,,
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective 1 / 24 Cutting stock problem 2 / 24 Problem description
More informationColumn Generation. MTech Seminar Report. Soumitra Pal Roll No: under the guidance of
Column Generation MTech Seminar Report by Soumitra Pal Roll No: 05305015 under the guidance of Prof. A. G. Ranade Computer Science and Engineering IIT-Bombay a Department of Computer Science and Engineering
More informationA Column Generation Based Heuristic for the Dial-A-Ride Problem
A Column Generation Based Heuristic for the Dial-A-Ride Problem Nastaran Rahmani 1, Boris Detienne 2,3, Ruslan Sadykov 3,2, François Vanderbeck 2,3 1 Kedge Business School, 680 Cours de la Libération,
More informationColumn Generation. ORLAB - Operations Research Laboratory. Stefano Gualandi. June 14, Politecnico di Milano, Italy
ORLAB - Operations Research Laboratory Politecnico di Milano, Italy June 14, 2011 Cutting Stock Problem (from wikipedia) Imagine that you work in a paper mill and you have a number of rolls of paper of
More informationAn Integer Cutting-Plane Procedure for the Dantzig-Wolfe Decomposition: Theory
An Integer Cutting-Plane Procedure for the Dantzig-Wolfe Decomposition: Theory by Troels Martin Range Discussion Papers on Business and Economics No. 10/2006 FURTHER INFORMATION Department of Business
More informationOn the exact solution of a large class of parallel machine scheduling problems
1 / 23 On the exact solution of a large class of parallel machine scheduling problems Teobaldo Bulhões 2 Ruslan Sadykov 1 Eduardo Uchoa 2 Anand Subramanian 3 1 Inria Bordeaux and Univ. Bordeaux, France
More informationInteger program reformulation for robust branch-and-cut-and-price
Integer program reformulation for robust branch-and-cut-and-price Marcus Poggi de Aragão Informática PUC-Rio Eduardo Uchoa Engenharia de Produção Universidade Federal Fluminense Outline of the talk Robust
More informationDecomposition-based Methods for Large-scale Discrete Optimization p.1
Decomposition-based Methods for Large-scale Discrete Optimization Matthew V Galati Ted K Ralphs Department of Industrial and Systems Engineering Lehigh University, Bethlehem, PA, USA Départment de Mathématiques
More informationDecomposition and Reformulation in Integer Programming
and Reformulation in Integer Programming Laurence A. WOLSEY 7/1/2008 / Aussois and Reformulation in Integer Programming Outline 1 Resource 2 and Reformulation in Integer Programming Outline Resource 1
More informationInterior-Point versus Simplex methods for Integer Programming Branch-and-Bound
Interior-Point versus Simplex methods for Integer Programming Branch-and-Bound Samir Elhedhli elhedhli@uwaterloo.ca Department of Management Sciences, University of Waterloo, Canada Page of 4 McMaster
More informationThe Fixed Charge Transportation Problem: A Strong Formulation Based On Lagrangian Decomposition and Column Generation
The Fixed Charge Transportation Problem: A Strong Formulation Based On Lagrangian Decomposition and Column Generation Yixin Zhao, Torbjörn Larsson and Department of Mathematics, Linköping University, Sweden
More informationDecomposition methods in optimization
Decomposition methods in optimization I Approach I: I Partition problem constraints into two groups: explicit and implicit I Approach II: I Partition decision variables into two groups: primary and secondary
More informationDecomposition Methods for Integer Programming
Decomposition Methods for Integer Programming J.M. Valério de Carvalho vc@dps.uminho.pt Departamento de Produção e Sistemas Escola de Engenharia, Universidade do Minho Portugal PhD Course Programa Doutoral
More informationLarge-scale optimization and decomposition methods: outline. Column Generation and Cutting Plane methods: a unified view
Large-scale optimization and decomposition methods: outline I Solution approaches for large-scaled problems: I Delayed column generation I Cutting plane methods (delayed constraint generation) 7 I Problems
More information1 Solution of a Large-Scale Traveling-Salesman Problem... 7 George B. Dantzig, Delbert R. Fulkerson, and Selmer M. Johnson
Part I The Early Years 1 Solution of a Large-Scale Traveling-Salesman Problem............ 7 George B. Dantzig, Delbert R. Fulkerson, and Selmer M. Johnson 2 The Hungarian Method for the Assignment Problem..............
More informationBenders Decomposition
Benders Decomposition Yuping Huang, Dr. Qipeng Phil Zheng Department of Industrial and Management Systems Engineering West Virginia University IENG 593G Nonlinear Programg, Spring 2012 Yuping Huang (IMSE@WVU)
More informationStrengthened Benders Cuts for Stochastic Integer Programs with Continuous Recourse
Strengthened Benders Cuts for Stochastic Integer Programs with Continuous Recourse Merve Bodur 1, Sanjeeb Dash 2, Otay Günlü 2, and James Luedte 3 1 Department of Mechanical and Industrial Engineering,
More informationNotes on Dantzig-Wolfe decomposition and column generation
Notes on Dantzig-Wolfe decomposition and column generation Mette Gamst November 11, 2010 1 Introduction This note introduces an exact solution method for mathematical programming problems. The method is
More information3.10 Lagrangian relaxation
3.10 Lagrangian relaxation Consider a generic ILP problem min {c t x : Ax b, Dx d, x Z n } with integer coefficients. Suppose Dx d are the complicating constraints. Often the linear relaxation and the
More informationGeneration and Representation of Piecewise Polyhedral Value Functions
Generation and Representation of Piecewise Polyhedral Value Functions Ted Ralphs 1 Joint work with Menal Güzelsoy 2 and Anahita Hassanzadeh 1 1 COR@L Lab, Department of Industrial and Systems Engineering,
More information3.10 Column generation method
3.10 Column generation method Many relevant decision-making problems can be formulated as ILP problems with a very large (exponential) number of variables. Examples: cutting stock, crew scheduling, vehicle
More informationA Benders Algorithm for Two-Stage Stochastic Optimization Problems With Mixed Integer Recourse
A Benders Algorithm for Two-Stage Stochastic Optimization Problems With Mixed Integer Recourse Ted Ralphs 1 Joint work with Menal Güzelsoy 2 and Anahita Hassanzadeh 1 1 COR@L Lab, Department of Industrial
More informationClassification of Dantzig-Wolfe Reformulations for MIP s
Classification of Dantzig-Wolfe Reformulations for MIP s Raf Jans Rotterdam School of Management HEC Montreal Workshop on Column Generation Aussois, June 2008 Outline and Motivation Dantzig-Wolfe reformulation
More informationLagrange Relaxation: Introduction and Applications
1 / 23 Lagrange Relaxation: Introduction and Applications Operations Research Anthony Papavasiliou 2 / 23 Contents 1 Context 2 Applications Application in Stochastic Programming Unit Commitment 3 / 23
More informationConvex Optimization. Dani Yogatama. School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA. February 12, 2014
Convex Optimization Dani Yogatama School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA February 12, 2014 Dani Yogatama (Carnegie Mellon University) Convex Optimization February 12,
More informationInteger Programming ISE 418. Lecture 16. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 16 Dr. Ted Ralphs ISE 418 Lecture 16 1 Reading for This Lecture Wolsey, Chapters 10 and 11 Nemhauser and Wolsey Sections II.3.1, II.3.6, II.3.7, II.5.4 CCZ Chapter 8
More informationA Hub Location Problem with Fully Interconnected Backbone and Access Networks
A Hub Location Problem with Fully Interconnected Backbone and Access Networks Tommy Thomadsen Informatics and Mathematical Modelling Technical University of Denmark 2800 Kgs. Lyngby Denmark tt@imm.dtu.dk
More informationIntroduction to Bin Packing Problems
Introduction to Bin Packing Problems Fabio Furini March 13, 2015 Outline Origins and applications Applications: Definition: Bin Packing Problem (BPP) Solution techniques for the BPP Heuristic Algorithms
More informationVehicle Routing and MIP
CORE, Université Catholique de Louvain 5th Porto Meeting on Mathematics for Industry, 11th April 2014 Contents: The Capacitated Vehicle Routing Problem Subproblems: Trees and the TSP CVRP Cutting Planes
More informationNetwork Flows. 6. Lagrangian Relaxation. Programming. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 6. Lagrangian Relaxation 6.3 Lagrangian Relaxation and Integer Programming Fall 2010 Instructor: Dr. Masoud Yaghini Integer Programming Outline Branch-and-Bound Technique
More information1. Gradient method. gradient method, first-order methods. quadratic bounds on convex functions. analysis of gradient method
L. Vandenberghe EE236C (Spring 2016) 1. Gradient method gradient method, first-order methods quadratic bounds on convex functions analysis of gradient method 1-1 Approximate course outline First-order
More informationApplications. Stephen J. Stoyan, Maged M. Dessouky*, and Xiaoqing Wang
Introduction to Large-Scale Linear Programming and Applications Stephen J. Stoyan, Maged M. Dessouky*, and Xiaoqing Wang Daniel J. Epstein Department of Industrial and Systems Engineering, University of
More information3.10 Column generation method
3.10 Column generation method Many relevant decision-making (discrete optimization) problems can be formulated as ILP problems with a very large (exponential) number of variables. Examples: cutting stock,
More informationLinear Programming Inverse Projection Theory Chapter 3
1 Linear Programming Inverse Projection Theory Chapter 3 University of Chicago Booth School of Business Kipp Martin September 26, 2017 2 Where We Are Headed We want to solve problems with special structure!
More informationSolving Dual Problems
Lecture 20 Solving Dual Problems We consider a constrained problem where, in addition to the constraint set X, there are also inequality and linear equality constraints. Specifically the minimization problem
More informationLecture 7: Lagrangian Relaxation and Duality Theory
Lecture 7: Lagrangian Relaxation and Duality Theory (3 units) Outline Lagrangian dual for linear IP Lagrangian dual for general IP Dual Search Lagrangian decomposition 1 / 23 Joseph Louis Lagrange Joseph
More informationLecture 10: Linear programming duality and sensitivity 0-0
Lecture 10: Linear programming duality and sensitivity 0-0 The canonical primal dual pair 1 A R m n, b R m, and c R n maximize z = c T x (1) subject to Ax b, x 0 n and minimize w = b T y (2) subject to
More informationLECTURE 25: REVIEW/EPILOGUE LECTURE OUTLINE
LECTURE 25: REVIEW/EPILOGUE LECTURE OUTLINE CONVEX ANALYSIS AND DUALITY Basic concepts of convex analysis Basic concepts of convex optimization Geometric duality framework - MC/MC Constrained optimization
More informationAn Integrated Column Generation and Lagrangian Relaxation for Flowshop Scheduling Problems
Proceedings of the 2009 IEEE International Conference on Systems, Man, and Cybernetics San Antonio, TX, USA - October 2009 An Integrated Column Generation and Lagrangian Relaxation for Flowshop Scheduling
More informationBilevel Integer Programming
Bilevel Integer Programming Ted Ralphs 1 Joint work with: Scott DeNegre 1, Menal Guzelsoy 1 COR@L Lab, Department of Industrial and Systems Engineering, Lehigh University Norfolk Southern Ralphs, et al.
More informationA branch-and-price algorithm for parallel machine scheduling using ZDDs and generic branching
A branch-and-price algorithm for parallel machine scheduling using ZDDs and generic branching Kowalczyk D, Leus R. KBI_1631 A branch-and-price algorithm for parallel machine scheduling using ZDDs and generic
More informationAn Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse
An Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse Yongjia Song, James Luedtke Virginia Commonwealth University, Richmond, VA, ysong3@vcu.edu University
More informationColumn Generation I. Teo Chung-Piaw (NUS)
Column Generation I Teo Chung-Piaw (NUS) 21 st February 2002 1 Outline Cutting Stock Problem Slide 1 Classical Integer Programming Formulation Set Covering Formulation Column Generation Approach Connection
More information8. Geometric problems
8. Geometric problems Convex Optimization Boyd & Vandenberghe extremal volume ellipsoids centering classification placement and facility location 8 Minimum volume ellipsoid around a set Löwner-John ellipsoid
More informationLecture 3: Lagrangian duality and algorithms for the Lagrangian dual problem
Lecture 3: Lagrangian duality and algorithms for the Lagrangian dual problem Michael Patriksson 0-0 The Relaxation Theorem 1 Problem: find f := infimum f(x), x subject to x S, (1a) (1b) where f : R n R
More information3.7 Cutting plane methods
3.7 Cutting plane methods Generic ILP problem min{ c t x : x X = {x Z n + : Ax b} } with m n matrix A and n 1 vector b of rationals. According to Meyer s theorem: There exists an ideal formulation: conv(x
More informationInteger Programming ISE 418. Lecture 8. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 8 Dr. Ted Ralphs ISE 418 Lecture 8 1 Reading for This Lecture Wolsey Chapter 2 Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Duality for Mixed-Integer
More informationBilevel Integer Linear Programming
Bilevel Integer Linear Programming TED RALPHS SCOTT DENEGRE ISE Department COR@L Lab Lehigh University ted@lehigh.edu MOPTA 2009, Lehigh University, 19 August 2009 Thanks: Work supported in part by the
More informationA PARALLEL INTERIOR POINT DECOMPOSITION ALGORITHM FOR BLOCK-ANGULAR SEMIDEFINITE PROGRAMS IN POLYNOMIAL OPTIMIZATION
A PARALLEL INTERIOR POINT DECOMPOSITION ALGORITHM FOR BLOCK-ANGULAR SEMIDEFINITE PROGRAMS IN POLYNOMIAL OPTIMIZATION Kartik K. Sivaramakrishnan Department of Mathematics North Carolina State University
More informationKNAPSACK PROBLEMS WITH SETUPS
7 e Conférence Francophone de MOdélisation et SIMulation - MOSIM 08 - du 31 mars au 2 avril 2008 - Paris - France Modélisation, Optimisation et Simulation des Systèmes : Communication, Coopération et Coordination
More informationUses of duality. Geoff Gordon & Ryan Tibshirani Optimization /
Uses of duality Geoff Gordon & Ryan Tibshirani Optimization 10-725 / 36-725 1 Remember conjugate functions Given f : R n R, the function is called its conjugate f (y) = max x R n yt x f(x) Conjugates appear
More informationDM545 Linear and Integer Programming. Lecture 13 Branch and Bound. Marco Chiarandini
DM545 Linear and Integer Programming Lecture 13 Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 2 Exam Tilladt Håndscanner/digital pen og ordbogsprogrammet
More informationSeparation, Inverse Optimization, and Decomposition. Some Observations. Ted Ralphs 1 Joint work with: Aykut Bulut 1
: Some Observations Ted Ralphs 1 Joint work with: Aykut Bulut 1 1 COR@L Lab, Department of Industrial and Systems Engineering, Lehigh University COLGEN 2016, Buzios, Brazil, 25 May 2016 What Is This Talk
More informationA Capacity Scaling Procedure for the Multi-Commodity Capacitated Network Design Problem. Ryutsu Keizai University Naoto KATAYAMA
A Capacity Scaling Procedure for the Multi-Commodity Capacitated Network Design Problem Ryutsu Keizai University Naoto KATAYAMA Problems 2006 1 Multi-Commodity Network Design Problem The basic model for
More informationSeparation, Inverse Optimization, and Decomposition. Some Observations. Ted Ralphs 1 Joint work with: Aykut Bulut 1
: Some Observations Ted Ralphs 1 Joint work with: Aykut Bulut 1 1 COR@L Lab, Department of Industrial and Systems Engineering, Lehigh University MOA 2016, Beijing, China, 27 June 2016 What Is This Talk
More informationIntroduction to Integer Linear Programming
Lecture 7/12/2006 p. 1/30 Introduction to Integer Linear Programming Leo Liberti, Ruslan Sadykov LIX, École Polytechnique liberti@lix.polytechnique.fr sadykov@lix.polytechnique.fr Lecture 7/12/2006 p.
More informationMulticommodity Flows and Column Generation
Lecture Notes Multicommodity Flows and Column Generation Marc Pfetsch Zuse Institute Berlin pfetsch@zib.de last change: 2/8/2006 Technische Universität Berlin Fakultät II, Institut für Mathematik WS 2006/07
More informationLagrangean relaxation
Lagrangean relaxation Giovanni Righini Corso di Complementi di Ricerca Operativa Joseph Louis de la Grange (Torino 1736 - Paris 1813) Relaxations Given a problem P, such as: minimize z P (x) s.t. x X P
More informationLinear Programming. Leo Liberti. LIX, École Polytechnique. Operations research courses / LP theory p.
Operations research courses / LP theory p. 1/47 Linear Programming Leo Liberti LIX, École Polytechnique liberti@lix.polytechnique.fr Operations research courses / LP theory p. 2/47 Contents LP formulations
More informationBig Data Analytics: Optimization and Randomization
Big Data Analytics: Optimization and Randomization Tianbao Yang Tutorial@ACML 2015 Hong Kong Department of Computer Science, The University of Iowa, IA, USA Nov. 20, 2015 Yang Tutorial for ACML 15 Nov.
More informationLinear and Integer Programming - ideas
Linear and Integer Programming - ideas Paweł Zieliński Institute of Mathematics and Computer Science, Wrocław University of Technology, Poland http://www.im.pwr.wroc.pl/ pziel/ Toulouse, France 2012 Literature
More informationSection Notes 8. Integer Programming II. Applied Math 121. Week of April 5, expand your knowledge of big M s and logical constraints.
Section Notes 8 Integer Programming II Applied Math 121 Week of April 5, 2010 Goals for the week understand IP relaxations be able to determine the relative strength of formulations understand the branch
More informationInterior Point Cutting Plane Methods in Integer Programming
Interior Point Cutting Plane Methods in Integer Programming by Joe Naoum-Sawaya A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy
More informationStochastic Integer Programming An Algorithmic Perspective
Stochastic Integer Programming An Algorithmic Perspective sahmed@isye.gatech.edu www.isye.gatech.edu/~sahmed School of Industrial & Systems Engineering 2 Outline Two-stage SIP Formulation Challenges Simple
More informationProjection in Logic, CP, and Optimization
Projection in Logic, CP, and Optimization John Hooker Carnegie Mellon University Workshop on Logic and Search Melbourne, 2017 Projection as a Unifying Concept Projection is a fundamental concept in logic,
More informationCutting Plane Separators in SCIP
Cutting Plane Separators in SCIP Kati Wolter Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies 1 / 36 General Cutting Plane Method MIP min{c T x : x X MIP }, X MIP := {x
More information0-1 Reformulations of the Network Loading Problem
0-1 Reformulations of the Network Loading Problem Antonio Frangioni 1 frangio@di.unipi.it Bernard Gendron 2 bernard@crt.umontreal.ca 1 Dipartimento di Informatica Università di Pisa Via Buonarroti, 2 56127
More information8. Geometric problems
8. Geometric problems Convex Optimization Boyd & Vandenberghe extremal volume ellipsoids centering classification placement and facility location 8 1 Minimum volume ellipsoid around a set Löwner-John ellipsoid
More informationStochastic Decomposition
IE 495 Lecture 18 Stochastic Decomposition Prof. Jeff Linderoth March 26, 2003 March 19, 2003 Stochastic Programming Lecture 17 Slide 1 Outline Review Monte Carlo Methods Interior Sampling Methods Stochastic
More informationAccelerating Benders Decomposition by Local Branching
Accelerating Benders Decomposition by Local Branching Walter Rei, Michel Gendreau Département d informatique et de recherche opérationnelle and Centre de recherche sur les transports, Université demontréal,
More informationELE539A: Optimization of Communication Systems Lecture 16: Pareto Optimization and Nonconvex Optimization
ELE539A: Optimization of Communication Systems Lecture 16: Pareto Optimization and Nonconvex Optimization Professor M. Chiang Electrical Engineering Department, Princeton University March 16, 2007 Lecture
More informationConvex Optimization Lecture 16
Convex Optimization Lecture 16 Today: Projected Gradient Descent Conditional Gradient Descent Stochastic Gradient Descent Random Coordinate Descent Recall: Gradient Descent (Steepest Descent w.r.t Euclidean
More informationSolution Methods for Stochastic Programs
Solution Methods for Stochastic Programs Huseyin Topaloglu School of Operations Research and Information Engineering Cornell University ht88@cornell.edu August 14, 2010 1 Outline Cutting plane methods
More informationAdvanced Linear Programming: The Exercises
Advanced Linear Programming: The Exercises The answers are sometimes not written out completely. 1.5 a) min c T x + d T y Ax + By b y = x (1) First reformulation, using z smallest number satisfying x z
More information15-780: LinearProgramming
15-780: LinearProgramming J. Zico Kolter February 1-3, 2016 1 Outline Introduction Some linear algebra review Linear programming Simplex algorithm Duality and dual simplex 2 Outline Introduction Some linear
More informationLocal Cuts and Two-Period Convex Hull Closures for Big-Bucket Lot-Sizing Problems
Local Cuts and Two-Period Convex Hull Closures for Big-Bucket Lot-Sizing Problems Kerem Akartunalı Dept. of Management Science, University of Strathclyde, Glasgow, G1 1QE, UK, kerem.akartunali@strath.ac.uk
More informationPartial Path Column Generation for the Elementary Shortest Path Problem with Resource Constraints
Partial Path Column Generation for the Elementary Shortest Path Problem with Resource Constraints Mads Kehlet Jepsen & Bjørn Petersen Department of Computer Science, University of Copenhagen Universitetsparken
More informationLogic-based Benders Decomposition
Logic-based Benders Decomposition A short Introduction Martin Riedler AC Retreat Contents 1 Introduction 2 Motivation 3 Further Notes MR Logic-based Benders Decomposition June 29 July 1 2 / 15 Basic idea
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationEE364a Review Session 5
EE364a Review Session 5 EE364a Review announcements: homeworks 1 and 2 graded homework 4 solutions (check solution to additional problem 1) scpd phone-in office hours: tuesdays 6-7pm (650-723-1156) 1 Complementary
More informationBilevel Integer Optimization: Theory and Algorithms
: Theory and Algorithms Ted Ralphs 1 Joint work with Sahar Tahernajad 1, Scott DeNegre 3, Menal Güzelsoy 2, Anahita Hassanzadeh 4 1 COR@L Lab, Department of Industrial and Systems Engineering, Lehigh University
More informationBenders Decomposition Methods for Structured Optimization, including Stochastic Optimization
Benders Decomposition Methods for Structured Optimization, including Stochastic Optimization Robert M. Freund April 29, 2004 c 2004 Massachusetts Institute of echnology. 1 1 Block Ladder Structure We consider
More informationCoordinate Descent and Ascent Methods
Coordinate Descent and Ascent Methods Julie Nutini Machine Learning Reading Group November 3 rd, 2015 1 / 22 Projected-Gradient Methods Motivation Rewrite non-smooth problem as smooth constrained problem:
More informationDual methods and ADMM. Barnabas Poczos & Ryan Tibshirani Convex Optimization /36-725
Dual methods and ADMM Barnabas Poczos & Ryan Tibshirani Convex Optimization 10-725/36-725 1 Given f : R n R, the function is called its conjugate Recall conjugate functions f (y) = max x R n yt x f(x)
More informationThe two-machine flowshop total completion time problem: A branch-and-bound based on network-flow formulation
The two-machine flowshop total completion time problem: A branch-and-bound based on network-flow formulation Boris Detienne 1, Ruslan Sadykov 1, Shunji Tanaka 2 1 : Team Inria RealOpt, University of Bordeaux,
More informationSection Notes 9. Midterm 2 Review. Applied Math / Engineering Sciences 121. Week of December 3, 2018
Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections
More information