Try if it works if you use rationals (-48/100) instead of floats (-0.48). On Monday, February 3, 2020 at 8:48:34 AM UTC-5, Murad Tukan wrote: > > Dear all, > > Computer info: > CPU: i7-6500U @ 2.50GHZ > Memory: 16.0GB > > Sage information: SageMath 9.0 (Windows, installed using the .exe file), > and even on Linux version. > > Platforms ran at: SageMath 9.0 Shell, SageMath 9.0 Notebook, From Python > script (and then running the command sage -python "script_name.py"). > > Description of problem: > > When running the following code: > > P = [[-0.48, 0.83], [-5.33, 9.23]] > Q = [[0.43, 0. ], [9.04, 0. ]] > Z = [[-0.35, -0.60], [-1.44, -2.49]] > > V = VoronoiDiagram(P) > R = [r for r in V.regions().items()] > s = R[0][1].hyperplane_arrangement() > > V = VoronoiDiagram(Q) > R = [r for r in V.regions().items()] > s = s.add_hyperplane(R[0][1].hyperplane_arrangement()) > > V = VoronoiDiagram(Z) > R = [r for r in V.regions().items()] > s = s.add_hyperplane(R[0][1].hyperplane_arrangement()) > > s.regions() > > Sage fails to compute the regions which is weird since it can compute the > vertices which means that there is even a bounded region (in my case). When > running the same code on SageMathCell <https://sagecell.sagemath.org/>, > It is computed smoothly while on my machine it is not! > > The full traceback I get is: > > --------------------------------------------------------------------------- > ValueError Traceback (most recent call last) > <ipython-input-6-4ba993d8025f> in <module>() > ----> 1 s.regions() > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/misc/cachefunc.pyx > in sage.misc.cachefunc.CachedMethodCallerNoArgs.__call__ > (build/cythonized/sage/misc/cachefunc.c:12712)() > 2308 if self.cache is None: > 2309 f = self.f > -> 2310 self.cache = f(self._instance) > 2311 return self.cache > 2312 > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/hyperplane_arrangement/arrangement.py > > in regions(self) > 1608 R = self.base_ring() > 1609 dim = self.dimension() > -> 1610 universe = Polyhedron(eqns=[[0] + [0] * dim], base_ring=R) > 1611 regions = [universe] > 1612 for hyperplane in self: > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/polyhedron/constructor.py > > in Polyhedron(vertices, rays, lines, ieqs, eqns, ambient_dim, base_ring, > minimize, verbose, backend) > 660 if got_Vrep: > 661 Vrep = [vertices, rays, lines] > --> 662 return parent(Vrep, Hrep, convert=convert, verbose=verbose) > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/structure/parent.pyx > in sage.structure.parent.Parent.__call__ > (build/cythonized/sage/structure/parent.c:9245)() > 900 return mor._call_(x) > 901 else: > --> 902 return mor._call_with_args(x, args, kwds) > 903 > 904 raise TypeError(_LazyString(_lazy_format, ("No conversion > defined from %s to %s", R, self), {})) > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/structure/coerce_maps.pyx > > in sage.structure.coerce_maps.DefaultConvertMap_unique._call_with_args > (build/cythonized/sage/structure/coerce_maps.c:5081)() > 178 print(type(C), C) > 179 print(type(C._element_constructor), > C._element_constructor) > --> 180 raise > 181 > 182 > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/structure/coerce_maps.pyx > > in sage.structure.coerce_maps.DefaultConvertMap_unique._call_with_args > (build/cythonized/sage/structure/coerce_maps.c:4969)() > 173 return C._element_constructor(x, *args) > 174 else: > --> 175 return C._element_constructor(x, *args, **kwds) > 176 except Exception: > 177 if print_warnings: > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/polyhedron/parent.py > > in _element_constructor_(self, *args, **kwds) > 523 if convert and Vrep: > 524 Vrep = [convert_base_ring(_) for _ in Vrep] > --> 525 return self.element_class(self, Vrep, Hrep, **kwds) > 526 if nargs == 1 and is_Polyhedron(args[0]): > 527 polyhedron = args[0] > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/polyhedron/backend_cdd.py > > in __init__(self, parent, Vrep, Hrep, **kwds) > 459 sage: TestSuite(p).run() > 460 """ > --> 461 Polyhedron_cdd.__init__(self, parent, Vrep, Hrep, **kwds) > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/polyhedron/base.py > > in __init__(self, parent, Vrep, Hrep, **kwds) > 125 elif Hrep is not None: > 126 ieqs, eqns = Hrep > --> 127 self._init_from_Hrepresentation(ieqs, eqns, **kwds) > 128 else: > 129 self._init_empty_polyhedron() > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/polyhedron/backend_cdd.py > > in _init_from_Hrepresentation(self, ieqs, eqns, verbose) > 114 s = cdd_Hrepresentation(self._cdd_type, ieqs, eqns) > 115 s = self._run_cdd(s, '--redcheck', verbose=verbose) > --> 116 s = self._run_cdd(s, '--repall', verbose=verbose) > 117 self._init_from_cdd_output(s) > 118 if not self.base_ring().is_exact(): > > /opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/geometry/polyhedron/backend_cdd.py > > in _run_cdd(self, cdd_input_string, cmdline_arg, verbose) > 163 if 'Error:' in ans + err: > 164 # cdd reports errors on stdout and misc information on > stderr > --> 165 raise ValueError(ans.strip()) > 166 return ans > 167 > > ValueError: *Input Error: Input format is not correct. > *Format: > begin > m n NumberType(real, rational or integer) > b -A > end > > How can i fix this problem!? > > > > >
-- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-support/98085cf6-b482-465d-98df-9b3da0d052a0%40googlegroups.com.
