Convex hulls of spheres and convex hulls of convex polytopes lying on parallel hyperplanes
|
|
- Francis Stephens
- 6 years ago
- Views:
Transcription
1 Convex hulls of spheres and convex hulls of convex polytopes lying on parallel hyperplanes Menelaos I. Karavelas joint work with Eleni Tzanaki University of Crete & FO.R.T.H. OrbiCG/ Workshop on Computational Geometry INRIA Sophia-Antipolis Méditerranée, December 9, 2010 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 1 / 42
2 1 Introduction OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 2 / 42
3 Problem setup Introduction Problem setup Previous work & problem history Our results We are given a set Σ of n spheres in E d, such that n i spheres have radius ρ i, where 1 i m and 1 m n. The radii are pairwise distinct, i.e., ρ i ρ j for i j. The dimension d is considered fixed. Problem What is the worst-case combinatorial complexity of the convex hull CH d (Σ) of Σ, when m is fixed? Was posed as an open problem by [Boissonnat & K. 2003] The problem is interesting only for odd dimensions. Throughout the talk d 3, and d odd. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 3 / 42
4 Problem setup Previous work & problem history Our results Previous/Related work Worst-case complexities For a point set P in E d the worst-case complexity of CH d (P ) is Θ(n d 2 ). Same for spheres with same radius. [Aurenhammer 1987]: The worst-case complexity of CH d (Σ) is O(n d 2 ). Reduction to power diagram in E d+1. [Boissonnat et al. 1996]: Showed that CH 3 (Σ) = Ω(n 2 ). [Boissonnat & K. 2003]: Showed that CH d (Σ) = Ω(n d 2 ) for all d 3. Correspondence between Möbius diagrams in E d 1 with additively weighted Voronoi cells in E d and convex hulls of spheres in E d (via inversions). Worst-case bound for Möbius diagrams in E d 1 is Θ(n d 2 ). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 4 / 42
5 Problem setup Previous work & problem history Our results Previous/Related work Worst-case optimal algorithms [Chazelle 1993]: The convex hull of n points in E d, d 2, can be computed in worst-case optimal O(n d 2 + n log n) time. Worst-case optimal algorithms already existed for d = 2, 3. [Boissonnat et al. 1996]: Presented a O(n d 2 + n log n) time algorithm: worst-case optimal only for even dimensions (at that point). Lifting map from spheres in E d to points in E d+1, using the radius as the last coordinate. [Boissonnat & K. 2003]: Due to the lower bound of Ω(n d 2 ) for CH d (Σ), the algorithm in [Boissonnat et al. 1996] is actually worst-case optimal for all d 2. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 5 / 42
6 Problem setup Previous work & problem history Our results Previous/Related work Output-sensitive algorithms Many output-sensitive algorithms for the convex hull of points. We mention those related to this presentation: [Seidel 1986]: Uses the notion of polytopal shellings to compute CH d (P ) in O(n 2 + f log n) time. [Matoušek & Schwarzkopf 1992]: Improved the n 2 part of Seidel s algorithm to get O(n 2 2/( d 2 +1)+ɛ + f log n). [Chan, Snoeyink & Yap 1997]: Works in 4D and has running time O((n + f) log 2 f). Many more algorithms for 2D and 3D, some of which are optimal (in the output-sensitive sense). Very few output-sensitive algorithms for spheres: [Nielsen & Yvinec 1998]: Computes CH 2 (Σ) in optimal O(n log f) time. [Boissonnat, Cérézo & Duquesne 1992]: Gift-wrapping algorithm for computing CH 3 (Σ); runs in O(nf) time. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 6 / 42
7 Our results Some definitions Problem setup Previous work & problem history Our results We have a set of spheres Σ consisting of n spheres, such that n i n spheres have radius ρ i, where 1 i m and m 2 and m is fixed. Definition We say that ρ λ dominates Σ if n λ = Θ(n). Definition We say that Σ is uniquely dominated if, for some λ, n λ = Θ(n), and n i = o(n) for all i λ. Definition We say that Σ is strongly dominated if, for some λ, n λ = Θ(n), and n i = O(1)) for all i λ. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 7 / 42
8 Problem setup Previous work & problem history Our results Our results Qualitative point-of-view We use the term generic worst-case complexity to refer to the worst-case complexity of CH d (Σ) where there is no restriction on the number of distinct radii in Σ. Result If Σ is dominated by at least two radii, the worst-case complexity of CH d (Σ) matches the generic worst-case complexity. Result If Σ is uniquely dominated, the worst-case complexity of CH d (Σ) is asymptotically smaller than the generic worst-case complexity. Result If Σ is strongly dominated, the worst-case complexity of CH d (Σ) matches the worst-case complexity of convex hulls of points. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 8 / 42
9 Problem setup Previous work & problem history Our results Our results Quantitative point-of-view Theorem The worst-case complexity of CH d (Σ), for m fixed, is Θ( 1 i j m n in d 2 j ). Result If Σ is dominated by at least two radii, the complexity of CH d (Σ) is Θ(n d 2 ). Result If Σ is uniquely dominated, the complexity of CH d (Σ) is o(n d 2 ). Result If Σ is strongly dominated, the complexity of CH d (Σ) is Θ(n d 2 ). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 9 / 42
10 Our results Methodology Problem setup Previous work & problem history Our results Theorem For the upper bound we reduce the sphere convex hull problem to the problem of computing the complexity of the convex hull of m d-polytopes lying on m parallel hyperplanes of E d+1. Let P be a set of m d-polytopes {P 1, P 2,..., P m } lying on m parallel hyperplanes of E d+1. The worst-case complexity of CH d+1 (P) is O( 1 i j m n in d 2 j ), where n i = f 0 (P i ), 1 i m. For the lower bound we first construct a set Σ of Θ(n 1 + n 2 ) spheres in E d such that the complexity of CH d (Σ) is Ω(n 1 n d n 2 n d 2 1 ). Then we generalize this construction for m 3. This construction also gives a matching lower bound for the parallel polytope convex hull problem. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 10 / 42
11 Definitions Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Polytope P: the convex hull of a set of points P Face of P: intersection of P with at least one supporting hyperplane k-face: k-dimensional face trivial face: the unique face of dimension d (P) proper faces: faces of dimension at most d 1 d-polytope: a polytope whose trivial face is d-dimensional vertices: 0-faces; edges: 1-faces facets: (d 1)-faces; ridges: (d 2)-faces; simplicial polytope: all proper faces are simplices f k (P): number of k-faces of P f-vector: (f 1 (P), f 0 (P),..., f d 1 (P)) f 1 (P) = 1 (empty set) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 11 / 42
12 Definitions (cont.) Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound For simplicial d-polytopes we can define the h-vector: (h 0 (P), h 1 (P),..., h d (P)), where h k (P) = k ( ) d i ( 1) k i f i 1 (P). d k i=0 h k (P): number of facets in a shelling of P whose restriction has size k. The elements of the f-vector determine the h-vector and vice versa. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 12 / 42
13 Dehn-Sommerville equations Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound The element of the f-vector are linearly dependent. They satisfy the so called Dehn-Sommerville equations: h k (P) = h d k (P), 0 k d. Important implication: if we know the face numbers f k (P), 0 k d 2 1, we can determine the remaining face numbers f k (P), d 2 k d 1. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 13 / 42
14 Simplicial vs. non-simplicial polytopes Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound For any non-simplicial d-polytope there exists a (nearby) simplicial polytope with the same number of vertices, and with at least as many faces as the non-simplicial polytope: Theorem ([Klee 1964],[McMullen 1970]) Let P be a d-polytope. 1 The d-polytope P we obtain by pulling a vertex of P has the same number of vertices with P, and f k (P) f k (P ) for all 1 k d 1. 2 The d-polytope P we obtain by successively pulling each of the vertices of P is simplicial, has the same number of vertices with P, and f k (P) f k (P ) for all 1 k d 1. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 14 / 42
15 Adaptation to parallel polytopes Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Lemma Let P = {P 1, P 2,..., P m } be a set of m 2 d-polytopes lying on m parallel hyperplanes Π 1, Π 2,..., Π m of E d+1, respectively, where Π j is above Π i for all j > i. Let P i be the vertex set of P i, 1 i m, P = P 1 P 2... P m, and P = CH d+1 (P ). The points in P can be perturbed in such a way that: 1 the points of P i remain in Π i, 1 i m, 2 all the faces of P, except possibly the facets P 1 and P m, are simplices, and, 3 f k (P) f k (P ) for all 1 k d, where P is the polytope we obtain after having perturbed the vertices of P in P. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 15 / 42
16 Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Adaptation to parallel polytopes Sketch of Proof First we pull the vertices of P 1 so that P 1 becomes simplicial (we consider P 1 embedded in E d ). Secondly we pull the vertices of P m so that P m becomes simplicial (we consider P m embedded in E d ). Finally, for each i, 2 i m 1, we pull the vertices of P i in P, so that all faces of P, except P 1 and P m, become simplicial. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 16 / 42
17 Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Adaptation to parallel polytopes Sketch of Proof First we pull the vertices of P 1 so that P 1 becomes simplicial (we consider P 1 embedded in E d ). Secondly we pull the vertices of P m so that P m becomes simplicial (we consider P m embedded in E d ). Finally, for each i, 2 i m 1, we pull the vertices of P i in P, so that all faces of P, except P 1 and P m, become simplicial. At each step above the number of vertices of the new polytope is the same as the old polytope. At each step above the number of faces of the new polytope is at least as big as the number of faces of the old polytope. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 16 / 42
18 The setting Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound P P4 F By means of the lemma, we need only consider polytopes P such that P is simplicial, except possibly for its two facets P 1 and P m. Π: hyperplane parallel and between Π m 1 and Π m. P3 F: set of faces of P with non-empty intersection with Π. L P2 For each k-face F F, at least one point comes from P m, whereas the remaining k points come from P 1,..., P m 1. P1 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 17 / 42
19 The easy upper bound Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Let A = (α 1..., α m), B = (β 1,..., β m) N m. Define: A = P m i=1 αi, and A B, iff α i β i, 1 i m. Theorem The number of k-faces of F is bounded from above as follows:! X my f 0(P i) f k (F), 1 k d. α i Proof. (0,...,0,1) A (k,...,k) A =k+1 A k-face F F is defined by k + 1 vertices of P, where at least one vertex comes from P m, whereas the remaining vertices are vertices of P 1,..., P m 1. Let α i be the number of vertices of F from P i. We have: 0 α i k for 1 i m 1, 1 α m k, and P m i=1 αi = k + 1. The maximum number of possible (α i 1)-faces of P i is `f 0 (P i ). Hence, the maximum possible number of k-faces of F is Q m i=1 α i we get the desired expression. i=1 α i `f0 (P i ). Summing over all possible values for the αi s OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 18 / 42
20 The easy upper bound Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Let A = (α 1..., α m), B = (β 1,..., β m) N m. Define: A = P m i=1 αi, and A B, iff α i β i, 1 i m. Theorem The number of k-faces of F is bounded from above as follows:! X my f 0(P i) f k (F), 1 k d. α i Corollary (0,...,0,1) A (k,...,k) A =k+1 Let n i = f 0(P i), 1 i m. The following asymptotic bounds hold: f k (F) = O(n k m m 1 X i=1 i=1 m 1 X n i + n m i=1 n k i ), 1 k d 2. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 18 / 42
21 Defining an auxilliary polytope Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Call C the polytopal complex with facets the facets of F. Call L the faces of P not in P m or F. Call L = C L. Call P m the boundary complex of P m. Define the point set Q so as to consist of the points of P m, L and two additional points y and z. y is below Π 1 and visible by the vertices of P 1 only. z is below Π m and visible by the vertices of P m only. Define Q = CH d+1 (Q); Q is simplicial (d + 1)-polytope. L is the link of y in Q. P m is the link of z in Q. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 19 / 42
22 Relations on faces Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound 1 For a vertex v Q, call S v its star. Then: f k (Q) = f k (F) + f k (S y ) + f k (S z ), 0 k d, where f 0 (F) = 0. 2 For the faces of S z and P m, we have for 0 k d: f k (S z ) = f k ( P m ) + f k 1 ( P m ) where f 1 ( P m ) = 1 and f d ( P m ) = 0. 3 For the faces of S y and L, we have for 0 k d: f k (S y ) = f k ( L)+f k 1 ( L) where f 1 ( L) = 1 and f d ( L) = 0. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 20 / 42
23 Relations on faces Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound 1 For a vertex v Q, call S v its star. Then: f k (Q) = f k (F) + f k (S y ) + f k (S z ), 0 k d, where f 0 (F) = 0. 2 For the faces of S z and P m, we have for 0 k d: f k (S z ) = f k ( P m ) + f k 1 ( P m ) = O(n d 2 m ) where f 1 ( P m ) = 1 and f d ( P m ) = 0. 3 For the faces of S y and L, we have for 0 k d: f k (S y ) = f k ( L)+f k 1 ( L) where f 1 ( L) = 1 and f d ( L) = 0. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 20 / 42
24 Relations on faces Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound 1 For a vertex v Q, call S v its star. Then: f k (Q) = f k (F) + f k (S y ) + f k (S z ), 0 k d, where f 0 (F) = 0. 2 For the faces of S z and P m, we have for 0 k d: f k (S z ) = f k ( P m ) + f k 1 ( P m ) = O(n d 2 m ) where f 1 ( P m ) = 1 and f d ( P m ) = 0. 3 For the faces of S y and L, we have for 0 k d: m 1 f k (S y ) = f k ( L)+f k 1 ( L) = O(( n i ) d 2 ) where f 1 ( L) = 1 and f d ( L) = 0. i=1 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 20 / 42
25 Relations on faces Introduction Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound 1 For a vertex v Q, call S v its star. Then: f k (Q) = f k (F) + f k (S y ) + f k (S z ), 0 k d, where f 0 (F) = 0. 2 For the faces of S z and P m, we have for 0 k d: f k (S z ) = f k ( P m ) + f k 1 ( P m ) = O(n d 2 m ) where f 1 ( P m ) = 1 and f d ( P m ) = 0. 3 For the faces of S y and L, we have for 0 k d: m 1 f k (S y ) = f k ( L)+f k 1 ( L) = O(( where f 1 ( L) = 1 and f d ( L) = 0. i=1 m 1 n i ) d 2 ) = O( i=1 n d 2 i ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 20 / 42
26 Using the Dehn-Sommerville equations Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Combining the bounds for f k (F), f k (S y ) and f k (S z ), for 0 k d 2, we get: f k (Q) = O(n d 2 m 1 m i=1 m 1 n i + n m i=1 n d 2 i ), 0 k d 2. Hence for the elements of the h-vector we have: k ( ) d + 1 i h k (Q) = ( 1) k i f i 1 (Q), 0 k d + 1 d + 1 k i=0 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 21 / 42
27 Using the Dehn-Sommerville equations Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Combining the bounds for f k (F), f k (S y ) and f k (S z ), for 0 k d 2, we get: f k (Q) = O(n d 2 m 1 m i=1 m 1 n i + n m i=1 n d 2 i ), 0 k d 2. Hence for the elements of the h-vector we have: k ( ) d + 1 i h k (Q) = ( 1) k i f i 1 (Q), 0 k d + 1 d + 1 k i=0 = O(n d 2 m 1 m i=1 m 1 n i + n m i=1 n d 2 i ), 0 k d+1 2 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 21 / 42
28 Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Using the Dehn-Sommerville equations (cont.) Using the Dehn-Sommerville equations we get for d+1 2 k d + 1: h k (Q) = h d+1 k (Q) = O(n d 2 m 1 m i=1 m 1 n i + n m i=1 Writing the f-vector in terms of the h-vector we have: k ( ) d + 1 i f k 1 (Q) = h i (Q), 0 k d k i i=0 n d 2 i ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 22 / 42
29 Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound Using the Dehn-Sommerville equations (cont.) Using the Dehn-Sommerville equations we get for d+1 2 k d + 1: h k (Q) = h d+1 k (Q) = O(n d 2 m 1 m i=1 m 1 n i + n m i=1 Writing the f-vector in terms of the h-vector we have: k ( ) d + 1 i f k 1 (Q) = h i (Q), 0 k d k i i=0 = O(n d 2 m 1 m i=1 m 1 n i + n m i=1 n d 2 i ) n d 2 i ), 0 k d. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 22 / 42
30 The complexity of CH d+1 (P) Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound P P4 F Theorem Let P = {P 1, P 2,..., P m } be a set of a fixed number of m 2 parallel d-polytopes, where d 3 and d is odd. The worst-case complexity of CH d+1 (P) is P3 L P2 O( 1 i j m n i n d 2 j ), P1 where n i = f 0 (P i ), 1 i m. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 23 / 42
31 Definitions & preliminaries Upper bound on number of faces cut by a hyperplane Inductive proof for upper bound The complexity of CH d+1 (P) Inductive proof P4 Proof. P3 L P P2 F Let T (m) be the worst-case complexity of CH d+1 (P). The set of faces of P is the disjoint union of L, F and the set of faces of P m. The faces in L are also faces of CH d+1 (P \ {P m}), which implies that the complexity of L is at most T (m 1). Therefore, we have: T (m) T (m 1) + O(n d 2 m ) T (m) c( + O(n d 2 m 1 X m i=1 X 1 i j m m 1 X n i + n m n i n d 2 j i=1 mx + i=1 n d 2 i ) n d 2 i ) P1 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 24 / 42
32 Context Introduction Upper bound Lower bound for two radii Lower bound for at least three radii Σ: set of n spheres in E d The radii of the spheres in Σ take m distinct values: ρ 1 < ρ 2 <... < ρ m CH d (Σ): convex hull of Σ Face of circularity l of CH d (Σ): maximal connected portion of the boundary of CH d (Σ), where the supporting hyperplanes are tangent to a given set of (d l) spheres n i : number of spheres in Σ with radius ρ i OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 25 / 42
33 Lifting map Introduction Upper bound Lower bound for two radii Lower bound for at least three radii Let H 0 be the hyperplane {x d+1 = 0} of E d+1. Map the sphere σ k = (c k, ρ k ) to the point p k = (c k, ρ k ) E d+1. Let P = {p 1, p 2,..., p n }. Call P = CH d+1 (P ). Let λ 0 be the lower halfcone in E d+1 with arbitrary apex, vertical axis, and apex angle equal to π 4. Let λ(p) be the translated copy of λ 0 with apex at p. Define the set Λ to be the set Λ = {λ(p 1 ), λ(p 2 ),..., λ(p n )} OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 26 / 42
34 Lifting map Introduction Upper bound Lower bound for two radii Lower bound for at least three radii Let H 0 be the hyperplane {x d+1 = 0} of E d+1. Map the sphere σ k = (c k, ρ k ) to the point p k = (c k, ρ k ) E d+1. Let P = {p 1, p 2,..., p n }. Call P = CH d+1 (P ). Let λ 0 be the lower halfcone in E d+1 with arbitrary apex, vertical axis, and apex angle equal to π 4. Let λ(p) be the translated copy of λ 0 with apex at p. Define the set Λ to be the set Λ = {λ(p 1 ), λ(p 2 ),..., λ(p n )} Fact λ(p k ) H 0 = σ k, 1 k n Fact CH d+1 (Λ) H 0 = CH d (Σ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 26 / 42
35 Upper bound Lower bound for two radii Lower bound for at least three radii Relation between faces of CH d+1 (Λ) and CH d (Σ) Let O be any point inside P. Theorem ([Boissonnat et al. 1996]) Any hyperplane of E d supporting CH d (Σ) is the intersection with H 0 of a unique hyperplane H of E d+1 satisfying the following three properties: 1. H supports P, 2. H is the translated copy of a hyperplane tangent to λ 0 along one of its generatrices, 3. H is above O. Conversely, let H be a hyperplane of E d+1 satisfying the above three properties. Its intersection with H 0 is a hyperplane of E d supporting CH d (Σ). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 27 / 42
36 Introduction Upper bound Lower bound for two radii Lower bound for at least three radii Π i : the hyperplane in E d+1 with eq. {x d+1 = ρ i } P i : the points of P in Π i n i : the cardinality of P i ˆP i : the vertex set of P i, i.e., CH d (P i ) = CH d ( ˆP i ) ˆn i : the cardinality of ˆP i ˆn i n i ˆP = ˆP 1 ˆP 2... ˆP m ˆP = CH d+1 ( ˆP ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 28 / 42
37 Introduction Upper bound Lower bound for two radii Lower bound for at least three radii Π i : the hyperplane in E d+1 with eq. {x d+1 = ρ i } P i : the points of P in Π i n i : the cardinality of P i ˆP i : the vertex set of P i, i.e., CH d (P i ) = CH d ( ˆP i ) ˆn i : the cardinality of ˆP i ˆn i n i ˆP = ˆP 1 ˆP 2... ˆP m ˆP = CH d+1 ( ˆP )! It is possible that P ˆP This can happen if P 1 ˆP 1 or P m ˆP m OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 28 / 42
38 Upper bound Introduction Upper bound Lower bound for two radii Lower bound for at least three radii There is an injection ϕ : CH d (Σ) P mapping a face of circularity (d l 1) of CH d (Σ) to a unique l-face of P. Points in P i \ ˆP i can never be points on a supporting hyperplane H of P of the theorem. The injection ϕ is in fact an injection from CH d (Σ) to ˆP ˆP is the convex hull of m parallel polytopes. Therefore the complexity of CH d (Σ) is O( ˆn iˆn d 2 j ) = O( 1 i j m 1 i j m n i n d 2 j ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 29 / 42
39 The trigonometric moment curve Upper bound Lower bound for two radii Lower bound for at least three radii For any even dimension 2δ we can define the so called trigonometric moment curve: γ2δ tr (t) = (cos t, sin t, cos 2t, sin 2t,..., cos δt, sin δt), t [0, π) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 30 / 42
40 The trigonometric moment curve Upper bound Lower bound for two radii Lower bound for at least three radii For any even dimension 2δ we can define the so called trigonometric moment curve: γ2δ tr (t) = (cos t, sin t, cos 2t, sin 2t,..., cos δt, sin δt), t [0, π) Fact For any set P of n points on γ tr 2δ (t), the convex hull CH 2δ (P ) is a polytope Q combinatorially equivalent to the cyclic polytope C 2δ (n) Fact f 2δ 1 (Q) = Θ(n δ ) Fact Points on γ tr δ (t) lie on the sphere of E2δ centered at the origin with radius δ. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 30 / 42
41 First step: the prism Upper bound Lower bound for two radii Lower bound for at least three radii We assume that the ambient space is E d, where d 3 and d odd. 1 Define H 1 = {x d = z 1 } and H 2 = {x d = z 2 }, z 2 > z 1 + 2(n 2 + 2) δ 2 Σ 1 : set of n points on γ tr d 1 (t) H 1, where for n 1 points t (0, π 2 ), whereas for the last point t ( π 2, π) 3 Σ 2 : vertical projection of Σ 1 on H 2. 4 Q i = CH d 1 (Σ i ), and = CH d (Σ 1 Σ 2 ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 31 / 42
42 First step: the prism Upper bound Lower bound for two radii Lower bound for at least three radii We assume that the ambient space is E d, where d 3 and d odd. 1 Define H 1 = {x d = z 1 } and H 2 = {x d = z 2 }, z 2 > z 1 + 2(n 2 + 2) δ 2 Σ 1 : set of n points on γ tr d 1 (t) H 1, where for n 1 points t (0, π 2 ), whereas for the last point t ( π 2, π) 3 Σ 2 : vertical projection of Σ 1 on H 2. 4 Q i = CH d 1 (Σ i ), and = CH d (Σ 1 Σ 2 ) Some facts: ➊ f d 2 (Q i ) = Θ(n d ) = Θ(n d 2 1 ) (# of facets of Q i ) ➋ consists of the bottom facet Q 1, the top facet Q 2, and Θ(n d 2 1 ) vertical facets OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 31 / 42
43 More facts Introduction Upper bound Lower bound for two radii Lower bound for at least three radii For a vertical facet F denote by ν F the unit normal vector and by F + and F the two open halfspaces bounded by F. Call vertical ridges, the ridges that are intersections of vertical facets of Call Y the hyperplane with unit normal vector ν = (1, 0,..., 0), and Y +, Y the two oriented halfspaces delimited by Y. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 32 / 42
44 More facts Introduction Upper bound Lower bound for two radii Lower bound for at least three radii For a vertical facet F denote by ν F the unit normal vector and by F + and F the two open halfspaces bounded by F. Call vertical ridges, the ridges that are intersections of vertical facets of Call Y the hyperplane with unit normal vector ν = (1, 0,..., 0), and Y +, Y the two oriented halfspaces delimited by Y. ➌ n 1 points of Σ 1 (or Σ 2 ) belong to Y +. ➍ 1 point of Σ 1 (or Σ 2 ) belongs to Y. ➎ The (d 2)-polytope Q i = Q i Y has at most n 1 vertices. ➏ The faces F i of Q i intersected by Y is at most O(n d ) = O(n d ). ➐ The number of vertices facets of in Y + is Θ(n d 2 1 ). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 32 / 42
45 Upper bound Lower bound for two radii Lower bound for at least three radii Second step: spheres with non-zero radius xd Y F Define a new set of n spheres Σ 3, where σ k = ((0, 0,..., 0, (2k + 1) δ), ρ) νf Choose ρ so that ➊ each σ k does not intersect any vertical ridge of ρ < δ ➋ each σ k intersects all vertical facets of ν x1 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 33 / 42
46 Upper bound Lower bound for two radii Lower bound for at least three radii Third step: perturd the spheres of Σ 3 View from the top xd Y Define the spheres σ k with radius ρ and center:! kx c ε k = c k + ν = c 2 l k +ε 2 12 «ν k l=0 F νf Choose ε so that ➊ each σ k does not intersect any vertical ridge of ρ < δ ➋ each σ k intersects all vertical facets of in Y + ➌ the (d 2)-sphere σ k σ k is contained in F for all vertical facets F of in Y + Call Σ 3 the set of perturbed spheres ν x1 Let Σ = Σ 1 Σ 2 Σ 3 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 34 / 42
47 Faces of circularity (d 1) Upper bound Lower bound for two radii Lower bound for at least three radii For each pair (σ k, F ), where 1 k n 2 and F a vertical facet of in Y +, we have a unique face of circularity (d 1) in CH d (Σ) CH d (Σ) has n 2 Θ(n d 2 1 ) faces of circularity (d 1) The complexity of CH d (Σ) is Ω(n 2 n d 2 1 ) WLOG n 2 n 1, which implies that the complexity of CH d (Σ) is Ω(n 2 n d n 1 n d 2 2 ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 35 / 42
48 Construction with m 3 radii Upper bound Lower bound for two radii Lower bound for at least three radii Let N 1 = m i=2 n i, N 2 = n 1. Construct the sets Σ 1, Σ 2 and Σ 3 as in the case m = 2, where Σ 1 and Σ 2 contain each N points and Σ 3 contains N spheres. Replace n i among the N 1 points of Σ 1, Σ 2 in Y + by a sphere of radius r i. Replace the point of Σ 1, Σ 2 in Y by a sphere of radius r 2. Choose r, where 0 < r 1, such that the following two conditions are satisfied: ➊ r = CH d (Σ 1 Σ 2 ) is combinatorially equivalent to 0 ➋ The two requirements for the spheres of Σ 3 should still be satisfied. CH d (Σ) has N 2 Θ(N d 2 1 ) faces of circularity (d 1) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 36 / 42
49 Results Introduction Upper bound Lower bound for two radii Lower bound for at least three radii The complexity of CH d (Σ) is Ω( 1 i j m n in d 2 j ) WLOG, assume that n 2 n 1 n i, 3 i m. Then: m n 1 ( i=2 n i ) d 2 n 1 n d m(m 1) ( 1 i j m n i n d 2 j ) OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 37 / 42
50 Results Introduction Upper bound Lower bound for two radii Lower bound for at least three radii The complexity of CH d (Σ) is Ω( 1 i j m n in d 2 j ) WLOG, assume that n 2 n 1 n i, 3 i m. Then: m n 1 ( i=2 n i ) d 2 n 1 n d m(m 1) ( 1 i j m n i n d 2 j ) Theorem Fix some odd d 3. There exists a set Σ of spheres in E d, consisting of n i spheres of radius ρ i, with ρ 1 < ρ 2 <... < ρ m and m 3 fixed, such that the complexity of CH d (Σ) is Ω( 1 i j m n in d 2 j ). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 37 / 42
51 Results Introduction Upper bound Lower bound for two radii Lower bound for at least three radii The complexity of CH d (Σ) is Ω( 1 i j m n in d 2 j ) WLOG, assume that n 2 n 1 n i, 3 i m. Then: m n 1 ( i=2 n i ) d 2 n 1 n d m(m 1) ( 1 i j m n i n d 2 j ) Corollary Let P = {P 1, P 2,..., P m } be a set of m d-polytopes, lying on m parallel hyperplanes of E d+1, with d 3 odd, and both d, m fixed. The worst-case complexity of CH d+1 (P) is Ω( 1 i j m n in d 2 j ), where n i = f 0 (P i ), 1 i m. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 37 / 42
52 Summary & extensions Open problems Minkowski sum of two convex d-polytopes P: n-vertex convex d-polytope Q: m-vertex convex d-polytope Embed P and Q in E d+1 and on the hyperplanes {x d+1 = 0} and {x d+1 = 1}, respectively. The Minkowski sum (1 λ)p λq, λ (0, 1) is combinatorially equivalent to the intersection of CH d+1 ({P, Q}) with the hyperplane {x d+1 = λ}. Corollary Let P and Q be two convex d-polytopes in E d, with n and m vertices, respectively, where d 3 and d odd. The complexity of the weighted Minkowski sum (1 λ)p λq, λ (0, 1), as well as the complexity of the Minkowski sum P Q, is Θ(mn d 2 + nm d 2 ). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 38 / 42
53 Summary & extensions Open problems Computing convex hulls of parallel polytopes We can use output-sensitive algorithms for computing CH d+1 (P). For d 5 the algorithm in [Seidel 1986] and the algorithm in [Matoušek & Schwarzkopf 1992] yield the same running time: O(( 1 i j m n i n d 2 j ) log n), n = m n i. i=1 OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 39 / 42
54 Summary & extensions Open problems Computing convex hulls of parallel polytopes We can use output-sensitive algorithms for computing CH d+1 (P). For d 5 the algorithm in [Seidel 1986] and the algorithm in [Matoušek & Schwarzkopf 1992] yield the same running time: O(( 1 i j m n i n d 2 j ) log n), n = m n i. For d = 3 the best two choices are the algorithm in [Matoušek & Schwarzkopf 1992], and that in [Chan, Snoeyink & Yap 1997]. The running time is in: O(min{n 4/3+ɛ +( n i n j ) log n, ( n i n j ) log 2 n}). 1 i j m i=1 1 i j m OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 39 / 42
55 Computing convex hulls of spheres Summary & extensions Open problems We can slightly modify the algorithm in [Boissonnat et al. 1996] to get a new algorithm adapted to the fact that the spheres have a constant number of radii. The algorithm first computes the convex hull for each subset of spheres having the same radius, and keeps only the spheres defining each convex hull. These spheres are then used to compute the convex hull of the entire sphere set. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 40 / 42
56 Computing convex hulls of spheres Summary & extensions Open problems We can slightly modify the algorithm in [Boissonnat et al. 1996] to get a new algorithm adapted to the fact that the spheres have a constant number of radii. The algorithm first computes the convex hull for each subset of spheres having the same radius, and keeps only the spheres defining each convex hull. These spheres are then used to compute the convex hull of the entire sphere set. Time complexity: O(n d 2 + n log n + T d+1 (n 1, n 2,..., n m )) where n i is the number of spheres of radius ρ i, n the total number of spheres, and T d+1 (n 1, n 2,..., n m ) stands for the time to compute the convex hull of m parallel convex d-polytopes in E d+1, where the i-th polytope has n i vertices (cf. previous slide). OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 40 / 42
57 Open problems Introduction Summary & extensions Open problems Worst-case optimal algorithms for all odd dimensions Refinement of worst-case complexity of additively weighted Voronoi cells, when the number of radii is fixed. Tight bound on the complexity of the Minkowski sum when we have at least three summands. Exact complexity for the Minkowski sum of two (or more) polytopes. OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 41 / 42
58 Open problems Introduction Summary & extensions Open problems Worst-case optimal algorithms for all odd dimensions Refinement of worst-case complexity of additively weighted Voronoi cells, when the number of radii is fixed. Tight bound on the complexity of the Minkowski sum when we have at least three summands. Exact complexity for the Minkowski sum of two (or more) polytopes. THANK YOU FOR YOUR ATTENTION OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 41 / 42
59 View from the top Introduction Summary & extensions Open problems Y x 1 F ν F ν Go back OrbiCG/Triangles Workshop on CG, December 9, 2010 Convex hulls of spheres/parallel polytopes 42 / 42
Shellings, the Euler-Poincaré Formula for Polytopes, Dehn-Sommerville Equations, the Upper Bound Theorem
Chapter 8 Shellings, the Euler-Poincaré Formula for Polytopes, Dehn-Sommerville Equations, the Upper Bound Theorem 8.1 Shellings The notion of shellability is motivated by the desire to give an inductive
More informationLecture 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 informationAutomorphism Groups of Cyclic Polytopes
8 Automorphism Groups of Cyclic Polytopes (Volker Kaibel and Arnold Waßmer ) It is probably well-known to most polytope theorists that the combinatorial automorphism group of a cyclic d-polytope with n
More informationVoronoi Diagrams, Delaunay Triangulations and Polytopes
Voronoi Diagrams, Delaunay Triangulations and Polytopes Jean-Daniel Boissonnat MPRI, Lecture 2 Computational Geometry Learning Voronoi, Delaunay & Polytopes MPRI, Lecture 2 1 / 43 Voronoi diagrams in nature
More informationLecture 3 Randomized Algorithms
Lecture 3 Randomized Algorithms Jean-Daniel Boissonnat Winter School on Computational Geometry and Topology University of Nice Sophia Antipolis January 23-27, 2017 Computational Geometry and Topology Randomized
More informationVoronoi diagram and Delaunay triangulation
Voronoi diagram and Delaunay triangulation Ioannis Emiris & Vissarion Fisikopoulos Dept. of Informatics & Telecommunications, University of Athens Computational Geometry, spring 2015 Outline 1 Voronoi
More informationVoronoi Diagrams and Delaunay Triangulations. O Rourke, Chapter 5
Voronoi Diagrams and Delaunay Triangulations O Rourke, Chapter 5 Outline Preliminaries Properties and Applications Computing the Delaunay Triangulation Preliminaries Given a function f: R 2 R, the tangent
More informationWe have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance.
Solid geometry We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance. First, note that everything we have proven for the
More informationSimplicial Complexes: Second Lecture
Simplicial Complexes: Second Lecture 4 Nov, 2010 1 Overview Today we have two main goals: Prove that every continuous map between triangulable spaces can be approximated by a simplicial map. To do this,
More informationLecture 2 - Introduction to Polytopes
Lecture 2 - Introduction to Polytopes Optimization and Approximation - ENS M1 Nicolas Bousquet 1 Reminder of Linear Algebra definitions Let x 1,..., x m be points in R n and λ 1,..., λ m be real numbers.
More informationApproximate Nearest Neighbor Problem: Improving Query Time CS468, 10/9/2006
Approximate Nearest Neighbor Problem: Improving Query Time CS468, 10/9/2006 Outline Reducing the constant from O ( ɛ d) to O ( ɛ (d 1)/2) in uery time Need to know ɛ ahead of time Preprocessing time and
More informationLecture 12 March 4th
Math 239: Discrete Mathematics for the Life Sciences Spring 2008 Lecture 12 March 4th Lecturer: Lior Pachter Scribe/ Editor: Wenjing Zheng/ Shaowei Lin 12.1 Alignment Polytopes Recall that the alignment
More informationOptimal detection of intersections between convex polyhedra
1 2 Optimal detection of intersections between convex polyhedra Luis Barba Stefan Langerman 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Abstract For a polyhedron P in R d, denote by P its combinatorial complexity,
More informationarxiv: v1 [math.co] 12 Aug 2018
CONVEX UNION REPRESENTABILITY AND CONVEX CODES R. AMZI JEFFS AND ISABELLA NOVIK arxiv:1808.03992v1 [math.co] 12 Aug 2018 Abstract. We introduce and investigate d-convex union representable complexes: the
More informationComplexity of Delaunay triangulation for points on lower-dimensional polyhedra
Complexity of Delaunay triangulation for points on lower-dimensional polyhedra Nina Amenta UC Davis Dominique Attali CNRS-LIS Grenoble Olivier Devillers INRIA Sophia-Antipolis Abstract We show that the
More informationComplexity of Delaunay Triangulation for Points on Lower-dimensional Polyhedra.
published in Proc. 18th ACM-SIAM Sympos. Discrete Algorithms (SODA) 2007, p.1106 1113. Complexity of Delaunay Triangulation for Points on Lower-dimensional Polyhedra. Nina Amenta UC Davis Dominique Attali
More informationLifting Transform, Voronoi, Delaunay, Convex Hulls
Lifting Transform, Voronoi, Delaunay, Convex Hulls Subhash Suri Department of Computer Science University of California Santa Barbara, CA 93106 1 Lifting Transform (A combination of Pless notes and my
More informationChapter 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 informationFigure 2.1: An example of a convex set and a nonconvex one.
Convex Hulls 2.1 Definitions 2 Convexity is the key to understanding and simplifying geometry, and the convex hull plays a role in geometry akin to the sorted order for a collection of numbers. So what
More informationLinear Programming in Small Dimensions
Linear Programming in Small Dimensions Lekcija 7 sergio.cabello@fmf.uni-lj.si FMF Univerza v Ljubljani Edited from slides by Antoine Vigneron Outline linear programming, motivation and definition one dimensional
More informationOptimal algorithms for constrained 1-center problems
Optimal algorithms for constrained 1-center problems Luis Barba 12, Prosenjit Bose 1, and Stefan Langerman 2 1 Carleton University, Ottawa, Canada, jit@scs.carleton.ca 2 Université Libre de Bruxelles,
More informationarxiv: v1 [cs.cc] 30 Jun 2017
On the Complexity of Polytopes in LI( Komei Fuuda May Szedlá July, 018 arxiv:170610114v1 [cscc] 30 Jun 017 Abstract In this paper we consider polytopes given by systems of n inequalities in d variables,
More informationPlanar 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 informationCrossing 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 informationOn the Hardness of Computing Intersection, Union and Minkowski Sum of Polytopes
On the Hardness of Computing Intersection, Union and Minkowski Sum of Polytopes Hans Raj Tiwary hansraj@cs.uni-sb.de FR Informatik Universität des Saarlandes D-66123 Saarbrücken, Germany Tel: +49 681 3023235
More informationLecture 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 informationA PROOF OF THE LOWER BOUND CONJECTURE FOR CONVEX POLYTOPES
PACIFIC JOURNAL OF MATHEMATICS Vol. 46, No. 2, 1973 A PROOF OF THE LOWER BOUND CONJECTURE FOR CONVEX POLYTOPES DAVID BARNETTE A d polytope is defined to be a cz-dimensional set that is the convex hull
More informationPebble 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 informationPolyhedral Voronoi Diagrams of Polyhedra in Three Dimensions
Polyhedral Voronoi Diagrams of Polyhedra in Three Dimensions Vladlen Koltun Micha Sharir April 21, 2003 Abstract We show that the complexity of the Voronoi diagram of a collection of disjoint polyhedra
More informationNumerical Optimization
Convex Sets Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on Let x 1, x 2 R n, x 1 x 2. Line and line segment Line passing through x 1 and x 2 : {y
More informationToric Cohomological Rigidity of Simple Convex Polytopes
Toric Cohomological Rigidity of Simple Convex Polytopes Dong Youp Suh (KAIST) The Second East Asian Conference on Algebraic Topology National University of Singapore December 15-19, 2008 1/ 28 This talk
More informationIn what follows, we will focus on Voronoi diagrams in Euclidean space. Later, we will generalize to other distance spaces.
Voronoi Diagrams 4 A city builds a set of post offices, and now needs to determine which houses will be served by which office. It would be wasteful for a postman to go out of their way to make a delivery
More informationCombinatorial Geometry & Topology arising in Game Theory and Optimization
Combinatorial Geometry & Topology arising in Game Theory and Optimization Jesús A. De Loera University of California, Davis LAST EPISODE... We discuss the content of the course... Convex Sets A set is
More informationLecture 0: Reivew of some basic material
Lecture 0: Reivew of some basic material September 12, 2018 1 Background material on the homotopy category We begin with the topological category TOP, whose objects are topological spaces and whose morphisms
More informationMonotone Paths in Geometric Triangulations
Monotone Paths in Geometric Triangulations Adrian Dumitrescu Ritankar Mandal Csaba D. Tóth November 19, 2017 Abstract (I) We prove that the (maximum) number of monotone paths in a geometric triangulation
More informationG 6i try. On the Number of Minimal 1-Steiner Trees* Discrete Comput Geom 12:29-34 (1994)
Discrete Comput Geom 12:29-34 (1994) G 6i try 9 1994 Springer-Verlag New York Inc. On the Number of Minimal 1-Steiner Trees* B. Aronov, 1 M. Bern, 2 and D. Eppstein 3 Computer Science Department, Polytechnic
More information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More informationLower Bounds for Geometric Diameter Problems
Lower Bounds for Geometric Diameter Problems p. 1/48 Lower Bounds for Geometric Diameter Problems Hervé Fournier University of Versailles St-Quentin en Yvelines Antoine Vigneron INRA Jouy-en-Josas Lower
More information2 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 informationHowever, 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 informationComputational Geometry
Computational Geometry 600.658 Convexity A set S is convex if for any two points p, q S the line segment pq S. S p S q Not convex Convex? Convexity A set S is convex if it is the intersection of (possibly
More informationDynamic Collision Detection
Distance Computation Between Non-Convex Polyhedra June 17, 2002 Applications Dynamic Collision Detection Applications Dynamic Collision Detection Evaluating Safety Tolerances Applications Dynamic Collision
More information3 No-Wait Job Shops with Variable Processing Times
3 No-Wait Job Shops with Variable Processing Times In this chapter we assume that, on top of the classical no-wait job shop setting, we are given a set of processing times for each operation. We may select
More informationFACES OF CONVEX SETS
FACES OF CONVEX SETS VERA ROSHCHINA Abstract. We remind the basic definitions of faces of convex sets and their basic properties. For more details see the classic references [1, 2] and [4] for polytopes.
More informationarxiv: v2 [math.co] 23 Jan 2018
CONNECTIVITY OF CUBICAL POLYTOPES HOA THI BUI, GUILLERMO PINEDA-VILLAVICENCIO, AND JULIEN UGON arxiv:1801.06747v2 [math.co] 23 Jan 2018 Abstract. A cubical polytope is a polytope with all its facets being
More informationCS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension
CS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension Antoine Vigneron King Abdullah University of Science and Technology November 7, 2012 Antoine Vigneron (KAUST) CS 372 Lecture
More informationHYPERPLANE SEPARABILITY AND CONVEXITY OF PROBABILISTIC POINT SETS
HYPERPLANE SEPARABILITY AND CONVEXITY OF PROBABILISTIC POINT SETS Martin Fink, John Hershberger, Nirman Kumar, and Subhash Suri Abstract. We describe an O(n d )-time algorithm for computing the exact probability
More informationREVIEW OF FUZZY SETS
REVIEW OF FUZZY SETS CONNER HANSEN 1. Introduction L. A. Zadeh s paper Fuzzy Sets* [1] introduces the concept of a fuzzy set, provides definitions for various fuzzy set operations, and proves several properties
More informationParallelohedra and the Voronoi Conjecture
Alexey Garber Bielefeld University November 28, 2012 Parallelohedra Denition Convex d-dimensional polytope P is called a parallelohedron if R d can be (face-to-face) tiled into parallel copies of P. Parallelohedra
More informationHausdorff Approximation of 3D Convex Polytopes
Hausdorff Approximation of 3D Convex Polytopes Mario A. Lopez Department of Mathematics University of Denver Denver, CO 80208, U.S.A. mlopez@cs.du.edu Shlomo Reisner Department of Mathematics University
More informationCOMP331/557. Chapter 2: The Geometry of Linear Programming. (Bertsimas & Tsitsiklis, Chapter 2)
COMP331/557 Chapter 2: The Geometry of Linear Programming (Bertsimas & Tsitsiklis, Chapter 2) 49 Polyhedra and Polytopes Definition 2.1. Let A 2 R m n and b 2 R m. a set {x 2 R n A x b} is called polyhedron
More informationarxiv: v2 [cs.cg] 24 Jul 2011
Ice-Creams and Wedge Graphs Eyal Ackerman Tsachik Gelander Rom Pinchasi Abstract arxiv:116.855v2 [cs.cg] 24 Jul 211 What is the minimum angle α > such that given any set of α-directional antennas (that
More informationLecture 5: Properties of convex sets
Lecture 5: Properties of convex sets Rajat Mittal IIT Kanpur This week we will see properties of convex sets. These properties make convex sets special and are the reason why convex optimization problems
More informationCAT(0)-spaces. Münster, June 22, 2004
CAT(0)-spaces Münster, June 22, 2004 CAT(0)-space is a term invented by Gromov. Also, called Hadamard space. Roughly, a space which is nonpositively curved and simply connected. C = Comparison or Cartan
More informationInternational Journal of Computational Geometry & Applications c World Scientific Publishing Company
International Journal of Computational Geometry & Applications c World Scientific Publishing Company FINDING SIMPLICES CONTAINING THE ORIGIN IN TWO AND THREE DIMENSIONS KHALED ELBASSIONI Max-Planck-Institut
More informationarxiv: v1 [cs.cg] 29 Dec 2018
Convex Polygons in Cartesian Products Jean-Lou De Carufel Adrian Dumitrescu Wouter Meulemans Tim Ophelders Claire Pennarun Csaba D. Tóth Sander Verdonschot arxiv:82.332v [cs.cg] 29 Dec 208 Abstract We
More informationVoronoi 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 informationA New Approach to Output-Sensitive Voronoi Diagrams and Delaunay Triangulations
A New Approach to Output-Sensitive Voronoi Diagrams and Delaunay Triangulations Don Sheehy INRIA Saclay, France soon: UConn Joint work with Gary Miller at CMU Voronoi Diagrams Voronoi Diagrams The Voronoi
More informationLecture 3. Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets. Tepper School of Business Carnegie Mellon University, Pittsburgh
Lecture 3 Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets Gérard Cornuéjols Tepper School of Business Carnegie Mellon University, Pittsburgh January 2016 Mixed Integer Linear Programming
More informationSimplicial Cells in Arrangements of Hyperplanes
Simplicial Cells in Arrangements of Hyperplanes Christoph Dätwyler 05.01.2013 This paper is a report written due to the authors presentation of a paper written by Shannon [1] in 1977. The presentation
More informationSimplicial 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 informationarxiv: v1 [cs.cg] 7 Oct 2017
A Proof of the Orbit Conjecture for Flipping Edge-Labelled Triangulations Anna Lubiw 1, Zuzana Masárová 2, and Uli Wagner 2 arxiv:1710.02741v1 [cs.cg] 7 Oct 2017 1 School of Computer Science, University
More informationThe Cut Locus and the Jordan Curve Theorem
The Cut Locus and the Jordan Curve Theorem Rich Schwartz November 19, 2015 1 Introduction A Jordan curve is a subset of R 2 which is homeomorphic to the circle, S 1. The famous Jordan Curve Theorem says
More informationAdvanced Algorithms Computational Geometry Prof. Karen Daniels. Fall, 2012
UMass Lowell Computer Science 91.504 Advanced Algorithms Computational Geometry Prof. Karen Daniels Fall, 2012 Voronoi Diagrams & Delaunay Triangulations O Rourke: Chapter 5 de Berg et al.: Chapters 7,
More informationIntegral Geometry and the Polynomial Hirsch Conjecture
Integral Geometry and the Polynomial Hirsch Conjecture Jonathan Kelner, MIT Partially based on joint work with Daniel Spielman Introduction n A lot of recent work on Polynomial Hirsch Conjecture has focused
More informationbe a polytope. has such a representation iff it contains the origin in its interior. For a generic, sort the inequalities so that
( Shelling (Bruggesser-Mani 1971) and Ranking Let be a polytope. has such a representation iff it contains the origin in its interior. For a generic, sort the inequalities so that. a ranking of vertices
More informationLinear Data Structures for Fast Ray-Shooting amidst Convex Polyhedra
Linear Data Structures for Fast Ray-Shooting amidst Convex Polyhedra Haim Kaplan, Natan Rubin, and Micha Sharir School of Computer Science, Tel Aviv University, Tel Aviv, Israel {haimk,rubinnat,michas}@post.tau.ac.il
More informationLecture 5: Simplicial Complex
Lecture 5: Simplicial Complex 2-Manifolds, Simplex and Simplicial Complex Scribed by: Lei Wang First part of this lecture finishes 2-Manifolds. Rest part of this lecture talks about simplicial complex.
More informationModeling and Analysis of Hybrid Systems
Modeling and Analysis of Hybrid Systems Convex polyhedra Prof. Dr. Erika Ábrahám Informatik 2 - LuFG Theory of Hybrid Systems RWTH Aachen University Szeged, Hungary, 27 September - 06 October 2017 Ábrahám
More informationModeling and Analysis of Hybrid Systems
Modeling and Analysis of Hybrid Systems 6. Convex polyhedra Prof. Dr. Erika Ábrahám Informatik 2 - LuFG Theory of Hybrid Systems RWTH Aachen University Szeged, Hungary, 27 September - 06 October 2017 Ábrahám
More informationGeometric Algorithms. 1. Objects
Geometric Algorithms 1. Definitions 2. Line-Segment properties 3. Inner Points of polygons 4. Simple polygons 5. Convex hulls 6. Closest pair of points 7. Diameter of a point set 1. Objects 2-dimensional
More informationEmpty-Ellipse Graphs
Empty-Ellipse Graphs Olivier Devillers Jeff Erickson Xavier Goaoc Submitted to SODA 2008: July 5, 2007 Abstract We define and study a geometric graph over points in the plane that captures the local behavior
More informationThe Medial Axis of the Union of Inner Voronoi Balls in the Plane
The Medial Axis of the Union of Inner Voronoi Balls in the Plane Joachim Giesen a, Balint Miklos b,, Mark Pauly b a Max-Planck Institut für Informatik, Saarbrücken, Germany b Applied Geometry Group, ETH
More informationBasics of Combinatorial Topology
Chapter 7 Basics of Combinatorial Topology 7.1 Simplicial and Polyhedral Complexes In order to study and manipulate complex shapes it is convenient to discretize these shapes and to view them as the union
More informationPolytopes Course Notes
Polytopes Course Notes Carl W. Lee Department of Mathematics University of Kentucky Lexington, KY 40506 lee@ms.uky.edu Fall 2013 i Contents 1 Polytopes 1 1.1 Convex Combinations and V-Polytopes.....................
More informationLINES AND FREE LINE SEGMENTS TANGENT TO ARBITRARY THREE-DIMENSIONAL CONVEX POLYHEDRA
LINES AND FREE LINE SEGMENTS TANGENT TO ARBITRARY THREE-DIMENSIONAL CONVEX POLYHEDRA HERVÉ BRÖNNIMANN, OLIVIER DEVILLERS, VIDA DUJMOVIĆ, HAZEL EVERETT, MARC GLISSE, XAVIER GOAOC, SYLVAIN LAZARD, HYEON-SUK
More informationPacific Journal of Mathematics
Pacific Journal of Mathematics SIMPLIFYING TRIANGULATIONS OF S 3 Aleksandar Mijatović Volume 208 No. 2 February 2003 PACIFIC JOURNAL OF MATHEMATICS Vol. 208, No. 2, 2003 SIMPLIFYING TRIANGULATIONS OF S
More informationCREPANT RESOLUTIONS OF GORENSTEIN TORIC SINGULARITIES AND UPPER BOUND THEOREM. Dimitrios I. Dais
Séminaires & Congrès 6, 2002, p. 187 192 CREPANT RESOLUTIONS OF GORENSTEIN TORIC SINGULARITIES AND UPPER BOUND THEOREM by Dimitrios I. Dais Abstract. A necessary condition for the existence of torus-equivariant
More informationRestricted-Orientation Convexity in Higher-Dimensional Spaces
Restricted-Orientation Convexity in Higher-Dimensional Spaces ABSTRACT Eugene Fink Derick Wood University of Waterloo, Waterloo, Ont, Canada N2L3G1 {efink, dwood}@violetwaterlooedu A restricted-oriented
More informationOn the Complexity of Many Faces in Arrangements of Pseudo-Segments and of Circles
On the Complexity of Many Faces in Arrangements of Pseudo-Segments and of Circles Pankaj K. Agarwal Boris Aronov Micha Sharir Abstract We obtain improved bounds on the complexity of m distinct faces in
More informationColoring. 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 informationChapter 4 Concepts from Geometry
Chapter 4 Concepts from Geometry An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Line Segments The line segment between two points and in R n is the set of points on the straight line joining
More informationDetecting intersections between convex polyhedra
Detecting intersections between convex polyhedra Luis Barba Stefan Langerman Abstract For a polyhedron P in R d, denote by P its combinatorial complexity, i.e., the number of faces, edges and vertices
More informationSIMION S TYPE B ASSOCIAHEDRON IS A PULLING TRIANGULATION OF THE LEGENDRE POLYTOPE
SIMION S TYPE B ASSOCIAHEDRON IS A PULLING TRIANGULATION OF THE LEGENDRE POLYTOPE RICHARD EHRENBORG, GÁBOR HETYEI AND MARGARET READDY Abstract. We show that the Simion type B associahedron is combinatorially
More informationOrthogonal Ham-Sandwich Theorem in R 3
Orthogonal Ham-Sandwich Theorem in R 3 Downloaded 11/24/17 to 37.44.201.8. Redistribution subject to SIAM license or copyright; see http://www.siam.org/journals/ojsa.php Abstract The ham-sandwich theorem
More informationLecture 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 informationOn the Complexity of Umbra and Penumbra
Author manuscript, published in "Computational Geometry 42, 8 (2009) 758--771" DOI : 10.1016/j.comgeo.2008.04.007 On the Complexity of Umbra and Penumbra J. Demouth O. Devillers H. Everett M. Glisse S.
More informationPlanarity. 1 Introduction. 2 Topological Results
Planarity 1 Introduction A notion of drawing a graph in the plane has led to some of the most deep results in graph theory. Vaguely speaking by a drawing or embedding of a graph G in the plane we mean
More informationAlgorithms for the Construction of Digital Convex Fuzzy Hulls
lgorithms for the Construction of Digital Convex Fuzzy Hulls Nebojša M. Ralević Faculty of Engineering, University of Novi Sad Trg Dositeja Obradovića 6, 000 Novi Sad, Serbia nralevic@uns.ns.ac.yu Lidija
More information4. Definition: topological space, open set, topology, trivial topology, discrete topology.
Topology Summary Note to the reader. If a statement is marked with [Not proved in the lecture], then the statement was stated but not proved in the lecture. Of course, you don t need to know the proof.
More informationBasic Measures for Imprecise Point Sets in R d
Basic Measures for Imprecise Point Sets in R d Heinrich Kruger Masters Thesis Game and Media Technology Department of Information and Computing Sciences Utrecht University September 2008 INF/SCR-08-07
More informationElementary Combinatorial Topology
Elementary Combinatorial Topology Frédéric Meunier Université Paris Est, CERMICS, Ecole des Ponts Paristech, 6-8 avenue Blaise Pascal, 77455 Marne-la-Vallée Cedex E-mail address: frederic.meunier@cermics.enpc.fr
More informationMath 734 Aug 22, Differential Geometry Fall 2002, USC
Math 734 Aug 22, 2002 1 Differential Geometry Fall 2002, USC Lecture Notes 1 1 Topological Manifolds The basic objects of study in this class are manifolds. Roughly speaking, these are objects which locally
More informationDense Point Sets Have Sparse Delaunay Triangulations or... But Not Too Nasty
Dense Point Sets Have Sparse Delaunay Triangulations or... But Not Too Nasty Jeff Erickson University of Illinois, Urbana-Champaign jeffe@cs.uiuc.edu http://www.cs.uiuc.edu/ jeffe Submitted to Discrete
More informationMath 5593 Linear Programming Lecture Notes
Math 5593 Linear Programming Lecture Notes Unit II: Theory & Foundations (Convex Analysis) University of Colorado Denver, Fall 2013 Topics 1 Convex Sets 1 1.1 Basic Properties (Luenberger-Ye Appendix B.1).........................
More informationPi at School. Arindama Singh Department of Mathematics Indian Institute of Technology Madras Chennai , India
Pi at School rindama Singh epartment of Mathematics Indian Institute of Technology Madras Chennai-600036, India Email: asingh@iitm.ac.in bstract: In this paper, an attempt has been made to define π by
More informationLecture 3: Convex sets
Lecture 3: Convex sets Rajat Mittal IIT Kanpur We denote the set of real numbers as R. Most of the time we will be working with space R n and its elements will be called vectors. Remember that a subspace
More informationConvexity: an introduction
Convexity: an introduction Geir Dahl CMA, Dept. of Mathematics and Dept. of Informatics University of Oslo 1 / 74 1. Introduction 1. Introduction what is convexity where does it arise main concepts and
More informationConvex Geometry arising in Optimization
Convex Geometry arising in Optimization Jesús A. De Loera University of California, Davis Berlin Mathematical School Summer 2015 WHAT IS THIS COURSE ABOUT? Combinatorial Convexity and Optimization PLAN
More informationConvex Sets. Pontus Giselsson
Convex Sets Pontus Giselsson 1 Today s lecture convex sets convex, affine, conical hulls closure, interior, relative interior, boundary, relative boundary separating and supporting hyperplane theorems
More information