Design of Global Supply Chains

Size: px
Start display at page:

Download "Design of Global Supply Chains"

Transcription

1 Design of Global Supply Chains Marc Goetschalckx Carlos Vidal Tjendera Santoso IERC Cleveland, Spring 2000 Marc Goetschalckx Overview Global Supply Chain Problem Bilinear Transfer Price Formulation Iterative Heuristic Global Optimization Procedure Computational Example Conclusions Marc Goetschalckx 1

2 Global Logistics Systems Models Domestic Plus Exchange Rates, Duties, Taxes Objective is Worldwide After-Tax Profit Maximization Decisions are Material Flows, Transportation Cost Allocations, and Transfer Prices Tax Rates and Profit Realization 34% 17% 12% 40% or more Marc Goetschalckx 2

3 Previous Research Nieckels (1976) NLP to solve TPP iteratively, local opt. Single commodity, no BOM Cohen et al. (1989) Dyn. NL MIP, solved iteratively, local opt. TP are markup Strong country tax reduction (feasible?) Implementation Cohen and Lee (1989) Previous Research Continued Arntzen et al. (1995) Comprehensive model for DEC No TP and taxes part of production costs Specialized MIP algorithms Canel and Khumawala (1997) TP fixed a priori at LB or UB MIP model Marc Goetschalckx 3

4 Previous Research Summary TP either Set a priori Determined iteratively, local optimum, and no bound Not included Global Transactions and Tax Rates a x,t,p b Country A Country B Tax Rate Net Income Before Tax Marc Goetschalckx 4

5 Impact of Transfer Prices on After Tax Profits A and B profit A loses, B profits A profits, B los es Transfe r price t Net income after tax A and B los e x = 20,000 x = 15,000 x = 12,000 x = 8,000 x = 6, The basic model Maximize global after tax profit: Subject to: After tax profit of internal suppliers + After tax profit at plants + After tax profit at distribution centers - Inventory costs Nonlinear expressions for the net income before tax Suppliers capacity (internal and external suppliers) Production capacity at plants Customer demand constraints Bill of materials (at plants) Balance constraints at DCs Minimum profit for internal suppliers, plants and DCs (optional) Bounds on decision variables Marc Goetschalckx 5

6 Before and After Tax Profit Objective and Constraint + max :(1 taxratek) ibtwfk ibtwfk subject to : 1 1 MPRICE w HANDC + TRCWM W w VPkp H TTWM klm + ( CSF ) SHIPFREQklm + SSFWkp l C( k) m T( k, l) p P E k TTWM klm wklmp lp klmp kp klm p klmp l C( k) m T( k, l) p P El l C( k) m T( k, l) p P Ek 1 tppldc ( 1 + DUTY ) + ( 1 propw ) TRCPW W x j M m T( j, k) p P( j) E j 1 + f FIXDCk = ibtwf k ibtwfk k W Ek jp jkp jkm jkm p jkmp The Basic Model (General Structure) Max s. to: T c x T T T P( x, t, p) c x + x A t + x B p = b ; r = 12,,..., m r Cx d 0 T t T 0 p 1 x 0, t 0 r + r r Marc Goetschalckx 6

7 Solution Methodology An optimization-based heuristic: Substitution of proportion variables Redefinition and substitution of TP variables Relaxation of nonlinear constraints Iterative procedure Global optimization Tightening of Dual Bound with Primal Heuristics Substitution of Transportation Cost Proportion Variables prosp W s = z i S, j M( i), m T( i, j) ijm r ijmr r R( i) R( j) ijm propw W x = z j M, k W, m T( j, k) jkm p jkmp p P( j) jkm Marc Goetschalckx 7

8 Transfer Prices Substitution and Constraints tpsupl s = y i S, j M( i), r R( i) R( j) ijr ijmr ijr m T(, i j) tppldc x = y j M, k W, p P( j) jkp jkmp jkp m T( j, k) y ij r ij r n = n+ 1 n+ 1 ijnmr ijn+ 1mr m T(, i j ) m T(, i j ) y n jk p jk p n = n+ 1 n+ 1 jknmp jkn+ 1mp m T( j, k ) m T( j, k ) n s x y y s x Transformed Formulation Max s. to: d v T 0 T T T T c x+ d v+ e y+ g z = f ; r = 1,2,..., m r r r r r Cx b l u Dx y Dx z Ex 0 T x F y = 0; q= 1,2,..., h (constraints to be relaxed) q x 0, y 0, v 0, z 0 Marc Goetschalckx 8

9 Optimization-Based Heuristic Solve RP(x, y, z): Upper Bound on P(x, t, p) and starting point Feasible? No Yes Optimal Solution to P(x, t, p) Solve P(x, t, z x) Solve P(x, t, z t) No Convergence criterion? Yes Local optimum of P(x, t, p) Starting Points for Heuristic Procedure Optimal TP (From Relaxation) Optimal Flows (From Relaxation) Tax Heuristic Lower Bound TP Upper Bound TP Middle Point TP Interval Zero Initial Flows Marc Goetschalckx 9

