4,683 research outputs found
Linear Expected Complexity for Directional and Multiplicative Voronoi Diagrams
While the standard unweighted Voronoi diagram in the plane has linear worst-case complexity, many of its natural generalizations do not. This paper considers two such previously studied generalizations, namely multiplicative and semi Voronoi diagrams. These diagrams both have quadratic worst-case complexity, though here we show that their expected complexity is linear for certain natural randomized inputs. Specifically, we argue that the expected complexity is linear for: (1) semi Voronoi diagrams when the visible direction is randomly sampled, and (2) for multiplicative diagrams when either weights are sampled from a constant-sized set, or the more challenging case when weights are arbitrary but locations are sampled from a square
Farthest-Polygon Voronoi Diagrams
Given a family of k disjoint connected polygonal sites in general position
and of total complexity n, we consider the farthest-site Voronoi diagram of
these sites, where the distance to a site is the distance to a closest point on
it. We show that the complexity of this diagram is O(n), and give an O(n log^3
n) time algorithm to compute it. We also prove a number of structural
properties of this diagram. In particular, a Voronoi region may consist of k-1
connected components, but if one component is bounded, then it is equal to the
entire region
A numerical algorithm for semi-discrete optimal transport in 3D
This paper introduces a numerical algorithm to compute the optimal
transport map between two measures and , where derives from a
density defined as a piecewise linear function (supported by a
tetrahedral mesh), and where is a sum of Dirac masses.
I first give an elementary presentation of some known results on optimal
transport and then observe a relation with another problem (optimal sampling).
This relation gives simple arguments to study the objective functions that
characterize both problems.
I then propose a practical algorithm to compute the optimal transport map
between a piecewise linear density and a sum of Dirac masses in 3D. In this
semi-discrete setting, Aurenhammer et.al [\emph{8th Symposium on Computational
Geometry conf. proc.}, ACM (1992)] showed that the optimal transport map is
determined by the weights of a power diagram. The optimal weights are computed
by minimizing a convex objective function with a quasi-Newton method. To
evaluate the value and gradient of this objective function, I propose an
efficient and robust algorithm, that computes at each iteration the
intersection between a power diagram and the tetrahedral mesh that defines the
measure .
The numerical algorithm is experimented and evaluated on several datasets,
with up to hundred thousands tetrahedra and one million Dirac masses.Comment: 23 pages, 14 figure
Discretization of the 3D Monge-Ampere operator, between Wide Stencils and Power Diagrams
We introduce a monotone (degenerate elliptic) discretization of the
Monge-Ampere operator, on domains discretized on cartesian grids. The scheme is
consistent provided the solution hessian condition number is uniformly bounded.
Our approach enjoys the simplicity of the Wide Stencil method, but
significantly improves its accuracy using ideas from discretizations of optimal
transport based on power diagrams. We establish the global convergence of a
damped Newton solver for the discrete system of equations. Numerical
experiments, in three dimensions, illustrate the scheme efficiency
- …