Open-pit mining problem and the BZ algorithm
|
|
- Albert Perkins
- 5 years ago
- Views:
Transcription
1 Open-pit mining problem and the BZ algorithm Eduardo Moreno (Universidad Adolfo Ibañez) Daniel Espinoza (Universidad de Chile) Marcos Goycoolea (Universidad Adolfo Ibañez) Gonzalo Muñoz (Ph.D. student Columbia Univ.)
2 Outline The Open-Pit Mine Scheduling Problem BZ algorithm The General Algorithm Template Uncertainty in Open-pite mining BZ Algorithm & Stochastic Programming
3 Open pit mining scheduling problem Which block should be extracted, when should they be extracted and how should they be processed First integer programming model due to Thys Johnson ( 69)
4 The problem Assumptions & Notation blocks B T periods Precedences G =(B, A) By decision variables x b,t Side constraints Dx apple d Discount rate ĉ b,t = t c b
5 The problem Time-Consistency constraints: x b,t x b,t+1 Precedence (wall-slope) constraints: x a,t x b,t Capacity (knapsack) constraints: b B a b (x b,t x b,t 1 ) C t Multi-period precedence constrained knapsack
6 The problem The by formulation (PCP by ) max cx x b,t apple x b,t+1, 8 t<t 1 x a,t apple x b,t, 8(a, b) 2 A, 8t Dx apple d x b,t 2 [0, 1]
7 Computational results (LP Relaxation) We count with several REAL instances of this problem. Most of them can t be solved using CPLEX 12.1 Intel(R) Xeon 2.33 GHz, 16 GB RAM. Mine Blocks Prec. Periods CPX1Time Tinkerbell PeterPan Pinochio Mowgli Hansel MEM! Ariel MEM! SnowWhite MEM! GreEel MEM! Simba CinderellaILimit MEM! Dorothy MEM! Rapunzel MEM! Beauty MEM! Aladdin MEM! Alice MEM! BigBadWolf MEM! Cinderella MEM!
8 Ad-hoc LP Algorithms Critical Multiplier Algorithm max cx st x P i apple x j 8(i, j) 2 A ai x i apple b x i 2 [0, 1] max(c a)x st x i apple x j 8(i, j) 2 A x i 2 [0, 1] Explicit construction of optimal multipliers BZ Algorithm
9 The BZ algorithm 1. Start with µ 1 = 0, C 0 1 = B, C 0 = {C 0 1} and k = Let y k be the solution of L(P, µ k 1 ), problem where the side-constraints are penalized with multipliers µ k 1. If k>1 and y k satisfies H k 1 x = 0 stop. 3. Define C k partition of B and impose the constraints H k x = 0 consistent of x i = x j, 8i, j 2 C k h. Also the solution yk must satisfy the constraints H k x = Solve P 2 k, obtained by adding the constraints H k x = 0 to the problem P, and call the solution x k with dual variables µ k associated to the sideconstraints. If µ k = µ k 1 stop. 5. Set k k + 1 and go to step 2.
10 BZ algorithm The algorithm L(P, µ k 1 ) Max-Closure Get duals and penalize with that duals Refine partitions P 2 k Small Problem
11 BZ algorithm P 2 0 Define Solve L(P,~0) Graphically and solve it get duals µ 0 (side constraints)
12 BZ algorithm Refine Solve L(P, the partition µ 0 )
13 BZ algorithm Define Refine Pthe 2 1 partition and solve it get duals µ 1 (side constraints)
14 BZ algorithm And on...
15 BZ algorithm BZ algorithm (+ a LOT of tricks) Mine Blocks Prec. Periods CPX1Time BZ1Time Tinkerbell PeterPan Pinochio Mowgli Hansel MEM! Ariel MEM! 7.89 SnowWhite MEM! GreEel MEM! Simba CinderellaILimit MEM! Dorothy MEM! Rapunzel MEM! Beauty MEM! Aladdin MEM! Alice MEM! BigBadWolf MEM! MEM! Cinderella MEM!
16 General algorithm template We consider the following problem (P ) max cx s.t Ax apple b Dx apple d Easy Hard
17 General algorithm template Feasible solution z 0 k =1 µ 0 = ~0 L(P, µ k 1 )
18 General algorithm template Feasible solution z 0 k =1 µ 0 = ~0 max cx µ k 1 (Dx d) s.t L(P, µ k 1 ) Ax apple b
19 General algorithm template Feasible solution z 0 k =1 µ 0 = ~0 L(P, µ k 1 ) w k
20 General algorithm template Feasible solution z 0 k =1 µ 0 = ~0 L(P, µ k 1 ) (P 2 k ) max cx s.t P 2 k Ax apple b Dx apple d H k x = h k
21 General algorithm template Feasible solution z 0 k =1 µ 0 = ~0 L(P, µ k 1 ) Get duals and penalize with that duals P 2 k
22 General algorithm template Clearly H k x = h k defines an affine subspace, h k =0 but for simplicity let s assume so the subspace is linear
23 General algorithm template This way we can rewrite the P 2 k problem as (P 2 k ) max c X i ix i! s.t A X i ix i! apple b D X i ix i! apple d where x i is a generator of {x H i k x =0}
24 Differences? BZ Dantzig-Wolfe In the master problem we move in an affine subspace that contains the generated point We have the freedom of choosing the affine subspace In the master problem we move in the convex hull of the generated extreme points We add one column (extreme point) at a time We solve the same pricing problem in both
25 Graphically Maybe DW would need to generate... before moving out of the upper lid
26 Graphically But BZ only needs... to describe the entire upper lid
27 Case Study Selective Mine, two metallurgical process Prices, Met. Recoveries, Mining, Process and Selling costs based on the usual numbers of the industry; Asymmetrical contribution of value by process: 10% of the reserves (Process A) give 45% of the cash flow of the project. Some comparison between both process (Process A / Process B): Process Cost 12:1
28 Inventory Reserves Total Tonnes
29 Inventory Reserves Metal
30 Inventory Reserves Cash Flow
31 Mine Schedule - Tonnage
32 Mill Feed - Tonnage
33 Mine Schedule Cumulated Metal Production
34 Mine Schedule Cash Flow
35 Ultimate Pit Whittle Pit Opt Pit
36 Mining Sequence Opt Pits Whittle Nested Pits
37 An important extension x b,d,t = 1 block b is extracted in time 1..t and destination 1..d 0 otherwise
38 Marvin Extracted/processed,weigth, Fixed, Extracted/processed'weigth' Model' " " " " " " Weigth(tons), " " " " " processed" extracted" restric5on_processed" Weight'(tons)' " " " " " restric0on_processed" restric0on_extracted" processed" " " restric5on_extracted" " " extracted" 0" 0" 0" 2" 4" 6" 8" 10" 12" 14" 16" 0" 2" 4" 6" 8" 10" 12" 14" 16" Period, Period' Grade&CU& Grade&AU& 0.7" 0.8" 0.6" 0.7" 0.5" 0.6" Grade& 0.4" 0.3" 0.2" model" fixed" Grade& 0.5" 0.4" 0.3" 0.2" model" fixed" 0.1" 0.1" 0" 0" 0" 2" 4" 6" 8" 10" 12" 14" 16" 0" 2" 4" 6" 8" 10" 12" 14" 16" Period& Period&
39 Delphos Titre,Extracted/Processed,weight, Fixed, Extracted/Processed'weight' Model' " " " " " " Weight(TONN), " " processed" extracted" Weight'(TONN)' " " processed" extracted" " " " " 0" 0" 5" 10" 15" 20" 25" 30" 35" 40" 45" 0" 0" 5" 10" 15" 20" 25" 30" 35" 40" 45" Period, Period' Average&grades&in&processed&blocks& Fixed& Average&grades&in&processed&blocks&& Model& 1.6" 1.8" 1.4" 1.6" average&grade& 1.2" 1" 0.8" 0.6" 0.4" 0.2" Cu" Au" As" average&grade& 1.4" 1.2" 1" 0.8" 0.6" 0.4" 0.2" Cu" Au" As" 0" 0" 5" 10" 15" 20" 25" 30" 35" 40" 45" 0" 0" 5" 10" 15" 20" 25" 30" 35" 40" 45" Period& Period&
40 NCL Extracted/processed'weight' Fixed' Extracted/Processed'weight' Model' " " " " Weight'(TONN)' " " " " " " processed" extracted" restric3on_processed" restric3on_extracted" Weight'(TONN)' " " " " " " restric.on_processed" restric.on_extracted" processed" extracted" 0" 0" 2" 4" 6" 8" 10" 0" 0" 2" 4" 6" 8" 10" Period' Period' Quan1ty'Su'in'processed'blocks' Mill$hours$for$processed$blocks$ " 9000" 8000" " 7000" Period'(TONN)' " " restric2on_su" model" fixed" Mill$(hours)$ 6000" 5000" 4000" 3000" restric2on_mill" model" fixed" 50000" 2000" 1000" 0" 0" 1" 2" 3" 4" 5" 6" 7" 8" 9" 10" Period' 0" 0" 1" 2" 3" 4" 5" 6" 7" 8" 9" 10" Period$
41 Uncertainty in Open-pit Mine Planning Grade Uncertainty Ore Prices x e b = ( 1 if b is extracted 0 if not x p b = ( 1 if b is processed 0 if not x p b apple xe b 8b Objective Function: (x e,x p )=(p p c p )x p c e x e ( b := grade of ore in block b)
42 Scenario approaches Multiple Blockmodels Opt. Multiple plans Extraction plan Multiple Blockmodels Opt. Extraction plan Most common in mining literature PH on ore prices Dimitrakopoulos et al. DP on grade uncertainty Risk minimization Robust programming
43 Our approach for ore prices uncertainty mean reversion process of ore prices (Ornstein Uhlenbeck process) (P ) max cx p p 2 U = {p p 9u, u 0 u apple 1,p p = p + u} A second-order conic constraint s.t Ax apple b Dx apple d Easy Hard
44 Other uses of BZ Algorithm Value-at-Risk VaR " ( ĉx) =min{t : F ĉx (t) 1 "} t Conditional Value-at-Risk (expected shortfall) CVaR " ( ĉx) =E( ĉx ĉx VaR " ( ĉx)) apple =min t t + 1 " E(( ĉx t)+ ) Convex function!
45 CVaR formulation probability of scenario i min CVaR " ( ĉx) s.t. Ax = b x 0 SAA min t + 1 " NX p i i i=1 s.t. Ax = b i c i x t 8i x, 0 scenario i CVaR " ( ĉx) =min t apple t + 1 " E(( ĉx t)+ )
46 Our idea x : A x = b min t + 1 " NX p i i Easy to solve! Upper bound i=1 s.t. c i x + t + i 0 8i NX 0 min t + 1 " p i i s.t. Ax = b i=1 x {P j } c i x + t + i 0 8i x, 0 min t + 1 " kx ˆp j ˆ j Smaller problem! N = {P j } k j=1 Lower bound ˆp j = X p i c j = 1ˆp X p i c i i2p j j i2p j j=1 s.t. Ax = b c j x + t +ˆ j 0 8j 2 {1,...,k} x, ˆ 0
47 Dual problems max y t b s.t. y t A + 1 " NX p i i c i apple 0 i=1 NX p i i = " i=1 0 apple i apple 1 8i x : A x = b Upper Lower bound N = {P j } k j=1 Lower Upper bound Additional constraints: i = i 0 8i, i 0 2 P j max s.t. max y t b 1 " NX i=1 p i i c i x NX p i i = " i=1 s.t. y t A + 1 " kx j=1 0 apple i apple 1 8i x kx j=1 ˆp j ˆj = " ˆp j ˆj c j apple 0 0 apple ˆj apple 1 8j 2 {1...K} (i) (i+1) = =1i apple i (i) =0i>i +1
48 Computational results for large N values for " =1%,10%,50%and75%. N =100, 000 N =1, 000, 000 N =10, 000, 000 Name time Iter. N time Iter. N time Iter. N adlittle afiro agg bandm beaconfd boeing bore3d brandy capri e etamacro f ganges grow grow grow israel kb lotfi perold pilot.ja
49 grow israel kb lotfi perold pilot.ja pilot pilot pilotnov recipe sc sc sc50a sc50b scagr scfxm scfxm scfxm share1b share2b stair standata standgub standmps stocfor tu vtp.base woodw
50 Thanks...
A primal-dual aggregation algorithm for minimizing Conditional Value-at-Risk in linear programs
Computational Optimization and Applications manuscript No. (will be inserted by the editor) A primal-dual aggregation algorithm for minimizing Conditional Value-at-Risk in linear programs Daniel Espinoza
More informationComputational Experience with Rigorous Error Bounds for the Netlib Linear Programming Library
Computational Experience with Rigorous Error Bounds for the Netlib Linear Programming Library Christian Keil (c.keil@tu-harburg.de) and Christian Jansson (jansson@tu-harburg.de) Hamburg University of Technology
More informationAN EFFICIENT APPROACH TO UPDATING SIMPLEX MULTIPLIERS IN THE SIMPLEX ALGORITHM
AN EFFICIENT APPROACH TO UPDATING SIMPLEX MULTIPLIERS IN THE SIMPLEX ALGORITHM JIAN-FENG HU AND PING-QI PAN Abstract. The simplex algorithm computes the simplex multipliers by solving a system (or two
More informationMcMaster University. Advanced Optimization Laboratory. Title: Computational Experience with Self-Regular Based Interior Point Methods
McMaster University Advanced Optimization Laboratory Title: Computational Experience with Self-Regular Based Interior Point Methods Authors: Guoqing Zhang, Jiming Peng, Tamás Terlaky, Lois Zhu AdvOl-Report
More informationA primal-dual aggregation algorithm for minimizing conditional value-at-risk in linear programs
Comput Optim Appl (2014) 59:617 638 DOI 10.1007/s10589-014-9692-6 A primal-dual aggregation algorithm for minimizing conditional value-at-risk in linear programs Daniel Espinoza Eduardo Moreno Received:
More informationResearch Article On the Simplex Algorithm Initializing
Abstract and Applied Analysis Volume 2012, Article ID 487870, 15 pages doi:10.1155/2012/487870 Research Article On the Simplex Algorithm Initializing Nebojša V. Stojković, 1 Predrag S. Stanimirović, 2
More informationPREDICTOR-CORRECTOR SMOOTHING METHODS FOR LINEAR PROGRAMS WITH A MORE FLEXIBLE UPDATE OF THE SMOOTHING PARAMETER 1
PREDICTOR-CORRECTOR SMOOTHING METHODS FOR LINEAR PROGRAMS WITH A MORE FLEXIBLE UPDATE OF THE SMOOTHING PARAMETER 1 Stephan Engelke and Christian Kanzow University of Hamburg Department of Mathematics Center
More informationHow good are projection methods for convex feasibility problems?
RAL-TR-2006-015 How good are projection methods for convex feasibility problems? N I M Gould September 21, 2006 Council for the Central Laboratory of the Research Councils c Council for the Central Laboratory
More informationOn Mehrotra-Type Predictor-Corrector Algorithms
On Mehrotra-Type Predictor-Corrector Algorithms M. Salahi, J. Peng, T. Terlaky October 10, 006 (Revised) Abstract In this paper we discuss the polynomiality of a feasible version of Mehrotra s predictor-corrector
More informationFinding an interior point in the optimal face of linear programs
Mathematical Programming 62 (1993) 497-515 497 North-Holland Finding an interior point in the optimal face of linear programs Sanjay Mehrotra* Department of lndustrial Engineering and Management Sciences,
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 informationA less conservative variant of Robust Optimization
A less conservative variant of Robust Optimization Ernst Roos, Dick den Hertog CentER and Department of Econometrics and Operations Research, Tilburg University e.j.roos@tilburguniversity.edu, d.denhertog@tilburguniversity.edu
More informationUsing Integer Programming for Strategic Underground and Open Pit-to-Underground Scheduling
Using Integer Programming for Strategic Underground and Open Pit-to-Underground Scheduling Barry King Advisor: Alexandra Newman Operations Research with Engineering PhD Program August 5-6, 2016 Cutoff
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 15, No. 3, pp. 635 653 c 2005 Society for Industrial and Applied Mathematics CONVERGENCE CONDITIONS AND KRYLOV SUBSPACE BASED CORRECTIONS FOR PRIMAL-DUAL INTERIOR-POINT METHOD SANJAY
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 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 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 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 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 informationSIMO CASE STUDY GEOVIA USER CONFERENCE PERTH CHRIS JAMES 31 MAY 2018
SIMO CASE STUDY GEOVIA USER CONFERENCE PERTH CHRIS JAMES 31 MAY 2018 2 [GEOVIA: SIMO] [May 2018] WHAT IS VALUE Important to understand potential levers that drive value Value What is it? NPV IRR Payback
More informationLIGHT ROBUSTNESS. Matteo Fischetti, Michele Monaci. DEI, University of Padova. 1st ARRIVAL Annual Workshop and Review Meeting, Utrecht, April 19, 2007
LIGHT ROBUSTNESS Matteo Fischetti, Michele Monaci DEI, University of Padova 1st ARRIVAL Annual Workshop and Review Meeting, Utrecht, April 19, 2007 Work supported by the Future and Emerging Technologies
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 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 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 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 informationTime Aggregation for Network Design to Meet Time-Constrained Demand
20th International Congress on Modelling and Simulation, Adelaide, Australia, 1 6 December 2013 www.mssanz.org.au/modsim2013 Time Aggregation for Network Design to Meet Time-Constrained Demand N. Boland
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 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 informationValid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints
Valid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints Nilay Noyan Andrzej Ruszczyński March 21, 2006 Abstract Stochastic dominance relations
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 informationIP Cut Homework from J and B Chapter 9: 14, 15, 16, 23, 24, You wish to solve the IP below with a cutting plane technique.
IP Cut Homework from J and B Chapter 9: 14, 15, 16, 23, 24, 31 14. You wish to solve the IP below with a cutting plane technique. Maximize 4x 1 + 2x 2 + x 3 subject to 14x 1 + 10x 2 + 11x 3 32 10x 1 +
More informationNonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control
Nonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control 19/4/2012 Lecture content Problem formulation and sample examples (ch 13.1) Theoretical background Graphical
More informationModeling Uncertainty in Linear Programs: Stochastic and Robust Linear Programming
Modeling Uncertainty in Linear Programs: Stochastic and Robust Programming DGA PhD Student - PhD Thesis EDF-INRIA 10 November 2011 and motivations In real life, Linear Programs are uncertain for several
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 informationA Stochastic-Oriented NLP Relaxation for Integer Programming
A Stochastic-Oriented NLP Relaxation for Integer Programming John Birge University of Chicago (With Mihai Anitescu (ANL/U of C), Cosmin Petra (ANL)) Motivation: The control of energy systems, particularly
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization
Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012
More informationMODIFICATION OF SIMPLEX METHOD AND ITS IMPLEMENTATION IN VISUAL BASIC. Nebojša V. Stojković, Predrag S. Stanimirović and Marko D.
MODIFICATION OF SIMPLEX METHOD AND ITS IMPLEMENTATION IN VISUAL BASIC Nebojša V Stojković, Predrag S Stanimirović and Marko D Petković Abstract We investigate the problem of finding the first basic solution
More informationOutline. Outline. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Scheduling CPM/PERT Resource Constrained Project Scheduling Model
Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING Lecture 3 and Mixed Integer Programg Marco Chiarandini 1. Resource Constrained Project Model 2. Mathematical Programg 2 Outline Outline 1. Resource Constrained
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 informationSection Notes 9. IP: Cutting Planes. Applied Math 121. Week of April 12, 2010
Section Notes 9 IP: Cutting Planes Applied Math 121 Week of April 12, 2010 Goals for the week understand what a strong formulations is. be familiar with the cutting planes algorithm and the types of cuts
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 informationLarge-scale Linear and Nonlinear Optimization in Quad Precision
Large-scale Linear and Nonlinear Optimization in Quad Precision Ding Ma and Michael Saunders MS&E and ICME, Stanford University US-Mexico Workshop on Optimization and its Applications Mérida, Yucatán,
More informationReformulation and Sampling to Solve a Stochastic Network Interdiction Problem
Network Interdiction Stochastic Network Interdiction and to Solve a Stochastic Network Interdiction Problem Udom Janjarassuk Jeff Linderoth ISE Department COR@L Lab Lehigh University jtl3@lehigh.edu informs
More informationStochastic Dual Dynamic Programming with CVaR Risk Constraints Applied to Hydrothermal Scheduling. ICSP 2013 Bergamo, July 8-12, 2012
Stochastic Dual Dynamic Programming with CVaR Risk Constraints Applied to Hydrothermal Scheduling Luiz Carlos da Costa Junior Mario V. F. Pereira Sérgio Granville Nora Campodónico Marcia Helena Costa Fampa
More informationWorkforce Scheduling. Outline DM87 SCHEDULING, TIMETABLING AND ROUTING. Outline. Workforce Scheduling. 1. Workforce Scheduling.
Outline DM87 SCHEDULING, TIMETABLING AND ROUTING Lecture 17 Workforce Scheduling 2. Crew Scheduling and Roering Marco Chiarandini DM87 Scheduling, Timetabling and Routing 2 Outline Workforce Scheduling
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 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 informationand to estimate the quality of feasible solutions I A new way to derive dual bounds:
Lagrangian Relaxations and Duality I Recall: I Relaxations provide dual bounds for the problem I So do feasible solutions of dual problems I Having tight dual bounds is important in algorithms (B&B), and
More informationOperations Research Lecture 6: Integer Programming
Operations Research Lecture 6: Integer Programming Notes taken by Kaiquan Xu@Business School, Nanjing University May 12th 2016 1 Integer programming (IP) formulations The integer programming (IP) is the
More informationStochastic Programming Models in Design OUTLINE
Stochastic Programming Models in Design John R. Birge University of Michigan OUTLINE Models General - Farming Structural design Design portfolio General Approximations Solutions Revisions Decision: European
More informationInteger Programming ISE 418. Lecture 12. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 12 Dr. Ted Ralphs ISE 418 Lecture 12 1 Reading for This Lecture Nemhauser and Wolsey Sections II.2.1 Wolsey Chapter 9 ISE 418 Lecture 12 2 Generating Stronger Valid
More information1. Introduction. We are concerned with the linear programming (LP) problem in the standard form. minimize subject to Ax = b, x 0,
A GENERALIZATION OF THE REVISED SIMPLEX ALGORITHM FOR LINEAR PROGRAMMING PING-QI PAN Abstract. Recently the concept of basis, which plays a fundamental role in the simplex methodology, were generalized
More informationDetermination of the transition point from open pit to underground mining
Determination of the transition point from open pit to underground mining Strategic Mine planning and Optimisation for Combination Mining Method J Chung, M Asad, E Topal, AK Ghosh Department of Mining
More informationScenario Grouping and Decomposition Algorithms for Chance-constrained Programs
Scenario Grouping and Decomposition Algorithms for Chance-constrained Programs Siqian Shen Dept. of Industrial and Operations Engineering University of Michigan Joint work with Yan Deng (UMich, Google)
More informationAdaptive Strategies for The Open-Pit Mine Optimal Scheduling Problem
Adaptive Strategies for The Open-Pit Mine Optimal Scheduling Problem Michel De Lara, Nelson Morales, Nathanaël Beeker arxiv:1706.08264v1 [math.oc] 26 Jun 2017 June 27, 2017 Abstract Within the mining discipline,
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 informationLecture 8. Strong Duality Results. September 22, 2008
Strong Duality Results September 22, 2008 Outline Lecture 8 Slater Condition and its Variations Convex Objective with Linear Inequality Constraints Quadratic Objective over Quadratic Constraints Representation
More informationIn the original knapsack problem, the value of the contents of the knapsack is maximized subject to a single capacity constraint, for example weight.
In the original knapsack problem, the value of the contents of the knapsack is maximized subject to a single capacity constraint, for example weight. In the multi-dimensional knapsack problem, additional
More information1 Column Generation and the Cutting Stock Problem
1 Column Generation and the Cutting Stock Problem In the linear programming approach to the traveling salesman problem we used the cutting plane approach. The cutting plane approach is appropriate when
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 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 informationDevelopment of the new MINLP Solver Decogo using SCIP - Status Report
Development of the new MINLP Solver Decogo using SCIP - Status Report Pavlo Muts with Norman Breitfeld, Vitali Gintner, Ivo Nowak SCIP Workshop 2018, Aachen Table of contents 1. Introduction 2. Automatic
More informationChapter 2: Linear Programming Basics. (Bertsimas & Tsitsiklis, Chapter 1)
Chapter 2: Linear Programming Basics (Bertsimas & Tsitsiklis, Chapter 1) 33 Example of a Linear Program Remarks. minimize 2x 1 x 2 + 4x 3 subject to x 1 + x 2 + x 4 2 3x 2 x 3 = 5 x 3 + x 4 3 x 1 0 x 3
More informationIntroduction to optimization and operations research
Introduction to optimization and operations research David Pisinger, Fall 2002 1 Smoked ham (Chvatal 1.6, adapted from Greene et al. (1957)) A meat packing plant produces 480 hams, 400 pork bellies, and
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 informationComputing with multi-row Gomory cuts
Computing with multi-row Gomory cuts Daniel G. Espinoza Departamento de Ingeniería Industrial, Universidad de Chile, Av. República 71, Santiago, 837-439, Chile Abstract Recent advances on the understanding
More informationTHE UNIVERSITY OF HONG KONG DEPARTMENT OF MATHEMATICS. Operations Research I
LN/MATH2901/CKC/MS/2008-09 THE UNIVERSITY OF HONG KONG DEPARTMENT OF MATHEMATICS Operations Research I Definition (Linear Programming) A linear programming (LP) problem is characterized by linear functions
More informationModeling with Integer Programming
Modeling with Integer Programg Laura Galli December 18, 2014 We can use 0-1 (binary) variables for a variety of purposes, such as: Modeling yes/no decisions Enforcing disjunctions Enforcing logical conditions
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 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 informationDecomposition Techniques in Mathematical Programming
Antonio J. Conejo Enrique Castillo Roberto Minguez Raquel Garcia-Bertrand Decomposition Techniques in Mathematical Programming Engineering and Science Applications Springer Contents Part I Motivation and
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 informationSummary of the simplex method
MVE165/MMG630, The simplex method; degeneracy; unbounded solutions; infeasibility; starting solutions; duality; interpretation Ann-Brith Strömberg 2012 03 16 Summary of the simplex method Optimality condition:
More informationEomhhhhEEEEEEEI EhhmhmhEmhhhEK
TEST PROBLEMS(U) STANFORD UNIV CA SYSTEMS OPTIMIZATION UNCLASSIFIEDb LAB I J LUSTIG UG 7 SOL-97-1± NOSS14-7-K-142 F/G 12/5 ML EomhhhhEEEEEEEI EhhmhmhEmhhhEK "V - is L S ystems in Optimization Laboratorv
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 informationInteger Programming Methods LNMB
Integer Programming Methods LNMB 2017 2018 Dion Gijswijt homepage.tudelft.nl/64a8q/intpm/ Dion Gijswijt Intro IntPM 2017-2018 1 / 24 Organisation Webpage: homepage.tudelft.nl/64a8q/intpm/ Book: Integer
More informationRO: Exercices Mixed Integer Programming
RO: Exercices Mixed Integer Programming N. Brauner Université Grenoble Alpes Exercice 1 : Knapsack A hiker wants to fill up his knapsack of capacity W = 6 maximizing the utility of the objects he takes.
More informationOn the Tightness of an LP Relaxation for Rational Optimization and its Applications
OPERATIONS RESEARCH Vol. 00, No. 0, Xxxxx 0000, pp. 000 000 issn 0030-364X eissn 526-5463 00 0000 000 INFORMS doi 0.287/xxxx.0000.0000 c 0000 INFORMS Authors are encouraged to submit new papers to INFORMS
More informationLagrangian Duality Theory
Lagrangian Duality Theory Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye Chapter 14.1-4 1 Recall Primal and Dual
More informationCSCI 1951-G Optimization Methods in Finance Part 01: Linear Programming
CSCI 1951-G Optimization Methods in Finance Part 01: Linear Programming January 26, 2018 1 / 38 Liability/asset cash-flow matching problem Recall the formulation of the problem: max w c 1 + p 1 e 1 = 150
More informationStabilized Branch-and-cut-and-price for the Generalized Assignment Problem
Stabilized Branch-and-cut-and-price for the Generalized Assignment Problem Alexandre Pigatti, Marcus Poggi de Aragão Departamento de Informática, PUC do Rio de Janeiro {apigatti, poggi}@inf.puc-rio.br
More informationQUALIFYING EXAM IN SYSTEMS ENGINEERING
QUALIFYING EXAM IN SYSTEMS ENGINEERING Written Exam: MAY 23, 2017, 9:00AM to 1:00PM, EMB 105 Oral Exam: May 25 or 26, 2017 Time/Location TBA (~1 hour per student) CLOSED BOOK, NO CHEAT SHEETS BASIC SCIENTIFIC
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 informationABSTRACT INTRODUCTION
SGS MINERALS SERVICES TECHNICAL BULLETIN 2010-1 JANUARY 2010 VISION FOR A RISK ADVERSE INTEGRATED GEOMETALLURGY FRAMEWORK GUILLERMO TURNER-SAAD, GLOBAL VICE PRESIDENT METALLURGY AND MINERALOGY, SGS MINERAL
More informationConvex optimization problems. Optimization problem in standard form
Convex optimization problems optimization problem in standard form convex optimization problems linear optimization quadratic optimization geometric programming quasiconvex optimization generalized inequality
More informationFinancial Optimization ISE 347/447. Lecture 21. Dr. Ted Ralphs
Financial Optimization ISE 347/447 Lecture 21 Dr. Ted Ralphs ISE 347/447 Lecture 21 1 Reading for This Lecture C&T Chapter 16 ISE 347/447 Lecture 21 2 Formalizing: Random Linear Optimization Consider the
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 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 informationMulti-Range Robust Optimization vs Stochastic Programming in Prioritizing Project Selection
Multi-Range Robust Optimization vs Stochastic Programming in Prioritizing Project Selection Ruken Düzgün Aurélie Thiele July 2012 Abstract This paper describes a multi-range robust optimization approach
More informationInteger Linear Programming Modeling
DM554/DM545 Linear and Lecture 9 Integer Linear Programming Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. Assignment Problem Knapsack Problem
More informationInteger Program Reformulation for Robust. Branch-and-Cut-and-Price Algorithms
Integer Program Reformulation for Robust Branch-and-Cut-and-Price Algorithms Marcus Poggi de Aragão 1, Eduardo Uchoa 2 1 Departamento de Informática, PUC-Rio, poggi@inf.puc-rio.br 2 Dep. de Engenharia
More informationSequential pairing of mixed integer inequalities
Sequential pairing of mixed integer inequalities Yongpei Guan, Shabbir Ahmed, George L. Nemhauser School of Industrial & Systems Engineering, Georgia Institute of Technology, 765 Ferst Drive, Atlanta,
More informationMEDIUM-TERM HYDROPOWER SCHEDULING BY STOCHASTIC DUAL DYNAMIC INTEGER PROGRAMMING: A CASE STUDY
MEDIUM-TERM HYDROPOWER SCHEDULING BY STOCHASTIC DUAL DYNAMIC INTEGER PROGRAMMING: A CASE STUDY Martin Hjelmeland NTNU Norwegian University of Science and Technology, Trondheim, Norway. martin.hjelmeland@ntnu.no
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 informationImproving Performance of The Interior Point Method by Preconditioning
Improving Performance of The Interior Point Method by Preconditioning Mid-Point Status Report Project by: Ken Ryals For: AMSC 663-664 Fall 27-Spring 28 6 December 27 Background / Refresher The IPM method
More informationStructured Problems and Algorithms
Integer and quadratic optimization problems Dept. of Engg. and Comp. Sci., Univ. of Cal., Davis Aug. 13, 2010 Table of contents Outline 1 2 3 Benefits of Structured Problems Optimization problems may become
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 informationReplacement: Material Increase to Tuzon Resource
Replacement: Material Increase to Tuzon Resource RNS Number : 1067C Hummingbird Resources PLC 12 March 2014 Hummingbird Resources plc ("Hummingbird Resources" or the "Company") (AIM: HUM) THE FOLLOWING
More informationOptimisation. 3/10/2010 Tibor Illés Optimisation
Optimisation Lectures 3 & 4: Linear Programming Problem Formulation Different forms of problems, elements of the simplex algorithm and sensitivity analysis Lecturer: Tibor Illés tibor.illes@strath.ac.uk
More informationDisjunctive Decomposition for Two-Stage Stochastic Mixed-Binary Programs with GUB Constraints
Disjunctive Decomposition for Two-Stage Stochastic Mixed-Binary Programs with GUB Constraints Brian Keller Booz Allen Hamilton, 134 National Business Parkway, Annapolis Junction, MD 20701, USA, keller
More information