10 Computational Test Case 50 Raw Material Suppliers 8 Plants, 10 Distribution Centers 80 Customers 35 Components, 12 Finished Products 3.1 Modes per Channel Variables, 2900 Constraints Heuristic Computation Results MEDIUM INSTANCES a No. Starting Point % gap from the upper bound Solution time (s) 3 Opt_Flows Opt_TP Heu_TP LB_TP UB_TP Mid_TP Init_Flows = No. of iterations b 4 Opt_Flows Opt_TP Heu_TP LB_TP UB_TP Mid_TP Init_Flows = Opt_Flows , Opt_TP , Heu_TP LB_TP , UB_TP , ,356 Mid_TP , ,177 Init_Flows = , ,028 Marc Goetschalckx 10

11 Global Optimization Procedure Specified Optimality Gap ε Ben-Tal (1994) Branch and Bound Method for Reducing Duality Gap f max g i g T i Acceleration Techniques Branching rule: largest TP interval Branching rule: interior transfer prices Computational Experiment Global Optimization Procedure No. Target Tolerance ε (%) MEDIUM INSTANCES % gap achieved Solution time (s) No. of branches b 3 N/A a (0.405) d c N/A a (0.568) d e e 5 N/A a (0.28) d e e Marc Goetschalckx 11

12 Profit Increases for Optimal Transfer Prices Transfer Price Heuristics Instance Middle Tax Lower Upper Point Rate Bound Bound Impact of Transfer Price on Corporate Profits NET INCOME AFTER TAX [$/YEAR] vs. TRANSFER PRICE FACTOR 24,000,000 23,500,000 23,000,000 22,500,000 Net Income After Tax [$/year] 22,000,000 21,500,000 21,000,000 20,500,000 20,000,000 What is the maximum transfer price factor that satisfies the operating profit constraints of all countries? 19,500,000 NIAT 19,000,000 18,500, Transfer Price Factor Marc Goetschalckx 12

13 Conclusions Transfer Price Formulation is Bilinear Iterative Heuristic Efficient, Case Dependent Gap Global Optimization Procedure A Priori Gap Efficient with Acceleration Techniques Significant Impact on After Tax Profits Supply Chain Modeling Challenges Multiple Periods Periodic demand Dynamic strategic systems Global Taxes and profit realization Local contents, duty drawback Stochastic Flexibility, robustness, risk, scenarios Marc Goetschalckx 13

14 Supply Chain Solution Algorithms Challenges Large Scale Models Non-Linear Models Stochastic Models Standard MIP Linear Algorithms Cannot Solve Very Large Cases NL-MIP or Stochastic Algorithms Only for Small Cases or Nonexistent Supply Chain Design Challenges Integrated models are large and complex Accommodate diversity of local characteristics Cost, flexibility, and responsiveness tradeoffs for performance measures Strategic design as a continuous effort Technology transfer to logistics professionals and students Marc Goetschalckx 14

15 From a Multicommodity Case......and Configuration by a Current Design Tool Marc Goetschalckx 15

16 To Design Tools for the Next Century Marc Goetschalckx 16

Capacity Planning with uncertainty in Industrial Gas Markets

Capacity Planning with uncertainty in Industrial Gas Markets Capacity Planning with uncertainty in Industrial Gas Markets A. Kandiraju, P. Garcia Herreros, E. Arslan, P. Misra, S. Mehta & I.E. Grossmann EWO meeting September, 2015 1 Motivation Industrial gas markets

More information

Overview of course. Introduction to Optimization, DIKU Monday 12 November David Pisinger

Overview of course. Introduction to Optimization, DIKU Monday 12 November David Pisinger Introduction to Optimization, DIKU 007-08 Monday November David Pisinger Lecture What is OR, linear models, standard form, slack form, simplex repetition, graphical interpretation, extreme points, basic

More information

- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs

- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs

More information

Designing the Distribution Network for an Integrated Supply Chain

Designing the Distribution Network for an Integrated Supply Chain Designing the Distribution Network for an Integrated Supply Chain Jia Shu and Jie Sun Abstract We consider an integrated distribution network design problem in which all the retailers face uncertain demand.

More information

Sensitivity Analysis and Duality in LP

Sensitivity Analysis and Duality in LP Sensitivity Analysis and Duality in LP Xiaoxi Li EMS & IAS, Wuhan University Oct. 13th, 2016 (week vi) Operations Research (Li, X.) Sensitivity Analysis and Duality in LP Oct. 13th, 2016 (week vi) 1 /

More information

Improvements to Benders' decomposition: systematic classification and performance comparison in a Transmission Expansion Planning problem

