Convex Hull of Points Lying on Lines in o(n log n) Time after Preprocessing

Motivated by the desire to cope with data imprecision, we study methods for taking advantage of preliminary information about point sets in order to speed up the computation of certain structures associated with them. In particular, we study the following problem: given a set L of n lines in the plane, we wish to preprocess L such that later, upon receiving a set P of n points, each of which lies on a distinct line of L, we can construct the convex hull of P efficiently. We show that in quadratic time and space it is possible to construct a data structure on L that enables us to compute the convex hull of any such point set P in O(n alpha(n) log* n) expected time. If we further assume that the points are "oblivious" with respect to the data structure, the running time improves to O(n alpha(n)). The analysis applies almost verbatim when L is a set of line-segments, and yields similar asymptotic bounds. We present several extensions, including a trade-off between space and query time and an output-sensitive algorithm. We also study the "dual problem" where we show how to efficiently compute the (<= k)-level of n lines in the plane, each of which lies on a distinct point (given in advance). We complement our results by Omega(n log n) lower bounds under the algebraic computation tree model for several related problems, including sorting a set of points (according to, say, their x-order), each of which lies on a given line known in advance. Therefore, the convex hull problem under our setting is easier than sorting, contrary to the "standard" convex hull and sorting problems, in which the two problems require Theta(n log n) steps in the worst case (under the algebraic computation tree model).Comment: 26 pages, 5 figures, 1 appendix; a preliminary version appeared at SoCG 201

Quasiconvex Programming

We define quasiconvex programming, a form of generalized linear programming in which one seeks the point minimizing the pointwise maximum of a collection of quasiconvex functions. We survey algorithms for solving quasiconvex programs either numerically or via generalizations of the dual simplex method from linear programming, and describe varied applications of this geometric optimization technique in meshing, scientific computation, information visualization, automated algorithm analysis, and robust statistics.Comment: 33 pages, 14 figure

Strengthening QC relaxations of optimal power flow problems by exploiting various coordinate changes

Motivated by the potential for improvements in electric power system economics, this dissertation studies the AC optimal power flow (AC OPF) problem. An AC OPF problem optimizes a specified objective function subject to constraints imposed by both the non-linear power flow equations and engineering limits. The difficulty of an AC OPF problem is strongly connected to its feasible space\u27s characteristics. This dissertation first investigates causes of nonconvexities in AC OPF problems. Understanding typical causes of nonconvexities is helpful for improving AC OPF solution methodologies. This dissertation next focuses on solution methods for AC OPF problems that are based on convex relaxations. The quadratic convex (QC) relaxation is one promising approach that constructs convex envelopes around the trigonometric and product terms in the polar representation of the power flow equations. This dissertation proposes several improvements to strengthen QC relaxations of OPF problems. The first group of improvements provides tighter envelopes for the trigonometric functions and product terms in the power flow equations. Methods for obtaining tighter envelopes includes implementing Meyer and Floudas envelopes that yield the convex hull of trilinear monomials. Furthermore, by leveraging a representation of line admittances in polar form, this dissertation proposes tighter envelopes for the trigonometric terms. Another proposed improvement exploits the ability to rotate the base power used in the per unit normalization in order to facilitate the application of tighter trigonometric envelopes. The second group of improvements propose additional constraints based on new variables that represent voltage magnitude differences between connected buses. Using \u27bound tightening\u27 techniques, the bounds on the voltage magnitude difference variables can be significantly tighter than the bounds on the voltage magnitudes themselves, so constraints based on voltage magnitude differences can improve the QC relaxation --Abstract, page iv

Curve Stabbing Depth: Data Depth for Plane Curves

Measures of data depth have been studied extensively for point data. Motivated by recent work on analysis, clustering, and identifying representative elements in sets of trajectories, we introduce {\em curve stabbing depth} to quantify how deeply a given curve $Q$ is located relative to a given set $\cal C$ of curves in $\mathbb{R}^2$. Curve stabbing depth evaluates the average number of elements of $\cal C$ stabbed by rays rooted along the length of $Q$. We describe an $O(n^3 + n^2 m\log^2m+nm^2\log^2 m)$-time algorithm for computing curve stabbing depth when $Q$ is an $m$-vertex polyline and $\cal C$ is a set of $n$ polylines, each with $O(m)$ vertices.Comment: Preprin

Classroom Examples of Robustness Problems in Geometric Computations

International audienceThe algorithms of computational geometry are designed for a machine model with exact real arithmetic. Substituting floating point arithmetic for the assumed real arithmetic may cause implementations to fail. Although this is well known, there is no comprehensive documentation of what can go wrong and why. In this extended abstract, we study a simple incremental algorithm for planar convex hulls and give examples which make the algorithm fail in all possible ways. We also show how to construct failure-examples semi-systematically and discuss the geometry of the floating point implementation of the orientation predicate. We hope that our work will be useful for teaching computational geometry. The full paper is available at http://hal.inria.fr/inria-00344310/. It contains further examples, more theory, and color pictures. We strongly recommend to read the full paper instead of this extended abstract