Relaxing Goodness is Still Good

Size: px
Start display at page:

Download "Relaxing Goodness is Still Good"

Transcription

1 Relaxing Goodness is Still Good Gordon Pace 1 Gerardo Schneider 2 1 Dept. of Computer Science and AI University of Malta 2 Dept. of Informatics University of Oslo ICTAC 08 September 1-3, Istanbul, Turkey Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

2 Reachability Analysis of GSPDIs Hybrid System: combines discrete and continuous dynamics Examples: thermostat, robot, chemical reaction Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

3 Reachability Analysis of GSPDIs Hybrid System: combines discrete and continuous dynamics Examples: thermostat, robot, chemical reaction e 3 R 4 e 2 R 2 e 4 e 5 R 5 R 3 R 7 x f? R 1 e 1 e 8 R 6 x 0 e 6 e 7 R 8 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

4 Outline 1 Polygonal Hybrid Systems (SPDIs) 2 Generalized Polygonal Hybrid Systems (GSPDIs) 3 Reachability Analysis of GSPDIs Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

5 Outline 1 Polygonal Hybrid Systems (SPDIs) 2 Generalized Polygonal Hybrid Systems (GSPDIs) 3 Reachability Analysis of GSPDIs Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

6 Polygonal Hybrid Systems (SPDIs) Preliminaries A constant differential inclusion (angle between vectors a and b): ẋ b a x b a R x Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

7 Polygonal Hybrid Systems (SPDIs) A finite partition of (a subset of) the plane into convex polygonal sets (regions) Dynamics given by the angle determined by two vectors: ẋ b a e 3 R 4 e 2 R 2 e 4 R 5 R 3 R 1 e 1 e 5 R 7 e 8 R 6 e 6 e 7 R 8 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

8 Polygonal Hybrid Systems (SPDIs) A finite partition of (a subset of) the plane into convex polygonal sets (regions) Dynamics given by the angle determined by two vectors: ẋ b a R 4 e 3 x f e 2 R 2 e 4 R 3 e 1 R 5 R 1 e 5 R 7 e 8 R 6 x 0 e 6 e 7 R 8 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

9 Polygonal Hybrid Systems (SPDIs) Goodness Goodness Assumption The dynamics of an SPDI only allows trajectories traversing any edge only in one direction Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

10 Polygonal Hybrid Systems (SPDIs) Goodness Goodness Assumption The dynamics of an SPDI only allows trajectories traversing any edge only in one direction e 4 e 3 e 5 e 5 b P b P e e 2 e 2 6 a a e 1 e 1 Good region e 6 e 4 Bad region e3 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

11 Polygonal Hybrid Systems (SPDIs) Goodness Goodness Assumption The dynamics of an SPDI only allows trajectories traversing any edge only in one direction entry only e 4 exit only e 3 exit only inout e 4 e 5 P e P 5e6 b a e e 2 6 e 1 e 1 Good region b a entry only Bad region e 2 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

12 Polygonal Hybrid Systems (SPDIs) Goodness Goodness Assumption The dynamics of an SPDI only allows trajectories traversing any edge only in one direction entry only e 4 exit only e 3 exit only inout e 4 e 5 P e P 5e6 b a e e 2 6 e 1 e 1 Good region b a entry only Bad region e 2 Theorem Under the goodness assumption, reachability for SPDIs is decidable Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

13 Outline 1 Polygonal Hybrid Systems (SPDIs) 2 Generalized Polygonal Hybrid Systems (GSPDIs) 3 Reachability Analysis of GSPDIs Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

14 GSPDI: Generalized SPDI Definition An SPDI without the goodness assumption is called a GSPDI e 3 R 4 e 2 e 4 R 2 e 1 R 5 R 3 R 7 R 1 e 5 e 8 R 6 e 6 e 7 R 8 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

15 Why Goodness is Good e e R R 2 e 4 e e R R 2 e 1 1 e 4 1 e e R e 1 R 3 R e 4 x 3 x f R 3 f xf R 5 R 1 R 5 R 1 R 5 R 1 R 7 R 7 R 7 R 2 e 1 e 5 e 8 e 5 e 8 e 5 e 8 R 6 x 0 e 6 e 7 R 8 R 6 x 0 e 6 e 7 R 8 R 6 x 0 e 6 e 7 R e 6 e 7 e 8 (e 1 e 2 e 3 e 4 e 5 e 6 e 7 e 8 ) 5 e 9 (e 6 e 7 e 8 e 1 e 2 e 3 e 4 e 5 ) 5 e 6 e 7 e 8 e 9 4 e 6 e 7 e 8 (e 1 e 2 e 3 e 4 e 5 e 6 e 7 e 8 ) e 9 r 1 s 1 r 2 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

16 Why Goodness is Good Lemma (Truncated Affine Multivalued Functions (TAMF)) Successors are positive TAMFs: Succ(x) = F({x} S) J where F(x) = [a 1 x + b 1, a 2 x + b 2 ] with a 1, a 2 > 0 Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

17 Why Goodness is Good Lemma (Truncated Affine Multivalued Functions (TAMF)) Successors are positive TAMFs: Succ(x) = F({x} S) J where F(x) = [a 1 x + b 1, a 2 x + b 2 ] with a 1, a 2 > 0 Theorem An edge-signature σ = e 1... e p can always be abstracted into types of signatures of the form σ A = r 1 s 1... r n s nr n+1, where r i is a sequence of pairwise different edges and all s i are disjoint simple cycle. There are finitely many type of signatures. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