Improvements to Benders' decomposition: systematic classification and performance comparison in a Transmission Expansion Planning problem Improvements to Benders' decomposition: systematic classification and performance comparison in a Transmission Expansion Planning problem Sara Lumbreras & Andrés Ramos July 2013 Agenda Motivation improvement

More information

Handout 1: Introduction to Dynamic Programming. 1 Dynamic Programming: Introduction and Examples

Handout 1: Introduction to Dynamic Programming. 1 Dynamic Programming: Introduction and Examples SEEM 3470: Dynamic Optimization and Applications 2013 14 Second Term Handout 1: Introduction to Dynamic Programming Instructor: Shiqian Ma January 6, 2014 Suggested Reading: Sections 1.1 1.5 of Chapter

More information

PART 4 INTEGER PROGRAMMING

PART 4 INTEGER PROGRAMMING PART 4 INTEGER PROGRAMMING 102 Read Chapters 11 and 12 in textbook 103 A capital budgeting problem We want to invest $19 000 Four investment opportunities which cannot be split (take it or leave it) 1.

More information

DESIGN OF OPTIMAL LINEAR SYSTEMS BY MULTIPLE OBJECTIVES

DESIGN OF OPTIMAL LINEAR SYSTEMS BY MULTIPLE OBJECTIVES etr Fiala ESIGN OF OTIMAL LINEAR SYSTEMS Y MULTILE OJECTIVES Abstract Traditional concepts of optimality focus on valuation of already given systems. A new concept of designing optimal systems is proposed.

More information

CHAPTER-3 MULTI-OBJECTIVE SUPPLY CHAIN NETWORK PROBLEM

CHAPTER-3 MULTI-OBJECTIVE SUPPLY CHAIN NETWORK PROBLEM CHAPTER-3 MULTI-OBJECTIVE SUPPLY CHAIN NETWORK PROBLEM 3.1 Introduction A supply chain consists of parties involved, directly or indirectly, in fulfilling customer s request. The supply chain includes

More information

The use of shadow price is an example of sensitivity analysis. Duality theory can be applied to do other kind of sensitivity analysis:

The use of shadow price is an example of sensitivity analysis. Duality theory can be applied to do other kind of sensitivity analysis: Sensitivity analysis The use of shadow price is an example of sensitivity analysis. Duality theory can be applied to do other kind of sensitivity analysis: Changing the coefficient of a nonbasic variable

More information

MEDIUM-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 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 information

System Planning Lecture 7, F7: Optimization

System 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 information

u = 50 u = 30 Generalized Maximum Flow Problem s u = 00 2 u = 00 u = 50 4 = 3=4 u = 00 t capacity u = 20 3 u = 20 = =2 5 u = 00 gain/loss factor Max o

u = 50 u = 30 Generalized Maximum Flow Problem s u = 00 2 u = 00 u = 50 4 = 3=4 u = 00 t capacity u = 20 3 u = 20 = =2 5 u = 00 gain/loss factor Max o Generalized Max Flows Kevin Wayne Cornell University www.orie.cornell.edu/~wayne = 3=4 40 v w 30 advisor: Eva Tardos u = 50 u = 30 Generalized Maximum Flow Problem s u = 00 2 u = 00 u = 50 4 = 3=4 u =

More information

Optimisation. 3/10/2010 Tibor Illés Optimisation

Optimisation. 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 information

Multiperiod Blend Scheduling Problem

Multiperiod Blend Scheduling Problem ExxonMobil Multiperiod Blend Scheduling Problem Juan Pablo Ruiz Ignacio E. Grossmann Department of Chemical Engineering Center for Advanced Process Decision-making University Pittsburgh, PA 15213 1 Motivation

More information

Duality in LPP Every LPP called the primal is associated with another LPP called dual. Either of the problems is primal with the other one as dual. The optimal solution of either problem reveals the information

More information

Column Generation. MTech Seminar Report. Soumitra Pal Roll No: under the guidance of

Column 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 information

Facility Location and Distribution System Planning. Thomas L. Magnanti

Facility Location and Distribution System Planning. Thomas L. Magnanti Facility Location and Distribution System Planning Thomas L. Magnanti Today s Agenda Why study facility location? Issues to be modeled Basic models Fixed charge problems Core uncapacitated and capacitated

More information

Nonconvex Generalized Benders Decomposition Paul I. Barton

Nonconvex Generalized Benders Decomposition Paul I. Barton Nonconvex Generalized Benders Decomposition Paul I. Barton Process Systems Engineering Laboratory Massachusetts Institute of Technology Motivation Two-stage Stochastic MINLPs & MILPs Natural Gas Production

More information

Outline. Outline. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Scheduling CPM/PERT Resource Constrained Project Scheduling Model

Outline. 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 information

MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands

MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands Axel Nyberg Åbo Aademi University Ignacio E. Grossmann Dept. of Chemical Engineering, Carnegie Mellon University,

More information

