Voronoi Graph -- Improved raycasting and integration schemes for high dimensional Voronoi diagrams

Alexander Sikorski,Martin Heida
2024-05-16
Abstract:The computation of Voronoi Diagrams, or their dual Delauney triangulations is difficult in high dimensions. In a recent publication Polianskii and Pokorny propose an iterative randomized algorithm facilitating the approximation of Voronoi tesselations in high dimensions. In this paper, we provide an improved vertex search method that is not only exact but even faster than the bisection method that was previously recommended. Building on this we also provide a depth-first graph-traversal algorithm which allows us to compute the entire Voronoi diagram. This enables us to compare the outcomes with those of classical algorithms like qHull, which we either match or marginally beat in terms of computation time. We furthermore show how the raycasting algorithm naturally lends to a Monte Carlo approximation for the volume and boundary integrals of the Voronoi cells, both of which are of importance for finite Volume methods. We compare the Monte-Carlo methods to the exact polygonal integration, as well as a hybrid approximation scheme.
Computational Geometry
What problem does this paper attempt to address?
This paper focuses on the optimization of Voronoi diagram computation in high-dimensional space, specifically the improvement of vertex search methods and ray shooting algorithms for Voronoi diagrams. The authors propose a new and more efficient vertex search algorithm that can terminate with just a few nearest neighbor searches, which is faster than the previously recommended binary search method. Based on this, they also develop a depth-first graph traversal algorithm to compute the entire Voronoi diagram and compare its computational time with traditional algorithms such as qHull, showing comparable or slightly superior performance. In addition, the paper introduces how to use ray shooting algorithms for Monte Carlo approximation, which is important for computing the volume and boundary integrals of Voronoi cells in finite volume methods. They compare these Monte Carlo methods with exact polygon integration and hybrid approximation schemes. The paper also discusses the performance of the algorithms, including computational complexity, pointing out that although they have the same Nlog(N) complexity as the Bowyer-Watson algorithm, previous methods require input nodes to be in general position, while the new algorithm can handle arbitrary point sets, including computing degenerate Voronoi diagrams. In conclusion, this paper aims to address the efficiency and accuracy issues in computing high-dimensional Voronoi diagrams, improving computational speed and accuracy through improved search strategies and integration methods.