Interior Point Methods and Column Generation
|
|
- Emory Nelson
- 6 years ago
- Views:
Transcription
1 School of Mathematics T H E U N I V E R S I T Y O H F R G E D I N B U Interior Point Methods and Column Generation Jacek Gondzio J.Gondzio@ed.ac.uk URL: Bromont, June
2 Outline Part 1: IPMs for Optimization IPM tricks: log barrier, central path polynomial complexity optimal partition Part 2: Warmstarting IPMs Part 3: Column Generation with IPM cutting stock problem vehicle routing problem with time windows Conclusions Bromont, June
3 Part 1: Interior Point Methods for Optimization Bromont, June
4 Interior Point Methods re-born in 1984 Narendra Karmarkar, AT&T Bell Labs Shocking mathematical concept: Take linear optimization problem and add nonlinear function to the objective. A step against common sense and centuries of mathematical practice: nonlinearize linear problem Bromont, June
5 Logarithmic barrier -ln x lnx j replaces the inequality x j 0. Observe that min e n j=1 lnx j max n j=1 x j 1 x The minimization of n j=1 lnx j is equivalent to the maximization of the product of distances from all hyperplanes defining the positive orthant: it prevents all x j from approaching zero. Bromont, June
6 LP Problem: min c T x s.t. Ax = b, x 0. LP Barrier Prob: min c T x µ n lnx j s.t. Ax = b. j=1 Lagrangian: L(x,y,µ) = c T x y T (Ax b) µ n j=1 Stationarity: x L(x,y,µ) = c A T y µx 1 e = 0 y L(x,y,µ) = Ax b = 0. Denote: s = µx 1 e, i.e. XSe = µe. lnx j, Bromont, June
7 Complementarity in the Interior Point Method The first order optimality conditions (FOC) Ax = b, A T y + s = c, XSe = µe, x,s 0, where X = diag{x j }, S = diag{s j } and e = (1,,1) R n. Analytic centre (µ-centre): a (unique) point (x(µ), y(µ), s(µ)), x(µ) > 0, s(µ) > 0 that satisfies FOC. The interior point method gradually reduces the complementarity products x j s j µ 0 j = 1,2,...,n. Bromont, June
8 Interior Point Methods Theory: convergence in O( n) or O(n) iterations Practice: convergence in O(log n) iterations Expected number of IPM iterations: Problem Dimension LP QP 1, , , ,000, ,000, ,000, ,000, but one iteration may be expensive! Bromont, June
9 Complementarity x j s j = 0 j = 1,2,...,n. Simplex Method guesses an optimal partition: For basic variables, s B = 0 and (x B ) j (s B ) j = 0 j B. For non-basic variables, x N = 0 hence (x N ) j (s N ) j = 0 j N. Interior Point Method uses ε-mathematics: Replace x j s j = 0 j = 1,2,...,n by x j s j = µ j = 1,2,...,n. Force convergence µ 0. Bromont, June
10 First Order Optimality Conditions Simplex Method: Ax = b A T y + s = c XSe = 0 x, s 0. Interior Point Method: Ax = b A T y + s = c XSe = µe x, s 0. s s s s x Basic: x > 0, s = 0 Nonbasic: x = 0, s > 0 x x "Basic": x > 0, s = 0 "Nonbasic": x = 0, s > 0 x G, IPMs 25 years later, EJOR 218 (2012), Bromont, June
11 Part 2: Warmstarting IPMs Bromont, June
12 A need to solve a sequence of similar problems column generation cutting plane methods subproblems in the block-angular LPs (Dantzig-Wolfe decomp., Benders decomp.) B&B, (and B&Cut, B&Cut&Price, etc) SQP any sequence of similar problems example: computing efficient frontier in Markowitz portfolio optimization Bromont, June
13 Warm Starts Which method should be used? Simplex Method, or Interior Point Method. When is the Simplex Method better? few indices change optimal partition B & B, adding one cut in CPM, etc. When is the Interior Point Method better? many indices change optimal partition adding many cuts in CPM, dealing with a general change of problem data, etc Conjecture: The more changes in the (large) problem the more attractive IPM-based warm starts are. Bromont, June
14 Difficulty of IPM Warm Starts Modified Problem Original Problem Bromont, June
15 Warm Starting in Mitchell, PhD Thesis, Cornell Univ Goffin & Vial et al., development of ACCPM G. & Sarkissian, development of PDCGM in 1995 G., Math. Prog. 83 (1998) G. & Vial, COAP 14 (1999) ACCPM Analytic Centre Cutting Plane Method PDCGM Primal-Dual Column Generation Method Bromont, June
16 Warmstarting Heuristic Idea: Start close to the (new) central path, not close to the (old) solution Modified Problem Original Problem G., Mathematical Programming 83 (1998) Bromont, June
17 Warm Start with µ-centres Old Problem: New Problem: min c T 0 x xt Q 0 x s.t. A 0 x = b 0, x 0, min c T x xt Qx s.t. Ax = b, x 0, We assume: c c 0, Q Q 0, A A 0, b b 0. Bromont, June
18 Warm Starting in Yildirim & Wright, SIOPT 12 (2002) G. & Grothey, SIOPT 13 (2003) Fliege, Maths of OR 31 (2006) Benson & Shanno, COAP 38 (2007) Benson & Shanno, COAP 40 (2008) G. & Grothey, SIOPT 19 (2008) John & Yildirim, COAP 41 (2008) Colombo, G. & Grothey, MP 127 (2011) Colombo & Grothey, follow-up reports in 09,10 Engau, Anjos & Vannelli, SIOPT 20 (2010) 1828 Benson & Mahanta, report in 2009 Ordonez & Waltz, report in 2009 Bromont, June
19 IPM Warmstarts: Theoretical Results Yildirim & Wright, SIOPT 12 (2002) G. & Grothey, SIOPT 13 (2003) Lemma. Let (x,y,s) N (γ 0 ) for problem (LP) then the full Newton step ( x, y, s) in the perturbed problem ( LP) is feasible and (x + x,y + y,s + s) Ñ (γ) provided that δ bc = ξ c 2 + A T (AA T ) 1 γ ξ b 2 P /γ µ, where P = I S 1 A T (AXS 1 A T ) 1 AX, ξ b = b Ax, ξc = c A T y s. Bromont, June
20 LOQO vs OOPS warmstarting NETLIB problems Benson & Shanno, COAP 38 (2007) G. & Grothey, SIOPT 19 (2008) Unblocking technique... Average savings: LOQO (B&S, 2007) % OOPS (G&G, 2008) % Bromont, June
21 Part 3: Primal-Dual Column Generation Method Joint work with two PhD students: Pablo Gonzalez-Brevis and Pedro Munari Bromont, June
22 Column Generation (CG) MP... RMP newrmp Bromont, June
23 Column Generation (CG) Consider an LP, called the master problem (MP): z := min j N c j λ j, s.t. j N a j λ j = b, λ j 0, j N. N is too big; The columns a j are implicit elements of A; We know how to generate them! Bromont, June
24 CG: Restricted master problem (RMP): N N z RMP := min j N c j λ j, s.t. j N a j λ j = b, λ j 0, j N. Optimal λ for the RMP feasible ˆλ for the MP; ˆλj = λj, j N, and ˆλj = 0 otherwise; Hence, z z RMP = UB (Upper Bound). How to know it is optimal? Call the oracle! Bromont, June
25 CG: Oracle: check the feasibility of the dual u; Reduced costs: s j = c j u T a j, j N; But the columns are not explicit and, hence, z SP := min{c j u T a j a j A}. (we reset z SP := 0, if z SP > 0); Lower Bound: LB = z RMP + κz SP z, where κ i N λ i, If z SP < 0, then new columns are generated; Otherwise, an optimal solution of the MP was found! Bromont, June
26 Appealing features of IPMs Use IPM to solve the RMP: no degeneracy issues Terminate RMP solution early: get stable dual solution ū Bromont, June
27 PDCGM Algorithm Parameters: ε max, D, δ, κ 1. set LB =, UB =, gap =, ε = 0.5; 2. while (gap > δ) do 3. find a well-centred ε-optimal ( λ,ũ) of the RMP; 4. UB = z RMP ; 5. call the oracle with the query point ũ; 6. LB = κ z SP + b T ũ; 7. gap = (UB LB)/(1 + UB ); 8. ε = min{ε max, gap/d}; 9. if ( z SP < 0) then add new columns into the RMP; 10. end(while) Bromont, June
28 CSP: Column Generation Formulation Gilmore and Gomory (1961) formulation: min s.t. p P p P λ p, a p λ p d, λ p 0 and integer, p P. Columns are cutting patterns; We do not need to enumerate all of them; They can be dynamically generated knapsack problem. Bromont, June
29 VRPTW: Column Generation Formulation Desrochers et al. (1992): min s.t. p P p P c p λ p a p λ p = 1, λ p {0,1}, p P. Columns are possible vehicle paths; The columns can be dynamically generated shortest path problem with resource constraints. Bromont, June
30 Computational experiments Solving LP relaxations Comparison of: Standard column generation (SCG): simplex-type methods of IBM/CPLEX v Primal dual column generation (PDCGM): interior point solver HOPDM. Analytic centre cutting plane (ACCPM): open-source solver OBOE/COIN-OR. Bromont, June
31 Cutting stock problem SCG PDCGM ACCPM Cuts Class iters time iters time iters time 10 Small Large Small Large Small Large instances: 178 small (m 199), 84 large (m 200) Bromont, June
32 CSP: Larger Instances: BPP-U09??? family SCG PDCGM ACCPM Instance m iters time iters time iters time U U U U U U U U Bromont, June
33 Vehicle routing problem with time windows SCG PDCGM ACCPM Cuts Class iters time iters time iters time 10 Small Medium Large Small Medium Large instances: 29 small (n = 25), 29 med (n = 50), 29 large (n = 100) Bromont, June
34 VRPTW: Larger Instances SCG PDCGM ACCPM Instance n iters time iters time iters time R C RC R C RC R C RC Bromont, June
35 Integer optimization Bromont, June
36 Integer VRPTW solved to optimality Branch-Price-and-Cut, Pedro Munari s PhD DLH08 I-PDCGM Problem cuts nodes time cuts nodes time C C RC RC R R DLH08: Desaulniers, Lessard & Hadjar, Transportation Science 42 (2008) Solomon, Operations Research 35 (1987) Homberger&Gehring, EJOR 162 (2005) Bromont, June
37 Conclusions A completely new perspective is needed to exploit the insight offered by IPMs in a number of combinatorial optimization applications: column generation cutting plane methods B & B, (and B & Cut, B & Cut & Price, etc) Warmstarting works well in the CG context: problems are re-optimized in 3-5 IPM iterations Bromont, June
38 References G., Gonzalez-Brevis, Munari, New developments in the primal-dual column generation technique, ERGO Tech Rep, Edinburgh, Munari, G., Using the primal-dual interior point algorithm within the branch-price-and-cut method, ERGO Tech Rep, Edinburgh, G., Gonzalez-Brevis, A new warmstarting strategy for the primal-dual column generation method, ERGO Tech Rep, Edinburgh, Bromont, June
39 Example: Cutting Stock Problem (CSP) A set N of large pieces of wood of length W is given. We need to cut them into smaller pieces. We need d j units of small piece j M of length w j. Minimize the number of units of large pieces of wood. Define binary variable y i which takes value 1 if i-th large piece of wood is cut and 0 if it is not used. Define integer variable x ij which determines the number of units of small piece of wood j M obtained by cutting the large piece i N. Bromont, June
40 Cutting Stock Problem (CSP) Kantorovich s formulation: min s.t. j M i N i N y i x ij d j j M, x ij w j Wy i i N, y i {0,1} i N, x ij 0 and integer j M, i N. LP relaxation gives very weak bound. Bromont, June
41 Vehicle Routing Problem with Time Windows A company delivers goods to customers i C. The company has vehicles k V and each of them starts at a depot, travels to several customers and returns to the depot. The visit of vehicle k to customer i needs to take place in a specific time window: a i s ik b i, where s ik is the time when vehicle k reaches customer i. Objective: Minimize the total cost of delivery. Define binary variable x ijk which takes value 1 if vehicle k travels from customer i to customer j (k V, i,j C) and takes value zero otherwise. Bromont, June
42 Vehicle Routing Problem with Time Windows Constraints: Exactly one vehicle leaves customer i: k V j N Vehicle capacity constraint: i C d i j N x ijk = 1, x ijk q, i C k V Each vehicle leaves the depot and returns to it: j N x 0jk = 1 and j N x i(n+1)k = 1, k V Bromont, June
43 VRPTW: Constraints (continued) Time-window constraint s ik + t ij M(1 x ijk ) s jk, i,j N, k V. Since x ijk is binary the above constraint has the following meaning: If x ijk = 1 (vehicle k travels from customer i to customer j) then s ik + t ij s jk that is, the arrival time of vehicle k to customer j is greater than or equal the sum of time when vehicle k arrives to customer i and the time t ij it takes to travel from i to j. Otherwise (if x ijk = 0) the presence of big M guarantees that the constraint is always inactive. Bromont, June
44 VRPTW min j N s.t. i C x 0jk = 1, k V k V d i j N i N i N j N i N x ihk j N j N c ij x ijk x ijk = 1, i C, x ijk q, k V, x i(n+1)k = 1, k V, x jhk = 0, h C, k V, s ik + t ij M(1 x ijk ) s jk, i,j N, k V, a i s ik b i, i N, k V, x ijk {0,1}, i,j N, k V. Bromont, June
A new Decomposition Algorithm for Multistage Stochastic Programs with Endogenous Uncertainties
A new Decomposition Algorithm for Multistage Stochastic Programs with Endogenous Uncertainties Vijay Gupta Ignacio E. Grossmann Department of Chemical Engineering Carnegie Mellon University, Pittsburgh
More informationA computational study of on-demand accuracy level decomposition for twostage stochastic programs
A computational study of on-demand accuracy level decomposition for twostage stochastic programs Christian Wolf, Universität Paderborn, Germany Csaba I. Fábián, Kecskemét College, Hungary Achim Koberstein,
More informationOptimizing Cyclist Parking in a Closed System
Optimizing Cyclist Parking in a Closed System Letu Qingge, Killian Smith Gianforte School of Computing, Montana State University, Bozeman, MT 59717, USA Abstract. In this paper, we consider the two different
More informationPrediction Market and Parimutuel Mechanism
Prediction Market and Parimutuel Mechanism Yinyu Ye MS&E and ICME Stanford University Joint work with Agrawal, Peters, So and Wang Math. of Ranking, AIM, 2 Outline World-Cup Betting Example Market for
More informationScheduling the Brazilian Soccer Championship. Celso C. Ribeiro* Sebastián Urrutia
Scheduling the Brazilian Soccer Championship Celso C. Ribeiro* Sebastián Urrutia Motivation Problem statement Solution approach Summary Phase 1: create all feasible Ps Phase 2: assign Ps to elite teams
More informationi) Linear programming
Problem statement Sailco Corporation must determine how many sailboats should be produced during each of the next four quarters (one quarter = three months). The demand during each of the next four quarters
More informationPolynomial DC decompositions
Polynomial DC decompositions Georgina Hall Princeton, ORFE Joint work with Amir Ali Ahmadi Princeton, ORFE 7/31/16 DIMACS Distance geometry workshop 1 Difference of convex (dc) programming Problems of
More informationProduct Decomposition in Supply Chain Planning
Mario R. Eden, Marianthi Ierapetritou and Gavin P. Towler (Editors) Proceedings of the 13 th International Symposium on Process Systems Engineering PSE 2018 July 1-5, 2018, San Diego, California, USA 2018
More informationLecture 10. Support Vector Machines (cont.)
Lecture 10. Support Vector Machines (cont.) COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Andrey Kan Copyright: University of Melbourne This lecture Soft margin SVM Intuition and problem
More informationCSE 3401: Intro to AI & LP Uninformed Search II
CSE 3401: Intro to AI & LP Uninformed Search II Required Readings: R & N Chapter 3, Sec. 1-4. 1 {Arad}, {Zerind, Timisoara, Sibiu}, {Zerind, Timisoara, Arad, Oradea, Fagaras, RimnicuVilcea }, {Zerind,
More informationDepth-bounded Discrepancy Search
Depth-bounded Discrepancy Search Toby Walsh* APES Group, Department of Computer Science University of Strathclyde, Glasgow Gl 1XL. Scotland tw@cs.strath.ac.uk Abstract Many search trees are impractically
More informationA IMPROVED VOGEL S APPROXIMATIO METHOD FOR THE TRA SPORTATIO PROBLEM. Serdar Korukoğlu 1 and Serkan Ballı 2.
Mathematical and Computational Applications, Vol. 16, No. 2, pp. 370-381, 2011. Association for Scientific Research A IMPROVED VOGEL S APPROXIMATIO METHOD FOR THE TRA SPORTATIO PROBLEM Serdar Korukoğlu
More informationBetter Search Improved Uninformed Search CIS 32
Better Search Improved Uninformed Search CIS 32 Functionally PROJECT 1: Lunar Lander Game - Demo + Concept - Open-Ended: No One Solution - Menu of Point Options - Get Started NOW!!! - Demo After Spring
More informationCSE 3402: Intro to Artificial Intelligence Uninformed Search II
CSE 3402: Intro to Artificial Intelligence Uninformed Search II Required Readings: Chapter 3, Sec. 1-4. 1 {Arad}, {Zerind, Timisoara, Sibiu}, {Zerind, Timisoara, Arad, Oradea, Fagaras, RimnicuVilcea },
More informationEE 364B: Wind Farm Layout Optimization via Sequential Convex Programming
EE 364B: Wind Farm Layout Optimization via Sequential Convex Programming Jinkyoo Park 1 Introduction In a wind farm, the wakes formed by upstream wind turbines decrease the power outputs of downstream
More informationDecomposition Techniques in Mathematical Programming
Decomposition Techniques in Mathematical Programming Antonio J. Conejo Roberto Mínguez Enrique Castillo Raquel García-Bertrand Decomposition Techniques in Mathematical Programming Engineering and Science
More informationOcean Fishing Fleet Scheduling Path Optimization Model Research. Based On Improved Ant Colony Algorithm
4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 205) Ocean Fishing Fleet Scheduling Path Optimization Model Research Based On Improved Ant Colony Algorithm Li Jia-lin2,
More informationSolving MINLPs with BARON. Mustafa Kılınç & Nick Sahinidis Department of Chemical Engineering Carnegie Mellon University
Solving MINLPs with BARON Mustafa Kılınç & Nick Sahinidis Department of Chemical Engineering Carnegie Mellon University MINLP 2014 Carnegie Mellon University June 4, 2014 MIXED-INTEGER NONLINEAR PROGRAMS
More informationThree New Methods to Find Initial Basic Feasible. Solution of Transportation Problems
Applied Mathematical Sciences, Vol. 11, 2017, no. 37, 1803-1814 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.75178 Three New Methods to Find Initial Basic Feasible Solution of Transportation
More informationMathematics of Pari-Mutuel Wagering
Millersville University of Pennsylvania April 17, 2014 Project Objectives Model the horse racing process to predict the outcome of a race. Use the win and exacta betting pools to estimate probabilities
More informationECE 697B (667) Spring 2003
ECE 667 - Synthesis & Verification - Lecture 2 ECE 697 (667) Spring 23 Synthesis and Verification of Digital Systems unctional Decomposition Slides adopted (with permission) from. Mishchenko, 23 Overview
More informationA chartered bus allocation problem
A chartered bus allocation problem Woo-Je Kim Department of Industrial Engineering, Daejin University, San11-1 Sundanri, Pochun, Kyungkido 487-711, Korea phone:+82-31-539-2001; fax:+82-31-539-2000; wjkim@road.daejin.ac.kr
More informationCS 4649/7649 Robot Intelligence: Planning
CS 4649/7649 Robot Intelligence: Planning Differential Kinematics, Probabilistic Roadmaps Sungmoon Joo School of Interactive Computing College of Computing Georgia Institute of Technology S. Joo (sungmoon.joo@cc.gatech.edu)
More informationTime Dependent Truck Routing and Driver Scheduling Problem with Hours of Service Regulations. Vidit Divyang Shah
Time Dependent Truck Routing and Driver Scheduling Problem with Hours of Service Regulations A Thesis Presented by Vidit Divyang Shah to The Department of Mechanical and Industrial Engineering In partial
More informationRobust planning and disruption management in roll-on roll-off liner shipping
Robust planning and disruption management in roll-on roll-off liner shipping Andreas Fischer a, Håkon Nokhart a, Henrik Olsen a, Kjetil Fagerholt a,b, Jørgen Glomvik Rakke c, Magnus Stålhane a a Norwegian
More informationconsist of friends, is open to all ages, and considers fair play of paramount importance. The matches are played without referees, since, according to
Indoor football scheduling Dries Goossens Frits Spieksma Abstract This paper deals with a real-life scheduling problem from an amateur indoor football league. The league consists of a number of divisions,
More informationSearch I. Tuomas Sandholm Carnegie Mellon University Computer Science Department. [Read Russell & Norvig Chapter 3]
Search I Tuomas Sandholm Carnegie Mellon University Computer Science Department [Read Russell & Norvig Chapter 3] Search I Goal-based agent (problem solving agent) Goal formulation (from preferences).
More informationAmmonia Synthesis with Aspen Plus V8.0
Ammonia Synthesis with Aspen Plus V8.0 Part 2 Closed Loop Simulation of Ammonia Synthesis 1. Lesson Objectives Review Aspen Plus convergence methods Build upon the open loop Ammonia Synthesis process simulation
More informationAGA Swiss McMahon Pairing Protocol Standards
AGA Swiss McMahon Pairing Protocol Standards Final Version 1: 2009-04-30 This document describes the Swiss McMahon pairing system used by the American Go Association (AGA). For questions related to user
More informationLecture 5. Optimisation. Regularisation
Lecture 5. Optimisation. Regularisation COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Andrey Kan Copyright: University of Melbourne Iterative optimisation Loss functions Coordinate
More informationBelow are the graphing equivalents of the above constraints.
Detailed Step- by- Step description of the Wilson Problem Part A. LP Formulation Step 1. Figure out the problem Although it may seem obvious, the first step in solving the Wilson Problem is to understand
More informationA SEMI-PRESSURE-DRIVEN APPROACH TO RELIABILITY ASSESSMENT OF WATER DISTRIBUTION NETWORKS
A SEMI-PRESSURE-DRIVEN APPROACH TO RELIABILITY ASSESSMENT OF WATER DISTRIBUTION NETWORKS S. S. OZGER PhD Student, Dept. of Civil and Envir. Engrg., Arizona State Univ., 85287, Tempe, AZ, US Phone: +1-480-965-3589
More informationCSC384: Introduction to Artificial Intelligence. Search
CSC384: Introduction to Artificial Intelligence Search Chapter 3 of R&N 3 rd edition is very useful reading. Chapter 4 of R&N 3 rd edition is worth reading for enrichment. We ll touch upon some of the
More information2 When Some or All Labels are Missing: The EM Algorithm
CS769 Spring Advanced Natural Language Processing The EM Algorithm Lecturer: Xiaojin Zhu jerryzhu@cs.wisc.edu Given labeled examples (x, y ),..., (x l, y l ), one can build a classifier. If in addition
More informationProblem Solving as Search - I
Problem Solving as Search - I Shobhanjana Kalita Dept. of Computer Science & Engineering Tezpur University Slides prepared from Artificial Intelligence A Modern approach by Russell & Norvig Problem-Solving
More informationMixture Models & EM. Nicholas Ruozzi University of Texas at Dallas. based on the slides of Vibhav Gogate
Mixture Models & EM Nicholas Ruozzi University of Texas at Dallas based on the slides of Vibhav Gogate Previously We looked at -means and hierarchical clustering as mechanisms for unsupervised learning
More information1.1 The size of the search space Modeling the problem Change over time Constraints... 21
Introduction : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 1 I What Are the Ages of My Three Sons? : : : : : : : : : : : : : : : : : 9 1 Why Are Some Problems Dicult to Solve? : : :
More informationThe Evolution of Transport Planning
The Evolution of Transport Planning On Proportionality and Uniqueness in Equilibrium Assignment Michael Florian Calin D. Morosan Background Several bush-based algorithms for computing equilibrium assignments
More informationAuthor s Name Name of the Paper Session. Positioning Committee. Marine Technology Society. DYNAMIC POSITIONING CONFERENCE September 18-19, 2001
Author s Name Name of the Paper Session PDynamic Positioning Committee Marine Technology Society DYNAMIC POSITIONING CONFERENCE September 18-19, 2001 POWER PLANT SESSION A New Concept for Fuel Tight DP
More informationSolving Problems by Searching chap3 1. Problem-Solving Agents
Chapter3 Solving Problems by Searching 20070315 chap3 1 Problem-Solving Agents Reflex agents cannot work well in those environments - state/action mapping too large - take too long to learn Problem-solving
More informationDesigning a Traffic Circle By David Bosworth For MATH 714
Designing a Traffic Circle By David Bosworth For MATH 714 Abstract We all have had the experience of sitting in a traffic jam, or we ve seen cars bunched up on a road for some no apparent reason. The use
More informationEE582 Physical Design Automation of VLSI Circuits and Systems
EE Prof. Dae Hyun Kim School of Electrical Engineering and Computer Science Washington State University Routing Grid Routing Grid Routing Grid Routing Grid Routing Grid Routing Lee s algorithm (Maze routing)
More informationCENG 466 Artificial Intelligence. Lecture 4 Solving Problems by Searching (II)
CENG 466 Artificial Intelligence Lecture 4 Solving Problems by Searching (II) Topics Search Categories Breadth First Search Uniform Cost Search Depth First Search Depth Limited Search Iterative Deepening
More informationStrategies for protecting supply chain networks against facility and transportation disruptions: an improved Benders decomposition approach
Ann Oper Res (2013) 210:125 163 DOI 10.1007/s10479-012-1146-x Strategies for protecting supply chain networks against facility and transportation disruptions: an improved Benders decomposition approach
More informationA Behavioral Theory. Teck H. Ho. & National University of Singapore. Joint Work with Noah Lim, Houston Juanjuan Zhang, MIT. July 2008 Teck H.
Designing Pricing Contracts for Homo Sapiens: A Behavioral Theory Teck H. Ho University of California, Berkeley & National University of Singapore Joint Work with Noah Lim, Houston Juanjuan Zhang, MIT
More informationCHAPTER 1 INTRODUCTION TO RELIABILITY
i CHAPTER 1 INTRODUCTION TO RELIABILITY ii CHAPTER-1 INTRODUCTION 1.1 Introduction: In the present scenario of global competition and liberalization, it is imperative that Indian industries become fully
More informationDynamic configuration of QC allocating problem based on multi-objective genetic algorithm
DOI.7/s8---7 Dynamic configuration of QC allocating problem based on multi-objective genetic algorithm ChengJi Liang MiaoMiao Li Bo Lu Tianyi Gu Jungbok Jo Yi Ding Received: vember / Accepted: January
More informationOPTIMAL FLOWSHOP SCHEDULING WITH DUE DATES AND PENALTY COSTS
J. Operation Research Soc. of Japan VoJ. 1, No., June 1971. 1971 The Operations Research Society of Japan OPTMAL FLOWSHOP SCHEDULNG WTH DUE DATES AND PENALTY COSTS JATNDER N.D. GUPTA Assistant Professor,
More informationChapter 20. Planning Accelerated Life Tests. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 20 Planning Accelerated Life Tests William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
More informationHydronic Systems Balance
Hydronic Systems Balance Balancing Is Misunderstood Balancing is application of fundamental hydronic system math Balance Adjustment of friction loss location Adjustment of pump to requirements By definition:
More informationThe Mobility and Safety of Walk-and-Ride Systems
The Mobility and Safety of Walk-and-Ride Systems By Dr. Hugh Medal Department of Industrial and Systems Engineering Mississippi State University (MSU) Phone: 662-325-3932 Email: hugh.medal@msstate.edu
More informationThe Orienteering Problem
The Orienteering Problem Bruce L. Golden, Larry Levy, and Rakesh Vohra* College of Business and Management, University of Maryland, College Park, Maryland 20742 Onenteenng is a sport in which start and
More informationTASK 4.2.1: "HAMILTON'S ROBOT" ==============================
TASK 4.2.1: "HAMILTON'S ROBOT" ============================== On a plane there are given N positions P1, P2,..., PN with integer coordinates (X1,Y1), (X2,Y2),..., (XN,YN). A robot should move through all
More informationUAID: Unconventional Arterial Intersection Design
UAID: Unconventional Arterial Intersection Design 1 Content Overview of Unconventional Arterial Intersection Design System Structure Continuous Flow Intersection Diverging Diamond Interchange Analysis
More informationEvacuation Time Minimization Model using Traffic Simulation and GIS Technology
Evacuation Time Minimization Model using Traffic Simulation and GIS Technology Daisik Danny Nam, Presenter Ph.D. Student Department of Civil and Environmental Engineering, University of California, Irvine
More informationTraffic circles. February 9, 2009
Traffic circles February 9, 2009 Abstract The use of a traffic circle is a relatively common means of controlling traffic in an intersection. Smaller Traffic circles can be especially effective in routing
More informationTIMETABLING IN SPORTS AND ENTERTAINMENT
TIMETABLING IN SPORTS AND ENTERTAINMENT Introduction More complicated applications of interval scheduling and timetabling: scheduling of games in tournaments scheduling of commercials on network television
More informationOn the convergence of fitting algorithms in computer vision
On the convergence of fitting algorithms in computer vision N. Chernov Department of Mathematics University of Alabama at Birmingham Birmingham, AL 35294 chernov@math.uab.edu Abstract We investigate several
More informationCS 4649/7649 Robot Intelligence: Planning
CS 4649/7649 Robot Intelligence: Planning Partially Observable MDP Sungmoon Joo School of Interactive Computing College of Computing Georgia Institute of Technology S. Joo (sungmoon.joo@cc.gatech.edu)
More informationThe system design must obey these constraints. The system is to have the minimum cost (capital plus operating) while meeting the constraints.
Computer Algorithms in Systems Engineering Spring 2010 Problem Set 6: Building ventilation design (dynamic programming) Due: 12 noon, Wednesday, April 21, 2010 Problem statement Buildings require exhaust
More informationif all agents follow RSS s interpretation then there will be zero accidents.
RSS Concept RSS - Mobileye SFF - Nvidia Safety Goal Guaranteeing that an agent will never be involved in an accident is impossible. Hence, our ultimate goal is to guarantee that an agent will be careful
More informationExcel Solver Case: Beach Town Lifeguard Scheduling
130 Gebauer/Matthews: MIS 213 Hands-on Tutorials and Cases, Spring 2015 Excel Solver Case: Beach Town Lifeguard Scheduling Purpose: Optimization under constraints. A. GETTING STARTED All Excel projects
More informationInverting a Batting Average - an Application of Continued Fractions (Preliminary Version)
Inverting a Batting Average - an Application of Continued Fractions (Preliminary Version) Allen Back August 30, 2000 Contents Introduction and Motivation 2 Continued Fractions 2 2. Comparison of Two Numbers
More informationRobust specification testing in regression: the FRESET test and autocorrelated disturbances
Robust specification testing in regression: the FRESET test and autocorrelated disturbances Linda F. DeBenedictis and David E. A. Giles * Policy and Research Division, Ministry of Human Resources, 614
More informationOptimization of a Golf Driver
Optimization of a Golf Driver By: Max Dreager Cole Snider Anthony Boyd Frank Rivera Final Report MAE 494: Design Optimization Due: May 10 Abstract In this study, the relationship between a golf ball and
More informationAryeh Rappaport Avinoam Meir. Schedule automation
Aryeh Rappaport Avinoam Meir Schedule automation Introduction In this project we tried to automate scheduling. Which is to help a student to pick the time he takes the courses so they want collide with
More informationSupport Vector Machines: Optimization of Decision Making. Christopher Katinas March 10, 2016
Support Vector Machines: Optimization of Decision Making Christopher Katinas March 10, 2016 Overview Background of Support Vector Machines Segregation Functions/Problem Statement Methodology Training/Testing
More informationMinimum Mean-Square Error (MMSE) and Linear MMSE (LMMSE) Estimation
Minimum Mean-Square Error (MMSE) and Linear MMSE (LMMSE) Estimation Outline: MMSE estimation, Linear MMSE (LMMSE) estimation, Geometric formulation of LMMSE estimation and orthogonality principle. Reading:
More informationLogistic Regression. Hongning Wang
Logistic Regression Hongning Wang CS@UVa Today s lecture Logistic regression model A discriminative classification model Two different perspectives to derive the model Parameter estimation CS@UVa CS 6501:
More informationArterial Traffic Analysis Actuated Signal Control
Arterial Traffic Analysis Actuated Signal Control Dr. Gang-Len Chang Professor and Director of Traffic Safety and Operations Lab. University of Maryland-College Park Actuated Signal Control Fully Actuated
More informationHIGH RESOLUTION DEPTH IMAGE RECOVERY ALGORITHM USING GRAYSCALE IMAGE.
HIGH RESOLUTION DEPTH IMAGE RECOVERY ALGORITHM USING GRAYSCALE IMAGE Kazunori Uruma 1, Katsumi Konishi 2, Tomohiro Takahashi 1 and Toshihiro Furukawa 1 1 Graduate School of Engineering, Tokyo University
More informationCycle Analysis and Construction of Protographs for QC LDPC Codes With Girth Larger Than 12 ½
Cycle Analysis and Construction of Protographs for QC LDPC Codes With Girth Larger Than 2 ½ Sunghwan Kim, Jong-Seon No School of Electrical Eng & Com Sci Seoul National University, Seoul, Korea Email:
More informationImperfectly Shared Randomness in Communication
Imperfectly Shared Randomness in Communication Madhu Sudan Harvard Joint work with Clément Canonne (Columbia), Venkatesan Guruswami (CMU) and Raghu Meka (UCLA). 11/16/2016 UofT: ISR in Communication 1
More informationINSTRUMENT INSTRUMENTAL ERROR (of full scale) INSTRUMENTAL RESOLUTION. Tutorial simulation. Tutorial simulation
Lab 1 Standing Waves on a String Learning Goals: To distinguish between traveling and standing waves To recognize how the wavelength of a standing wave is measured To recognize the necessary conditions
More informationOptimization and Search. Jim Tørresen Optimization and Search
Optimization and Search INF3490 - Biologically inspired computing Lecture 1: Marsland chapter 9.1, 9.4-9.6 2017 Optimization and Search Jim Tørresen 2 Optimization and Search Methods (selection) Optimization
More informationA Fair Target Score Calculation Method for Reduced-Over One day and T20 International Cricket Matches
A Fair Target Score Calculation Method for Reduced-Over One day and T20 International Cricket Matches Rohan de Silva, PhD. Abstract In one day internationals and T20 cricket games, the par score is defined
More informationUninformed Search (Ch )
1 Uninformed Search (Ch. 3-3.4) 7 Small examples 8-Queens: how to fit 8 queens on a 8x8 board so no 2 queens can capture each other Two ways to model this: Incremental = each action is to add a queen to
More informationDecision Trees. Nicholas Ruozzi University of Texas at Dallas. Based on the slides of Vibhav Gogate and David Sontag
Decision Trees Nicholas Ruozzi University of Texas at Dallas Based on the slides of Vibhav Gogate and David Sontag Announcements Course TA: Hao Xiong Office hours: Friday 2pm-4pm in ECSS2.104A1 First homework
More informationDESIGN AND ANALYSIS OF ALGORITHMS (DAA 2017)
DESIGN AND ANALYSIS OF ALGORITHMS (DAA 2017) Veli Mäkinen 12/05/2017 1 COURSE STRUCTURE 7 weeks: video lecture -> demo lecture -> study group -> exercise Video lecture: Overview, main concepts, algorithm
More informationUninformed search methods
Lecture 3 Uninformed search methods Milos Hauskrecht milos@cs.pitt.edu 5329 Sennott Square Announcements Homework 1 Access through the course web page http://www.cs.pitt.edu/~milos/courses/cs2710/ Two
More informationCalculation of Trail Usage from Counter Data
1. Introduction 1 Calculation of Trail Usage from Counter Data 1/17/17 Stephen Martin, Ph.D. Automatic counters are used on trails to measure how many people are using the trail. A fundamental question
More informationOptimal Weather Routing Using Ensemble Weather Forecasts
Optimal Weather Routing Using Ensemble Weather Forecasts Asher Treby Department of Engineering Science University of Auckland New Zealand Abstract In the United States and the United Kingdom it is commonplace
More informationCS 4649/7649 Robot Intelligence: Planning
CS 4649/7649 Robot Intelligence: Planning Roadmap Approaches Sungmoon Joo School of Interactive Computing College of Computing Georgia Institute of Technology S. Joo (sungmoon.joo@cc.gatech.edu) 1 *Slides
More informationBézier Curves and Splines
CS-C3100 Computer Graphics Bézier Curves and Splines Majority of slides from Frédo Durand vectorportal.com CS-C3100 Fall 2016 Lehtinen Before We Begin Anything on your mind concerning Assignment 1? Linux
More informationUsing GPOPS-II to optimize sum of squared torques of a double pendulum as a prosthesis leg. Abstract
Using GPOPS-II to optimize sum of squared torques of a double pendulum as a prosthesis leg Abstract Milad Zarei MCE 593 Prosthesis Design & Control A two-dimensional, two links pendulum is developed to
More informationIAC-06-D4.1.2 CORRELATIONS BETWEEN CEV AND PLANETARY SURFACE SYSTEMS ARCHITECTURE PLANNING Larry Bell
IAC-06-D4.1.2 CORRELATIONS BETWEEN CEV AND PLANETARY SURFACE SYSTEMS ARCHITECTURE PLANNING Larry Bell Sasakawa International Center for Space Architecture (SICSA), University of Houston, USA e-mail: lbell@uh.edu
More informationOPTIMAL TRAJECTORY GENERATION OF COMPASS-GAIT BIPED BASED ON PASSIVE DYNAMIC WALKING
OPTIMAL TRAJECTORY GENERATION OF COMPASS-GAIT BIPED BASED ON PASSIVE DYNAMIC WALKING Minseung Kim Dept. of Computer Science Illinois Institute of Technology 3201 S. State St. Box 2082 Chicago IL 60616
More informationUninformed search methods II.
CS 1571 Introduction to AI Lecture 5 Uninformed search methods II. Milos Hauskrecht milos@cs.pitt.edu 5329 Sennott Square Uninformed methods Uninformed search methods use only information available in
More informationOnline Companion to Using Simulation to Help Manage the Pace of Play in Golf
Online Companion to Using Simulation to Help Manage the Pace of Play in Golf MoonSoo Choi Industrial Engineering and Operations Research, Columbia University, New York, NY, USA {moonsoo.choi@columbia.edu}
More informationCVEN Computer Applications in Engineering and Construction. Programming Assignment #4 Analysis of Wave Data Using Root-Finding Methods
CVEN 30-501 Computer Applications in Engineering and Construction Programming Assignment #4 Analysis of Wave Data Using Root-Finding Methods Date distributed: 9/30/016 Date due: 10/14/016 at 3:00 PM (electronic
More informationRobust Task Execution: Procedural and Model-based. Outline. Desiderata: Robust Task-level Execution
Robust Task Execution: Procedural and Model-based Mission Goals and Environment Constraints Temporal Planner Temporal Network Solver Projective Task Expansion Initial Conditions Temporal Plan Dynamic Scheduling
More informationUninformed Search (Ch )
1 Uninformed Search (Ch. 3-3.4) 3 Terminology review State: a representation of a possible configuration of our problem Action: -how our agent interacts with the problem -can be different depending on
More informationLesson 22: Average Rate of Change
Student Outcomes Students know how to compute the average rate of change in the height of water level when water is poured into a conical container at a constant rate. MP.1 Lesson Notes This lesson focuses
More informationDetermination of the Design Load for Structural Safety Assessment against Gas Explosion in Offshore Topside
Determination of the Design Load for Structural Safety Assessment against Gas Explosion in Offshore Topside Migyeong Kim a, Gyusung Kim a, *, Jongjin Jung a and Wooseung Sim a a Advanced Technology Institute,
More informationA Point-Based Algorithm to Generate Final Table States of Football Tournaments
International Journal of Computing Academic Research (IJCAR) ISSN 2305-9184, Volume 7, Number 3 (June 2018), pp.38-42 MEACSE Publications http://www.meacse.org/ijcar A Point-Based Algorithm to Generate
More informationApplying Bi-objective Shortest Path Methods to Model Cycle Route-choice
Applying Bi-objective Shortest Path Methods to Model Cycle Route-choice Chris Van Houtte, Judith Y. T. Wang, and Matthias Ehrgott September 30, 2009 Outline Commuter Cyclists Motivation Choice Set Route
More informationChapter 4 (sort of): How Universal Are Turing Machines? CS105: Great Insights in Computer Science
Chapter 4 (sort of): How Universal Are Turing Machines? CS105: Great Insights in Computer Science Some Probability Theory If event has probability p, 1/p tries before it happens (on average). If n distinct
More informationYAN GU. Assistant Professor, University of Massachusetts Lowell. Frederick N. Andrews Fellowship, Graduate School, Purdue University ( )
YAN GU Assistant Professor, University of Massachusetts Lowell CONTACT INFORMATION 31 University Avenue Cumnock 4E Lowell, MA 01854 yan_gu@uml.edu 765-421-5092 http://www.locomotionandcontrolslab.com RESEARCH
More informationA Branch and Cut Search for the Deterministic Optimization of the Thermal Unit Commitment Problem. Part II: Computational Results
1 A Branch and Cut Search for the Deterministic Optimization of the Thermal Unit Commitment Problem. Part II: Computational Results Marian G. Marcovecchio, Augusto Q. Novais, Ignacio E. Grossmann Abstract--
More informationA Game-Theoretic Intelligent Agent for the Board Game Football Strategy
A Game-Theoretic Intelligent Agent for the Board Game Football Strategy Sean McCulloch Ohio Wesleyan University stmccull@owu.edu Abstract While much work has been done in the development of intelligent
More information