Advances in Convex Optimization: Theory, Algorithms, and Applications

Advances in Convex Optimization: Theory, Algorithms, and Applications Advances in Convex Optimization: Theory, Algorithms, and Applications Stephen Boyd Electrical Engineering Department Stanford University (joint work with Lieven Vandenberghe, UCLA) ISIT 02 ISIT 02 Lausanne

More information

Resource Constrained Project Scheduling Linear and Integer Programming (1)

Resource 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 information

MULTI-PERIOD MULTI-DIMENSIONAL KNAPSACK PROBLEM AND ITS APPLICATION TO AVAILABLE-TO-PROMISE

MULTI-PERIOD MULTI-DIMENSIONAL KNAPSACK PROBLEM AND ITS APPLICATION TO AVAILABLE-TO-PROMISE MULTI-PERIOD MULTI-DIMENSIONAL KNAPSACK PROBLEM AND ITS APPLICATION TO AVAILABLE-TO-PROMISE Hoong Chuin LAU and Min Kwang LIM School of Computing National University of Singapore, 3 Science Drive 2, Singapore

More information

A three-level MILP model for generation and transmission expansion planning

A three-level MILP model for generation and transmission expansion planning A three-level MILP model for generation and transmission expansion planning David Pozo Cámara (UCLM) Enzo E. Sauma Santís (PUC) Javier Contreras Sanz (UCLM) Contents 1. Introduction 2. Aims and contributions

More information

Introduction to Bin Packing Problems

Introduction 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 information

Notes on Dantzig-Wolfe decomposition and column generation

Notes 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 information

Agenda today. Introduction to prescriptive modeling. Linear optimization models through three examples: Beyond linear optimization

Agenda today. Introduction to prescriptive modeling. Linear optimization models through three examples: Beyond linear optimization Agenda today Introduction to prescriptive modeling Linear optimization models through three examples: 1 Production and inventory optimization 2 Distribution system design 3 Stochastic optimization Beyond

More information

Linear Programming: Chapter 5 Duality

Linear Programming: Chapter 5 Duality Linear Programming: Chapter 5 Duality Robert J. Vanderbei September 30, 2010 Slides last edited on October 5, 2010 Operations Research and Financial Engineering Princeton University Princeton, NJ 08544

More information

A Lagrangian relaxation method for solving choice-based mixed linear optimization models that integrate supply and demand interactions

A Lagrangian relaxation method for solving choice-based mixed linear optimization models that integrate supply and demand interactions A Lagrangian relaxation method for solving choice-based mixed linear optimization models that integrate supply and demand interactions Meritxell Pacheco Shadi Sharif Azadeh, Michel Bierlaire, Bernard Gendron

More information

Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS. Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano

Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS. Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano ... Our contribution PIPS-PSBB*: Multi-level parallelism for Stochastic

More information

Pareto Efficiency in Robust Optimization

Pareto Efficiency in Robust Optimization Pareto Efficiency in Robust Optimization Dan Iancu Graduate School of Business Stanford University joint work with Nikolaos Trichakis (HBS) 1/26 Classical Robust Optimization Typical linear optimization

More information

Fundamentals of Operations Research. Prof. G. Srinivasan. Indian Institute of Technology Madras. Lecture No. # 15

Fundamentals of Operations Research. Prof. G. Srinivasan. Indian Institute of Technology Madras. Lecture No. # 15 Fundamentals of Operations Research Prof. G. Srinivasan Indian Institute of Technology Madras Lecture No. # 15 Transportation Problem - Other Issues Assignment Problem - Introduction In the last lecture

More information

Department of Agricultural Economics. PhD Qualifier Examination. May 2009

Department of Agricultural Economics. PhD Qualifier Examination. May 2009 Department of Agricultural Economics PhD Qualifier Examination May 009 Instructions: The exam consists of six questions. You must answer all questions. If you need an assumption to complete a question,

More information

Linear Programming. H. R. Alvarez A., Ph. D. 1

Linear Programming. H. R. Alvarez A., Ph. D. 1 Linear Programming H. R. Alvarez A., Ph. D. 1 Introduction It is a mathematical technique that allows the selection of the best course of action defining a program of feasible actions. The objective of

More information

Integer programming: an introduction. Alessandro Astolfi

Integer programming: an introduction. Alessandro Astolfi Integer programming: an introduction Alessandro Astolfi Outline Introduction Examples Methods for solving ILP Optimization on graphs LP problems with integer solutions Summary Introduction Integer programming

More information

Introduction to LP. Types of Linear Programming. There are five common types of decisions in which LP may play a role

Introduction to LP. Types of Linear Programming. There are five common types of decisions in which LP may play a role Linear Programming RK Jana Lecture Outline Introduction to Linear Programming (LP) Historical Perspective Model Formulation Graphical Solution Method Simplex Method Introduction to LP Continued Today many

More information

Some Results in Duality

