ASKSAGE: Sage Q&A Forum - RSS feedhttps://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Thu, 14 Dec 2023 10:45:15 +0100Is there a problem in Voronoi cell computation ?https://ask.sagemath.org/question/74836/is-there-a-problem-in-voronoi-cell-computation/Dear community,
I was wondering if there is a problem in my following computation of voronoi cell ?
Here is the code I am using :
from sage.modules.free_module_integer import IntegerLattice
l = [7, 0, -1, -2, -1, -2, 7, -2, 0, 0, -2, 0, 7, -2, 0, -1, -2, -1, 7, 0, -1, -1, 0, -2, 7]
M = matrix(5, 5, l)
print (IntegerLattice(M).voronoi_cell().volume() == M.det() )
If I am not mistaken, the last print should be true which is not the case.
Am I doing something wrong here ?alesavouThu, 14 Dec 2023 10:45:15 +0100https://ask.sagemath.org/question/74836/VoronoiDiagram returns empty regionshttps://ask.sagemath.org/question/49749/voronoidiagram-returns-empty-regions/Dear all,
Given the set of 2 points
P = [[-2687.19000000000, -2088.53000000000], [-2686.81000000000, -2084.19000000000]]
I would like to compute the voronoi diagram of $P$. When doing so i saw that the regions of the voronoi are:
{P(-2687.19000000000, -2088.53000000000): The empty polyhedron in RDF^0,
P(-2686.81000000000, -2084.19000000000): A 1-dimensional polyhedron in RDF^2 defined as the convex hull of 1 vertex and 1 ray}
In my code, I do need that those regions are not empty! Why is this happening and how can i fix it?
Please advise and thanks in advance.
P.s. Note that when I make the points a little I get that the voronoi is not empty.
Full Code:
P = [[-2687.19, -2088.53], [-2686.81, -2084.19]]
V = VoronoiDiagram(P)
R = V.regions()muradtukSun, 02 Feb 2020 17:43:31 +0100https://ask.sagemath.org/question/49749/