A Scheme for Integrated Optimization
|
|
- Roger Kelley
- 5 years ago
- Views:
Transcription
1 A Scheme for Integrated Optimization John Hooker ZIB November 2009 Slide 1
2 Outline Overview of integrated methods A taste of the underlying theory Eamples, with results from SIMPL Proposal for SCIP/SIMPL cooperation Slide 2
3 Overview of integrated methods Integration principles. Search-infer-rela Classical solution methods. Evolution of SIMPL Slide 3
4 Integration Principles Integrate MILP, constraint programming, global optimization in a unified approach. Slide 4
5 Integration Principles Integrate MILP, constraint programming, global optimization in a unified approach. Low-level integration with high-level modeling. Slide 5
6 Integration Principles Integrate MILP, constraint programming, global optimization in a unified approach. Low-level integration with high-level modeling. Succinct modeling with meta-constraints. Model communicates problem structure to the solver. Slide 6
7 Integration Principles Integrate MILP, constraint programming, global optimization in a unified approach. Low-level integration with high-level modeling. Succinct modeling with meta-constraints. Model communicates problem structure to the solver. General search-infer-rela solution algorithm. Enumerate problem restrictions. Branching or logic-based Benders. Underlying search/inference and search/relaation dualities. Slide 7
8 Integration Principles Integrate MILP, constraint programming, global optimization in a unified approach. Low-level integration with high-level modeling. Succinct modeling with meta-constraints. Model communicates problem structure to the solver. General search-infer-rela solution algorithm. Enumerate problem restrictions. Branching or logic-based Benders. Underlying search/inference and search/relaation dualities Constraint-based control. Filtering, relaation, branching. Slide 8
9 Search-Infer-Rela Search: Enumerate problem restrictions. Branching tree nodes, Benders subproblems, local search neighborhoods, etc. Infer: Deduce constraints from current restriction Nogoods, cutting planes, filtering, etc. Rela: Solve relaation of current restriction LP, Lagrangean, domain store, Benders master, etc. Slide 9
10 Classical solution methods CP solver Search: Branching Inference: Flltering Relaation: Domain store MILP solver Search: Branching Inference: Cutting planes, presolve, reduced cost variable fiing Relaation: LP Benders Search: Enumerate subproblems. Inference: Benders cuts Relaation: Master problem Slide 10
11 Classical solution methods Global optimization Search: Enumerate boes Inference: Domain reduction, dual-based variable bounding Relaation: Conveification SAT Search: Branching Inference: Conflict clauses Relaation: Same as restriction Local search Search: Enumerate neighborhoods. Inference: Tabu list, etc. Relaation: Same as restriction Slide 11
12 Interaction Inference Strengthens Fies variables Reduces domains Adds IP cuts Adds conflict clause Adds Benders cuts Shrinks bo Creates neighborhood Guides Separating cut Relaation Slide 12 Activates IP cut Fied variable Filtering/propagation Reduced domain Benders cut Subproblem dual Restriction Defines Identifies net branch Fractional variable Nonsingleton domain Violated constraint Defines subproblem Solution of master Defines neighborhood Center on previous solution
13 Evolution of SIMPL Basic Framework JNH. Logic-based methods for optimization, CP 1994 JNH, M. A. Osorio. Mied logical/linear programming, Discrete Applied Mathematics JHH, G. Ottosson, E. S. Thorsteinsson, H.-J. Kim. On integrating CP and LP for combinatorial optimization. AAAI G. Ottosson, E. Thorsteinsson, JNH. Mied global constraints and inference in hybrid IP-CLP solvers. CP 1999 JNH, G. Ottosson, E. S. Thorsteinsson, H.-J. Kim. A scheme for unifying optimization and constraint satisfaction methods. Knowledge Engineering Review, Slide 13
14 Evolution of SIMPL JNH. A framework for integrating solution methods. ICS I. Aron, JNH, T. H. Yunes. SIMPL: A system for integrating optimization techniques. CPAIOR T. Yunes, I. Aron, JNH. An integrated solver for optimization problems. Operations Research, to appear. Slide 14
15 Evolution of SIMPL Search JNH. Unifying local and ehaustive search. ENC JNH, A search-infer-and-rela framework for integrating solution methods. CPAIOR 2005 Modeling JNH, H.-J. Kim, G. Ottosson. A declarative modeling framework that integrates solution methods. Annals of OR JNH. Hybrid modeling. Ten Years of CPAIOR, to appear. Slide 15
16 Evolution of SIMPL Theory and background JNH. Logic-Based Methods for Optimization JNH. Logic, optimization and constraint programming. INFORMS Journal on Computing, A. Bockmayr, JNH. Constraint programming. Handbook of Discrete Optimization, JNH. Operations research methods in constraint programming, Handbook of Constraint Programming, J. N. Hooker. Duality in optimization and constraint satisfaction. CPAIOR 2006 JNH. Integrated Methods for Optimization Slide 16
17 Evolution of SIMPL Nogood-based search (Logic-based Benders) JNH, H. Yan. Logic circuit verification by Benders decomposition. CP JNH. Logic-Based Methods for Optimization E. Thorsteinsson. Branch and check: A hybrid framework for integrating MIP and CLP. CP JNH, G. Ottosson. Logic-based Benders decomposition. Mathematical Programming, JNH. A hybrid method for planning and scheduling. Constraints JNH. Planning and scheduling to minimize tardiness. CP JNH. Planning and scheduling by logic-based Benders decomposition. Operations Research, Slide 17
18 Evolution of SIMPL Relaation methods I. E. Grossmann, JNH, R. Raman, H. Yan. Logic cuts for processing networks with fied charges. Computers and OR, E. Thorsteinsson and G. Ottosson. Linear relaations and reduced-cost based propagation of continuous variable subscripts. Annals of OR, S. Bollapragada, O. Ghattas, JNH. Optimal design of truss structures by mied logical and linear programming. Operations Research JNH. Conve programming methods for global optimization, COCOS 2003 L. Genc-Kaya and JNH. The circuit polytope Slide 18
19 A taste of the underlying theory Inference duality. LP, Lagrangean, surrogate duals Nogood-based search Logic-based and classical Benders Eample: SAT DPL DPL + conflict clauses Partial order dynamic backtracking Slide 19
20 Inference duality All optimization duals are inference duals Also relaation duals Solution of inference dual is proof of optimality Primal in co-np when dual is in NP Provides nogoods Including Benders cuts, conflict clauses Nogood = conditions under which proof is still valid Directs the search All search is nogood-based search Net restriction defined by nogood Postoptimality analysis Result of altering premises of proof Slide 20
21 Inference dual Primal min f ( ) S Feasible set ma Dual v P S v f ( ) P P Family of proofs (inference method) Follows using proof P Dual is defined relative to an inference method. Strong duality applies if inference method is complete. Slide 21
22 Eample: LP dual Primal min c A b 0 ma Dual v A b P c v 0 P P = ma λb λa c λ 0 Nonnegative linear combination + domination 0 λa λb dominates c v A b c v iff for some λ 0 λa c and λb v Slide 22 Inference method is complete (assuming feasibility) due to Farkas Lemma. So we have strong duality (assuming feasibility).
23 Eample: Lagrangean dual Primal min f ( ) g( ) 0 S ma v P g( ) b f ( ) v P P Dual λ 0 Nonnegative linear combination + domination { λ } = ma min ( ) ( ) f g S S λg( ) 0 dominates f ( ) v 0 g( ) 0 f ( ) v iff for some λ 0 λg() f() v for all S That is, v f() λg() for all S Or v min { f ( ) λg( ) } S Slide 23 Inference method is incomplete
24 Eample: Surrogate dual Primal min f ( ) g( ) 0 S ma v P g( ) b f ( ) v P P Dual λ 0 S { f λg } = ma min ( ) ( ) 0 Nonnegative linear combination + implication S λg( ) 0 implies f ( ) v g( ) 0 f ( ) v iff for some λ 0 Any S with λg() 0 satisfies f() v So, min {f() λg() 0, S} v Slide 24 Inference method is incomplete
25 Nogood-based search All search is nogood-based search Each solution eamined generates a nogood. Net solution must satisfy current nogood set. Nogoods are derived by solving inference dual of the subproblem. Subproblem is normally defined by fiing variables to current values in the search. Slide 25
26 Eample: Logic-based Benders Partition variables,y and search over values of Subproblem results from fiing min f (, y ) min f (, y ) (, y ) S (, y ) S Let proof P be solution of subproblem dual for = Let B(P,) be lower bound obtained by P for given. Add Benders cut v B(P,) to master problem: Solve master problem for net min v Benders cuts Slide 26
27 Eample: Classical Benders Partition variables,y and search over values of min f ( ) + cy g( ) + Ay b D, 0 y Subproblem results from fiing k min f ( ) + cy k Ay b g( ) y 0 ( λ) Let proof λ be solution of subproblem dual for Let B(λ,) = f() + λ(b g()) be bound obtained by λ for given. Add Benders cut v B(λ,) to master problem: Solve master problem for net = min v Benders cuts Slide 27
28 Eample: SAT Solve SAT by chronological backtracking + unit clause rule = DPL. Chronological = fied branching order. To get nogood, solve subproblem at current node. Solve with unit clause rule Nogood identifies branches that create infeasibility. Simplest scheme: nogood rules out path to current leaf node. Process nogood set with parallel resolution Nogood set is a relaation of the problem. Solve relaation without branching Select solution with preference for 0 Slide 28
29 DPL with chronological backtracking 1 = 0 2 = 0 3 = 0 5 = 0 4 = 0 Branch to here. Solve subproblem with unit clause, which proves infeasibility. ( 1, 5 ) = (0,, 0) creates the infeasibility. SARA 07 Slide 29
30 DPL with chronological backtracking 1 = 0 2 = 0 3 = 0 5 = 0 4 = 0 Branch to here. Solve subproblem with unit clause, which proves infeasibility. ( 1, 5 ) = (0,, 0) creates the infeasibility Generate nogood. SARA 07 Slide 30
31 DPL with chronological backtracking 1 = 0 Consists of processed nogoods 3 = 0 2 = 0 k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) k k = 0 5 = 0 Conflict clause appears as nogood induced by solution of R k SARA 07 Slide 31
32 DPL with chronological backtracking 1 = 0 Consists of processed nogoods 3 = 0 2 = 0 k Relaation R Solution of R Nogoods k 0 (0,0,0,0,0, ) (0,0,0,0,1, ) k 5 = 0 4 = 0 5 = 1 Go to solution that solves relaation, with priority to 0 SARA 07 Slide 32
33 DPL with chronological backtracking 1 = 0 Consists of processed nogoods 3 = 0 2 = 0 k Relaation R Solution of R Nogoods k 0 (0,0,0,0,0, ) (0,0,0,0,1, ) k 5 = 0 4 = 0 5 = Process nogood set with parallel resolution parallel-absorbs SARA 07 Slide 33
34 DPL with chronological backtracking 1 = 0 Consists of processed nogoods 5 = 0 4 = 0 3 = 0 2 = 0 5 = 1 k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (0,0,0,0,1, ) k k Process nogood set with parallel resolution parallel-absorbs SARA 07 Slide 34
35 DPL with chronological backtracking 1 = 0 Consists of processed nogoods 5 = 0 4 = 0 3 = 0 2 = 0 k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (0,0,0,0,1, ) 2 (0,0,0,1,0, ) k Solve relaation again, continue. So backtracking is nogood-based search with parallel resolution k SARA 07 Slide 35
36 Eample: SAT + conflict clauses Nogoods = conflict clauses. Nogoods rule out only branches that play a role in unit clause refutation. Slide 36
37 DPL with conflict clauses 1 = 0 2 = 0 3 = 0 4 = 0 = 0 5 Branch to here. Unit clause rule proves infeasibility. ( 1, 5 ) = (0,0) is only premise of unit clause proof. SARA 07 Slide 37
38 DPL with conflict clauses 1 = 0 3 = 0 2 = 0 k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) k k = 0 5 = 0 Conflict clause appears as nogood induced by solution of R k. 1 5 SARA 07 Slide 38
39 DPL with conflict clauses 1 = 0 Consists of processed nogoods 3 = 0 2 = 0 k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) k 1 (0,0,0,0,1, ) k 4 = 0 5 = 0 5 = SARA 07 Slide 39
40 DPL with conflict clauses 1 = 0 Consists of processed nogoods 4 = 0 3 = 0 2 = 0 5 = 0 5 = k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (0,0,0,0,1, ) k parallel-resolve to yield 1 2 parallel-absorbs k SARA 07 Slide 40
41 DPL with conflict clauses 1 = 0 2 = 0 2 = 1 4 = 0 3 = 0 5 = 0 = SARA 07 Slide k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (0,0,0,0,1, ) 2 (0,1,,,, ) k parallel-resolve to yield k
42 DPL with conflict clauses 1 = 0 1 = 1 4 = 0 3 = 0 SARA 07 Slide 42 2 = 0 5 = 0 5 = k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (0,0,0,0,1, ) 2 (0,1,,,, ) 3 (1,,,,, ) 4 k Search terminates k
43 Eample: SAT + partial order dynamic backtracking Solve relaation by selecting a solution that conforms to nogoods. Conform = takes opposite sign than in nogoods. More freedom than in branching. Slide 43
44 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) k k 5 1 Arbitrarily choose one variable to be last SARA 07 Slide 44
45 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) k k 5 1 Other variables are penultimate Arbitrarily choose one variable to be last SARA 07 Slide 45
46 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (1,,,,0, ) k k 5 1 Since 5 is penultimate in at least one nogood, it must conform to nogoods. It must take value opposite its sign in the nogoods. 5 will have the same sign in all nogoods where it is penultimate. SARA 07 Slide 46 This allows more freedom than chronological backtracking.
47 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (1,,,,0, ) 2 k k Choice of last variable is arbitrary but must be consistent with partial order implied by previous choices. Since 5 is penultimate in at least one nogood, it must conform to nogoods. It must take value opposite its sign in the nogoods. 5 will have the same sign in all nogoods where it is penultimate. SARA 07 Slide 47 This allows more freedom than chronological backtracking.
48 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (1,,,,0, ) k k Choice of last variable is arbitrary but must be consistent with partial order implied by previous choices Parallel-resolve to yield 5 SARA 07 Slide 48 Since 5 is penultimate in at least one nogood, it must conform to nogoods. It must take value opposite its sign in the nogoods. 5 will have the same sign in all nogoods where it is penultimate. This allows more freedom than chronological backtracking.
49 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (1,,,,0, ) 2 (,0,,,1, ) k k 5 does not parallel-resolve with because 5 is not last in both clauses 5 2 SARA 07 Slide 49
50 Partial Order Dynamic Backtracking k Relaation R Solution of R Nogoods 0 (0,0,0,0,0, ) 1 (1,,,,0, ) 2 (,0,,,1, ) (,1,,,1, ) 4 k k 2 Search terminates Must conform SARA 07 Slide 50
51 Eamples, with results from SIMPL Production planning. Semicontinuous piecewise linear functions Product configuration. Variable indices Machine scheduling. Logic-based Benders Truss structure design. Global optimization. Slide 51
52 Production Planning Maimize profit, which is a piecewise linear function of output. ma f ( ) i i i i C i Each f i is a piecewise linear semicontinuous function Slide 52
53 Production Planning Semicontinuous piecewise linear function f() Slide 53
54 Production Planning MILP model SOS2 branching not useful min i i ik = λ L + µ U, all i λ + µ = 1, all i 0 λ y, all i, k 0 µ y, all i, k { } ik ik ik ik i ik ik ik ik k k k y ik y ik ik ik ik λ c + µ d C ik ik ik = 1, all i 0,1, all i, k Slide 54
55 Production Planning MILP model SOS2 branching not useful = 1 if i is in interval k min i i ik = λ L + µ U, all i λ + µ = 1, all i 0 λ y, all i, k 0 µ y, all i, k { } ik ik ik ik i ik ik ik ik k k k y ik y ik ik ik ik λ c + µ d C ik ik ik = 1, all i 0,1, all i, k Slide 55
56 Production Planning Integrated model ma i i i u C i ( u L U c d ) piecewise,,,,,, all i i i i i i i Metaconstraint (global constraint in CP) Original hand-coded method: Ottosson, Thorsteinsson and JNH Slide 56
57 Production Planning Semicontinuous piecewise linear function f() Tight linear relaation Slide 57
58 Production Planning Semicontinuous piecewise linear function f() Tighter relaation after branching Slide 58 Value of in solution of current linear relaation
59 Production Planning SIMPL model Slide 59
60 Production Planning SIMPL model Recognized as a linear system. Slide 60
61 Production Planning SIMPL model Is its own LP relaation. CP relaation propagates bounds. Slide 61
62 Production Planning SIMPL model Piecewise linear metaconstraint. Slide 62
63 Production Planning SIMPL model LP relaation is conve hull. CP relaation propagates bounds. Slide 63
64 Production Planning SIMPL model Branch-and-bound search. Dive to leaf node from node with best lower bound. Slide 64
65 Production Planning SIMPL model Branch on piecewise constraint with greatest violation. Slide 65
66 Production Planning Computational Results (seconds) Hand-coded integrated method was comparable to CPLEX 9 No. Products MILP CPLEX 9 MILP CPLEX 11 SIMPL Slide 66
67 Production Planning CPLEX has become orders of magnitude faster, but still slower than SIMPL No. Products MILP CPLEX 9 MILP CPLEX 11 SIMPL Slide 67
68 SIMPL s advantage grows with the problem size No. Products MILP CPLEX 9 Seconds MILP CPLEX 11 SIMPL Slide 68
69 SIMPL s advantage grows with the problem size No. Products MILP CPLEX 9 Seconds Nodes MILP CPLEX 11 SIMPL Slide 69 No. Products MILP CPLEX 9 MILP CPLEX 11 SIMPL , , , , , , ,732 1,297,
70 Product configuration Choose what type of each component, and how many Memory Memory Memory Memory Memory Memory Personal computer Power supply Disk drive Disk drive Disk drive Disk drive Disk drive Power supply Power supply Power supply Slide 70
71 Product configuration Unit cost of producing attribute Amount of attribute produced (< 0 if consumed): memory, heat, power, weight, etc. MILP model min v = A q, all L v U, all q M, all i, k k ik ik tk ik i ik ik c v =1, all i Quantity of component i installed Slide 71
72 Product configuration Unit cost of producing attribute Amount of attribute produced (< 0 if consumed): memory, heat, power, weight, etc. MILP model min v = A q, all L v U, all q M, all i, k k ik ik tk ik i ik ik c v =1, all i Amount of attribute produced by type k of component i Quantity of component i installed Slide 72
73 Product configuration Unit cost of producing attribute Amount of attribute produced (< 0 if consumed): memory, heat, power, weight, etc. MILP model min v = A q, all L v U, all q M, all i, k k ik ik tk ik i ik ik c v =1, all i Amount of attribute produced by type k of component i Quantity of component i installed 1 if type k of component i is used Slide 73
74 Product configuration Unit cost of producing attribute Amount of attribute produced (< 0 if consumed): memory, heat, power, weight, etc. min c v v = A q, all iti i ik L v U, all Integrated model Quantity of component i installed Original hand-coded method: Thorsteinsson and Ottosson Slide 74
75 Product configuration Unit cost of producing attribute Amount of attribute produced by type t i of component i Amount of attribute produced (< 0 if consumed): memory, heat, power, weight, etc. min c v v = A q, all iti i ik L v U, all Integrated model Quantity of component i installed Original hand-coded method: Thorsteinsson and Ottosson Slide 75
76 Product configuration Unit cost of producing attribute Amount of attribute produced by type t i of component i Amount of attribute produced (< 0 if consumed): memory, heat, power, weight, etc. min c v v = A q, all iti i ik L v U, all t i is a variable inde Integrated model Original hand-coded method: Thorsteinsson and Ottosson Slide 76
77 Product configuration Linear inequality metaconstraint min c v v = q A, all i it ik L v U, all i Original hand-coded method: Thorsteinsson and Ottosson Slide 77
78 Product configuration Indeed linear metaconstraint min c v v = q A, all i it ik L v U, all i Original hand-coded method: Thorsteinsson and Ottosson Slide 78
79 Product configuration Propagation min c v v = q A, all i it ik L v U, all i This is propagated in the usual way Slide 79
80 Product configuration Propagation v z, all = i i ( ) element t,( q, A,, q A ), z, all i, i i i1 i in i min c v v = q A, all i it ik L v U, all i This is rewritten as This is propagated in the usual way Slide 80
81 Product configuration Propagation v z, all = i i ( ) element t,( q, A,, q A ), z, all i, i i i1 i in i This is propagated by (a) using specialized filters for element constraints of this form Slide 81
82 Product configuration Propagation v z, all = i i ( ) element t,( q, A,, q A ), z, all i, i i i1 i in i Slide 82 This is propagated by (a) using specialized filters for element constraints of this form, (b) adding knapsack cuts for the valid inequalities: i i ma k D t i min k D t i { } A q v, all ik { } i A q v, all ik i and (c) propagating the knapsack cuts. [ v, v ] is current domain of v
83 Product configuration Relaation min c v v = q A, all i it ik L v U, all i This is relaed as v v v Slide 83
84 Product configuration Relaation v z, all = i i ( ) element t,( q, A,, q A ), z, all i, i i i1 i in i min c v v = q A, all i it ik L v U, all i This is relaed by relaing this and adding the knapsack cuts. This is relaed as v v v Slide 84
85 Product configuration Relaation v z, all = i i ( ) element t,( q, A,, q A ), z, all i, i i i1 i in i This is relaed by writing each element constraint as a disunction of linear systems and writing a conve hull relaation of the disunction: z = A q, q = q i ik ik i ik k D k D ti ti Slide 85
86 Product configuration SIMPL model Recognized as indeed linear system Slide 86
87 Product configuration SIMPL model LP relaation is conve hull of disunction. CP relaation propagates bounds. Slide 87
88 Product configuration SIMPL model Generate knapsack cuts from associated valid inequalities. Slide 88
89 Product configuration SIMPL model Logical constraint on quantities Slide 89
90 Product configuration SIMPL model First branch on violated logical constraint on q i variables Slide 90
91 Product configuration SIMPL model Then branch on most violated t i in-domain constraint. Violated when domain of t i is not a singleton, or two or more associated q ik s are positive. Slide 91
92 Product configuration SIMPL model Then branch on least violated q i in-domain constraint. Create three branches: q i = nearest integer q i, q i < q i, q i > q i Slide 92
93 Product configuration SIMPL model Then branch on most violated logical constraint on t i variables (omitted) Slide 93
94 Product configuration SIMPL model Reduced-cost variable fiing for q i s Slide 94
95 Product configuration MILP (CPLEX 11) SIMPL Computational results SIMPL matches hand-coded integrated method, which was orders of magnitude faster than CPLEX. Again, CPLEX has become much faster, now somewhat faster than SIMPL. Slide 95 Nodes Sec. Nodes Sec
96 Machine scheduling Assign obs to machines, and schedule the machines assigned to each machine within time windows. The obective is to minimize processing cost. Slide 96
97 Machine scheduling Job Data Eample Assign 5 obs to 2 machines. Schedule obs assigned to each machine without overlap. Machine A Machine B Slide 97
98 Machine scheduling MILP continuoustime model from Jain & Grossmann min i i i r s d p, all i i i i i = 1, all y + y 1, all > y + y + + 2, all >, i i i i y + y + 1, all >, i i i s s + p M(1 y ), all i i i i c { } { } p ma d min r, all i i i { } { } 0,1, y 0,1, all Slide 98
99 Machine scheduling MILP continuoustime model from Jain & Grossmann = 1 of ob precedes ob start time min i i i r s d p, all i i i i i y + y 1, all > y + y + + 2, all >, i i i i y + y + 1, all >, i i i s s + p M(1 y ), all i i i i c = 1, all { } { } p ma d min r, all i i i { } { } = 1 if ob assigned to machine i 0,1, y 0,1, all Slide 99
100 Machine scheduling Integrated model min c r s d p, all Start time of ob ( = i = ) Time windows disunctive ( s i),( p i), all i Jobs cannot overlap Machine assigned to ob Original hand-coded method: Jain and Grossman Slide 100
101 Machine scheduling Integrated approach Assign the obs in the master problem, to be solved by MILP. Schedule the obs in the subproblem, to be solved by CP. Original hand-coded method: Jain and Grossman Slide 101
102 Machine scheduling Integrated approach Assign the obs in the master problem, to be solved by MILP. Schedule the obs in the subproblem, to be solved by CP. The subproblem decouples into a separate scheduling problem on each machine. In this problem, the subproblem is a feasibility problem. Original hand-coded method: Jain and Grossman Slide 102
103 Machine scheduling Integrated model min c r s d p, all ( = i = ) Indeed linear metaconstraint disunctive ( s i ),( p i ), all i Slide 103
104 Machine scheduling Integrated model min c r s d p, all ( = i = ) Indeed linear metaconstraint disunctive ( s i ),( p i ), all i Disunctive scheduling metaconstraint Slide 104
105 Machine scheduling Integrated model min M M s + p, all r s d p, all Start time of ob ( = i = ) Time windows disunctive ( s i ),( p i ), all i Jobs cannot overlap For a fied assignment the subproblem on each machine i is Slide 105 min M M s + p, all with = i r s d p, all with = i ( s = i pi = i ) disunctive ( ),( )
106 Machine scheduling Logic-based Benders approach Suppose we assign obs 1,2,3,5 to machine A in iteration k. We can prove that there is no feasible schedule. Edge finding derives infeasibility by reasoning only with obs 2,3,5. So these obs alone create infeasibility. So we have a Benders cut ( = = = A) Slide 106
107 Machine scheduling Logic-based Benders approach We want the master problem to be an MILP, which is good for assignment problems. So we write the Benders cut ( A) = = = Using 0-1 variables: A A A = 1 if ob 5 is assigned to machine A Slide 107
108 Machine scheduling The master problem is a relaation, formulated as an MILP: min 5 = 1 5 = 1 i A B A2 A3 A5 i p p c i A B { 0,1} i 9, etc. 9, etc Constraints derived from time windows Benders cut from machine A Slide 108
109 Machine scheduling The master problem is a relaation, formulated as an MILP: min 5 = 1 5 = 1 A2 A3 A5 i i p A p B c i A B { 0,1} i 9, etc. 9, etc Constraints derived from time windows Benders cut from machine A Benders cuts have been developed for min makespan and min tardiness (subproblem is an optimization problem) Slide 109 Also for cumulative scheduling.
110 Machine scheduling SIMPL model Machine assignment constraint Slide 110
111 Machine scheduling SIMPL model MILP relaation of the constraint (which is the constraint itself) goes into master problem Slide 111
112 Machine scheduling SIMPL model Definition of i variables for MILP master problem Slide 112
113 Machine scheduling SIMPL model CP-based propagation Slide 113
114 Machine scheduling SIMPL model Time window constraints recognized as indeed linear system Slide 114
115 Machine scheduling SIMPL model MILP formulation goes into master problem CP-based propagation Slide 115
116 Machine scheduling SIMPL model Disunctive scheduling constraint written as special case of cumulative scheduling constraint (resource consumption = 1, capacity = 1) Slide 116
117 Machine scheduling SIMPL model The CP problem goes into the Benders subproblem. A relaation of the constraint goes into the master Slide 117
118 Machine scheduling SIMPL model Generate logic-based Benders cuts when the subproblem is infeasible Slide 118
119 Machine scheduling SIMPL model Benders-based search, where problem restrictions are Benders subproblems and problem relaations are master problems. Slide 119
120 Machine scheduling Computational results Long processing times Jobs Machines MILP (CPLEX 11) Nodes Sec. SIMPL Benders Iter. Cuts Sec. Slide , , , ,309 10, SIMPL results are similar to original hand-coded results.
121 Machine scheduling Computational results Short processing times Jobs Machines MILP (CPLEX 11) Nodes Sec. SIMPL Benders Iter. Cuts Sec , > 27.5 mil. > 48 hr > 5.4 mil. > 19 hr* Slide 121 *out of memory
122 Machine scheduling Benders cut for minimum makespan, with a cumulative scheduling subproblem M M p (1 ) + ma d min d i { } { } * i i i Ji Ji J Jobs currently assigned to machine i Minimum makespan on machine i for obs currently assigned Slide 122
123 Logic-based Benders Decomposition In general, Benders cuts are obtained by solving the inference dual of the subproblem. The dual solution is a proof of optimality. LP dual is a special case, where the proof is encoded by dual multipliers. Slide 123
124 Logic-based Benders Decomposition In general, Benders cuts are obtained by solving the inference dual of the subproblem. The dual solution is a proof of optimality. LP dual is a special case, where the proof is encoded by dual multipliers. The Benders cut states conditions on the master problem variables under which the proof remains valid. Classical Benders cut is a special case. Slide 124
125 Logic-based Benders Decomposition In general, Benders cuts are obtained by solving the inference dual of the subproblem. The dual solution is a proof of optimality. LP dual is a special case, where the proof is encoded by dual multipliers. The Benders cut states conditions on the master problem variables under which the proof remains valid. Classical Benders cut is a special case. LP, surrogate, Lagrangean, and superadditive duals are special cases of inference duality and relaation duality. Whence the prevalence of relaation and inference dualities in problem solving. Slide 125
126 Truss Structure Design Select size of each bar (possibly zero) to support the load while minimizing weight. 10-bar cantilever truss 0 deg. freedom 2 deg. freedom Slide 126 Total 8 degrees of freedom Load
127 Truss Structure Design Notation v i = elongation of bar s i = force along bar h = length of bar i i d = node displacement p = load along d.f. A i = cross-sectional area of bar Slide 127
128 Truss Structure Design min i i i } Minimize total weight s. t. cos θ s = p, all i h A i i cos θ d = v, alli i i } Equilibrium } Compatibility nonlinear Ei Av i i = s i, all i h i L U v v v, alli L U d d d, all \/ i i i ( A = A ) k i ik } Hooke s law } Elongation bounds } Displacement bounds } Logical disunction Area must be one of several discrete values A ik Slide 128 Constraints can be imposed for multiple loading conditions
129 Truss Structure Design Introducing new variables linearizes the problem but makes it much larger. MILP model Slide 129 min h A y i i ik ik k s. t. cos θ s = p, all i E h i i i i cos θ d = v, alli k i ik k A v = s, alli ik ik i L U v v v, alli i i i L U d d d, all k y ik = 1, all i 0-1 variables indicating size of bar i Elongation variable disaggregated by bar size Hooke s law becomes linear
130 Truss Structure Design Integrated approach Use the original model (don t introduce new variables) Branch by splitting the range of areas A i Generate a quasi-relaation, which is linear and much smaller than MILP model. Original hand-coded method: Bollapragada, Ghattas, and JNH Slide 130
131 Truss Structure Design Theorem (JNH 2005) If g(,y) is semihomogeneous in R n and concave in scalar y, then the following is a quasi-relaation of g(,y) 0: g y g y 1 2 (, L ) + (, U ) 0 L 1 U α α L 2 (1 α ) (1 α ) = U Slide 131
132 Truss Structure Design Theorem (JNH 2005) If g(,y) is semihomogeneous in R n and concave in scalar y, then the following is a quasi-relaation of g(,y) 0: g y g y 1 2 (, L ) + (, U ) 0 L 1 U α α L 2 (1 α ) (1 α ) = U Its optimal value is a lower bound on the optimal value of the original problem, if cost is a function of alone. Slide 132
133 Truss Structure Design Theorem (JNH 2005) If g(,y) is semihomogeneous in R n and concave in scalar y, then the following is a quasi-relaation of g(,y) 0: g y g y 1 2 (, L ) + (, U ) 0 L 1 U α α L 2 (1 α ) (1 α ) = U g( α, y ) αg(, y ) for all, y and α [0,1] g(0, y ) = 0 for all y Slide 133
134 Truss Structure Design Theorem (JNH 2005) If g(,y) is semihomogeneous in R n and concave in scalar y, then the following is a quasi-relaation of g(,y) 0: g y g y 1 2 (, L ) + (, U ) 0 L 1 U α α L 2 (1 α ) (1 α ) = U Bounds on y Slide 134
135 Truss Structure Design Theorem (JNH 2005) If g(,y) is semihomogeneous in R n and concave in scalar y, then the following is a quasi-relaation of g(,y) 0: g y g y 1 2 (, L ) + (, U ) 0 L 1 U α α L 2 (1 α ) (1 α ) = U Bounds on Slide 135
136 Truss Structure Design Theorem (JNH 2005) If g(,y) is semihomogeneous in R n and concave in scalar y, then the following is a quasi-relaation of g(,y) 0: g y g y 1 2 (, L ) + (, U ) 0 L 1 U α α L 2 (1 α ) (1 α ) = U Ei Av i i = s i hi Slide 136 has the form g(,y) = 0 with g semihomogenous in because we can write it with = (A i,s i ), y = v i. Ei Av i i s i hi = 0
137 Truss Structure Design So we have a quasi-relaation of the truss problem: Slide 137 L U min h [ A y + A (1 y )] i i i i i i s. t. cos θ s = p, all i i i cos θ d = v + v, alli i i 0 i1 E ( A v + A v ) = s, all i h i L U i i 0 i i1 i i L U i i i 0 i i, all L U i i i1 i i L U, all v y v v y i v (1 y ) v v (1 y ), alli d d d 0 y 1, all i i Hooke s law is linearized Elongation bounds split into 2 sets of bounds
138 Truss Structure Design Logic cuts v i0 and v i1 must have same sign in a feasible solution. If not, we branch by adding logic cuts v, v 0, v, v 0 i 0 i1 i 0 i1 Slide 138
139 Truss Structure Design In general, we can have a metaconstraint to represent the semihomogeneous constraint g(,y) 0 and generate a quasirelaation. Since a bilinear constraint y = α is always semihomogeneous, we will use a bilinear metaconstraint with a quasi-relaation option. Slide 139
140 Truss Structure Design SIMPL model Recognized as linear systems Slide 140
141 Truss Structure Design SIMPL model Recognized as bilinear system Slide 141
142 Truss Structure Design SIMPL model Generate quasirelaation for semihomogenous function Slide 142
143 Truss Structure Design SIMPL model Branch first on violated logic cuts for quasirelaation Slide 143
144 Truss Structure Design SIMPL model Then branch on A i in-domain constraint. Violated when A i is not one of the discrete bar sizes. Take upper branch first. Slide 144
145 Truss Structure Design 10-bar cantilever truss Load Slide 145
146 Truss Structure Design Computational results (seconds) Hand-coded integrated method No. bars Loads BARON CPLEX 11 Hand coded SIMPL * * Slide 146 *plus displacement bounds
147 Truss Structure Design 25-bar problem Slide 147
148 Truss Structure Design 72-bar problem Slide 148
149 Truss Structure Design Computational results (seconds) Hand-coded integrated method No. bars Loads BARON CPLEX 11 Hand coded SIMPL , , , > 24 hr* > 24 hr* > 24 hr* > 24 hr** > 24 hr*** Slide 149 * no feasible solution found ** best feasible solution has cost 32,748 *** best feasible solution has cost 32,700
150 Current Version of SIMPL To download: Click the link to SIMPL on John Hooker s website. See readme file for complete instructions. Download eecutable and associated files Operational on GNU/Linu only Requires subsidiary solvers CPLEX (version 9, 10, or 11) Eclipse (any version or later), free download Download problem instances Including all reported in this talk. Slide 150
151 Proposal for SCIP/SIMPL cooperation A unique opportunity What SCIP/SIMPL can offer Primary tasks Slide 151
152 A unique opportunity Net generation of general-purpose solvers High-level modeling that integrates CP, MIP, GO, LS. Efficient code Eploits program micro-structure through meta-constraints Non-commercial SCIP brings efficient algorithmic tools Advanced MIP technology. SIMPL brings eperience in integrated problem-solving High-level modeling Slide 152
153 SCIP offers highly efficient noncommercial LP solver SIMPL currently uses CPLEX. well-developed MIP code. SIMPL has rudimentary MIP features. many cutting planes SIMPL implements few. constraint handles For user-supplied constraints Slide 153
154 SIMPL offers high-level modeling language SCIP uses ZIMPL language. built-in integration of CP, MIP, GO. SCIP leaves this largely to user. some CP filters (needs more) SCIP leaves it to user to write constraint handles. nogood-based search (logic-based Benders) Not available in SCIP Slide 154
155 Primary tasks Implement library of metaconstraints Design high-level modeling language or user interface Filters Relaations / conveification / rules for pooling relaations Cutting planes / disunctive MIP formulations Constraint-based branching choices Implement single branching scheme for CP & MIP Combine propagation and relaation Implement other search schemes Nogood-based search Continuous global optimization Local search? Perhaps integrate these with CP/MIP branching Slide 155
An Integrated Approach to Truss Structure Design
Slide 1 An Integrated Approach to Truss Structure Design J. N. Hooker Tallys Yunes CPAIOR Workshop on Hybrid Methods for Nonlinear Combinatorial Problems Bologna, June 2010 How to Solve Nonlinear Combinatorial
More informationIntegrating Solution Methods. through Duality. US-Mexico Workshop on Optimization. Oaxaca, January John Hooker
Integrating Solution Methods through Duality John Hooker US-Meico Workshop on Optimization Oaaca, January 2011 Zapotec Duality mask Late classical period Represents life/death, day/night, heat/cold See
More informationA Framework for Integrating Optimization and Constraint Programming
A Framework for Integrating Optimization and Constraint Programming John Hooker Carnegie Mellon University SARA 2007 SARA 07 Slide Underlying theme Model vs. solution method. How can we match the model
More informationLessons from MIP Search. John Hooker Carnegie Mellon University November 2009
Lessons from MIP Search John Hooker Carnegie Mellon University November 2009 Outline MIP search The main ideas Duality and nogoods From MIP to AI (and back) Binary decision diagrams From MIP to constraint
More informationOperations Research Methods in Constraint Programming
Operations Research Methods in Constraint Programming John Hooker Carnegie Mellon University Prepared for Lloret de Mar, Spain June 2007 2007 Slide 1 CP and OR Have Complementary Strengths CP: Inference
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 informationCombining Optimization and Constraint Programming
Combining Optimization and Constraint Programming John Hooker Carnegie Mellon University GE Research Center 7 September 2007 GE 7 Sep 07 Slide 1 Optimization and Constraint Programming Optimization: focus
More informationConsistency as Projection
Consistency as Projection John Hooker Carnegie Mellon University INFORMS 2015, Philadelphia USA Consistency as Projection Reconceive consistency in constraint programming as a form of projection. For eample,
More informationA Framework for Integrating Exact and Heuristic Optimization Methods
A Framework for Integrating Eact and Heuristic Optimization Methods John Hooker Carnegie Mellon University Matheuristics 2012 Eact and Heuristic Methods Generally regarded as very different. Eact methods
More informationProjection, Inference, and Consistency
Projection, Inference, and Consistency John Hooker Carnegie Mellon University IJCAI 2016, New York City A high-level presentation. Don t worry about the details. 2 Projection as a Unifying Concept Projection
More informationProjection, Consistency, and George Boole
Projection, Consistency, and George Boole John Hooker Carnegie Mellon University CP 2015, Cork, Ireland Projection as a Unifying Concept Projection underlies both optimization and logical inference. Optimization
More informationA Search-Infer-and-Relax Framework for Integrating Solution Methods
Carnegie Mellon University Research Showcase @ CMU Tepper School of Business 2005 A Search-Infer-and-Relax Framework for Integrating Solution Methods John N. Hooker Carnegie Mellon University, john@hooker.tepper.cmu.edu
More informationAn Integrated Solver for Optimization Problems
WORKING PAPER WPS-MAS-08-01 Department of Management Science School of Business Administration, University of Miami Coral Gables, FL 33124-8237 Created: December 2005 Last update: July 2008 An Integrated
More informationLogic, Optimization and Data Analytics
Logic, Optimization and Data Analytics John Hooker Carnegie Mellon University United Technologies Research Center, Cork, Ireland August 2015 Thesis Logic and optimization have an underlying unity. Ideas
More informationDuality in Optimization and Constraint Satisfaction
Duality in Optimization and Constraint Satisfaction J. N. Hooker Carnegie Mellon University, Pittsburgh, USA john@hooker.tepper.cmu.edu Abstract. We show that various duals that occur in optimization and
More informationLogic-Based Benders Decomposition
Logic-Based Benders Decomposition J. N. Hooker Graduate School of Industrial Administration Carnegie Mellon University, Pittsburgh, PA 15213 USA G. Ottosson Department of Information Technology Computing
More informationScheduling Home Hospice Care with Logic-Based Benders Decomposition
Scheduling Home Hospice Care with Logic-Based Benders Decomposition John Hooker Carnegie Mellon University Joint work with Aliza Heching Ryo Kimura Compassionate Care Hospice CMU Lehigh University October
More informationSingle-Facility Scheduling by Logic-Based Benders Decomposition
Single-Facility Scheduling by Logic-Based Benders Decomposition Elvin Coban J. N. Hooker Tepper School of Business Carnegie Mellon University ecoban@andrew.cmu.edu john@hooker.tepper.cmu.edu September
More informationAn Integrated Method for Planning and Scheduling to Minimize Tardiness
An Integrated Method for Planning and Scheduling to Minimize Tardiness J N Hooker Carnegie Mellon University john@hookerteppercmuedu Abstract We combine mixed integer linear programming (MILP) and constraint
More informationHow to Relax. CP 2008 Slide 1. John Hooker Carnegie Mellon University September 2008
How to Relax Slide 1 John Hooker Carnegie Mellon University September 2008 Two ways to relax Relax your mind and body. Relax your problem formulations. Slide 2 Relaxing a problem Feasible set of original
More informationDecision Diagrams: Tutorial
Decision Diagrams: Tutorial John Hooker Carnegie Mellon University CP Summer School Cork, Ireland, June 2016 Decision Diagrams Used in computer science and AI for decades Logic circuit design Product configuration
More informationGlobal Optimization of Truss. Structure Design INFORMS J. N. Hooker. Tallys Yunes. Slide 1
Slde 1 Global Optmzaton of Truss Structure Desgn J. N. Hooker Tallys Yunes INFORMS 2010 Truss Structure Desgn Select sze of each bar (possbly zero) to support the load whle mnmzng weght. Bar szes are dscrete.
More informationBenders decomposition [7, 17] uses a problem-solving strategy that can be generalized to a larger context. It assigns some of the variables trial valu
Logic-Based Benders Decomposition Λ J. N. Hooker Graduate School of Industrial Administration Carnegie Mellon University, Pittsburgh, PA 15213 USA G. Ottosson Department of Information Technology Computing
More informationTutorial: Operations Research in Constraint Programming
Tutorial: Operations Research in Constraint Programming John Hooker Carnegie Mellon University May 2009 Revised June 2009 May 2009 Slide 1 Why Integrate OR and CP? Complementary strengths Computational
More informationInteger Programming and Branch and Bound
Courtesy of Sommer Gentry. Used with permission. Integer Programming and Branch and Bound Sommer Gentry November 4 th, 003 Adapted from slides by Eric Feron and Brian Williams, 6.40, 00. Integer Programming
More informationPlanning and Scheduling by Logic-Based Benders Decomposition
Planning and Scheduling by Logic-Based Benders Decomposition J. N. Hooker Carnegie Mellon University john@hooker.tepper.cmu.edu July 2004, Revised October 2005 Abstract We combine mixed integer linear
More informationInteger Programming, Constraint Programming, and their Combination
Integer Programming, Constraint Programming, and their Combination Alexander Bockmayr Freie Universität Berlin & DFG Research Center Matheon Eindhoven, 27 January 2006 Discrete Optimization General framework
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 informationDecision Diagrams for Discrete Optimization
Decision Diagrams for Discrete Optimization Willem Jan van Hoeve Tepper School of Business Carnegie Mellon University www.andrew.cmu.edu/user/vanhoeve/mdd/ Acknowledgments: David Bergman, Andre Cire, Samid
More informationIntegrating CP and Mathematical Programming
Integrating CP and Mathematical Programming John Hooker Carnegie Mellon University June 2011 June 2011 Slide 1 Why Integrate CP and MP? Complementary strengths Computational advantages Outline of the Tutorial
More informationOptimization in Process Systems Engineering
Optimization in Process Systems Engineering M.Sc. Jan Kronqvist Process Design & Systems Engineering Laboratory Faculty of Science and Engineering Åbo Akademi University Most optimization problems in production
More informationRobust Scheduling with Logic-Based Benders Decomposition
Robust Scheduling with Logic-Based Benders Decomposition Elvin Çoban and Aliza Heching and J N Hooker and Alan Scheller-Wolf Abstract We study project scheduling at a large IT services delivery center
More informationInteger Programming Part II
Be the first in your neighborhood to master this delightful little algorithm. Integer Programming Part II The Branch and Bound Algorithm learn about fathoming, bounding, branching, pruning, and much more!
More informationLift-and-Project cuts: an efficient solution method for mixed-integer programs
Lift-and-Project cuts: an efficient solution method for mied-integer programs Sebastian Ceria Graduate School of Business and Computational Optimization Research Center http://www.columbia.edu/~sc44 Columbia
More informationOptimization Methods in Logic
Optimization Methods in Logic John Hooker Carnegie Mellon University February 2003, Revised December 2008 1 Numerical Semantics for Logic Optimization can make at least two contributions to boolean logic.
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 informationStochastic Decision Diagrams
Stochastic Decision Diagrams John Hooker CORS/INFORMS Montréal June 2015 Objective Relaxed decision diagrams provide an generalpurpose method for discrete optimization. When the problem has a dynamic programming
More informationDecision Diagram Relaxations for Integer Programming
Decision Diagram Relaxations for Integer Programming Christian Tjandraatmadja April, 2018 Tepper School of Business Carnegie Mellon University Submitted to the Tepper School of Business in Partial Fulfillment
More informationAlternative Methods for Obtaining. Optimization Bounds. AFOSR Program Review, April Carnegie Mellon University. Grant FA
Alternative Methods for Obtaining Optimization Bounds J. N. Hooker Carnegie Mellon University AFOSR Program Review, April 2012 Grant FA9550-11-1-0180 Integrating OR and CP/AI Early support by AFOSR First
More informationOperations Research as Empirical Modeling
Operations Research as Empirical Modeling John Hooker Carnegie Mellon University June 2009 June 2009 Slide 1 OR as Modeling Current practice: research in OR is primarily about algorithms. Computational
More informationConstraint Programming Overview based on Examples
DM841 Discrete Optimization Part I Lecture 2 Constraint Programming Overview based on Examples Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. An
More informationSolving Mixed-Integer Nonlinear Programs
Solving Mixed-Integer Nonlinear Programs (with SCIP) Ambros M. Gleixner Zuse Institute Berlin MATHEON Berlin Mathematical School 5th Porto Meeting on Mathematics for Industry, April 10 11, 2014, Porto
More informationComputational Integer Programming. Lecture 2: Modeling and Formulation. Dr. Ted Ralphs
Computational Integer Programming Lecture 2: Modeling and Formulation Dr. Ted Ralphs Computational MILP Lecture 2 1 Reading for This Lecture N&W Sections I.1.1-I.1.6 Wolsey Chapter 1 CCZ Chapter 2 Computational
More informationRecent Developments in Integrated Methods for Optimization
Recent Developments in Integrated Methods for Optimization John Hooker Carnegie Mellon University August 2012 Premise Constraint programming and mathematical programming can work together. There is underlying
More informationOrbital Shrinking: A New Tool for Hybrid MIP/CP Methods
Orbital Shrinking: A New Tool for Hybrid MIP/CP Methods Domenico Salvagnin DEI, University of Padova salvagni@dei.unipd.it Abstract. Orbital shrinking is a newly developed technique in the MIP community
More informationCutting Planes in SCIP
Cutting Planes in SCIP Kati Wolter Zuse-Institute Berlin Department Optimization Berlin, 6th June 2007 Outline 1 Cutting Planes in SCIP 2 Cutting Planes for the 0-1 Knapsack Problem 2.1 Cover Cuts 2.2
More informationMixed Integer Programming Solvers: from Where to Where. Andrea Lodi University of Bologna, Italy
Mixed Integer Programming Solvers: from Where to Where Andrea Lodi University of Bologna, Italy andrea.lodi@unibo.it November 30, 2011 @ Explanatory Workshop on Locational Analysis, Sevilla A. Lodi, MIP
More informationarxiv: v1 [cs.cc] 5 Dec 2018
Consistency for 0 1 Programming Danial Davarnia 1 and J. N. Hooker 2 1 Iowa state University davarnia@iastate.edu 2 Carnegie Mellon University jh38@andrew.cmu.edu arxiv:1812.02215v1 [cs.cc] 5 Dec 2018
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 informationScheduling with Constraint Programming. Job Shop Cumulative Job Shop
Scheduling with Constraint Programming Job Shop Cumulative Job Shop CP vs MIP: Task Sequencing We need to sequence a set of tasks on a machine Each task i has a specific fixed processing time p i Each
More informationRecoverable Robustness in Scheduling Problems
Master Thesis Computing Science Recoverable Robustness in Scheduling Problems Author: J.M.J. Stoef (3470997) J.M.J.Stoef@uu.nl Supervisors: dr. J.A. Hoogeveen J.A.Hoogeveen@uu.nl dr. ir. J.M. van den Akker
More informationECE Optimization for wireless networks Final. minimize f o (x) s.t. Ax = b,
ECE 788 - Optimization for wireless networks Final Please provide clear and complete answers. PART I: Questions - Q.. Discuss an iterative algorithm that converges to the solution of the problem minimize
More informationInteger vs. constraint programming. IP vs. CP: Language
Discrete Math for Bioinformatics WS 0/, by A. Bockmayr/K. Reinert,. Januar 0, 0:6 00 Integer vs. constraint programming Practical Problem Solving Model building: Language Model solving: Algorithms IP vs.
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 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 informationarxiv:cs/ v2 [cs.dm] 21 Aug 2001
Solving Assembly Line Balancing Problems by Combining IP and CP Alexander Bockmayr and Nicolai Pisaruk arxiv:cs/0106002v2 [cs.dm] 21 Aug 2001 Université Henri Poincaré, LORIA B.P. 239, F-54506 Vandœuvre-lès-Nancy,
More informationInexact Solution of NLP Subproblems in MINLP
Ineact Solution of NLP Subproblems in MINLP M. Li L. N. Vicente April 4, 2011 Abstract In the contet of conve mied-integer nonlinear programming (MINLP, we investigate how the outer approimation method
More information3.4 Relaxations and bounds
3.4 Relaxations and bounds Consider a generic Discrete Optimization problem z = min{c(x) : x X} with an optimal solution x X. In general, the algorithms generate not only a decreasing sequence of upper
More informationProbabilistic Graphical Models
Probabilistic Graphical Models David Sontag New York University Lecture 6, March 7, 2013 David Sontag (NYU) Graphical Models Lecture 6, March 7, 2013 1 / 25 Today s lecture 1 Dual decomposition 2 MAP inference
More informationMINLP: Theory, Algorithms, Applications: Lecture 3, Basics of Algorothms
MINLP: Theory, Algorithms, Applications: Lecture 3, Basics of Algorothms Jeff Linderoth Industrial and Systems Engineering University of Wisconsin-Madison Jonas Schweiger Friedrich-Alexander-Universität
More informationDecision Diagrams for Sequencing and Scheduling
Decision Diagrams for Sequencing and Scheduling Willem-Jan van Hoeve Tepper School of Business Carnegie Mellon University www.andrew.cmu.edu/user/vanhoeve/mdd/ Plan What can MDDs do for Combinatorial Optimization?
More informationMDD-based Postoptimality Analysis for Mixed-integer Programs
MDD-based Postoptimality Analysis for Mixed-integer Programs John Hooker, Ryo Kimura Carnegie Mellon University Thiago Serra Mitsubishi Electric Research Laboratories Symposium on Decision Diagrams for
More informationThe Separation Problem for Binary Decision Diagrams
The Separation Problem for Binary Decision Diagrams J. N. Hooker Joint work with André Ciré Carnegie Mellon University ISAIM 2014 Separation Problem in Optimization Given a relaxation of an optimization
More informationMVE165/MMG630, Applied Optimization Lecture 6 Integer linear programming: models and applications; complexity. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming: models and applications; complexity Ann-Brith Strömberg 2011 04 01 Modelling with integer variables (Ch. 13.1) Variables Linear programming (LP) uses continuous
More informationLecture 23 Branch-and-Bound Algorithm. November 3, 2009
Branch-and-Bound Algorithm November 3, 2009 Outline Lecture 23 Modeling aspect: Either-Or requirement Special ILPs: Totally unimodular matrices Branch-and-Bound Algorithm Underlying idea Terminology Formal
More informationThe core of solving constraint problems using Constraint Programming (CP), with emphasis on:
What is it about? l Theory The core of solving constraint problems using Constraint Programming (CP), with emphasis on: l Modeling l Solving: Local consistency and propagation; Backtracking search + heuristics.
More informationCritical Reading of Optimization Methods for Logical Inference [1]
Critical Reading of Optimization Methods for Logical Inference [1] Undergraduate Research Internship Department of Management Sciences Fall 2007 Supervisor: Dr. Miguel Anjos UNIVERSITY OF WATERLOO Rajesh
More informationof a bimatrix game David Avis McGill University Gabriel Rosenberg Yale University Rahul Savani University of Warwick
Finding all Nash equilibria of a bimatrix game David Avis McGill University Gabriel Rosenberg Yale University Rahul Savani University of Warwick Bernhard von Stengel London School of Economics Nash equilibria
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 informationIndicator Constraints in Mixed-Integer Programming
Indicator Constraints in Mixed-Integer Programming Andrea Lodi University of Bologna, Italy - andrea.lodi@unibo.it Amaya Nogales-Gómez, Universidad de Sevilla, Spain Pietro Belotti, FICO, UK Matteo Fischetti,
More informationResource Constrained Project Scheduling Linear and Integer Programming (1)
DM204, 2010 SCHEDULING, TIMETABLING AND ROUTING Lecture 3 Resource Constrained Project Linear and Integer Programming (1) Marco Chiarandini Department of Mathematics & Computer Science University of Southern
More informationInteger Linear Programming (ILP)
Integer Linear Programming (ILP) Zdeněk Hanzálek, Přemysl Šůcha hanzalek@fel.cvut.cz CTU in Prague March 8, 2017 Z. Hanzálek (CTU) Integer Linear Programming (ILP) March 8, 2017 1 / 43 Table of contents
More informationChapter 11. Approximation Algorithms. Slides by Kevin Wayne Pearson-Addison Wesley. All rights reserved.
Chapter 11 Approximation Algorithms Slides by Kevin Wayne. Copyright @ 2005 Pearson-Addison Wesley. All rights reserved. 1 P and NP P: The family of problems that can be solved quickly in polynomial time.
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 information4/12/2011. Chapter 8. NP and Computational Intractability. Directed Hamiltonian Cycle. Traveling Salesman Problem. Directed Hamiltonian Cycle
Directed Hamiltonian Cycle Chapter 8 NP and Computational Intractability Claim. G has a Hamiltonian cycle iff G' does. Pf. Suppose G has a directed Hamiltonian cycle Γ. Then G' has an undirected Hamiltonian
More informationR u t c o r Research R e p o r t. A Method to Schedule Both Transportation and Production at the Same Time in a Special FMS.
R u t c o r Research R e p o r t A Method to Schedule Both Transportation and Production at the Same Time in a Special FMS Navid Hashemian a Béla Vizvári b RRR 12-08, September 2008 RUTCOR Rutgers Center
More informationChapter 11. Approximation Algorithms. Slides by Kevin Wayne Pearson-Addison Wesley. All rights reserved.
Chapter 11 Approximation Algorithms Slides by Kevin Wayne. Copyright @ 2005 Pearson-Addison Wesley. All rights reserved. 1 Approximation Algorithms Q. Suppose I need to solve an NP-hard problem. What should
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 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 informationMVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms
MVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms Ann-Brith Strömberg 2017 04 07 Lecture 8 Linear and integer optimization with applications
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 informationThe two-dimensional bin-packing problem is the problem of orthogonally packing a given set of rectangles
INFORMS Journal on Computing Vol. 19, No. 1, Winter 2007, pp. 36 51 issn 1091-9856 eissn 1526-5528 07 1901 0036 informs doi 10.1287/ijoc.1060.0181 2007 INFORMS Using Decomposition Techniques and Constraint
More informationComputational Mixed-Integer Programming
Computational Mixed-Integer Programming Ambros Gleixner and the SCIP team Zuse Institute Berlin gleixner@zib.de SCIP Optimization Suite http://scip.zib.de Theory and Practice of Satisfiability Solving
More informationMultivalued Decision Diagrams. Postoptimality Analysis Using. J. N. Hooker. Tarik Hadzic. Cork Constraint Computation Centre
Postoptimality Analysis Using Multivalued Decision Diagrams Tarik Hadzic Cork Constraint Computation Centre J. N. Hooker Carnegie Mellon University London School of Economics June 2008 Postoptimality Analysis
More informationSystem Planning Lecture 7, F7: Optimization
System Planning 04 Lecture 7, F7: Optimization System Planning 04 Lecture 7, F7: Optimization Course goals Appendi A Content: Generally about optimization Formulate optimization problems Linear Programming
More informationStochastic Integer Programming
IE 495 Lecture 20 Stochastic Integer Programming Prof. Jeff Linderoth April 14, 2003 April 14, 2002 Stochastic Programming Lecture 20 Slide 1 Outline Stochastic Integer Programming Integer LShaped Method
More informationInteger Programming Duality
Integer Programming Duality M. Guzelsoy T. K. Ralphs July, 2010 1 Introduction This article describes what is known about duality for integer programs. It is perhaps surprising that many of the results
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 informationDiscrete (and Continuous) Optimization WI4 131
Discrete (and Continuous) Optimization WI4 131 Kees Roos Technische Universiteit Delft Faculteit Electrotechniek, Wiskunde en Informatica Afdeling Informatie, Systemen en Algoritmiek e-mail: C.Roos@ewi.tudelft.nl
More informationSolving LP and MIP Models with Piecewise Linear Objective Functions
Solving LP and MIP Models with Piecewise Linear Obective Functions Zonghao Gu Gurobi Optimization Inc. Columbus, July 23, 2014 Overview } Introduction } Piecewise linear (PWL) function Convex and convex
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 informationInteger Linear Programming
Integer Linear Programming Solution : cutting planes and Branch and Bound Hugues Talbot Laboratoire CVN April 13, 2018 IP Resolution Gomory s cutting planes Solution branch-and-bound General method Resolution
More informationFrom structures to heuristics to global solvers
From structures to heuristics to global solvers Timo Berthold Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies OR2013, 04/Sep/13, Rotterdam Outline From structures to
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 informationGestion de la production. Book: Linear Programming, Vasek Chvatal, McGill University, W.H. Freeman and Company, New York, USA
Gestion de la production Book: Linear Programming, Vasek Chvatal, McGill University, W.H. Freeman and Company, New York, USA 1 Contents 1 Integer Linear Programming 3 1.1 Definitions and notations......................................
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 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 informationMixed Integer Non Linear Programming
Mixed Integer Non Linear Programming Claudia D Ambrosio CNRS Research Scientist CNRS & LIX, École Polytechnique MPRO PMA 2016-2017 Outline What is a MINLP? Dealing with nonconvexities Global Optimization
More informationis called an integer programming (IP) problem. model is called a mixed integer programming (MIP)
INTEGER PROGRAMMING Integer Programming g In many problems the decision variables must have integer values. Example: assign people, machines, and vehicles to activities in integer quantities. If this is
More informationIntSat: From SAT to Integer Linear Programming
IntSat: From SAT to Integer Linear Programming CPAIOR 2015 (invited talk) Robert Nieuwenhuis Barcelogic.com - Computer Science Department BarcelonaTech (UPC) 1 Proposed travel arrangements (next time):
More information