Approximating Median Points in a Convex Polygon

We develop two simple and efficient approximation algorithms for the continuous $k$-medians problems, where we seek to find the optimal location of $k$ facilities among a continuum of client points in a convex polygon $C$ with $n$ vertices in a way that the total (average) Euclidean distance between clients and their nearest facility is minimized. Both algorithms run in $\mathcal{O}(n + k + k \log n)$ time. Our algorithms produce solutions within a factor of 2.002 of optimality. In addition, our simulation results applied to the convex hulls of the State of Massachusetts and the Town of Brookline, MA show that our algorithms generally perform within a range of 5\% to 22\% of optimality in practice

Data depth and floating body

Little known relations of the renown concept of the halfspace depth for multivariate data with notions from convex and affine geometry are discussed. Halfspace depth may be regarded as a measure of symmetry for random vectors. As such, the depth stands as a generalization of a measure of symmetry for convex sets, well studied in geometry. Under a mild assumption, the upper level sets of the halfspace depth coincide with the convex floating bodies used in the definition of the affine surface area for convex bodies in Euclidean spaces. These connections enable us to partially resolve some persistent open problems regarding theoretical properties of the depth

Automatic, computer aided geometric design of free-knot, regression splines

A new algorithm for Computer Aided Geometric Design of least squares (LS) splines with variable knots, named GeDS, is presented. It is based on interpreting functional spline regression as a parametric B-spline curve, and on using the shape preserving property of its control polygon. The GeDS algorithm includes two major stages. For the first stage, an automatic adaptive, knot location algorithm is developed. By adding knots, one at a time, it sequentially "breaks" a straight line segment into pieces in order to construct a linear LS B-spline fit, which captures the "shape" of the data. A stopping rule is applied which avoids both over and under fitting and selects the number of knots for the second stage of GeDS, in which smoother, higher order (quadratic, cubic, etc.) fits are generated. The knots appropriate for the second stage are determined, according to a new knot location method, called the averaging method. It approximately preserves the linear precision property of B-spline curves and allows the attachment of smooth higher order LS B-spline fits to a control polygon, so that the shape of the linear polygon of stage one is followed. The GeDS method produces simultaneously linear, quadratic, cubic (and possibly higher order) spline fits with one and the same number of B-spline regression functions. The GeDS algorithm is very fast, since no deterministic or stochastic knot insertion/deletion and relocation search strategies are involved, neither in the first nor the second stage. Extensive numerical examples are provided, illustrating the performance of GeDS and the quality of the resulting LS spline fits. The GeDS procedure is compared with other existing variable knot spline methods and smoothing techniques, such as SARS, HAS, MDL, AGS methods and is shown to produce models with fewer parameters but with similar goodness of fit characteristics, and visual quality

Geometric partitioning algorithms for fair division of geographic resources

University of Minnesota Ph.D. dissertation. July 2014. Major: Industrial and Systems Engineering. Advisor: John Gunnar Carlsson. 1 computer file (PDF): vi, 140 pages, appendices p. 129-140.This dissertation focuses on a fundamental but under-researched problem: how does one divide a piece of territory into smaller pieces in an efficient way? In particular, we are interested in \emph{map segmentation problem} of partitioning a geographic region into smaller subregions for allocating resources or distributing a workload among multiple agents. This work would result in useful solutions for a variety of fundamental problems, ranging from congressional districting, facility location, and supply chain management to air traffic control and vehicle routing. In a typical map segmentation problem, we are given a geographic region $R$, a probability density function defined on $R$ (representing, say population density, distribution of a natural resource, or locations of clients) and a set of points in $R$ (representing, say service facilities or vehicle depots). We seek a \emph{partition} of $R$ that is a collection of disjoint sub-regions $\{R_1, . . . , R_n\}$ such that $\bigcup_i R_i = R$, that optimizes some objective function while satisfying a shape condition. As examples of shape conditions, we may require that all sub-regions be compact, convex, star convex, simply connected (not having holes), connected, or merely measurable.Such problems are difficult because the search space is infinite-dimensional (since we are designing boundaries between sub-regions) and because the shape conditions are generally difficult to enforce using standard optimization methods. There are also many interesting variants and extensions to this problem. It is often the case that the optimal partition for a problem changes over time as new information about the region is collected. In that case, we have an \emph{online} problem and we must re-draw the sub-region boundaries as time progresses. In addition, we often prefer to construct these sub-regions in a \emph{decentralized} fashion: that is, the sub-region assigned to agent $i$ should be computable using only local information to agent $i$ (such as nearby neighbors or information about its surroundings), and the optimal boundary between two sub-regions should be computable using only knowledge available to those two agents.This dissertation is an attempt to design geometric algorithms aiming to solve the above mentioned problems keeping in view the various design constraints. We describe the drawbacks of the current approach to solving map segmentation problems, its ineffectiveness to impose geometric shape conditions and its limited utility in solving the online version of the problem. Using an intrinsically interdisciplinary approach, combining elements from variational calculus, computational geometry, geometric probability theory, and vector space optimization, we present an approach where we formulate the problems geometrically and then use a fast geometric algorithm to solve them. We demonstrate our success by solving problems having a particular choice of objective function and enforcing certain shape conditions. In fact, it turns out that such methods actually give useful insights and algorithms into classical location problems such as the continuous $k$-medians problem, where the aim is to find optimal locations for facilities. We use a map segmentation technique to present a constant factor approximation algorithm to solve the continuous $k$-medians problem in a convex polygon. We conclude this thesis by describing how we intend to build on this success and develop algorithms to solve larger classes of these problems