18 Why Goodness is Good Lemma (Truncated Affine Multivalued Functions (TAMF)) Successors are positive TAMFs: Succ(x) = F({x} S) J where F(x) = [a 1 x + b 1, a 2 x + b 2 ] with a 1, a 2 > 0 Theorem An edge-signature σ = e 1... e p can always be abstracted into types of signatures of the form σ A = r 1 s 1... r n s nr n+1, where r i is a sequence of pairwise different edges and all s i are disjoint simple cycle. There are finitely many type of signatures. Many proofs (decidability, soundess, completeness) depend on the goodness assumption Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

19 Problems when Relaxing Goodness Successors are not longer guaranteed to be positive TAMFs Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

20 Problems when Relaxing Goodness Successors are not longer guaranteed to be positive TAMFs Proper inout edges and bounces: e e I 5 I 4 I 1 e I 4 e e e I 3 I 3 I 1 I 0 I 2 I 2 e I 0 e Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

21 Problems when Relaxing Goodness Finiteness argument for types of signature is broken for GSPDIs b c Type of signature: d (abcd) (dcba) (abcd) a a d Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

22 Problems when Relaxing Goodness Finiteness argument for types of signature is broken for GSPDIs c b Type of signature: d (abcd) (dcba) (abcd) a a d Decidability of Reachability for GSPDIs Reduce GSPDI reachability to SPDI reachability (not possible: [SAC 08]) A decidability proof extending that of SPDI Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

23 Outline 1 Polygonal Hybrid Systems (SPDIs) 2 Generalized Polygonal Hybrid Systems (GSPDIs) 3 Reachability Analysis of GSPDIs Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

24 Decidability for GSPDI 1 Prove it is enough to consider trajectories without self-crossing 2 Redefine successors and prove them to be positive TAMFs 3 Duplicate inout edges (e and e 1 ) s.t. each edge is traversed in one direction only 4 Prove that signatures not containing bounces (sub-sequences ee 1 ) behave as for SPDIs 5 Show that we can treat signatures containing bounces 6 Prove soundness and termination Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

25 1. It is enough to consider trajectories without self-crossing Idem as for SPDIs Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

26 2. Successors are positive TAMFs Lemma In GSPDIs, successors can always be written as positive TAMFs. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

27 2. Successors are positive TAMFs Lemma In GSPDIs, successors can always be written as positive TAMFs. We apply a transformation to successors of regions containing inout edges e e b b a a e e Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

28 3. Duplicate inout edges (e and e 1 ) Out Out a b R 2 In In Out Out a b R 2 In In Out Out a b R 2 In In Out In In? R 1 b a In Out Out Out In In e 1 1 e 1 R 1 b a In Out Out Out In In R 1 b a In Out Out In Out In Out In Out (a) (b) (c) Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

29 4. Signatures without bounces behaves as for SPDIs Lemma For any two edges e 0 and e 1, Succ e0 e 1 is always a positive TAMF, whenever e 1 e 1 0. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

30 5. We can treat signatures containing bounces Flip Lemma Let Flip[l, u] = [1 u, 1 l]. Then Succ ee 1 = Flip Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

31 5. We can treat signatures containing bounces Flip Lemma Let Flip[l, u] = [1 u, 1 l]. Then Succ ee 1 = Flip Lemma Composing Flip with an inverted TAMF gives a positive TAMF and an inverted TAMF if we compose it with a positive TAMF. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

32 5. We can treat signatures containing bounces Signatures with even number of bounces Corollary Any signature with an even number of bounces has its behaviour characterised by a positive TAMF. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

33 5. We can treat signatures containing bounces Signatures with even number of bounces Corollary Any signature with an even number of bounces has its behaviour characterised by a positive TAMF. Lemma Given a simple cycle σ containing an even number of bounces, its iterated behaviour can be calculated as for SPDIs. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

34 5. We can treat signatures containing bounces Signatures with odd number of bounces Lemma Given a simple cycle s with an odd number of bounces, we can calculate the limit of its iterated behaviour without iterating. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

35 6. Soundness, Completeness and Termination Algorithm 1 Pre-processing: Redefine successors Transform H: e and e 1 as different edges 2 Generate the finite set of types of signatures Σ = {σ 0,..., σ n } Simple cycles are all distinct 3 Apply Reach σi (x 0, x f ) for each σ i Σ 4 Reach(H, x 0, x f ) = Yes iff for some σ i Σ, Reach σi (x 0, x f ) = Yes. Theorem Reach(H, x 0, x f ) is a sound and complete algorithm calculating GSPDI reachability. The algorithm always terminates. Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

36 Final Remarks This presentation A (constructive) proof that reachability is decidable for GSPDIs Reuse of reachability algorithm for SPDI Acceleration of simple cycles Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

37 Final Remarks This presentation A (constructive) proof that reachability is decidable for GSPDIs Reuse of reachability algorithm for SPDI Acceleration of simple cycles Current Work Implementation of the algorithm (Hallstein A. Hansen) Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

38 Final Remarks This presentation A (constructive) proof that reachability is decidable for GSPDIs Reuse of reachability algorithm for SPDI Acceleration of simple cycles Current Work Implementation of the algorithm (Hallstein A. Hansen) Future Work Application to analysis of nonlinear differential equations Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

39 Motivation Use of SPDIs for approximating nonlinear differential equations Example Pendulum with friction coefficient k, mass M, pendulum length R and gravitational constant g. Behaviour: ẋ = y and ẏ = ky g sin(x) MR 2 R Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

40 Motivation Use of SPDIs for approximating nonlinear differential equations Example Pendulum with friction coefficient k, mass M, pendulum length R and gravitational constant g. Behaviour: ẋ = y and ẏ = ky g sin(x) MR 2 R Triangulation of the plane: Huge number of regions Need to reduce the complexity... without too much overhead Relax Goodness: GSPDI Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

