Re: [sage-release] Sage 8.7.beta5 released
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 > > 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 =
Re: [sage-release] Sage 8.7.beta5 released
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). Vincent Le 25/02/2019 à 15:39, E. Madison Bray a écrit : On Mon, Feb 25, 2019 at 3:36 PM E. Madison Bray 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
Re: [sage-release] Sage 8.7.beta5 released
On Mon, Feb 25, 2019 at 3:36 PM E. Madison Bray 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
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 (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 +
Re: [sage-release] Sage 8.7.beta5 released
On Ubuntu 14.04 I am getting a strange error when trying to build the
docs. Tried `make doc-clean` to no avail:
[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 "", line 1, in
[dochtml] [categorie] File "", line 2, in
[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
[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]