Some Results in Duality Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 30, 1493-1501 Some Results in Duality Vanita Ben Dhagat and Savita Tiwari Jai Narain college of Technology Bairasia Road, Bhopal M.P., India vanita1_dhagat@yahoo.co.in

More information

Linear Programming Applications. Transportation Problem

Linear Programming Applications. Transportation Problem Linear Programming Applications Transportation Problem 1 Introduction Transportation problem is a special problem of its own structure. Planning model that allocates resources, machines, materials, capital

More information

3.10 Lagrangian relaxation

3.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 information

Bilinear Programming: Applications in the Supply Chain Management

Bilinear Programming: Applications in the Supply Chain Management Bilinear Programming: Applications in the Supply Chain Management Artyom G. Nahapetyan Center for Applied Optimization Industrial and Systems Engineering Department University of Florida Gainesville, Florida

More information

The multi-period incremental service facility location problem

The multi-period incremental service facility location problem Computers & Operations Research ( ) www.elsevier.com/locate/cor The multi-period incremental service facility location problem Maria Albareda-Sambola a,, Elena Fernández a, Yolanda Hinojosa b, Justo Puerto

More information

Section 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 Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections

More information

Lecture 11 October 7, 2013

Lecture 11 October 7, 2013 CS 4: Advanced Algorithms Fall 03 Prof. Jelani Nelson Lecture October 7, 03 Scribe: David Ding Overview In the last lecture we talked about set cover: Sets S,..., S m {,..., n}. S has cost c S. Goal: Cover

More information

Feasibility Pump Heuristics for Column Generation Approaches

Feasibility 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 information

MAT016: Optimization

MAT016: Optimization MAT016: Optimization M.El Ghami e-mail: melghami@ii.uib.no URL: http://www.ii.uib.no/ melghami/ March 29, 2011 Outline for today The Simplex method in matrix notation Managing a production facility The

More information

i.e., into a monomial, using the Arithmetic-Geometric Mean Inequality, the result will be a posynomial approximation!

i.e., into a monomial, using the Arithmetic-Geometric Mean Inequality, the result will be a posynomial approximation! Dennis L. Bricker Dept of Mechanical & Industrial Engineering The University of Iowa i.e., 1 1 1 Minimize X X X subject to XX 4 X 1 0.5X 1 Minimize X X X X 1X X s.t. 4 1 1 1 1 4X X 1 1 1 1 0.5X X X 1 1

More information

Branch-and-Price-and-Cut for the Split Delivery Vehicle Routing Problem with Time Windows

Branch-and-Price-and-Cut for the Split Delivery Vehicle Routing Problem with Time Windows Branch-and-Price-and-Cut for the Split Delivery Vehicle Routing Problem with Time Windows Guy Desaulniers École Polytechnique de Montréal and GERAD Column Generation 2008 Aussois, France Outline Introduction

More information

3.10 Column generation method

3.10 Column generation method 3.10 Column generation method Many relevant decision-making (discrete optimization) problems can be formulated as ILP problems with a very large (exponential) number of variables. Examples: cutting stock,

More information

Public Economics Ben Heijdra Chapter 9: Introduction to Normative Public Economics