41 Motivation Use of SPDIs for approximating nonlinear differential equations Gerardo Schneider (UiO) Relaxing Goodness is Still Good ICTAC / 27

GSPeeDI a Verification Tool for Generalized Polygonal Hybrid Systems

GSPeeDI a Verification Tool for Generalized Polygonal Hybrid Systems GSPeeDI a Verification Tool for Generalized Polygonal Hybrid Systems Hallstein A. Hansen 1 and Gerardo Schneider 2 1 Buskerud University College, Kongsberg, Norway Hallstein.Asheim.Hansen@hibu.no 2 Dept.

More information

Generell Topologi. Richard Williamson. May 6, 2013

Generell Topologi. Richard Williamson. May 6, 2013 Generell Topologi Richard Williamson May 6, 2013 1 8 Thursday 7th February 8.1 Using connectedness to distinguish between topological spaces I Proposition 8.1. Let (, O ) and (Y, O Y ) be topological spaces.

More information

Automatic synthesis of switching controllers for linear hybrid systems: Reachability control

Automatic synthesis of switching controllers for linear hybrid systems: Reachability control Automatic synthesis of switching controllers for linear hybrid systems: Reachability control Massimo Benerecetti and Marco Faella Università di Napoli Federico II, Italy Abstract. We consider the problem

More information

Cell-Like Maps (Lecture 5)

Cell-Like Maps (Lecture 5) Cell-Like Maps (Lecture 5) September 15, 2014 In the last two lectures, we discussed the notion of a simple homotopy equivalences between finite CW complexes. A priori, the question of whether or not a

More information

Lecture 15: The subspace topology, Closed sets

Lecture 15: The subspace topology, Closed sets Lecture 15: The subspace topology, Closed sets 1 The Subspace Topology Definition 1.1. Let (X, T) be a topological space with topology T. subset of X, the collection If Y is a T Y = {Y U U T} is a topology

More information

Introduction to Rational Billiards II. Talk by John Smillie. August 21, 2007

Introduction to Rational Billiards II. Talk by John Smillie. August 21, 2007 Introduction to Rational Billiards II Talk by John Smillie August 21, 2007 Translation surfaces and their singularities Last time we described the Zemlyakov-Katok construction for billiards on a triangular

More information

Lecture notes for Topology MMA100

