Piecewise-Linear Approximation: Multidimensional
|
|
- Lawrence Bradford
- 5 years ago
- Views:
Transcription
1 OptIntro 1 / 38 Piecewise-Linear Approximation: Multidimensional Eduardo Camponogara Department of Automation and Systems Engineering Federal University of Santa Catarina October 2016
2 OptIntro 2 / 38 Summary Motivation CC Model DCC Model SOS2 Model
3 OptIntro 3 / 38 Motivation Summary Motivation CC Model DCC Model SOS2 Model
4 OptIntro 4 / 38 Motivation Well Production (Lift-gas, Well-head Pressure).
5 OptIntro 5 / 38 Motivation Pressure Loss in Flowlines
6 OptIntro 5 / 38 Motivation Pressure Loss in Flowlines
7 OptIntro 6 / 38 Motivation Multidimensional Models Let us consider functions with multidimensional domains: f : D R, in which D R n is the domain. An issue that comes up is how to partition the domain. Let P be a set of polyhedrons defining a partitioning of D: P P P = D, int(p) P =, P, P P The piecewise-linear models can be applied depending on the structure of the partitioning P (geometry/topology).
8 OptIntro 6 / 38 Motivation Multidimensional Models Let us consider functions with multidimensional domains: f : D R, in which D R n is the domain. An issue that comes up is how to partition the domain. Let P be a set of polyhedrons defining a partitioning of D: P P P = D, int(p) P =, P, P P The piecewise-linear models can be applied depending on the structure of the partitioning P (geometry/topology).
9 OptIntro 7 / 38 Motivation Models and Partitioning These models can be applied to any convex partitioning of the domain: CC, DCC, DLog e Multiple Choice. Can be applied only to a J-1 triangulation of the domain: Log. Can be applied only to regular lattices: SOS2.
10 OptIntro 7 / 38 Motivation Models and Partitioning These models can be applied to any convex partitioning of the domain: CC, DCC, DLog e Multiple Choice. Can be applied only to a J-1 triangulation of the domain: Log. Can be applied only to regular lattices: SOS2.
11 OptIntro 7 / 38 Motivation Models and Partitioning These models can be applied to any convex partitioning of the domain: CC, DCC, DLog e Multiple Choice. Can be applied only to a J-1 triangulation of the domain: Log. Can be applied only to regular lattices: SOS2.
12 OptIntro 8 / 38 Motivation Models and Partitioning
13 OptIntro 9 / 38 CC Model Summary Motivation CC Model DCC Model SOS2 Model
14 OptIntro 10 / 38 CC Model CC Model Definitions: P is the set of polyhedra defining the domain. V is the set of breakpoints. V (P) V is the set of vertices of polyhedron P P. Variables: λ v [0, 1] is a weighting variable associated with vertex v V. y P {0, 1} assumes value 1 if x P, P P, otherwise assumes value 0.
15 OptIntro 10 / 38 CC Model CC Model Definitions: P is the set of polyhedra defining the domain. V is the set of breakpoints. V (P) V is the set of vertices of polyhedron P P. Variables: λ v [0, 1] is a weighting variable associated with vertex v V. y P {0, 1} assumes value 1 if x P, P P, otherwise assumes value 0.
16 OptIntro 11 / 38 CC Model CC Model f = f (v)λ v, x = vλ v, v V v V 1 = λ v, 1 = y P, v V P P λ v y P, P P(v) λ v 0, v V, y P {0, 1}, P P in which: P(v) = {P P : v V (P)} is the set of polyhedra that contain breakpoint v.
17 OptIntro 12 / 38 CC Model Sample Function Consider the nonlinear function: f (x, y) = e y + e x2 x 2 + xy 2y y x 1 2
18 OptIntro 12 / 38 CC Model Sample Function Consider the nonlinear function: f (x, y) = e y + e x2 x 2 + xy 2y y x 1 2
19 OptIntro 13 / 38 CC Model Example: CC Model Vertex set: V = { 2, 1, 0, 1, 2} { 2, 1, 0, 1, 2}} Polyhedra set: P = { 2, 1, 0, 1} { 2, 1, 0, 1}} in which, the index (i, j) P corresponds to the polyhedron P i,j = {(x, y) R 2 : i x i + 1, j y j + 1}
20 OptIntro 13 / 38 CC Model Example: CC Model Vertex set: V = { 2, 1, 0, 1, 2} { 2, 1, 0, 1, 2}} Polyhedra set: P = { 2, 1, 0, 1} { 2, 1, 0, 1}} in which, the index (i, j) P corresponds to the polyhedron P i,j = {(x, y) R 2 : i x i + 1, j y j + 1}
21 OptIntro 14 / 38 CC Model Example: CC Model f = v V f (v)λ v = f ( 2, 2)λ ( 2, 2) + f ( 2, 1)λ ( 2, 1) + f ( 2, 0)λ ( 2,0). + f ( 2, 1)λ ( 2,1) + f ( 2, 2)λ ( 2,2) + f (2, 2)λ (2, 2) + f (2, 1)λ (2, 1) + f (2, 0)λ (2,0) + f (2, 1)λ (2,1) + f (2, 2)λ (2,2)
22 OptIntro 15 / 38 CC Model Example: CC Model x = v V v xλ v = ( 2)λ ( 2, 2) + ( 2)λ ( 2, 1) + ( 2)λ ( 2,0) + ( 2)λ ( 2,1) + ( 2)λ ( 2,2). + (2)λ (2, 2) + (2)λ (2, 1) + (2)λ (2,0) + (2)λ (2,1) + (2)λ (2,2)
23 OptIntro 16 / 38 CC Model Example: CC Model y = v V v y λ v = ( 2)λ ( 2, 2) + ( 1)λ ( 2, 1) + (0)λ ( 2,0) + (1)λ ( 2,1) + (2)λ ( 2,2). + ( 2)λ (2, 2) + ( 1)λ (2, 1) + (0)λ (2,0) + (1)λ (2,1) + (2)λ (2,2)
24 OptIntro 17 / 38 CC Model Example: CC Model 1 = v V λ v = λ ( 2, 2) + λ ( 2, 1) + λ ( 2,0) + λ ( 2,1) + λ ( 2,2) + λ ( 1, 2) + λ ( 1, 1) + λ ( 1,0) + λ ( 1,1) + λ ( 1,2). + λ (2, 2) + λ (2, 1) + λ (2,0) + λ (2,1) + λ (2,2)
25 OptIntro 18 / 38 CC Model Example: CC Model 1 = P P y P = y ( 2, 2) + y ( 2, 1) + y ( 2,0) + y ( 2,1) + y ( 1, 2) + y ( 1, 1) + y ( 1,0) + y ( 1,1). + y (2, 2) + y (2, 1) + y (2,0) + y (2,1)
26 OptIntro 19 / 38 CC Model Example: CC Model The constraint family λ v y P, P P(v) becomes λ ( 2, 2) y ( 2, 2) λ ( 2, 1) y ( 2, 2) + y ( 2, 1) λ ( 2,0) y ( 2, 1) + y ( 2,0) λ ( 2,1) y ( 2,0) + y ( 2,1) λ ( 2,2) y ( 2,1).. λ ( 1, 2) y ( 2, 2) + y ( 2, 1) λ ( 1, 1) y ( 2, 2) + y ( 2, 1) + y 1, 2 + y ( 1, 1) λ ( 1,0) y ( 2, 1) + y ( 1, 1) + y ( 2,0) + y ( 1,0) λ ( 1,1) y ( 2,0) + y ( 1,0) + y ( 2,1) + y ( 1,1) λ ( 1,2) y ( 2,1) + y ( 1,1)..
27 OptIntro 20 / 38 DCC Model Summary Motivation CC Model DCC Model SOS2 Model
28 OptIntro 21 / 38 DCC Model Disaggregated Convex Combination (DCC) Model f = f (v)λ P,v, P P v V(P) y P = λ P,v, P P, v V(P) 1 = P P y P, x = P P vλ P,v, v V(P) λ P,v 0, P P, v V(P), y P {0, 1}, P P where: λ P,v is associated to vertex v of polyhedron P (breakpoint).
29 OptIntro 22 / 38 DCC Model Example: Function
30 OptIntro 23 / 38 DCC Model Example: Domain y 2 P( 2,1) P( 1,1) P(0,1) P(1,1) x 1 P( 2, 2) P( 1, 2) P( 0, 2) P(1, 2) 2
31 OptIntro 24 / 38 DCC Model Example: DCC Model f = P P v V(P) f (v)λ P,v, = f ( 2, 2)λ P( 2, 2) ( 2, 2) +f ( 1, 2)λ P( 2, 2) ( 1, 2) +f ( 1, 2)λ P( 1, 2) ( 1, 2) +f (0, 2)λ P( 1, 2) (0, 2). +f (1, 1)λ P(1,1) (1,1) + f (1, 2)λ P(1,1) (1,2) + f ( 2, 1)λ P( 2, 2) ( 2, 1) + f ( 1, 1)λ P( 2, 2) ( 1, 1) + f ( 1, 1)λ P( 1, 2) ( 1, 1) + f (0, 1)λ P( 1, 2) (0, 1) +f (2, 1)λ P(1,1) (2,1) + f (2, 2)λ P(1,1) (2,2)
32 OptIntro 25 / 38 DCC Model Example: DCC Model (x, y) = P P v V(P) vλ P,v, = ( 2, 2)λ P( 2, 2) ( 2, 2) +( 1, 2)λ P( 2, 2) ( 1, 2) +( 1, 2)λ P( 1, 2) ( 1, 2) +(0, 2)λ P( 1, 2) (0, 2). +(1, 1)λ P(1,1) (1,1) + (1, 2)λ P(1,1) (1,2) + ( 2, 1)λ P( 2, 2) ( 2, 1) + ( 1, 1)λ P( 2, 2) ( 1, 1) + ( 1, 1)λ P( 1, 2) ( 1, 1) + (0, 1)λ P( 1, 2) (0, 1) +(2, 1)λ P(1,1) (2,1) + (2, 2)λ P(1,1) (2,2)
33 OptIntro 26 / 38 DCC Model Example: DCC Model The constraint family y P = λ P,v, P P, v V(P) becomes y P( 2, 2) = λ P( 2, 2) ( 2, 2) y P( 1, 2) = λ P( 1, 2) ( 1, 2).. + λ P( 2, 2) ( 2, 1) + λ P( 1, 2) ( 1, 1) + λ P( 2, 2) ( 1, 2) + λ P( 1, 2) (0, 2) y P(1,1) = λ P(1,1) (1,1) + λ P(1,1) (1,2) + λ P(1,1) (2,1) + λ P(1,1) (2,2) + λ P( 2, 2) ( 1, 1) + λ P( 1, 2) (0, 1)
34 OptIntro 27 / 38 DCC Model Example: DCC Model The constraint family 1 = P P y P, becomes 1 = y P( 2, 2) + y P( 1, 2) + y P(0, 2) y P(1,1)
35 OptIntro 28 / 38 SOS2 Model Summary Motivation CC Model DCC Model SOS2 Model
36 OptIntro 29 / 38 SOS2 Model SOS2 Model Objectives: Piecewise-linear approximation of the function f : R 2 R. Definitions: Breakpoint set in x : X = {x 1, x 2,..., x n}. Breakpoint set in y : Y = {y 1, y 2,..., y n}. Given the values of the function f for (x, y) X Y. Polyhedron set given by: P = {[x i, x i+1 ] [y j, y j+1 ] : i = 1,..., n 1, j = 1,..., n 1}.
37 OptIntro 29 / 38 SOS2 Model SOS2 Model Objectives: Piecewise-linear approximation of the function f : R 2 R. Definitions: Breakpoint set in x : X = {x 1, x 2,..., x n}. Breakpoint set in y : Y = {y 1, y 2,..., y n}. Given the values of the function f for (x, y) X Y. Polyhedron set given by: P = {[x i, x i+1 ] [y j, y j+1 ] : i = 1,..., n 1, j = 1,..., n 1}.
38 OptIntro 30 / 38 SOS2 Model SOS2 Model f = f (x, y)λ x,y, x X y Y x = xλ x,y, x X y Y y = yλ x,y, x X y Y 1 = λ x,y, λ x,y 0, x X, y Y x X y Y η x = λ x,y, x X, y Y η y = λ x,y, y Y, x X (η x) x X is SOS2 (η y) y Y is SOS2
39 OptIntro 31 / 38 SOS2 Model Example: Function
40 OptIntro 32 / 38 SOS2 Model Example: Domain y x 1 2
41 OptIntro 33 / 38 SOS2 Model Example: SOS2 Model Function approximation becomes f = f (x, y)λ x,y x X y Y = +f ( 2, 2)λ ( 2, 2) + f ( 2, 1)λ ( 2, 1) + f ( 2, 0)λ ( 2,0) +f ( 2, 1)λ ( 2,1) + f ( 2, 2)λ ( 2,2). +f (2, 2)λ (2, 2) + f (2, 1)λ (2, 1) + f (2, 0)λ (2,0) +f (2, 1)λ (2,1) + f (2, 2)λ (2,2)
42 OptIntro 34 / 38 SOS2 Model Example: SOS2 Model Domain variables (x, y) = x X (x, y)λ x,y y Y = +( 2, 2)λ ( 2, 2) + ( 2, 1)λ ( 2, 1) + ( 2, 0)λ ( 2,0) +( 2, 1)λ ( 2,1) + ( 2, 2)λ ( 2,2). +(2, 2)λ (2, 2) + (2, 1)λ (2, 1) + (2, 0)λ (2,0) +(2, 1)λ (2,1) + (2, 2)λ (2,2)
43 OptIntro 35 / 38 SOS2 Model Example: SOS2 Model 1 = x X y Y λ x,y = λ ( 2, 2) + λ ( 2, 1) + λ ( 2,0) + λ ( 2,1) + λ ( 2,2). +λ (2, 2) + λ (2, 1) + λ (2,0) + λ (2,1) + λ (2,2)
44 OptIntro 36 / 38 SOS2 Model SOS2 Model The constraint family η x = y Y λ x,y, x X becomes η x 2 = λ ( 2, 2) + λ ( 2, 1) + λ ( 2,0) + λ ( 2,1) + λ ( 2,2) η x 1 = λ ( 1, 2) + λ ( 1, 1) + λ ( 1,0) + λ ( 1,1) + λ ( 1,2). =. η x 2 = λ (2, 2) + λ (2, 1) + λ (2,0) + λ (2,1) + λ (2,2)
45 OptIntro 37 / 38 SOS2 Model SOS2 Model The constraint family η y = x X λ x,y, y Y becomes η y 2 = λ ( 2, 2) + λ ( 1, 2) + λ (0, 2) + λ (1, 2) + λ (2, 2) η y 1 = λ ( 2, 1) + λ ( 1, 1) + λ (0, 1) + λ (1, 1) + λ (2, 1). =. η y 2 = λ ( 2,2) + λ ( 1,2) + λ (0,2) + λ (1,2) + λ (2,2)
46 OptIntro 38 / 38 SOS2 Model Piecewise-Linear Approximation: Multidimensional End! Thank you for your attention.
Integer Programming: Cutting Planes
OptIntro 1 / 39 Integer Programming: Cutting Planes Eduardo Camponogara Department of Automation and Systems Engineering Federal University of Santa Catarina October 2016 OptIntro 2 / 39 Summary Introduction
More informationMixed Integer Programming Models for Non-Separable Piecewise Linear Cost Functions
Mixed Integer Programming Models for Non-Separable Piecewise Linear Cost Functions Juan Pablo Vielma H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Joint
More informationLecture 10 (Submodular function)
Discrete Methods in Informatics January 16, 2006 Lecture 10 (Submodular function) Lecturer: Satoru Iwata Scribe: Masaru Iwasa and Yukie Nagai Submodular functions are the functions that frequently appear
More informationarxiv: v1 [math.mg] 23 Oct 2014
arxiv:1410.6442v1 [math.mg] 23 Oct 2014 Viviani s Theorem and its Extensions Revisited; Canghareeb, Sanghareeb and a New Definition of the Ellipse October 24, 2014 A two dimensional outsider called Canghareeb
More informationLOCAL CHARACTERIZATION OF POLYHEDRAL SPACES
LOCAL CHARACTERIZATION OF POLYHEDRAL SPACES NINA LEBEDEVA AND ANTON PETRUNIN Abstract. We show that a compact length space is polyhedral if a small spherical neighborhood of any point is conic. arxiv:1402.6670v2
More informationDecidability of consistency of function and derivative information for a triangle and a convex quadrilateral
Decidability of consistency of function and derivative information for a triangle and a convex quadrilateral Abbas Edalat Department of Computing Imperial College London Abstract Given a triangle in the
More informationMarch 2002, December Introduction. We investigate the facial structure of the convex hull of the mixed integer knapsack set
ON THE FACETS OF THE MIXED INTEGER KNAPSACK POLYHEDRON ALPER ATAMTÜRK Abstract. We study the mixed integer knapsack polyhedron, that is, the convex hull of the mixed integer set defined by an arbitrary
More informationLecture Simplex Issues: Number of Pivots. ORIE 6300 Mathematical Programming I October 9, 2014
ORIE 6300 Mathematical Programming I October 9, 2014 Lecturer: David P. Williamson Lecture 14 Scribe: Calvin Wylie 1 Simplex Issues: Number of Pivots Question: How many pivots does the simplex algorithm
More information- 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 informationActive sets, steepest descent, and smooth approximation of functions
Active sets, steepest descent, and smooth approximation of functions Dmitriy Drusvyatskiy School of ORIE, Cornell University Joint work with Alex D. Ioffe (Technion), Martin Larsson (EPFL), and Adrian
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationDiscrete Differential Geometry. Discrete Laplace-Beltrami operator and discrete conformal mappings.
Discrete differential geometry. Discrete Laplace-Beltrami operator and discrete conformal mappings. Technische Universität Berlin Geometric Methods in Classical and Quantum Lattice Systems, Caputh, September
More informationLinear programs, convex polyhedra, extreme points
MVE165/MMG631 Extreme points of convex polyhedra; reformulations; basic feasible solutions; the simplex method Ann-Brith Strömberg 2015 03 27 Linear programs, convex polyhedra, extreme points A linear
More informationMIXED INTEGER PROGRAMMING APPROACHES FOR NONLINEAR AND STOCHASTIC PROGRAMMING
MIXED INTEGER PROGRAMMING APPROACHES FOR NONLINEAR AND STOCHASTIC PROGRAMMING A Thesis Presented to The Academic Faculty by Juan Pablo Vielma Centeno In Partial Fulfillment of the Requirements for the
More informationAmoebas and Instantons
Joint Meeting of Pacific Region Particle Physics Communities (2006, 10/29-11/03, Waikiki) Amoebas and Instantons Toshio Nakatsu Talk based on Takashi Maeda and T.N., Amoebas and Instantons, hep-th/0601233.
More informationIntroduction to Mathematical Programming IE406. Lecture 3. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 3 Dr. Ted Ralphs IE406 Lecture 3 1 Reading for This Lecture Bertsimas 2.1-2.2 IE406 Lecture 3 2 From Last Time Recall the Two Crude Petroleum example.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationComputable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization
Computable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization Heinrich Voss voss@tuhh.de Joint work with Kemal Yildiztekin Hamburg University of Technology Institute
More informationA note on : A Superior Representation Method for Piecewise Linear Functions
A note on : A Superior Representation Method for Piecewise Linear Functions Juan Pablo Vielma Business Analytics and Mathematical Sciences Department, IBM T. J. Watson Research Center, Yorktown Heights,
More informationR-cross-sections of the semigroup of order-preserving transformations of a finite chain
R-cross-sections of the semigroup of order-preserving transformations of a finite chain Eugenija A. Bondar Ural Federal University AAA9+NSAC 07 Novi Sad, Serbia, June -8 Bondar Е.А. R-cross-sections of
More informationOn the Volume of Unbounded Polyhedra in the Hyperbolic Space
Beiträge zur Algebra und eometrie Contributions to Algebra and eometry Volume 44 (2003), No. 1, 145-154. On the Volume of Unbounded Polyhedra in the Hyperbolic Space S. Kántor Institute of Mathematics
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 informationHard-Core Model on Random Graphs
Hard-Core Model on Random Graphs Antar Bandyopadhyay Theoretical Statistics and Mathematics Unit Seminar Theoretical Statistics and Mathematics Unit Indian Statistical Institute, New Delhi Centre New Delhi,
More informationSplit Rank of Triangle and Quadrilateral Inequalities
Split Rank of Triangle and Quadrilateral Inequalities Santanu Dey 1 Quentin Louveaux 2 June 4, 2009 Abstract A simple relaxation of two rows of a simplex tableau is a mixed integer set consisting of two
More informationCOMBINATORICS OF HYPERGEOMETRIC FUNCTIONS ASSOCIATED WITH POSITIVE ROOTS
COMBINATORICS OF HYPERGEOMETRIC FUNCTIONS ASSOCIATED WITH POSITIVE ROOTS Israel M. Gelfand Department of Mathematics Rutgers University New Brunswick, NJ 08903, U.S.A. E-mail: igelfand@math.rutgers.edu
More informationA Branch-and-Refine Method for Nonconvex Mixed-Integer Optimization
A Branch-and-Refine Method for Nonconvex Mixed-Integer Optimization Sven Leyffer 2 Annick Sartenaer 1 Emilie Wanufelle 1 1 University of Namur, Belgium 2 Argonne National Laboratory, USA IMA Workshop,
More informationMa/CS 6a Class 24: More Generating Functions
Ma/CS 6a Class 24: More Generating Functions Road signs from http://www.math.ucla.edu/~pak/ By Adam Sheffer Reminder: Generating Functions Given an infinite sequence of numbers a 0, a 1, a 2, the generating
More informationUSAC Colloquium. Bending Polyhedra. Andrejs Treibergs. September 4, Figure 1: A Rigid Polyhedron. University of Utah
USAC Colloquium Bending Polyhedra Andrejs Treibergs University of Utah September 4, 2013 Figure 1: A Rigid Polyhedron. 2. USAC Lecture: Bending Polyhedra The URL for these Beamer Slides: BendingPolyhedra
More informationAffine Geometry and the Discrete Legendre Transfrom
Affine Geometry and the Discrete Legendre Transfrom Andrey Novoseltsev Department of Mathematics University of Washington April 25, 2008 Andrey Novoseltsev (UW) Affine Geometry April 25, 2008 1 / 24 Outline
More informationPentahedral Volume, Chaos, and Quantum Gravity
Pentahedral Volume, Chaos, and Quantum Gravity Hal Haggard May 30, 2012 Volume Polyhedral Volume (Bianchi, Doná and Speziale): ˆV Pol = The volume of a quantum polyhedron Outline 1 Pentahedral Volume 2
More informationQR Decomposition. When solving an overdetermined system by projection (or a least squares solution) often the following method is used:
(In practice not Gram-Schmidt, but another process Householder Transformations are used.) QR Decomposition When solving an overdetermined system by projection (or a least squares solution) often the following
More informationAffine Geometry and Discrete Legendre Transform
Affine Geometry and Discrete Legendre Transform Andrey Novoseltsev April 24, 2008 Abstract The combinatorial duality used in Gross-Siebert approach to mirror symmetry is the discrete Legendre transform
More informationDiscrete Newton s Algorithm for Parametric Submodular Function Minimization
Discrete Newton s Algorithm for Parametric Submodular Function Minimization Michel X. Goemans, Swati Gupta, and Patrick Jaillet Massachusetts Institute of Technology, USA, goemans@math.mit.edu, swatig@mit.edu,
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 informationEquivariant Perturbation in Gomory and Johnson s Infinite Group Problem. I. The One-Dimensional Case
Equivariant Perturbation in Gomory and Johnson s Infinite Group Problem. I. The One-Dimensional Case Amitabh Basu Robert Hildebrand Matthias Köppe March 25, 204 Abstract We give an algorithm for testing
More informationHW6 Solutions. Micha l Dereziński. March 20, 2015
HW6 Solutions Micha l Dereziński March 20, 2015 1 Exercise 5.5 (a) The PDA accepts odd-length strings whose middle symbol is a and whose other letters are as and bs. Its diagram is below. b, Z 0 /XZ 0
More informationInteger programming, Barvinok s counting algorithm and Gomory relaxations
Integer programming, Barvinok s counting algorithm and Gomory relaxations Jean B. Lasserre LAAS-CNRS, Toulouse, France Abstract We propose an algorithm based on Barvinok s counting algorithm for P max{c
More informationOn the relation between concavity cuts and the surrogate dual for convex maximization problems
On the relation between concavity cuts and the surrogate dual for convex imization problems Marco Locatelli Dipartimento di Ingegneria Informatica, Università di Parma Via G.P. Usberti, 181/A, 43124 Parma,
More informationCombinatorial Reciprocity Theorems
Combinatorial Reciprocity Theorems Matthias Beck San Francisco State University math.sfsu.edu/beck based on joint work with Raman Sanyal Universität Frankfurt JCCA 2018 Sendai Thomas Zaslavsky Binghamton
More information8. Geometric problems
8. Geometric problems Convex Optimization Boyd & Vandenberghe extremal volume ellipsoids centering classification placement and facility location 8 Minimum volume ellipsoid around a set Löwner-John ellipsoid
More informationCombinatorial Harmonic Coordinates
Combinatorial Harmonic Coordinates Saar Hersonsky Fractals 5 - Cornell University. Saar Hersonsky (UGA) June 14, 2014 1 / 15 Perspective Perspective Motivating questions Saar Hersonsky (UGA) June 14, 2014
More informationAntiferromagnetic Potts models and random colorings
Antiferromagnetic Potts models and random colorings of planar graphs. joint with A.D. Sokal (New York) and R. Kotecký (Prague) October 9, 0 Gibbs measures Let G = (V, E) be a finite graph and let S be
More information4 Packing T-joins and T-cuts
4 Packing T-joins and T-cuts Introduction Graft: A graft consists of a connected graph G = (V, E) with a distinguished subset T V where T is even. T-cut: A T -cut of G is an edge-cut C which separates
More informationBarycentric coordinates for Lagrange interpolation over lattices on a simplex
Barycentric coordinates for Lagrange interpolation over lattices on a simplex Gašper Jaklič gasper.jaklic@fmf.uni-lj.si, Jernej Kozak jernej.kozak@fmf.uni-lj.si, Marjeta Krajnc marjetka.krajnc@fmf.uni-lj.si,
More informationEquivalence of the Combinatorial Definition (Lecture 11)
Equivalence of the Combinatorial Definition (Lecture 11) September 26, 2014 Our goal in this lecture is to complete the proof of our first main theorem by proving the following: Theorem 1. The map of simplicial
More informationTropical Elliptic Curves and their j-invariant
and their j-invariant (joint work with Eric Katz and Hannah Markwig) Elliptic Tropical Thomas Markwig Technische Universität Kaiserslautern 15th February, 2008 Elliptic Tropical The tropical j-invariant
More informationDevelopment of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks
Development of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks E. Wanufelle 1 S. Leyffer 2 A. Sartenaer 1 Ph. Toint 1 1 FUNDP, University of Namur 2
More informationCorner Polyhedron and Intersection Cuts
Corner Polyhedron and Intersection Cuts Michele Conforti 1,5, Gérard Cornuéjols 2,4 Giacomo Zambelli 3,5 August 2010 Revised March 2011 Abstract Four decades ago, Gomory introduced the corner polyhedron
More informationConvex Optimization and Modeling
Convex Optimization and Modeling Convex Optimization Fourth lecture, 05.05.2010 Jun.-Prof. Matthias Hein Reminder from last time Convex functions: first-order condition: f(y) f(x) + f x,y x, second-order
More informationarxiv: v2 [math.ag] 24 Jun 2015
TRIANGULATIONS OF MONOTONE FAMILIES I: TWO-DIMENSIONAL FAMILIES arxiv:1402.0460v2 [math.ag] 24 Jun 2015 SAUGATA BASU, ANDREI GABRIELOV, AND NICOLAI VOROBJOV Abstract. Let K R n be a compact definable set
More informationLecture 2: Cubulations
Lecture 2: Cubulations Pocsets CAT(0) Cube Complexes Σ - locally finite, finite width pocset An ultrafilter on Σ is a subset α Σ satisfying Choice: For every A Σ, A α or A* α (not both) Consistency: A
More informationChange in the State of the Art of (Mixed) Integer Programming
Change in the State of the Art of (Mixed) Integer Programming 1977 Vancouver Advanced Research Institute 24 papers 16 reports 1 paper computational, 4 small instances Report on computational aspects: only
More informationVariational methods for lattice systems
Variational methods for lattice systems Andrea Braides Università di Roma Tor Vergata and Mathematical Institute, Oxford Atomistic to Continuum Modelling Workshop Oxford Solid Mechanics November 29 2013
More informationRectilinear Polygon Exchange Transformations via Cut-and-Project Methods
Rectilinear Polygon Exchange Transformations via Cut-and-Project Methods Ren Yi (joint work with Ian Alevy and Richard Kenyon) Zero Entropy System, CIRM 1/39 Outline 1. Polygon Exchange Transformations
More informationFinding parking when not commuting. Jon McCammond U.C. Santa Barbara
Finding parking when not commuting PSfrag replacements 7 8 1 2 6 3 5 {{1, 4, 5}, {2, 3}, {6, 8}, {7}} 4 Jon McCammond U.C. Santa Barbara 1 A common structure The goal of this talk will be to introduce
More informationFinding parking when not commuting. Jon McCammond U.C. Santa Barbara
Finding parking when not commuting 7 6 8 1 2 3 5 4 {{1,4,5}, {2,3}, {6,8}, {7}} Jon McCammond U.C. Santa Barbara 1 A common structure The main goal of this talk will be to introduce you to a mathematical
More informationTHE BANG-BANG PRINCIPLE FOR THE GOURSAT-DARBOUX PROBLEM*
THE BANG-BANG PRINCIPLE FOR THE GOURSAT-DARBOUX PROBLEM* DARIUSZ IDCZAK Abstract. In the paper, the bang-bang principle for a control system connected with a system of linear nonautonomous partial differential
More informationMinimax Problems. Daniel P. Palomar. Hong Kong University of Science and Technolgy (HKUST)
Mini Problems Daniel P. Palomar Hong Kong University of Science and Technolgy (HKUST) ELEC547 - Convex Optimization Fall 2009-10, HKUST, Hong Kong Outline of Lecture Introduction Matrix games Bilinear
More informationA notion of Total Dual Integrality for Convex, Semidefinite and Extended Formulations
A notion of for Convex, Semidefinite and Extended Formulations Marcel de Carli Silva Levent Tunçel April 26, 2018 A vector in R n is integral if each of its components is an integer, A vector in R n is
More informationStructure of Valid Inequalities for Mixed Integer Conic Programs
Structure of Valid Inequalities for Mixed Integer Conic Programs Fatma Kılınç-Karzan Tepper School of Business Carnegie Mellon University 18 th Combinatorial Optimization Workshop Aussois, France January
More informationAdvanced Mixed Integer Programming (MIP) Formulation Techniques
Advanced Mixed Integer Programming (MIP) Formulation Techniques Juan Pablo Vielma Massachusetts Institute of Technology Center for Nonlinear Studies, Los Alamos National Laboratory. Los Alamos, New Mexico,
More informationOperations that preserve the covering property of the lifting region
Operations that preserve the covering property of the lifting region Amitabh Basu and Joe Paat June 23, 2015 Abstract We contribute to the theory for minimal liftings of cut-generating functions. In particular,
More informationTopics Covered in Math 115
Topics Covered in Math 115 Basic Concepts Integer Exponents Use bases and exponents. Evaluate exponential expressions. Apply the product, quotient, and power rules. Polynomial Expressions Perform addition
More informationChapter 2. Vectors and Vector Spaces
2.1. Operations on Vectors 1 Chapter 2. Vectors and Vector Spaces Section 2.1. Operations on Vectors Note. In this section, we define several arithmetic operations on vectors (especially, vector addition
More informationMulti-Row Cuts in Integer Programming. Tepper School of Business Carnegie Mellon University, Pittsburgh
Multi-Row Cuts in Integer Programming Gérard Cornuéjols Tepper School o Business Carnegie Mellon University, Pittsburgh March 2011 Mixed Integer Linear Programming min s.t. cx Ax = b x j Z or j = 1,...,
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationTropicalizations of Positive Parts of Cluster Algebras The conjectures of Fock and Goncharov David Speyer
Tropicalizations of Positive Parts of Cluster Algebras The conjectures of Fock and Goncharov David Speyer This talk is based on arxiv:math/0311245, section 4. We saw before that tropicalizations look like
More informationCalibrated Surrogate Losses
EECS 598: Statistical Learning Theory, Winter 2014 Topic 14 Calibrated Surrogate Losses Lecturer: Clayton Scott Scribe: Efrén Cruz Cortés Disclaimer: These notes have not been subjected to the usual scrutiny
More informationThe Degree of Central Curve in Quadratic Programming
The in Quadratic Programming Mathematics Department San Francisco State University 15 October 2014 Quadratic Programs minimize 1 2 x t Qx + x t c subject to Ax = b x 0 where Q is n n positive definite
More informationCSC 412 (Lecture 4): Undirected Graphical Models
CSC 412 (Lecture 4): Undirected Graphical Models Raquel Urtasun University of Toronto Feb 2, 2016 R Urtasun (UofT) CSC 412 Feb 2, 2016 1 / 37 Today Undirected Graphical Models: Semantics of the graph:
More informationSOME SPECIAL KLEINIAN GROUPS AND THEIR ORBIFOLDS
Proyecciones Vol. 21, N o 1, pp. 21-50, May 2002. Universidad Católica del Norte Antofagasta - Chile SOME SPECIAL KLEINIAN GROUPS AND THEIR ORBIFOLDS RUBÉN HIDALGO Universidad Técnica Federico Santa María
More informationStrong Formulations for Convex Functions over Non-Convex Domains
Strong Formulations for Convex Functions over Non-Convex Domains Daniel Bienstock and Alexander Michalka University JMM 2013 Introduction Generic Problem: min Q(x), s.t. x F, Introduction Generic Problem:
More informationConvex Optimization. Ofer Meshi. Lecture 6: Lower Bounds Constrained Optimization
Convex Optimization Ofer Meshi Lecture 6: Lower Bounds Constrained Optimization Lower Bounds Some upper bounds: #iter μ 2 M #iter 2 M #iter L L μ 2 Oracle/ops GD κ log 1/ε M x # ε L # x # L # ε # με f
More informationLinear Algebra Review: Linear Independence. IE418 Integer Programming. Linear Algebra Review: Subspaces. Linear Algebra Review: Affine Independence
Linear Algebra Review: Linear Independence IE418: Integer Programming Department of Industrial and Systems Engineering Lehigh University 21st March 2005 A finite collection of vectors x 1,..., x k R n
More informationThe partial-fractions method for counting solutions to integral linear systems
The partial-fractions method for counting solutions to integral linear systems Matthias Beck, MSRI www.msri.org/people/members/matthias/ arxiv: math.co/0309332 Vector partition functions A an (m d)-integral
More informationK-Lists. Anindya Sen and Corrie Scalisi. June 11, University of California, Santa Cruz. Anindya Sen and Corrie Scalisi (UCSC) K-Lists 1 / 25
K-Lists Anindya Sen and Corrie Scalisi University of California, Santa Cruz June 11, 2007 Anindya Sen and Corrie Scalisi (UCSC) K-Lists 1 / 25 Outline 1 Introduction 2 Experts 3 Noise-Free Case Deterministic
More information8. Geometric problems
8. Geometric problems Convex Optimization Boyd & Vandenberghe extremal volume ellipsoids centering classification placement and facility location 8 1 Minimum volume ellipsoid around a set Löwner-John ellipsoid
More informationOptimizing Convex Functions over Non-Convex Domains
Optimizing Convex Functions over Non-Convex Domains Daniel Bienstock and Alexander Michalka University Berlin 2012 Introduction Generic Problem: min Q(x), s.t. x F, Introduction Generic Problem: min Q(x),
More informationOn the Relative Strength of Split, Triangle and Quadrilateral Cuts
On the Relative Strength of Split, Triangle and Quadrilateral Cuts Amitabh Basu Pierre Bonami Gérard Cornuéjols François Margot Abstract Integer programs defined by two equations with two free integer
More informationCave Model. Jim Emery 9/10/ A Complete Polyhedron Model of the Cave 4. 5 Clipping the Model With A Plane to Model Water Level 5
Cave Model Jim Emery 9/10/2010 Contents 1 Model 2 2 Basic Hexahedron 2 3 Transition Hexahedron 4 4 A Complete Polyhedron Model of the Cave 4 5 Clipping the Model With A Plane to Model Water Level 5 6 Curl,
More informationarxiv:math/ v1 [math.co] 3 Sep 2000
arxiv:math/0009026v1 [math.co] 3 Sep 2000 Max Min Representation of Piecewise Linear Functions Sergei Ovchinnikov Mathematics Department San Francisco State University San Francisco, CA 94132 sergei@sfsu.edu
More informationLecture 1 Introduction
L. Vandenberghe EE236A (Fall 2013-14) Lecture 1 Introduction course overview linear optimization examples history approximate syllabus basic definitions linear optimization in vector and matrix notation
More informationOn the Relative Strength of Split, Triangle and Quadrilateral Cuts
On the Relative Strength of Split, Triangle and Quadrilateral Cuts Amitabh Basu Tepper School of Business, Carnegie Mellon University, Pittsburgh, PA 53 abasu@andrew.cmu.edu Pierre Bonami LIF, Faculté
More informationA Principled Approach to Mixed Integer/Linear Problem Formulation
A Principled Approach to Mixed Integer/Linear Problem Formulation J N Hooker September 9, 2008 Abstract We view mixed integer/linear problem formulation as a process of identifying disjunctive and knapsack
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationIntroduction to Integer Linear Programming
Lecture 7/12/2006 p. 1/30 Introduction to Integer Linear Programming Leo Liberti, Ruslan Sadykov LIX, École Polytechnique liberti@lix.polytechnique.fr sadykov@lix.polytechnique.fr Lecture 7/12/2006 p.
More information7 Matrix Operations. 7.0 Matrix Multiplication + 3 = 3 = 4
7 Matrix Operations Copyright 017, Gregory G. Smith 9 October 017 The product of two matrices is a sophisticated operations with a wide range of applications. In this chapter, we defined this binary operation,
More informationBilevel Integer Linear Programming
Bilevel Integer Linear Programming TED RALPHS SCOTT DENEGRE ISE Department COR@L Lab Lehigh University ted@lehigh.edu MOPTA 2009, Lehigh University, 19 August 2009 Thanks: Work supported in part by the
More informationBuilding a finite element program in MATLAB Linear elements in 1d and 2d
Building a finite element program in MATLAB Linear elements in 1d and 2d D. Peschka TU Berlin Supplemental material for the course Numerische Mathematik 2 für Ingenieure at the Technical University Berlin,
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationA Geometric Study of the Split Decomposition
A Geometric Study of the Split Decomposition Hiroshi HIRAI Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan hirai@kurims.kyoto-u.ac.jp January 2006 Abstract This paper
More informationCORRECT ANISOTROPIC AND HETEROGENEOUS TRANSMISSIBILITIES CALCULATION IN CORNER-POINT GRIDS
Mecánica Computacional Vol XXIX, págs. 8717-8724 (artículo completo) Eduardo Dvorkin, Marcela Goldschmit, Mario Storti (Eds.) uenos ires, rgentina, 15-18 Noviembre 2010 CORRECT NISOTROPIC ND HETEROGENEOUS
More informationTropical Varieties. Jan Verschelde
Tropical Varieties Jan Verschelde University of Illinois at Chicago Department of Mathematics, Statistics, and Computer Science http://www.math.uic.edu/ jan jan@math.uic.edu Graduate Computational Algebraic
More informationObserver dependent background geometries arxiv:
Observer dependent background geometries arxiv:1403.4005 Manuel Hohmann Laboratory of Theoretical Physics Physics Institute University of Tartu DPG-Tagung Berlin Session MP 4 18. März 2014 Manuel Hohmann
More informationinforms DOI /moor.xxxx.xxxx
MATHEMATICS OF OPERATIONS RESEARCH Vol. 00, No. 0, Xxxxxx 20xx, pp. xxx xxx ISSN 0364-765X EISSN 1526-5471 xx 0000 0xxx informs DOI 10.1287/moor.xxxx.xxxx c 20xx INFORMS Polynomial-Time Approximation Schemes
More informationInteger 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 informationM-convex and L-convex Fuctions over the Real Space Two Conjugate Classes of Combinatorial Convex Functions 1
M-convex and L-convex Fuctions over the Real Space Two Conjugate Classes of Combinatorial Convex Functions 1 Kazuo MUROTA 2 and Akiyoshi SHIOURA 3 July 24, 2002 Abstract By extracting combinatorial structures
More informationMATH8808: ALGEBRAIC TOPOLOGY
MATH8808: ALGEBRAIC TOPOLOGY DAWEI CHEN Contents 1. Underlying Geometric Notions 2 1.1. Homotopy 2 1.2. Cell Complexes 3 1.3. Operations on Cell Complexes 3 1.4. Criteria for Homotopy Equivalence 4 1.5.
More informationOn Strong Unimodality of Multivariate Discrete Distributions
R u t c o r Research R e p o r t On Strong Unimodality of Multivariate Discrete Distributions Ersoy Subasi a Mine Subasi b András Prékopa c RRR 47-2004, DECEMBER, 2004 RUTCOR Rutgers Center for Operations
More informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More information