Public Economics Ben Heijdra Chapter 9: Introduction to Normative Public Economics Public Economics: Chapter 9 1 Public Economics Ben Heijdra Chapter 9: Introduction to Normative Public Economics Objectives of this chapter Public Economics: Chapter 9 2 Read Atkinson & Stiglitz (1980,

More information

Robust technological and emission trajectories for long-term stabilization targets with an energyenvironment

Robust technological and emission trajectories for long-term stabilization targets with an energyenvironment Robust technological and emission trajectories for long-term stabilization targets with an energyenvironment model Nicolas Claire Tchung-Ming Stéphane (IFPEN) Bahn Olivier (HEC Montréal) Delage Erick (HEC

More information

3.4 Relaxations and bounds

3.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 information

Deterministic Operations Research, ME 366Q and ORI 391 Chapter 2: Homework #2 Solutions

Deterministic Operations Research, ME 366Q and ORI 391 Chapter 2: Homework #2 Solutions Deterministic Operations Research, ME 366Q and ORI 391 Chapter 2: Homework #2 Solutions 11. Consider the following linear program. Maximize z = 6x 1 + 3x 2 subject to x 1 + 2x 2 2x 1 + x 2 20 x 1 x 2 x

More information

Solving LP and MIP Models with Piecewise Linear Objective Functions

Solving 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 information

Reducing the Cost of Demand Uncertainty through Accurate Response to Early Sales

Reducing the Cost of Demand Uncertainty through Accurate Response to Early Sales Reducing the Cost of Demand Uncertainty through Accurate Response to Early Sales Fisher and Raman 1996 Operations Research,, vol. 44 (1), 87-99 Presentation by Hongmin Li This summary presentation is based

More information

A Hub Location Problem with Fully Interconnected Backbone and Access Networks

A Hub Location Problem with Fully Interconnected Backbone and Access Networks A Hub Location Problem with Fully Interconnected Backbone and Access Networks Tommy Thomadsen Informatics and Mathematical Modelling Technical University of Denmark 2800 Kgs. Lyngby Denmark tt@imm.dtu.dk

More information

Benders Decomposition Methods for Structured Optimization, including Stochastic Optimization

Benders Decomposition Methods for Structured Optimization, including Stochastic Optimization Benders Decomposition Methods for Structured Optimization, including Stochastic Optimization Robert M. Freund May 2, 2001 Block Ladder Structure Basic Model minimize x;y c T x + f T y s:t: Ax = b Bx +

More information

A Priori Route Evaluation for the Lateral Transhipment Problem (ARELTP) with Piecewise Linear Profits

A Priori Route Evaluation for the Lateral Transhipment Problem (ARELTP) with Piecewise Linear Profits 1 / 47 A Priori Route Evaluation for the Lateral Transhipment Problem (ARELTP) with Piecewise Linear Profits Martin Romauch 1 Richard Hartl 1 Thibaut Vidal 2 1 University of Vienna 2 PUC-Rio, Rio de Janeiro,

More information

Routing. Topics: 6.976/ESD.937 1

Routing. Topics: 6.976/ESD.937 1 Routing Topics: Definition Architecture for routing data plane algorithm Current routing algorithm control plane algorithm Optimal routing algorithm known algorithms and implementation issues new solution

More information

Duality Methods in Portfolio Allocation with Transaction Constraints and Uncertainty

Duality Methods in Portfolio Allocation with Transaction Constraints and Uncertainty Duality Methods in Portfolio Allocation with Transaction Constraints and Uncertainty Christopher W. Miller Department of Mathematics University of California, Berkeley January 9, 2014 Project Overview

More information

Recoverable Robustness in Scheduling Problems

Recoverable 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 information

Robust Inventory Routing under Demand Uncertainty

Robust Inventory Routing under Demand Uncertainty TRANSPORTATION SCIENCE Vol. 00, No. 0, Xxxxx 0000, pp. 000 000 issn0041-1655 eissn1526-5447 00 0000 0001 INFORMS doi 10.1287/xxxx.0000.0000 c 0000 INFORMS Robust Inventory Routing under Demand Uncertainty

More information

Benders Decomposition

Benders 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 information

Module 8 Linear Programming. CS 886 Sequential Decision Making and Reinforcement Learning University of Waterloo

Module 8 Linear Programming. CS 886 Sequential Decision Making and Reinforcement Learning University of Waterloo Module 8 Linear Programming CS 886 Sequential Decision Making and Reinforcement Learning University of Waterloo Policy Optimization Value and policy iteration Iterative algorithms that implicitly solve

More information

Integrated schedule planning with supply-demand interactions for a new generation of aircrafts

Integrated schedule planning with supply-demand interactions for a new generation of aircrafts Integrated schedule planning with supply-demand interactions for a new generation of aircrafts Bilge Atasoy, Matteo Salani and Michel Bierlaire Abstract We present an integrated schedule planning model

More information

56:270 Final Exam - May

56:270  Final Exam - May @ @ 56:270 Linear Programming @ @ Final Exam - May 4, 1989 @ @ @ @ @ @ @ @ @ @ @ @ @ @ Select any 7 of the 9 problems below: (1.) ANALYSIS OF MPSX OUTPUT: Please refer to the attached materials on the

More information

A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem

A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem Dimitris J. Bertsimas Dan A. Iancu Pablo A. Parrilo Sloan School of Management and Operations Research Center,

More information

OPTIMIZATION. joint course with. Ottimizzazione Discreta and Complementi di R.O. Edoardo Amaldi. DEIB Politecnico di Milano

OPTIMIZATION. joint course with. Ottimizzazione Discreta and Complementi di R.O. Edoardo Amaldi. DEIB Politecnico di Milano OPTIMIZATION joint course with Ottimizzazione Discreta and Complementi di R.O. Edoardo Amaldi DEIB Politecnico di Milano edoardo.amaldi@polimi.it Website: http://home.deib.polimi.it/amaldi/opt-15-16.shtml

More information

BBM402-Lecture 20: LP Duality

BBM402-Lecture 20: LP Duality BBM402-Lecture 20: LP Duality Lecturer: Lale Özkahya Resources for the presentation: https://courses.engr.illinois.edu/cs473/fa2016/lectures.html An easy LP? which is compact form for max cx subject to

More information

Mixed Integer Linear Programming Formulation for Chance Constrained Mathematical Programs with Equilibrium Constraints

Mixed Integer Linear Programming Formulation for Chance Constrained Mathematical Programs with Equilibrium Constraints Mixed Integer Linear Programming Formulation for Chance Constrained Mathematical Programs with Equilibrium Constraints ayed A. adat and Lingling Fan University of outh Florida, email: linglingfan@usf.edu

More information

Dynamic supply chain design with inventory

Dynamic supply chain design with inventory Computers & Operations Research ( ) www.elsevier.com/locate/cor Dynamic supply chain design with inventory Y. Hinoosa a,, J. Kalcsics b, S. Nickel c, J. Puerto d, S. Velten e a Departamento de Economía

More information

Part III: Traveling salesman problems

Part III: Traveling salesman problems Transportation Logistics Part III: Traveling salesman problems c R.F. Hartl, S.N. Parragh 1/282 Motivation Motivation Why do we study the TSP? c R.F. Hartl, S.N. Parragh 2/282 Motivation Motivation Why

More information

Algorithms and Theory of Computation. Lecture 13: Linear Programming (2)

Algorithms and Theory of Computation. Lecture 13: Linear Programming (2) Algorithms and Theory of Computation Lecture 13: Linear Programming (2) Xiaohui Bei MAS 714 September 25, 2018 Nanyang Technological University MAS 714 September 25, 2018 1 / 15 LP Duality Primal problem

More information

Exercises - Linear Programming

Exercises - Linear Programming Chapter 38 Exercises - Linear Programming By Sariel Har-Peled, December 10, 2007 1 Version: 1.0 This chapter include problems that are related to linear programming. 38.1 Miscellaneous Exercise 38.1.1

More information

Integer Programming ISE 418. Lecture 8. Dr. Ted Ralphs

Integer 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 information

Stochastic Optimization

Stochastic Optimization Chapter 27 Page 1 Stochastic Optimization Operations research has been particularly successful in two areas of decision analysis: (i) optimization of problems involving many variables when the outcome

More information

SINGLE OBJECTIVE RISK- BASED TRANSMISSION EXPANSION

SINGLE OBJECTIVE RISK- BASED TRANSMISSION EXPANSION Vol.2, Issue.1, Jan-Feb 2012 pp-424-430 ISSN: 2249-6645 SINGLE OBJECTIVE RISK- BASED TRANSMISSION EXPANSION V.Sumadeepthi 1, K.Sarada 2 1 (Student, Department of Electrical and Electronics Engineering,

More information

maxz = 3x 1 +4x 2 2x 1 +x 2 6 2x 1 +3x 2 9 x 1,x 2

maxz = 3x 1 +4x 2 2x 1 +x 2 6 2x 1 +3x 2 9 x 1,x 2 ex-5.-5. Foundations of Operations Research Prof. E. Amaldi 5. Branch-and-Bound Given the integer linear program maxz = x +x x +x 6 x +x 9 x,x integer solve it via the Branch-and-Bound method (solving

More information

3.10 Column generation method

3.10 Column generation method 3.10 Column generation method Many relevant decision-making problems can be formulated as ILP problems with a very large (exponential) number of variables. Examples: cutting stock, crew scheduling, vehicle

More information

Duality Theory, Optimality Conditions

Duality Theory, Optimality Conditions 5.1 Duality Theory, Optimality Conditions Katta G. Murty, IOE 510, LP, U. Of Michigan, Ann Arbor We only consider single objective LPs here. Concept of duality not defined for multiobjective LPs. Every

More information

Benders Decomposition for the Uncapacitated Multicommodity Network Design Problem

Benders Decomposition for the Uncapacitated Multicommodity Network Design Problem Benders Decomposition for the Uncapacitated Multicommodity Network Design Problem 1 Carlos Armando Zetina, 1 Ivan Contreras, 2 Jean-François Cordeau 1 Concordia University and CIRRELT, Montréal, Canada

More information

Lecture 1. Stochastic Optimization: Introduction. January 8, 2018

Lecture 1. Stochastic Optimization: Introduction. January 8, 2018 Lecture 1 Stochastic Optimization: Introduction January 8, 2018 Optimization Concerned with mininmization/maximization of mathematical functions Often subject to constraints Euler (1707-1783): Nothing

More information

Large-scale optimization and decomposition methods: outline. Column Generation and Cutting Plane methods: a unified view

Large-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 information

Lecture - 02 Rules for Pinch Design Method (PEM) - Part 02

Lecture - 02 Rules for Pinch Design Method (PEM) - Part 02 Process Integration Prof. Bikash Mohanty Department of Chemical Engineering Indian Institute of Technology, Roorkee Module - 05 Pinch Design Method for HEN synthesis Lecture - 02 Rules for Pinch Design

More information

Pedro Munari - COA 2017, February 10th, University of Edinburgh, Scotland, UK 2

Pedro Munari - COA 2017, February 10th, University of Edinburgh, Scotland, UK 2 Pedro Munari [munari@dep.ufscar.br] - COA 2017, February 10th, University of Edinburgh, Scotland, UK 2 Outline Vehicle routing problem; How interior point methods can help; Interior point branch-price-and-cut:

More information

Lecture 9: Dantzig-Wolfe Decomposition

Lecture 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 information

The CPLEX Library: Mixed Integer Programming

The CPLEX Library: Mixed Integer Programming The CPLEX Library: Mixed Programming Ed Rothberg, ILOG, Inc. 1 The Diet Problem Revisited Nutritional values Bob considered the following foods: Food Serving Size Energy (kcal) Protein (g) Calcium (mg)

More information

The Transportation Problem

The Transportation Problem CHAPTER 12 The Transportation Problem Basic Concepts 1. Transportation Problem: BASIC CONCEPTS AND FORMULA This type of problem deals with optimization of transportation cost in a distribution scenario

More information

Lagrangian Duality. Richard Lusby. Department of Management Engineering Technical University of Denmark

Lagrangian Duality. Richard Lusby. Department of Management Engineering Technical University of Denmark Lagrangian Duality Richard Lusby Department of Management Engineering Technical University of Denmark Today s Topics (jg Lagrange Multipliers Lagrangian Relaxation Lagrangian Duality R Lusby (42111) Lagrangian

More information

From structures to heuristics to global solvers

From 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 information

5 Integer Linear Programming (ILP) E. Amaldi Foundations of Operations Research Politecnico di Milano 1

5 Integer Linear Programming (ILP) E. Amaldi Foundations of Operations Research Politecnico di Milano 1 5 Integer Linear Programming (ILP) E. Amaldi Foundations of Operations Research Politecnico di Milano 1 Definition: An Integer Linear Programming problem is an optimization problem of the form (ILP) min

More information

Decision Models Lecture 5 1. Lecture 5. Foreign-Currency Trading Integer Programming Plant-location example Summary and Preparation for next class

Decision Models Lecture 5 1. Lecture 5. Foreign-Currency Trading Integer Programming Plant-location example Summary and Preparation for next class Decision Models Lecture 5 1 Lecture 5 Foreign-Currency Trading Integer Programming Plant-location example Summary and Preparation for next class Foreign Exchange (FX) Markets Decision Models Lecture 5

More information

Almost Robust Optimization with Binary Variables

Almost Robust Optimization with Binary Variables Almost Robust Optimization with Binary Variables Opher Baron, Oded Berman, Mohammad M. Fazel-Zarandi Rotman School of Management, University of Toronto, Toronto, Ontario M5S 3E6, Canada, Opher.Baron@Rotman.utoronto.ca,

More information

Computations with Disjunctive Cuts for Two-Stage Stochastic Mixed 0-1 Integer Programs

Computations with Disjunctive Cuts for Two-Stage Stochastic Mixed 0-1 Integer Programs Computations with Disjunctive Cuts for Two-Stage Stochastic Mixed 0-1 Integer Programs Lewis Ntaimo and Matthew W. Tanner Department of Industrial and Systems Engineering, Texas A&M University, 3131 TAMU,

More information

Network Flows. 6. Lagrangian Relaxation. Programming. Fall 2010 Instructor: Dr. Masoud Yaghini

Network 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 information

Decentralized Control of Stochastic Systems

Decentralized Control of Stochastic Systems Decentralized Control of Stochastic Systems Sanjay Lall Stanford University CDC-ECC Workshop, December 11, 2005 2 S. Lall, Stanford 2005.12.11.02 Decentralized Control G 1 G 2 G 3 G 4 G 5 y 1 u 1 y 2 u

More information

A LAGRANGIAN RELAXATION FOR CAPACITATED SINGLE ALLOCATION P-HUB MEDIAN PROBLEM WITH MULTIPLE CAPACITY LEVELS

A LAGRANGIAN RELAXATION FOR CAPACITATED SINGLE ALLOCATION P-HUB MEDIAN PROBLEM WITH MULTIPLE CAPACITY LEVELS A LAGRANGIAN RELAXATION FOR CAPACITATED SINGLE ALLOCATION P-HUB MEDIAN PROBLEM WITH MULTIPLE CAPACITY LEVELS Ching-Jung Ting Department of Industrial Engineering and Management, Yuan Ze University Kuo-Rui

More information

18 hours nodes, first feasible 3.7% gap Time: 92 days!! LP relaxation at root node: Branch and bound

18 hours nodes, first feasible 3.7% gap Time: 92 days!! LP relaxation at root node: Branch and bound The MIP Landscape 1 Example 1: LP still can be HARD SGM: Schedule Generation Model Example 157323 1: LP rows, still can 182812 be HARD columns, 6348437 nzs LP relaxation at root node: 18 hours Branch and

More information

Copositive Programming and Combinatorial Optimization

Copositive Programming and Combinatorial Optimization Copositive Programming and Combinatorial Optimization Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria joint work with I.M. Bomze (Wien) and F. Jarre (Düsseldorf) IMA

More information