Non-Discriminatory Service Robot Placement Using Geometric Median

Service robots are becoming increasingly common, and businesses are adopting their use at an increasingly rapid rate in order to reduce costs and provide efficiencies in performing mundane tasks. However, very little research has been performed in order to understand and address ethical concerns regarding their deployment and use. One such concern is how one can ensure placement of a service robot such that is does not discriminate either in favor of or against individuals. This research explores techniques that can be used to provide a quantitative methodology to ensure fairness in terms of service robot placement such that discrimination does not occur. These techniques include the development and further enhancement of a heuristic hill climbing algorithm used to approximate the Geometric Median (GM). This algorithm is then benchmarked against Weiszfeld’s Algorithm, a well-known algorithm commonly used to solve the GM problem. iii These two algorithms are then visualized using Dynamics Explorer, an open source software tool, to create 2d maps of the dynamics of their convergence rates along with maps of F(), the “sum of the Euclidean distances” function underlying the calculations used by both GM approximation algorithms. The heuristic hill climbing algorithm is also extended to handle obstacles being introduced into the service robot’s workspace. It is further shown that as the size of ξ approaches ∞+, the Geometric Median converges to the centroid, given certain assumptions, such as the target points being evenly distributed in the plane

Isometric Embeddings in Trees and Their Use in Distance Problems

International audienceWe present powerful techniques for computing the diameter, all the eccentricities, and other related distance problems on some geometric graph classes, by exploiting their "tree-likeness" properties. We illustrate the usefulness of our approach as follows: (1) We propose a subquadratic-time algorithm for computing all eccentricities on partial cubes of bounded lattice dimension and isometric dimension O(n^{0.5−ε}). This is one of the first positive results achieved for the diameter problem on a subclass of partial cubes beyond median graphs. (2) Then, we obtain almost linear-time algorithms for computing all eccentricities in some classes of face-regular plane graphs, including benzenoid systems, with applications to chemistry. Previously, only a linear-time algorithm for computing the diameter and the center was known (and an O(n^{5/3})-time algorithm for computing all the eccentricities). (3) We also present an almost linear-time algorithm for computing the eccentricities in a polygon graph with an additive one-sided error of at most 2. (4) Finally, on any cube-free median graph, we can compute its absolute center in almost linear time. Independently from this work, Bergé and Habib have recently presented a linear-time algorithm for computing all eccentricities in this graph class (LAGOS'21), which also implies a linear-time algorithm for the absolute center problem. Our strategy here consists in exploiting the existence of some embeddings of these graphs in either a system or a product of trees, or in a single tree but where each vertex of the graph is embedded in a subset of nodes. While this may look like a natural idea, the way it can be done efficiently, which is our main technical contribution in the paper, is surprisingly intricate

Computing with Liquid Crystal Fingers: Models of geometric and logical computation

When a voltage is applied across a thin layer of cholesteric liquid crystal, fingers of cholesteric alignment can form and propagate in the layer. In computer simulation, based on experimental laboratory results, we demonstrate that these cholesteric fingers can solve selected problems of computational geometry, logic and arithmetics. We show that branching fingers approximate a planar Voronoi diagram, and non-branching fingers produce a convex subdivision of concave polygons. We also provide a detailed blue-print and simulation of a one-bit half-adder functioning on the principles of collision-based computing, where the implementation is via collision of liquid crystal fingers with obstacles and other fingers.Comment: submitted Sept 201