324 research outputs found
Kinetic and Dynamic Delaunay tetrahedralizations in three dimensions
We describe the implementation of algorithms to construct and maintain
three-dimensional dynamic Delaunay triangulations with kinetic vertices using a
three-simplex data structure. The code is capable of constructing the geometric
dual, the Voronoi or Dirichlet tessellation. Initially, a given list of points
is triangulated. Time evolution of the triangulation is not only governed by
kinetic vertices but also by a changing number of vertices. We use
three-dimensional simplex flip algorithms, a stochastic visibility walk
algorithm for point location and in addition, we propose a new simple method of
deleting vertices from an existing three-dimensional Delaunay triangulation
while maintaining the Delaunay property. The dual Dirichlet tessellation can be
used to solve differential equations on an irregular grid, to define partitions
in cell tissue simulations, for collision detection etc.Comment: 29 pg (preprint), 12 figures, 1 table Title changed (mainly
nomenclature), referee suggestions included, typos corrected, bibliography
update
Regular triangulations of dynamic sets of points
The Delaunay triangulations of a set of points are a class of
triangulations which play an important role in a variety of
different disciplines of science. Regular triangulations are a
generalization of Delaunay triangulations that maintain both their
relationship with convex hulls and with Voronoi diagrams. In regular
triangulations, a real value, its weight, is assigned to each point.
In this paper a simple data structure is presented that allows
regular triangulations of sets of points to be dynamically updated,
that is, new points can be incrementally inserted in the set and old
points can be deleted from it. The algorithms we propose for
insertion and deletion are based on a geometrical interpretation of
the history data structure in one more dimension and use lifted
flips as the unique topological operation. This results in rather
simple and efficient algorithms. The algorithms have been
implemented and experimental results are given.Postprint (published version
On Deletion in Delaunay Triangulation
This paper presents how the space of spheres and shelling may be used to
delete a point from a -dimensional triangulation efficiently. In dimension
two, if k is the degree of the deleted vertex, the complexity is O(k log k),
but we notice that this number only applies to low cost operations, while time
consuming computations are only done a linear number of times.
This algorithm may be viewed as a variation of Heller's algorithm, which is
popular in the geographic information system community. Unfortunately, Heller
algorithm is false, as explained in this paper.Comment: 15 pages 5 figures. in Proc. 15th Annu. ACM Sympos. Comput. Geom.,
181--188, 199
Three-dimensional alpha shapes
Frequently, data in scientific computing is in its abstract form a finite
point set in space, and it is sometimes useful or required to compute what one
might call the ``shape'' of the set. For that purpose, this paper introduces
the formal notion of the family of -shapes of a finite point set in
\Real^3. Each shape is a well-defined polytope, derived from the Delaunay
triangulation of the point set, with a parameter \alpha \in \Real controlling
the desired level of detail. An algorithm is presented that constructs the
entire family of shapes for a given set of size in time , worst
case. A robust implementation of the algorithm is discussed and several
applications in the area of scientific computing are mentioned.Comment: 32 page
Gap Processing for Adaptive Maximal Poisson-Disk Sampling
In this paper, we study the generation of maximal Poisson-disk sets with
varying radii. First, we present a geometric analysis of gaps in such disk
sets. This analysis is the basis for maximal and adaptive sampling in Euclidean
space and on manifolds. Second, we propose efficient algorithms and data
structures to detect gaps and update gaps when disks are inserted, deleted,
moved, or have their radius changed. We build on the concepts of the regular
triangulation and the power diagram. Third, we will show how our analysis can
make a contribution to the state-of-the-art in surface remeshing.Comment: 16 pages. ACM Transactions on Graphics, 201
A Systematic Review of Algorithms with Linear-time Behaviour to Generate Delaunay and Voronoi Tessellations
Triangulations and tetrahedrizations are important geometrical discretization procedures applied to several areas, such as the reconstruction of surfaces and data visualization. Delaunay and Voronoi tessellations are discretization structures of domains with desirable geometrical properties. In this work, a systematic review of algorithms with linear-time behaviour to generate 2D/3D Delaunay and/or Voronoi tessellations is presented
One machine, one minute, three billion tetrahedra
This paper presents a new scalable parallelization scheme to generate the 3D
Delaunay triangulation of a given set of points. Our first contribution is an
efficient serial implementation of the incremental Delaunay insertion
algorithm. A simple dedicated data structure, an efficient sorting of the
points and the optimization of the insertion algorithm have permitted to
accelerate reference implementations by a factor three. Our second contribution
is a multi-threaded version of the Delaunay kernel that is able to concurrently
insert vertices. Moore curve coordinates are used to partition the point set,
avoiding heavy synchronization overheads. Conflicts are managed by modifying
the partitions with a simple rescaling of the space-filling curve. The
performances of our implementation have been measured on three different
processors, an Intel core-i7, an Intel Xeon Phi and an AMD EPYC, on which we
have been able to compute 3 billion tetrahedra in 53 seconds. This corresponds
to a generation rate of over 55 million tetrahedra per second. We finally show
how this very efficient parallel Delaunay triangulation can be integrated in a
Delaunay refinement mesh generator which takes as input the triangulated
surface boundary of the volume to mesh
- …