jeffrey-hokanson
commented
4 years ago This is currently implemented in voronoi_vertex and the old voronoi_vertex has been moved to voronoi_vertex_sample.
However, it is still necessary to fix this code so that it can handle equality constraints in the domain. Once this is done, I will close the issue.
See documentation from scipy: https://docs.scipy.org/doc/scipy/reference/generated/scipy.spatial.HalfspaceIntersection.html#scipy.spatial.HalfspaceIntersection
Idea: The structure of LinIneqDomain yields a series of halfplanes that define the space. We can think of the equidistant constraints for Voronoi diagrams as similarly defining halfplanes. Thus enumerating through the points x, we can construct the bounded Voronoi diagram. Likely this is not feasible for more than a few dimensions.