Parallel Techniques for Paths, Visibility, and Related Problems (Thesis)

Total Page:16

File Type:pdf, Size:1020Kb

Parallel Techniques for Paths, Visibility, and Related Problems (Thesis) Purdue University Purdue e-Pubs Department of Computer Science Technical Reports Department of Computer Science 1992 Parallel Techniques for Paths, Visibility, and Related Problems (Thesis) Ziyi Danny Chen Report Number: 92-051 Chen, Ziyi Danny, "Parallel Techniques for Paths, Visibility, and Related Problems (Thesis)" (1992). Department of Computer Science Technical Reports. Paper 972. https://docs.lib.purdue.edu/cstech/972 This document has been made available through Purdue e-Pubs, a service of the Purdue University Libraries. Please contact [email protected] for additional information. PARALLEL TECHNIQUES FOR PATHS, VISIBILITY, AND RELATED PROBLEMS Ziyi Danny Chen CSD·TR-92·0S1 August 1992 PARALLEL TECHNIQUES FOR PATHS, VISIBILITY, AND RELATED PROBLEMS A Thesis Submitted to the Faculty of Purdue University by Ziyi Danny Chen In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy August 1992 11 To the memory of my father III ACKNOWLEDGMENTS There are many people that I would like to thank. First, I would like to express my greatest appreciation to my major professor Mikhail J. Atallah. He has always given me the freedom to develop my own research interest and at the same time provided all the support, guidance, and encouragement I needed. Learning how to do research from him and discussing problems with him have become some of my most joyful experience at Purdue. Furthermore, I am greatly indebted to Professor Greg N. Frederickson and Professor Susanne E. Hambrusch, for serving as my committee members and for giving me many needed suggestions and helps during my study at Purdue. I also like to thank Professor Christoph M. Hoffmann for serving as my committee member. I would like to thank the people in the CS department with whom I had discussed research problems, especially Professor Chandrajit Bajaj, Professor Chia-Hoang Lee, Scott McFaddin, Jyh-Jong Tsay, and Tony Tu. I would also like to thank all the people who have been my friends and helped me during my stay at Purdue. In particular, I thank Peter Ng, Thomas Tsang, Rose Ng, and Hiu Leung for many very needed helps during my early stage at Purdue and for many warm dinners and exciting badminton games, my office mates JiaXun Yu and Ching-Shoei Chiang for all the interesting chats and for giving me various helps, my office mate Susan Stratton for tolerating my frequent interruptions and kindly answering all my English questions, and Jindong Chen, Jung-Hong Chuang, Peggy Fahl, D.-J. Guan, Ajay Gupta, Hsin Pan, Bonita Rais, and many Chinese students for their friendship. I would certainly like to thank my "American parents in law," Dick and Gladys, for their love. I enjoy my weekend visits to their house, which always make me feel like being at home. IV I can never thank enough my father and mother. They gave me the most love and support. My father's love for education was my first motivation to pursue a Ph.D. degree. His inspiration, expectation, and confidence have been and will always be one of the greatest sources of support to me. My mother is great; she devoted all her time to raising and educating her eight children. lowe a great deal to my grandfather YueYe, whom I never have a chance to meet. It was his lifetime hard work that provided for my education for years. Many of my thanks go to my sister HuaYing; she and I shared the same dream of getting a higher education. She gave me the best support she could come up with. I thank all my brothers and sisters for their encouragement to my study. I cannot forget it was my uncle GengOn and aunt May who created the opportunity for me to study in the United States and provided a place for me to stay during my first three years in San Francisco. Finally, I am deeply grateful to my dear wife Xiaobo, for her love and patience, for her understanding and appreciation of my work, and for driving many times on the icy and endless road between Detroit and West Lafayette to bring her support and encouragement to me. This research was partially supported by the Office of Naval Research under Grants N00014-84-K-0502 and N00014-86-K-0689, the Air Force Office of Scientific Research under Grant AFOSR-90-0107, the National Science Foundation under Grant DCR­ 8451393, and the National Library of Medicine under Grant R01-LM05118. v TABLE OF CONTENTS Page LIST OF FIGURES Vll ABSTRACT ..... IX 1. INTRODUCTION . 1 1.1 The Parallel Computational Model 3 1.2 Some Basic Parallel Operations and Techniques 6 1.3 Previous Work 9 1.4 Our Thesis. .................... 12 2. RECTILINEAR SHORTEST PATHS WITH RECTANGULAR OBSTACLES 15 2.1 Preliminaries . 18 2.2 Computing a Staircase Separator 26 2.3 Other Building Blocks ...... 30 2.4 Computing the Lengths Matrix D p When IPI = O(IRI) . 34 2.5 Path Lengths between Arbitrary Points . 36 2.5.1 Some Useful Observations ... 36 2.5.2 The B(P)-to-VR Path Lengths .. 44 2.5.3 The VR-to- VR Path Lengths ... 45 2.5.4 Path Lengths with Arbitrary Query Points 51 2.6 Path Lengths When !PI >> IRI ....... 52 2.7 Computing the Actual Paths. ........... 55 2.8 A Note on the Sequential Time Complexity ... 58 3. VISIBILITY OF A SIMPLE POLYGONAL CHAIN FROM A POINT 61 3.1 Preliminaries .... ........... 63 3.2 An Overview of the Algorithm. ..... 66 3.3 Visibility Chains and Their Intersections 68 3.3.1 Simple Geometric Facts ..... 68 3.3.2 Simple Computational Observations. 73 VI Page 3.3.3 The Relative Positions of Ai and VIS(Ci) ..... 74 3.3.4 Computing the Portions of VIS(Ci) Hidden by Ai 78 3.4 The EREW PRAM Implementation. 84 3.5 Applications......................... 86 4. SOLVING GEOMETRIC PROBLEMS ON THE EREW PRAM. 87 4.1 Notation and Basic Algorithm Structure ... 90 4.2 Rank Trees and Parallel Operations ...... 91 4.3 Computing the Convex Hull of Sorted Points. 94 4.4 Triangulating a Trapezoidally Decomposed Polygon 97 4.5 Other Geometric Algorithms . 99 5. DETECTING WEAK VISIBILITY OF A SIMPLE POLYGON 102 5.1 Preliminaries ....... 105 5.2 Some Useful Observations .. ............... 107 5.3 Detecting the Weak Visibility of a Star-Shaped Polygon. 113 5.3.1 Phase 1.A ................ 114 5.3.2 Phase 1.B ...................... 129 5.3.3 Computing the Union of Bad Intervals ...... 131 5.4 Checking the Weak Visibility of a Polygon from an Edge 132 5.4.1 The Preprocess 133 5.4.2 The Basic Idea 134 5.4.3 Phase 2.A ... 138 5.4.4 Phase 2.B ... 144 5.5 Detecting the Weak Visibility of a Simple Polygon. 152 5.5.1 The Reduction by Sack and Suri 152 5.5.2 Our Reduction Procedure .. ........ 154 5.6 Applications...................... 162 5.6.1 Computing Shortest Paths in a Weakly Visible Polygon. 163 5.6.2 Triangulating a Weakly Visible Polygon . 166 5.6.3 Solving the One-Cruising-Guard Problem. 167 5.6.4 Other Applications 170 6. CONCLUSION . 172 BIBLIOGRAPHY .................." 176 VITA 187 VB LIST OF FIGURES 2.1 A rectilinear shortest path between v and w. 16 2.2 Illustrating MAXNE(R') and MAXsw(R'). .. 20 2.3 Illustrating Env(R') and the circular ordering on Bound(Env(R')). 22 2.4 Illustrating B(Q). 23 2.5 Illustrating Monge and non-Monge matrices of path lengths. 24 2.6 Illustrating NE(p) and WS(p). .. 27 2.7 Illustrating the algorithm for Sep. 29 2.8 Illustrating array Vert. 31 2.9 Illustrating Lemma 2.8. 33 2.10 Illustrating the proof of Theorem 2.2. 35 2.11 Illustrating U, U', W, and W' at a node v of T. 37 2.12 Illustrating the proof of Lemma 2.14. 42 2.13 Illustrating Lemma 2.15. ........ 44 2.14 Illustrating the computation of (a) mode 1, and (b) mode 2. 49 2.15 Illustrating the partition of Bound(P). ... 53 3.1 An example of the visibility chain VIS(P). 62 3.2 Illustrating the definitions. ..... 65 3.3 Illustrating the proof of Lemma 3.1. 70 3.4 Illustrating the two possible two-intersection cases. 72 VIll Figure Page 3.5 Illustrating the proofs of Lemmas 3.4 and 3.5. 81 4.1 The convex hull for a set S of planar points. 94 4.2 A triangulation of a monotone polygon. ... 97 5.1 The weakly visible edges of P are ell e2, and e4. 103 5.2 Illustrating HE; and HE;. E.g., HEi = {e3' e4} and HE;; = {es}. 105 5.3 Illustrating the proof of Lemma 5.5. (a) Vj is in Q. (b) Vj is not in Q.. 112 5.4 The four possible cases of IICP(Cj)l. 117 5.5 Illustrating the proof of Lemma 5.6. 119 5.6 Partitioning P into PI' P2 , and P3 • 133 5.7 Illustrating Lemma 5.10. ...... 135 5.8 Illustrating the proof of Lemma 5.12. 137 5.9 Illustrating the proof of Lemma 5.14. 141 5.10 Illustrating the structure of tree Te. 143 5.11 Illustrating Lemma 5.15. ....... 145 5.12 Partitioning Pinto Lp , Mp , and R p . 153 5.13 Illustrating the case where the procedure fails to find p*. 160 5.14 Partitioning P based on edge e. .... 163 5.15 Computing the shortest path SP(p, q). 164 5.16 Illustrating the two types of characteristic intervals. 169 IX ABSTRACT Chen, Ziyi Danny. Ph.D., Purdue University, August 1992.
Recommended publications
  • 30 POLYGONS Joseph O’Rourke, Subhash Suri, and Csaba D
    30 POLYGONS Joseph O'Rourke, Subhash Suri, and Csaba D. T´oth INTRODUCTION Polygons are one of the fundamental building blocks in geometric modeling, and they are used to represent a wide variety of shapes and figures in computer graph- ics, vision, pattern recognition, robotics, and other computational fields. By a polygon we mean a region of the plane enclosed by a simple cycle of straight line segments; a simple cycle means that nonadjacent segments do not intersect and two adjacent segments intersect only at their common endpoint. This chapter de- scribes a collection of results on polygons with both combinatorial and algorithmic flavors. After classifying polygons in the opening section, Section 30.2 looks at sim- ple polygonizations, Section 30.3 covers polygon decomposition, and Section 30.4 polygon intersection. Sections 30.5 addresses polygon containment problems and Section 30.6 touches upon a few miscellaneous problems and results. 30.1 POLYGON CLASSIFICATION Polygons can be classified in several different ways depending on their domain of application. In chip-masking applications, for instance, the most commonly used polygons have their sides parallel to the coordinate axes. GLOSSARY Simple polygon: A closed region of the plane enclosed by a simple cycle of straight line segments. Convex polygon: The line segment joining any two points of the polygon lies within the polygon. Monotone polygon: Any line orthogonal to the direction of monotonicity inter- sects the polygon in a single connected piece. Star-shaped polygon: The entire polygon is visible from some point inside the polygon. Orthogonal polygon: A polygon with sides parallel to the (orthogonal) coordi- nate axes.
    [Show full text]
  • Minimum Rectilinear Polygons for Given Angle Sequences ∗
    Minimum Rectilinear Polygons for Given Angle Sequences ∗ William S. Evansyz Krzysztof Fleszarx Philipp Kindermann{k Noushin Saeediy∗∗ Chan-Su Shinyy Alexander Wolff {zz Abstract A rectilinear polygon is a polygon whose edges are axis-aligned. Walking counterclock- wise on the boundary of such a polygon yields a sequence of left turns and right turns. The number of left turns always equals the number of right turns plus 4. It is known that any such sequence can be realized by a rectilinear polygon. In this paper, we consider the problem of finding realizations that minimize the perimeter or the area of the polygon or the area of the bounding box of the polygon. We show that all three problems are NP-hard in general. This answers an open question of Patrignani [CGTA 2001], who showed that it is NP-hard to minimize the area of the bounding box of an orthogonal drawing of a given planar graph. We also show that realizing polylines with minimum bounding box area is NP-hard. Then we consider the special cases of x-monotone and xy-monotone rectilinear polygons. For these, we can optimize the three objectives efficiently. 1 Introduction In this paper, we consider the problem of computing, for a given rectilinear angle sequence, a \small" rectilinear polygon that realizes the sequence. A rectilinear angle sequence S is a sequence of left (+90◦) turns, denoted by L, and right (−90◦) turns, denoted by R. We n write S = (s1; : : : ; sn) 2 fL; Rg , where n is the length of S. As we consider only rectilinear angle sequences, we usually drop the term \rectilinear." A polygon P realizes an angle sequence S if there is a counterclockwise (ccw) walk along the boundary of P such that the turns at the vertices of P , encountered during the walk, form the sequence S.
    [Show full text]
  • Optimal Uniformly Monotone Partitioning of Polygons with Holes
    Computer-Aided Design 44 (2012) 1235–1252 Contents lists available at SciVerse ScienceDirect Computer-Aided Design journal homepage: www.elsevier.com/locate/cad Optimal uniformly monotone partitioning of polygons with holes Xiangzhi Wei a, Ajay Joneja a,∗, David M. Mount b a Department of Industrial Engineering and Logistics Management, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong b Department of Computer Science and Institute for Advanced Computer Studies, University of Maryland, College Park, MD, USA article info a b s t r a c t Article history: Polygon partitioning is an important problem in computational geometry with a long history. In this Received 15 November 2011 paper we consider the problem of partitioning a polygon with holes into a minimum number of uniformly Accepted 11 June 2012 monotone components allowing arbitrary Steiner points. We call this the MUMC problem. We show that, given a polygon with n vertices and h holes and a scan direction, the MUMC problem relative to this Keywords: direction can be solved in time O.n log nCh log3 h/. Our algorithm produces a compressed representation Polygon subdivision of the subdivision of size O.n/, from which it is possible to extract either the entire decomposition or just Uniformly monotone subdivision the boundary of any desired component, in time proportional to the output size. When the scan direction Simple polygons 3 Network flow is not given, the problem can be solved in time O.K.n log n C h log h//, where K is the number of edges Planar graphs in the polygon's visibility graph.
    [Show full text]
  • Triangulating Polygons
    TRIANGULATING POLYGONS Vera Sacrist´an Computational Geometry Facultat d'Inform`aticade Barcelona Universitat Politcnica de Catalunya TRIANGULATING POLYGONS Computational Geometry, Facultat d'Inform`aticade Barcelona, UPC TRIANGULATING POLYGONS A polygon triangulation is the decomposition of a polygon into triangles. This is done by inserting internal diagonals. An internal diagonal is any segment... • connecting two vertices of the polygon and • completely enclosed in the polygon. Computational Geometry, Facultat d'Inform`aticade Barcelona, UPC TRIANGULATING POLYGONS A polygon triangulation is the decomposition of a polygon into triangles. This is done by inserting internal diagonals. An internal diagonal is any segment... • connecting two vertices of the polygon and • completely enclosed in the polygon. Computational Geometry, Facultat d'Inform`aticade Barcelona, UPC TRIANGULATING POLYGONS A polygon triangulation is the decomposition of a polygon into triangles. This is done by inserting internal diagonals. An internal diagonal is any segment... • connecting two vertices of the polygon and • completely enclosed in the polygon. Computational Geometry, Facultat d'Inform`aticade Barcelona, UPC TRIANGULATING POLYGONS A polygon triangulation is the decomposition of a polygon into triangles. This is done by inserting internal diagonals. An internal diagonal is any segment... • connecting two vertices of the polygon and • completely enclosed in the polygon. Computational Geometry, Facultat d'Inform`aticade Barcelona, UPC TRIANGULATING POLYGONS A polygon triangulation is the decomposition of a polygon into triangles. This is done by inserting internal diagonals. An internal diagonal is any segment... • connecting two vertices of the polygon and • completely enclosed in the polygon. Computational Geometry, Facultat d'Inform`aticade Barcelona, UPC TRIANGULATING POLYGONS A polygon triangulation is the decomposition of a polygon into triangles.
    [Show full text]
  • Multi-Goal Path Planning for Cooperative Sensing
    CZECH TECHNICAL UNIVERSITY IN PRAGUE FACULTY OF ELECTRICAL ENGINEERING Doctoral Thesis February, 2010 Jan Faigl CZECH TECHNICAL UNIVERSITY IN PRAGUE FACULTY OF ELECTRICAL ENGINEERING DEPARTMENT OF CYBERNETICS GERSTNER LABORATORY MULTI-GOAL PATH PLANNING FOR COOPERATIVE SENSING Doctoral Thesis Jan Faigl Prague, February, 2010 Ph.D. Programme: Electrical Engineering and Information Technology Branch of study: Artificial Intelligence and Biocybernetics Supervisor: Ing. Libor Pˇreuˇcil,CSc. Abstract The thesis deals with the multi-goal path planning problem. The problem is studied in the context of the inspection task of cooperating robots in a search and rescue mission. Two models of the sensing are considered: continuous and discrete. An environment is a priori known and it is represented by a polygonal domain in both approaches. The con- tinuous sensing assumes relatively inexpensive cost of sensing in comparison to the cost of motion and can be formulated as the Watchman Route Problem (WRP). The discrete sensing leads to the decoupled approach that is motivated by problems where the cost of sensing is dominant. The decoupled approach consists of two problems: the problem of finding minimal set of sensing locations, and the multi-goal path planning problem to find a path visiting the found locations. The set of sensing locations can be found as a solution of the Art Gallery Problem (AGP) or sensor placement problem if additional vis- ibility constraints have to be considered, e.g. visibility range, incident angle. The multi- goal path planning problem can be formulated as the Traveling Salesman Problem (TSP) in which paths between cities have to be traversable by the robot.
    [Show full text]
  • CS 473G Algorithmsx
    CS 473g Algorithms¬ Sariel Har-Peled December 7, 2009 ­This work is licensed under the Creative Commons Attribution-Noncommercial 3.0 License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc/3.0/ or send a letter to Creative Commons, 171 Second Street, Suite 300, San Francisco, California, 94105, USA. 2 Contents Preface 11 I NP Completeness 13 1 NP Completeness I 15 1.1 Introduction . 15 1.2 Complexity classes . 17 1.2.1 Reductions . 18 1.3 More NP-Complete problems . 19 1.3.1 3SAT ..................................... 19 1.4 Bibliographical Notes . 21 2 NP Completeness II 23 2.1 Max-Clique . 23 2.2 Independent Set . 25 2.3 Vertex Cover . 25 2.4 Graph Coloring . 26 3 NP Completeness III 29 3.1 Hamiltonian Cycle . 29 3.2 Traveling Salesman Problem . 30 3.3 Subset Sum . 31 3.4 3 dimensional Matching (3DM) . 32 3.5 Partition . 33 4 Dynamic programming 35 4.1 Basic Idea - Partition Number . 35 4.1.1 Memoization: . 37 4.2 Fibonacci numbers . 38 4.3 Edit Distance . 39 5 Dynamic programming II - The Recursion Strikes Back 43 5.1 Optimal Search Trees . 43 5.2 Optimal Triangulations . 45 3 5.3 Matrix Multiplication . 46 5.4 Longest Ascending Subsequence . 47 5.5 Pattern Matching . 47 6 Approximation algorithms 49 6.1 Greedy algorithms and approximation algorithms . 49 6.1.1 Alternative algorithm – two for the price of one . 51 6.2 Traveling Salesman Person . 51 6.2.1 TSP with the triangle inequality . 52 6.2.1.1 A 2-approximation .
    [Show full text]
  • Fast Geometric Algorithms
    FAST GEOMETRICALGORITHMS by Mark T. Noga Dissertation submitted to the Faculty of the Yi rgin ia Polytechnic Institute and State University in partial fulfillment of the requirements for the degree of DOCTOR OF PH I LOSO PHY in Computer Science and Applications APPROVED: D.C.S. Allison D. P. Rosel le R.M. Haralick R.W. Ehrich J .W. Roach January, 1984 Blacksburg, Virginia To Naida ii WORDS OF THANKS "I' II note you in my book of memory." - William Shakespeare, King Henry VI This thesis describes the culmination of five years of prolonged study at Virginia Polytechnic Institute. During that time, have benefitted from the close companionship of many people from all different parts of the United States and the world. Without your concern, understanding, encouragement, wisdom, and prodding, completion of this document would have been an all but impossible task. To all my friends I would Ii ke to say "thank you." Some of you gave more than what one would normally expect out of friendship and curiosity, and it is an honor for me to acknowledge these individually (and not in any special order). Occasionally you meet an individual who has the ability to take several seemingly unrelated results and somehow combine them to form what is often a colorful and unique solution to a particular problem. Doug Smith is one of these rare individuals. We spent many hours discussing, amongst other things, why art is pleasing to the eye (the entropy of art), or what features programming languages would contain in the year 2050 (or even if they would exist in a form equivalent to today's languages).
    [Show full text]
  • Algorithm for Triangulating a Simple Polygon
    On the Agenda Computational Geometry The Art Gallery Problem Polygon Triangulation Chapter 3 Polygons and Triangulation 13. 23. Art Gallery Problem Observations Given a simple polygon P, say The entire interior of a convex that two points p and q can see polygon is visible from any each other if the open segment pq R interior point. lies entirely within P. p q A point guards a region R ⊆ P if p convex sees all q∈R. Given a polygon P, what is the p non star-shaped minimal number of guards required A star-shaped polygon requires to guard P, and what are their only one guard located in its locations ? kernel . (actually this is the definition of a star-shape star-shaped polygon) 33. 43. Art Gallery Problem – Easy Upper Diagonals in Polygons Bound A diagonal of a polygon P is a line segment connecting two non-adjacent vertices, which lies entirely within P. Theorem: Any polygon with n>3 vertices n-2 guards suffice: has a diagonal, which may be found in Subdivide the polygon into n-2 O( n) time. triangles (triangulation) u Proof: Find the leftmost vertex v. Connect (we will show that this is the correct number of triangles) its two neighbors u and w. If this is not a diagonal there are other vertices inside v Place one guard in each triangle. the triangle uvw. Connect v to the vertex v’ v’ furthest from the (line containing the) w Theorem : Any simple planar polygon segment uw . with n vertices has a triangulation of size n-2.
    [Show full text]
  • Computational Geometry Fall 2002
    CMSC 754: Computational Geometry Fall 2002 http://www.cs.umd.edu/~mount/754/ Instructor: Dave Mount. Office: AVW 3209. Email: [email protected]. Office phone: (301) 405–2704. Office hours: Mon 3:00–4:00, Wed 3:30–4:30. I am also available immediately after class for questions. Feel free to send me email if you cannot make these times to set up another time. Teaching Assistant: Justin Wan. Office: AVW 1151. Email: [email protected]. Office hours: Mon 2:30–3:30, Thu 5:00–6:00. Class Time: Tue, Thur 2:00-3:15. Room: CSI 1122. Course Objectives: This is an introductory course to computational geometry and its appli- cations. We will discuss techniques needed in designing and analyzing efficient algorithms for problems in geometry, including convex hulls, geometric intersections, Voronoi diagrams, Delaunay triangulations, geometric data structures, and motion planning. Text: (Required) Computational Geometry: Algorithms and Applications (2nd Edition), M. de Berg, M. van Kreveld, M. Overmars, O. Schwarzkopf, Springer-Verlag, 2000. Prerequisites: CMSC 451, or its equivalent. • Knowledge of basic algorithm design techniques, such as divide-and-conquer, greedy algorithms, dynamic programming, • Knowledge of basic analysis techniques such as asymptotic notation, solving summations and recurrences, basic probability theory. • Knowledge of basic data structures (e.g. balanced binary trees and heaps). Course Work: There will be roughly 5 homework assignments. Normally, late homeworks will not be accepted. If you think that you have a legitimate excuse for handing in a homework late, contact me before the due date. There will be a midterm and comprehensive final exam.
    [Show full text]