Lecture notes for Topology MMA100 Lecture notes for Topology MMA100 J A S, S-11 1 Simplicial Complexes 1.1 Affine independence A collection of points v 0, v 1,..., v n in some Euclidean space R N are affinely independent if the (affine

More information

Principles of AI Planning. Principles of AI Planning. 7.1 How to obtain a heuristic. 7.2 Relaxed planning tasks. 7.1 How to obtain a heuristic

Principles of AI Planning. Principles of AI Planning. 7.1 How to obtain a heuristic. 7.2 Relaxed planning tasks. 7.1 How to obtain a heuristic Principles of AI Planning June 8th, 2010 7. Planning as search: relaxed planning tasks Principles of AI Planning 7. Planning as search: relaxed planning tasks Malte Helmert and Bernhard Nebel 7.1 How to

More information

Integer Programming Theory

Integer Programming Theory Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x

More information

Foundations of Computer Science Spring Mathematical Preliminaries

Foundations of Computer Science Spring Mathematical Preliminaries Foundations of Computer Science Spring 2017 Equivalence Relation, Recursive Definition, and Mathematical Induction Mathematical Preliminaries Mohammad Ashiqur Rahman Department of Computer Science College

More information

Week 7 Convex Hulls in 3D

Week 7 Convex Hulls in 3D 1 Week 7 Convex Hulls in 3D 2 Polyhedra A polyhedron is the natural generalization of a 2D polygon to 3D 3 Closed Polyhedral Surface A closed polyhedral surface is a finite set of interior disjoint polygons

More information

Lecture 1. 1 Notation

Lecture 1. 1 Notation Lecture 1 (The material on mathematical logic is covered in the textbook starting with Chapter 5; however, for the first few lectures, I will be providing some required background topics and will not be

More information

On the number of distinct directions of planes determined by n points in R 3

On the number of distinct directions of planes determined by n points in R 3 On the number of distinct directions of planes determined by n points in R 3 Rom Pinchasi August 27, 2007 Abstract We show that any set of n points in R 3, that is not contained in a plane, determines

More information

Rec 8 Notes. 1 Divide-and-Conquer Convex Hull. 1.1 Algorithm Idea. 1.2 Algorithm. Here s another random Divide-And-Conquer algorithm.

Rec 8 Notes. 1 Divide-and-Conquer Convex Hull. 1.1 Algorithm Idea. 1.2 Algorithm. Here s another random Divide-And-Conquer algorithm. Divide-and-Conquer Convex Hull Here s another random Divide-And-Conquer algorithm. Recall, possibly from CS 2, that the convex hull H(Q) of a set of points Q is the smallest convex set containing Q (that

More information

Pebble Sets in Convex Polygons

Pebble Sets in Convex Polygons 2 1 Pebble Sets in Convex Polygons Kevin Iga, Randall Maddox June 15, 2005 Abstract Lukács and András posed the problem of showing the existence of a set of n 2 points in the interior of a convex n-gon

More information

CLASSIFICATION OF SURFACES

CLASSIFICATION OF SURFACES CLASSIFICATION OF SURFACES JUSTIN HUANG Abstract. We will classify compact, connected surfaces into three classes: the sphere, the connected sum of tori, and the connected sum of projective planes. Contents

More information

However, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example).

However, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example). 98 CHAPTER 3. PROPERTIES OF CONVEX SETS: A GLIMPSE 3.2 Separation Theorems It seems intuitively rather obvious that if A and B are two nonempty disjoint convex sets in A 2, then there is a line, H, separating

More information

Simplified Walking: A New Way to Generate Flexible Biped Patterns

Simplified Walking: A New Way to Generate Flexible Biped Patterns 1 Simplified Walking: A New Way to Generate Flexible Biped Patterns Jinsu Liu 1, Xiaoping Chen 1 and Manuela Veloso 2 1 Computer Science Department, University of Science and Technology of China, Hefei,

More information

Humanoid Robotics. Path Planning and Walking. Maren Bennewitz

Humanoid Robotics. Path Planning and Walking. Maren Bennewitz Humanoid Robotics Path Planning and Walking Maren Bennewitz 1 Introduction Given the robot s pose in a model of the environment Compute a path to a target location First: 2D path in a 2D grid map representation

More information

Lecture 9: Reachability

Lecture 9: Reachability Lecture 9: Reachability Outline of Lecture Reachability General Transition Systems Algorithms for Reachability Safety through Reachability Backward Reachability Algorithm Given hybrid automaton H : set

More information

Optimization of a two-link Robotic Manipulator

Optimization of a two-link Robotic Manipulator Optimization of a two-link Robotic Manipulator Zachary Renwick, Yalım Yıldırım April 22, 2016 Abstract Although robots are used in many processes in research and industry, they are generally not customized

More information

Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams

Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24

More information

Graph Algorithms. Tours in Graphs. Graph Algorithms

Graph Algorithms. Tours in Graphs. Graph Algorithms Graph Algorithms Tours in Graphs Graph Algorithms Special Paths and Cycles in Graphs Euler Path: A path that traverses all the edges of the graph exactly once. Euler Cycle: A cycle that traverses all the

More information

T. Background material: Topology

T. Background material: Topology MATH41071/MATH61071 Algebraic topology Autumn Semester 2017 2018 T. Background material: Topology For convenience this is an overview of basic topological ideas which will be used in the course. This material

More information

Adjacent: Two distinct vertices u, v are adjacent if there is an edge with ends u, v. In this case we let uv denote such an edge.

Adjacent: Two distinct vertices u, v are adjacent if there is an edge with ends u, v. In this case we let uv denote such an edge. 1 Graph Basics What is a graph? Graph: a graph G consists of a set of vertices, denoted V (G), a set of edges, denoted E(G), and a relation called incidence so that each edge is incident with either one

More information

Math 454 Final Exam, Fall 2005

Math 454 Final Exam, Fall 2005 c IIT Dept. Applied Mathematics, December 12, 2005 1 PRINT Last name: Signature: First name: Student ID: Math 454 Final Exam, Fall 2005 I. Examples, Counterexamples and short answer. (6 2 ea.) Do not give

More information

Mathematical and Algorithmic Foundations Linear Programming and Matchings

Mathematical and Algorithmic Foundations Linear Programming and Matchings Adavnced Algorithms Lectures Mathematical and Algorithmic Foundations Linear Programming and Matchings Paul G. Spirakis Department of Computer Science University of Patras and Liverpool Paul G. Spirakis

More information

Parameterization of triangular meshes

Parameterization of triangular meshes Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to

More information

ON THE EMPTY CONVEX PARTITION OF A FINITE SET IN THE PLANE**

ON THE EMPTY CONVEX PARTITION OF A FINITE SET IN THE PLANE** Chin. Ann. of Math. 23B:(2002),87-9. ON THE EMPTY CONVEX PARTITION OF A FINITE SET IN THE PLANE** XU Changqing* DING Ren* Abstract The authors discuss the partition of a finite set of points in the plane

More information

arxiv: v1 [math.co] 27 Feb 2015

arxiv: v1 [math.co] 27 Feb 2015 Mode Poset Probability Polytopes Guido Montúfar 1 and Johannes Rauh 2 arxiv:1503.00572v1 [math.co] 27 Feb 2015 1 Max Planck Institute for Mathematics in the Sciences, Inselstraße 22, 04103 Leipzig, Germany,

More information

FLIP GRAPHS, YOKE GRAPHS AND DIAMETER

FLIP GRAPHS, YOKE GRAPHS AND DIAMETER FLIP GRAPHS, AND Roy H. Jennings Bar-Ilan University, Israel 9th Sèminaire Lotharingien de Combinatoire 0- September, 0 Bertinoro, Italy FLIP GRAPHS, AND FLIP GRAPH TRIANGULATIONS PERMUTATIONS TREES KNOWN

More information

Revisiting the bijection between planar maps and well labeled trees

Revisiting the bijection between planar maps and well labeled trees Revisiting the bijection between planar maps and well labeled trees Daniel Cosmin Porumbel September 1, 2014 Abstract The bijection between planar graphs and well labeled trees was published by Cori and

More information

Chapter 8. Voronoi Diagrams. 8.1 Post Oce Problem

Chapter 8. Voronoi Diagrams. 8.1 Post Oce Problem Chapter 8 Voronoi Diagrams 8.1 Post Oce Problem Suppose there are n post oces p 1,... p n in a city. Someone who is located at a position q within the city would like to know which post oce is closest

More information

THE INSULATION SEQUENCE OF A GRAPH

THE INSULATION SEQUENCE OF A GRAPH THE INSULATION SEQUENCE OF A GRAPH ELENA GRIGORESCU Abstract. In a graph G, a k-insulated set S is a subset of the vertices of G such that every vertex in S is adjacent to at most k vertices in S, and

More information

Delaunay Triangulations

Delaunay Triangulations Delaunay Triangulations (slides mostly by Glenn Eguchi) Motivation: Terrains Set of data points A R 2 Height ƒ(p) defined at each point p in A How can we most naturally approximate height of points not

More information

Jordan Curves. A curve is a subset of IR 2 of the form

Jordan Curves. A curve is a subset of IR 2 of the form Jordan Curves A curve is a subset of IR 2 of the form α = {γ(x) : x [0,1]}, where γ : [0,1] IR 2 is a continuous mapping from the closed interval [0,1] to the plane. γ(0) and γ(1) are called the endpoints

More information

Coloring. Radhika Gupta. Problem 1. What is the chromatic number of the arc graph of a polygonal disc of N sides?

Coloring. Radhika Gupta. Problem 1. What is the chromatic number of the arc graph of a polygonal disc of N sides? Coloring Radhika Gupta 1 Coloring of A N Let A N be the arc graph of a polygonal disc with N sides, N > 4 Problem 1 What is the chromatic number of the arc graph of a polygonal disc of N sides? Or we would

More information

RELATIVELY OPTIMAL CONTROL: THE STATIC SOLUTION

RELATIVELY OPTIMAL CONTROL: THE STATIC SOLUTION RELATIVELY OPTIMAL CONTROL: THE STATIC SOLUTION Franco Blanchini,1 Felice Andrea Pellegrino Dipartimento di Matematica e Informatica Università di Udine via delle Scienze, 208 33100, Udine, Italy blanchini@uniud.it,

More information

STANLEY S SIMPLICIAL POSET CONJECTURE, AFTER M. MASUDA

STANLEY S SIMPLICIAL POSET CONJECTURE, AFTER M. MASUDA Communications in Algebra, 34: 1049 1053, 2006 Copyright Taylor & Francis Group, LLC ISSN: 0092-7872 print/1532-4125 online DOI: 10.1080/00927870500442005 STANLEY S SIMPLICIAL POSET CONJECTURE, AFTER M.

More information

The Borsuk-Ulam theorem- A Combinatorial Proof

The Borsuk-Ulam theorem- A Combinatorial Proof The Borsuk-Ulam theorem- A Combinatorial Proof Shreejit Bandyopadhyay April 14, 2015 1 Introduction The Borsuk-Ulam theorem is perhaps among the results in algebraic topology having the greatest number

More information

CLAN: A Tool for Contract Analysis and Conflict Discovery

CLAN: A Tool for Contract Analysis and Conflict Discovery CLAN: A Tool for Contract Analysis and Conflict Discovery Stephen Fenech 1, Gordon J. Pace 1, and Gerardo Schneider 2 1 Dept. of Computer Science, University of Malta, Malta 2 Dept. of Informatics, University

More information

9.5 Equivalence Relations

9.5 Equivalence Relations 9.5 Equivalence Relations You know from your early study of fractions that each fraction has many equivalent forms. For example, 2, 2 4, 3 6, 2, 3 6, 5 30,... are all different ways to represent the same

More information

Robotics. Chapter 25-b. Chapter 25-b 1

Robotics. Chapter 25-b. Chapter 25-b 1 Robotics Chapter 25-b Chapter 25-b 1 Particle Filtering Particle filtering uses a population of particles (each particle is a state estimate) to localize a robot s position. This is called Monte Carlo

More information

Betweenness and the Crossbar Theorem. Lemma: Let A, B, and C be distinct points. If A*B*C, then neither A*C*B nor B*A*C.

Betweenness and the Crossbar Theorem. Lemma: Let A, B, and C be distinct points. If A*B*C, then neither A*C*B nor B*A*C. Betweenness and the Crossbar Theorem Lemma: Let A, B, and C be distinct points. If A*B*C, then neither A*C*B nor B*A*C. Suppose that both A*B*C and A*C*B. Thus AB+BC =AC, and AC +CB = AB. From this we

More information

Crossing Families. Abstract

Crossing Families. Abstract Crossing Families Boris Aronov 1, Paul Erdős 2, Wayne Goddard 3, Daniel J. Kleitman 3, Michael Klugerman 3, János Pach 2,4, Leonard J. Schulman 3 Abstract Given a set of points in the plane, a crossing

More information

Softwaretechnik. Program verification. Albert-Ludwigs-Universität Freiburg. June 28, Softwaretechnik June 28, / 24

Softwaretechnik. Program verification. Albert-Ludwigs-Universität Freiburg. June 28, Softwaretechnik June 28, / 24 Softwaretechnik Program verification Albert-Ludwigs-Universität Freiburg June 28, 2012 Softwaretechnik June 28, 2012 1 / 24 Road Map Program verification Automatic program verification Programs with loops

More information

The Probabilistic Method

The Probabilistic Method The Probabilistic Method Po-Shen Loh June 2010 1 Warm-up 1. (Russia 1996/4 In the Duma there are 1600 delegates, who have formed 16000 committees of 80 persons each. Prove that one can find two committees

More information

Foundations of Computing

Foundations of Computing Foundations of Computing Darmstadt University of Technology Dept. Computer Science Winter Term 2005 / 2006 Copyright c 2004 by Matthias Müller-Hannemann and Karsten Weihe All rights reserved http://www.algo.informatik.tu-darmstadt.de/

More information

Stable sets, corner polyhedra and the Chvátal closure

Stable sets, corner polyhedra and the Chvátal closure Stable sets, corner polyhedra and the Chvátal closure Manoel Campêlo Departamento de Estatística e Matemática Aplicada, Universidade Federal do Ceará, Brazil, mcampelo@lia.ufc.br. Gérard Cornuéjols Tepper

More information

Network flows and Menger s theorem

Network flows and Menger s theorem Network flows and Menger s theorem Recall... Theorem (max flow, min cut strong duality). Let G be a network. The maximum value of a flow equals the minimum capacity of a cut. We prove this strong duality

More information

Discrete Mathematics Lecture 4. Harper Langston New York University

Discrete Mathematics Lecture 4. Harper Langston New York University Discrete Mathematics Lecture 4 Harper Langston New York University Sequences Sequence is a set of (usually infinite number of) ordered elements: a 1, a 2,, a n, Each individual element a k is called a

More information

1 Undirected Vertex Geography UVG

1 Undirected Vertex Geography UVG Geography Start with a chip sitting on a vertex v of a graph or digraph G. A move consists of moving the chip to a neighbouring vertex. In edge geography, moving the chip from x to y deletes the edge (x,

More information

The Topology of Bendless Orthogonal Three-Dimensional Graph Drawing. David Eppstein Computer Science Dept. Univ. of California, Irvine

The Topology of Bendless Orthogonal Three-Dimensional Graph Drawing. David Eppstein Computer Science Dept. Univ. of California, Irvine The Topology of Bendless Orthogonal Three-Dimensional Graph Drawing David Eppstein Computer Science Dept. Univ. of California, Irvine Graph drawing: visual display of symbolic information Vertices and

More information

Lecture 17: Continuous Functions

Lecture 17: Continuous Functions Lecture 17: Continuous Functions 1 Continuous Functions Let (X, T X ) and (Y, T Y ) be topological spaces. Definition 1.1 (Continuous Function). A function f : X Y is said to be continuous if the inverse

More information

Visibilty: Finding the Staircase Kernel in Orthogonal Polygons

Visibilty: Finding the Staircase Kernel in Orthogonal Polygons American Journal of Computational and Applied Mathematics 2012, 2(2): 17-24 DOI: 10.5923/j.ajcam.20120202.04 Visibilty: Finding the Staircase Kernel in Orthogonal Polygons Stefan A. Pape, Tzvetalin S.

More information

Planar Graphs. 1 Graphs and maps. 1.1 Planarity and duality

Planar Graphs. 1 Graphs and maps. 1.1 Planarity and duality Planar Graphs In the first half of this book, we consider mostly planar graphs and their geometric representations, mostly in the plane. We start with a survey of basic results on planar graphs. This chapter

More information

Lecture 3: Graphs and flows

Lecture 3: Graphs and flows Chapter 3 Lecture 3: Graphs and flows Graphs: a useful combinatorial structure. Definitions: graph, directed and undirected graph, edge as ordered pair, path, cycle, connected graph, strongly connected

More information

The Graphs of Triangulations of Polygons

The Graphs of Triangulations of Polygons The Graphs of Triangulations of Polygons Matthew O Meara Research Experience for Undergraduates Summer 006 Basic Considerations Let Γ(n) be the graph with vertices being the labeled planar triangulation

More information

Number Theory and Graph Theory

Number Theory and Graph Theory 1 Number Theory and Graph Theory Chapter 6 Basic concepts and definitions of graph theory By A. Satyanarayana Reddy Department of Mathematics Shiv Nadar University Uttar Pradesh, India E-mail: satya8118@gmail.com

More information

Matching Theory. Figure 1: Is this graph bipartite?

Matching Theory. Figure 1: Is this graph bipartite? Matching Theory 1 Introduction A matching M of a graph is a subset of E such that no two edges in M share a vertex; edges which have this property are called independent edges. A matching M is said to

More information

CHAPTER 8. Copyright Cengage Learning. All rights reserved.

CHAPTER 8. Copyright Cengage Learning. All rights reserved. CHAPTER 8 RELATIONS Copyright Cengage Learning. All rights reserved. SECTION 8.3 Equivalence Relations Copyright Cengage Learning. All rights reserved. The Relation Induced by a Partition 3 The Relation

More information

Error correction guarantees

Error correction guarantees Error correction guarantees Drawback of asymptotic analyses Valid only as long as the incoming messages are independent. (independence assumption) The messages are independent for l iterations only if

More information

Pantographic polygons

Pantographic polygons 203 Pantographic polygons John Miller and Emanuel Strzelecki Abstract Necessary and sufficient conditions are given for a polygon to be pantographic. The property is held by all regular polygons and by

More information

SFU CMPT Lecture: Week 8

SFU CMPT Lecture: Week 8 SFU CMPT-307 2008-2 1 Lecture: Week 8 SFU CMPT-307 2008-2 Lecture: Week 8 Ján Maňuch E-mail: jmanuch@sfu.ca Lecture on June 24, 2008, 5.30pm-8.20pm SFU CMPT-307 2008-2 2 Lecture: Week 8 Universal hashing

More information

Bayesian Networks, Winter Yoav Haimovitch & Ariel Raviv

Bayesian Networks, Winter Yoav Haimovitch & Ariel Raviv Bayesian Networks, Winter 2009-2010 Yoav Haimovitch & Ariel Raviv 1 Chordal Graph Warm up Theorem 7 Perfect Vertex Elimination Scheme Maximal cliques Tree Bibliography M.C.Golumbic Algorithmic Graph Theory

More information

Hypergraphs With a Unique Perfect Matching

Hypergraphs With a Unique Perfect Matching Hypergraphs With a Unique Perfect Matching Aaron Spindel Under the direction of Dr. John S. Caughman February 26, 2012 Introduction This presentation discusses the paper On the maximum number of edges

More information

Lecture 1 Discrete Geometric Structures

Lecture 1 Discrete Geometric Structures Lecture 1 Discrete Geometric Structures Jean-Daniel Boissonnat Winter School on Computational Geometry and Topology University of Nice Sophia Antipolis January 23-27, 2017 Computational Geometry and Topology

More information

Chapter 3 Numerical Methods

Chapter 3 Numerical Methods Chapter 3 Numerical Methods Part 1 3.1 Linearization and Optimization of Functions of Vectors 1 Problem Notation 2 Outline 3.1.1 Linearization 3.1.2 Optimization of Objective Functions 3.1.3 Constrained

More information

Voronoi diagrams Delaunay Triangulations. Pierre Alliez Inria

Voronoi diagrams Delaunay Triangulations. Pierre Alliez Inria Voronoi diagrams Delaunay Triangulations Pierre Alliez Inria Voronoi Diagram Voronoi Diagram Voronoi Diagram The collection of the non-empty Voronoi regions and their faces, together with their incidence

More information

Convex Optimization CMU-10725

Convex Optimization CMU-10725 Convex Optimization CMU-10725 Conjugate Direction Methods Barnabás Póczos & Ryan Tibshirani Conjugate Direction Methods 2 Books to Read David G. Luenberger, Yinyu Ye: Linear and Nonlinear Programming Nesterov:

More information

CSE 20 DISCRETE MATH WINTER

CSE 20 DISCRETE MATH WINTER CSE 20 DISCRETE MATH WINTER 2016 http://cseweb.ucsd.edu/classes/wi16/cse20-ab/ Today's learning goals Explain the steps in a proof by (strong) mathematical induction Use (strong) mathematical induction

More information

Point-Set Topology II

Point-Set Topology II Point-Set Topology II Charles Staats September 14, 2010 1 More on Quotients Universal Property of Quotients. Let X be a topological space with equivalence relation. Suppose that f : X Y is continuous and

More information

2 Geometry Solutions

2 Geometry Solutions 2 Geometry Solutions jacques@ucsd.edu Here is give problems and solutions in increasing order of difficulty. 2.1 Easier problems Problem 1. What is the minimum number of hyperplanar slices to make a d-dimensional

More information

Delaunay Triangulations. Presented by Glenn Eguchi Computational Geometry October 11, 2001

Delaunay Triangulations. Presented by Glenn Eguchi Computational Geometry October 11, 2001 Delaunay Triangulations Presented by Glenn Eguchi 6.838 Computational Geometry October 11, 2001 Motivation: Terrains Set of data points A R 2 Height ƒ(p) defined at each point p in A How can we most naturally

More information

A NOTE ON PROPER MAPS

A NOTE ON PROPER MAPS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 51, Number 1, August 1975 A NOTE ON PROPER MAPS CHUNG-WU HO1 ABSTRACT. The authot establishes some necessary and sufficient conditions on a Hausdorff

More information

1. Relations 2. Equivalence relations 3. Modular arithmetics. ! Think of relations as directed graphs! xry means there in an edge x!

1. Relations 2. Equivalence relations 3. Modular arithmetics. ! Think of relations as directed graphs! xry means there in an edge x! 2 Today s Topics: CSE 20: Discrete Mathematics for Computer Science Prof. Miles Jones 1. Relations 2. 3. Modular arithmetics 3 4 Relations are graphs! Think of relations as directed graphs! xry means there

More information

Computational Geometry

Computational Geometry Motivation Motivation Polygons and visibility Visibility in polygons Triangulation Proof of the Art gallery theorem Two points in a simple polygon can see each other if their connecting line segment is

More information

QED Q: Why is it called the triangle inequality? A: Analogue with euclidean distance in the plane: picture Defn: Minimum Distance of a code C:

QED Q: Why is it called the triangle inequality? A: Analogue with euclidean distance in the plane: picture Defn: Minimum Distance of a code C: Lecture 3: Lecture notes posted online each week. Recall Defn Hamming distance: for words x = x 1... x n, y = y 1... y n of the same length over the same alphabet, d(x, y) = {1 i n : x i y i } i.e., d(x,

More information

6 Randomized rounding of semidefinite programs

6 Randomized rounding of semidefinite programs 6 Randomized rounding of semidefinite programs We now turn to a new tool which gives substantially improved performance guarantees for some problems We now show how nonlinear programming relaxations can

More information

Simplicial Hyperbolic Surfaces

Simplicial Hyperbolic Surfaces Simplicial Hyperbolic Surfaces Talk by Ken Bromberg August 21, 2007 1-Lipschitz Surfaces- In this lecture we will discuss geometrically meaningful ways of mapping a surface S into a hyperbolic manifold

More information

Polygon Triangulation. (slides partially by Daniel Vlasic )

Polygon Triangulation. (slides partially by Daniel Vlasic ) Polygon Triangulation (slides partially by Daniel Vlasic ) Triangulation: Definition Triangulation of a simple polygon P: decomposition of P into triangles by a maximal set of non-intersecting diagonals

More information

CMPUT 412 Motion Control Wheeled robots. Csaba Szepesvári University of Alberta

CMPUT 412 Motion Control Wheeled robots. Csaba Szepesvári University of Alberta CMPUT 412 Motion Control Wheeled robots Csaba Szepesvári University of Alberta 1 Motion Control (wheeled robots) Requirements Kinematic/dynamic model of the robot Model of the interaction between the wheel

More information

The Touring Polygons Problem (TPP)

The Touring Polygons Problem (TPP) The Touring Polygons Problem (TPP) [Dror-Efrat-Lubiw-M]: Given a sequence of k polygons in the plane, a start point s, and a target point, t, we seek a shortest path that starts at s, visits in order each

More information

Paths, Flowers and Vertex Cover

Paths, Flowers and Vertex Cover Paths, Flowers and Vertex Cover Venkatesh Raman M. S. Ramanujan Saket Saurabh Abstract It is well known that in a bipartite (and more generally in a König) graph, the size of the minimum vertex cover is

More information

Convex Optimization Lecture 2

Convex Optimization Lecture 2 Convex Optimization Lecture 2 Today: Convex Analysis Center-of-mass Algorithm 1 Convex Analysis Convex Sets Definition: A set C R n is convex if for all x, y C and all 0 λ 1, λx + (1 λ)y C Operations that

More information

2009 HMMT Team Round. Writing proofs. Misha Lavrov. ARML Practice 3/2/2014

2009 HMMT Team Round. Writing proofs. Misha Lavrov. ARML Practice 3/2/2014 Writing proofs Misha Lavrov ARML Practice 3/2/2014 Warm-up / Review 1 (From my research) If x n = 2 1 x n 1 for n 2, solve for x n in terms of x 1. (For a more concrete problem, set x 1 = 2.) 2 (From this

More information

THREE LECTURES ON BASIC TOPOLOGY. 1. Basic notions.

THREE LECTURES ON BASIC TOPOLOGY. 1. Basic notions. THREE LECTURES ON BASIC TOPOLOGY PHILIP FOTH 1. Basic notions. Let X be a set. To make a topological space out of X, one must specify a collection T of subsets of X, which are said to be open subsets of

More information

5. THE ISOPERIMETRIC PROBLEM

5. THE ISOPERIMETRIC PROBLEM Math 501 - Differential Geometry Herman Gluck March 1, 2012 5. THE ISOPERIMETRIC PROBLEM Theorem. Let C be a simple closed curve in the plane with length L and bounding a region of area A. Then L 2 4 A,

More information

Lesson 19: The Graph of a Linear Equation in Two Variables is a Line

Lesson 19: The Graph of a Linear Equation in Two Variables is a Line Lesson 19: The Graph of a Linear Equation in Two Variables is a Line Classwork Exercises Theorem: The graph of a linear equation y = mx + b is a non-vertical line with slope m and passing through (0, b),

More information

PCP and Hardness of Approximation

PCP and Hardness of Approximation PCP and Hardness of Approximation January 30, 2009 Our goal herein is to define and prove basic concepts regarding hardness of approximation. We will state but obviously not prove a PCP theorem as a starting

More information

Final Test in MAT 410: Introduction to Topology Answers to the Test Questions

Final Test in MAT 410: Introduction to Topology Answers to the Test Questions Final Test in MAT 410: Introduction to Topology Answers to the Test Questions Stefan Kohl Question 1: Give the definition of a topological space. (3 credits) A topological space (X, τ) is a pair consisting

More information

1 Matching in Non-Bipartite Graphs

1 Matching in Non-Bipartite Graphs CS 369P: Polyhedral techniques in combinatorial optimization Instructor: Jan Vondrák Lecture date: September 30, 2010 Scribe: David Tobin 1 Matching in Non-Bipartite Graphs There are several differences

More information

One-to-One Piecewise Linear Mappings Over Triangulations

One-to-One Piecewise Linear Mappings Over Triangulations One-to-One Piecewise Linear Mappings Over Triangulations Michael S. Floater Presentation by Elif Tosun Geometric Modeling Fall 02 Outline Introduction Triangulations Convex Combination Mappings Proof:

More information

Vesa Halava Tero Harju. Walks on Borders of Polygons

Vesa Halava Tero Harju. Walks on Borders of Polygons Vesa Halava Tero Harju Walks on Borders of Polygons TUCS Technical Report No 698, June 2005 Walks on Borders of Polygons Vesa Halava Tero Harju Department of Mathematics and TUCS - Turku Centre for Computer

More information

Math 210 Manifold III, Spring 2018 Euler Characteristics of Surfaces Hirotaka Tamanoi

Math 210 Manifold III, Spring 2018 Euler Characteristics of Surfaces Hirotaka Tamanoi Math 210 Manifold III, Spring 2018 Euler Characteristics of Surfaces Hirotaka Tamanoi 1. Euler Characteristic of Surfaces Leonhard Euler noticed that the number v of vertices, the number e of edges and

More information

arxiv: v1 [math.na] 20 Sep 2016

arxiv: v1 [math.na] 20 Sep 2016 arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.

More information

Let and be a differentiable function. Let Then be the level surface given by

Let and be a differentiable function. Let Then be the level surface given by Module 12 : Total differential, Tangent planes and normals Lecture 35 : Tangent plane and normal [Section 35.1] > Objectives In this section you will learn the following : The notion tangent plane to a

More information

Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves

Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves Csaba Attila Vigh, Dr. Michael Bader Department of Informatics, TU München JASS 2006, course 2: Numerical Simulation:

More information

ON INDEX EXPECTATION AND CURVATURE FOR NETWORKS

ON INDEX EXPECTATION AND CURVATURE FOR NETWORKS ON INDEX EXPECTATION AND CURVATURE FOR NETWORKS OLIVER KNILL Abstract. We prove that the expectation value of the index function i f (x) over a probability space of injective function f on any finite simple

More information

Evolutionary Computation. Chao Lan

Evolutionary Computation. Chao Lan Evolutionary Computation Chao Lan Outline Introduction Genetic Algorithm Evolutionary Strategy Genetic Programming Introduction Evolutionary strategy can jointly optimize multiple variables. - e.g., max

More information