On Mon, Feb 25, 2019 at 4:04 PM Vincent Delecroix <20100.delecr...@gmail.com> wrote: > > I can not reproduce the error. What is the configuration of > your machine? > > The comparison routine in QQbar that is called > (cmp_elements_with_same_minpoly) assumes that the elements > a and b have the same minimal polynomial p. But the error message > is weird as it does not seem to be the case. The only roots of p > that looks like a and b are the following pair of conjugates > > 2.917983368234029?e15 - 1.630474617080222?e18*I > 2.917983368234029?e15 + 1.630474617080222?e18*I > > But the imaginary part is non-trivial! It seems that a and b > that appears in the error are here the *real part* of some roots > but not roots. > > It would be helpful if you could provide more data (ie via pdb or > adding prints in qqbar.py). More precisely at lines 3367 (qqbar.py) > the values and types of > * sd._left._descr > * sd._right._descr > * sd._op > and also at lines 3372 the values of the polynomials > * lp > * rp > > (the datatype of an element in QQbar is mostly a tree and > in the current situation sd._left and sd._right are the > two subtrees) > > Ideally, it should be possible to reconstruct the element that > triggers the error... you can try to print a sage_input(self) > in the __bool__ method (line 3367).
Thank you for taking the time to look into this, and I apologize for that because I just found that it was a false alarm: I apparently somehow (despite just having rebuilt it) wound up with a partially-broken PARI, at least I think. After a manual `rm -f local/lib/*pari*` and rebuilding pari, cypari, and their dependents everything is fixed now. Of course, that could indicate a bug in the build system, but I'm not sure now how I would reproduce the situation. This arose in part because I was playing around with https://trac.sagemath.org/ticket/27267 > Le 25/02/2019 à 15:39, E. Madison Bray a écrit : > > On Mon, Feb 25, 2019 at 3:36 PM E. Madison Bray <erik.m.b...@gmail.com> > > wrote: > >> > >> On Ubuntu 14.04 I am getting a strange error when trying to build the > >> docs. Tried `make doc-clean` to no avail: > > > > Also, when I try running the tests for > > sage.categories.finite_coxeter_groups I get the same failure at: > > > > sage -t src/sage/categories/finite_coxeter_groups.py > > ********************************************************************** > > File "src/sage/categories/finite_coxeter_groups.py", line 725, in > > sage.categories.finite_coxeter_groups.FiniteCoxeterGroups.ParentMethods.permutahedron > > Failed example: > > W.permutahedron() > > > > Followed with a long traceback, the relevant portion being: > > > > W.permutahedron() > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/categories/finite_coxeter_groups.py", > > line 765, in permutahedron > > return Polyhedron(vertices=vertices, base_ring=base_ring) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/constructor.py", > > line 600, in Polyhedron > > return parent(Vrep, Hrep, convert=convert, verbose=verbose) > > File "sage/structure/parent.pyx", line 901, in > > sage.structure.parent.Parent.__call__ > > (build/cythonized/sage/structure/parent.c:9205) > > return mor._call_with_args(x, args, kwds) > > File "sage/structure/coerce_maps.pyx", line 181, in > > sage.structure.coerce_maps.DefaultConvertMap_unique._call_with_args > > (build/cythonized/sage/structure/coerce_maps.c:5076) > > raise > > File "sage/structure/coerce_maps.pyx", line 176, in > > sage.structure.coerce_maps.DefaultConvertMap_unique._call_with_args > > (build/cythonized/sage/structure/coerce_maps.c:4964) > > return C._element_constructor(x, *args, **kwds) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/parent.py", > > line 517, in _element_constructor_ > > return self.element_class(self, Vrep, Hrep, **kwds) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/backend_field.py", > > line 177, in __init__ > > super(Polyhedron_field, self).__init__(parent, Vrep, Hrep, **kwds) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/base.py", > > line 124, in __init__ > > self._init_from_Vrepresentation(vertices, rays, lines, **kwds) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/backend_field.py", > > line 208, in _init_from_Vrepresentation > > H = Vrep2Hrep(self.base_ring(), self.ambient_dim(), vertices, > > rays, lines) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description_inhomogeneous.py", > > line 446, in __init__ > > DD = Algorithm(A).run() > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > > line 764, in run > > DD.add_inequality(a) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > > line 715, in add_inequality > > if not self.are_adjacent(rp, rn): > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > > line 447, in are_adjacent > > Z = self.zero_set(r1).intersection(self.zero_set(r2)) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > > line 375, in zero_set > > t.update(self.A[i] for i in range(n,len(self.A)) if > > self.A[i].inner_product(ray) == self.zero) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > > line 375, in <genexpr> > > t.update(self.A[i] for i in range(n,len(self.A)) if > > self.A[i].inner_product(ray) == self.zero) > > File "sage/structure/element.pyx", line 1091, in > > sage.structure.element.Element.__richcmp__ > > (build/cythonized/sage/structure/element.c:9923) > > return (<Element>self)._richcmp_(other, op) > > File "sage/structure/element.pyx", line 1095, in > > sage.structure.element.Element._richcmp_ > > (build/cythonized/sage/structure/element.c:10030) > > cpdef _richcmp_(left, right, int op): > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/rings/qqbar.py", > > line 4641, in _richcmp_ > > return bool(self) == (op == op_NE) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/rings/qqbar.py", > > line 3376, in __bool__ > > c = cmp_elements_with_same_minpoly(left, right, left.minpoly()) > > File > > "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/rings/qqbar.py", > > line 2225, in cmp_elements_with_same_minpoly > > raise RuntimeError('a = {}\nb = {}\np = {}'.format(a, b, p)) > > RuntimeError: a = 2.917983368234029?e15 > > b = 2.917983368234029?e15 > > p = x^6 - 5835966736468056*x^5 + > > 2658455991569831813923974927379870544*x^4 + > > 2658455991569831814199401327502424704*x^3 + > > 2658455991569831809001982851260652032*x^2 + 253512790118400*x + 4096 > > > > Any ideas? > > > > > > > >> [dochtml] [categorie] > >> /home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/categories/finite_coxeter_groups.py:docstring > >> of > >> sage.categories.finite_coxeter_groups.FiniteCoxeterGroups.ParentMethods.permutahedron:50: > >> WARNING: Exception occurred in plotting finite_coxeter_groups-1 > >> [dochtml] [categorie] from > >> /home/embray/src/sagemath/sage/src/doc/en/reference/categories/sage/categories/finite_coxeter_groups.rst: > >> [dochtml] [categorie] Traceback (most recent call last): > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/matplotlib/sphinxext/plot_directive.py", > >> line 524, in run_code > >> [dochtml] [categorie] six.exec_(code, ns) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/six.py", > >> line 709, in exec_ > >> [dochtml] [categorie] exec("""exec _code_ in _globs_, _locs_""") > >> [dochtml] [categorie] File "<string>", line 1, in <module> > >> [dochtml] [categorie] File "<string>", line 2, in <module> > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/categories/finite_coxeter_groups.py", > >> line 765, in permutahedron > >> [dochtml] [categorie] return Polyhedron(vertices=vertices, > >> base_ring=base_ring) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/constructor.py", > >> line 600, in Polyhedron > >> [dochtml] [categorie] return parent(Vrep, Hrep, convert=convert, > >> verbose=verbose) > >> [dochtml] [categorie] File "sage/structure/parent.pyx", line 901, in > >> sage.structure.parent.Parent.__call__ > >> (build/cythonized/sage/structure/parent.c:9205) > >> [dochtml] [categorie] return mor._call_with_args(x, args, kwds) > >> [dochtml] [categorie] File "sage/structure/coerce_maps.pyx", line > >> 181, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_with_args > >> (build/cythonized/sage/structure/coerce_maps.c:5076) > >> [dochtml] [categorie] raise > >> [dochtml] [categorie] File "sage/structure/coerce_maps.pyx", line > >> 176, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_with_args > >> (build/cythonized/sage/structure/coerce_maps.c:4964) > >> [dochtml] [categorie] return C._element_constructor(x, *args, **kwds) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/parent.py", > >> line 517, in _element_constructor_ > >> [dochtml] [categorie] return self.element_class(self, Vrep, Hrep, > >> **kwds) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/backend_field.py", > >> line 177, in __init__ > >> [dochtml] [categorie] super(Polyhedron_field, > >> self).__init__(parent, Vrep, Hrep, **kwds) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/base.py", > >> line 124, in __init__ > >> [dochtml] [categorie] self._init_from_Vrepresentation(vertices, > >> rays, lines, **kwds) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/backend_field.py", > >> line 208, in _init_from_Vrepresentation > >> [dochtml] [categorie] H = Vrep2Hrep(self.base_ring(), > >> self.ambient_dim(), vertices, rays, lines) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description_inhomogeneous.py", > >> line 446, in __init__ > >> [dochtml] [categorie] DD = Algorithm(A).run() > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > >> line 764, in run > >> [dochtml] [categorie] DD.add_inequality(a) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > >> line 715, in add_inequality > >> [dochtml] [categorie] if not self.are_adjacent(rp, rn): > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > >> line 447, in are_adjacent > >> [dochtml] [categorie] Z = > >> self.zero_set(r1).intersection(self.zero_set(r2)) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > >> line 375, in zero_set > >> [dochtml] [categorie] t.update(self.A[i] for i in > >> range(n,len(self.A)) if self.A[i].inner_product(ray) == self.zero) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/geometry/polyhedron/double_description.py", > >> line 375, in <genexpr> > >> [dochtml] [categorie] t.update(self.A[i] for i in > >> range(n,len(self.A)) if self.A[i].inner_product(ray) == self.zero) > >> [dochtml] [categorie] File "sage/structure/element.pyx", line 1091, > >> in sage.structure.element.Element.__richcmp__ > >> (build/cythonized/sage/structure/element.c:9923) > >> [dochtml] [categorie] return (<Element>self)._richcmp_(other, op) > >> [dochtml] [categorie] File "sage/structure/element.pyx", line 1095, > >> in sage.structure.element.Element._richcmp_ > >> (build/cythonized/sage/structure/element.c:10030) > >> [dochtml] [categorie] cpdef _richcmp_(left, right, int op): > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/rings/qqbar.py", > >> line 4641, in _richcmp_ > >> [dochtml] [categorie] return bool(self) == (op == op_NE) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/rings/qqbar.py", > >> line 3376, in __bool__ > >> [dochtml] [categorie] c = cmp_elements_with_same_minpoly(left, > >> right, left.minpoly()) > >> [dochtml] [categorie] File > >> "/home/embray/src/sagemath/sage/local/lib/python2.7/site-packages/sage/rings/qqbar.py", > >> line 2225, in cmp_elements_with_same_minpoly > >> [dochtml] [categorie] raise RuntimeError('a = {}\nb = {}\np = > >> {}'.format(a, b, p)) > >> [dochtml] [categorie] RuntimeError: a = 2.917983368234029?e15 > >> [dochtml] [categorie] b = 2.917983368234029?e15 > >> [dochtml] [categorie] p = x^6 - 5835966736468056*x^5 + > >> 2658455991569831813923974927379870544*x^4 + > >> 2658455991569831814199401327502424704*x^3 + > >> 2658455991569831809001982851260652032*x^2 + 253512790118400*x + 4096 > >> > >> On Mon, Feb 25, 2019 at 7:21 AM 'Justin C. Walker' via sage-release > >> <sage-release@googlegroups.com> wrote: > >>> > >>> > >>> > >>>> On Feb 24, 2019, at 02:39 , Volker Braun <vbraun.n...@gmail.com> wrote: > >>>> > >>>> As always, you can get the latest beta version from the "develop" git > >>>> branch. Alternatively, the self-contained source tarball is at > >>>> http://www.sagemath.org/download-latest.html > >>> > >>> Built w/o problems on three macOS platforms from a fresh clone of the > >>> develop tree. > >>> > >>> 10.11.6 (mid-2015 MBP, Quad-core Core i7): all tests passed! > >>> > >>> 10.13.6 (2017 iMac Pro, 18-core Xeon W): all tests passed! > >>> However, repeating the “loop” testing of polynomials.pyx showed the same > >>> intermittent failures as before (this failure occurred in 4 of 8 tests). > >>> > >>> 10.14.2 (2017 MBP, Quad-core Core i7): one testing failure: > >>> sage -t --long --warn-long 77.8 > >>> src/sage/numerical/backends/logging_backend.py # Killed due to abort > >>> > >>> When run by itself the test passed. I then reran the full test cycle, > >>> with no failures (other than a batch of timeouts, due to my being unable > >>> to keep my MBP from going to sleep at intervals). > >>> > >>> Log for the logging_bsckend failure attached. > >>> > >>> -- > >>> You received this message because you are subscribed to the Google Groups > >>> "sage-release" group. > >>> To unsubscribe from this group and stop receiving emails from it, send an > >>> email to sage-release+unsubscr...@googlegroups.com. > >>> To post to this group, send email to sage-release@googlegroups.com. > >>> Visit this group at https://groups.google.com/group/sage-release. > >>> For more options, visit https://groups.google.com/d/optout. > >>> > >>> > >>> -- > >>> Justin C. Walker, Curmudgeon-At-Large > >>> Institute for the Enhancement of the Director's Income > >>> -------- > >>> Experience is what you get > >>> when you don't get what you want. > >>> -------- > >>> > >>> > >>> > >>> -- > >>> You received this message because you are subscribed to the Google Groups > >>> "sage-release" group. > >>> To unsubscribe from this group and stop receiving emails from it, send an > >>> email to sage-release+unsubscr...@googlegroups.com. > >>> To post to this group, send email to sage-release@googlegroups.com. > >>> Visit this group at https://groups.google.com/group/sage-release. > >>> For more options, visit https://groups.google.com/d/optout. > > > > -- > You received this message because you are subscribed to the Google Groups > "sage-release" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to sage-release+unsubscr...@googlegroups.com. > To post to this group, send email to sage-release@googlegroups.com. > Visit this group at https://groups.google.com/group/sage-release. > For more options, visit https://groups.google.com/d/optout. -- You received this message because you are subscribed to the Google Groups "sage-release" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-release+unsubscr...@googlegroups.com. To post to this group, send email to sage-release@googlegroups.com. Visit this group at https://groups.google.com/group/sage-release. For more options, visit https://groups.google.com/d/optout.
