#13054: Pernicious bug for algebraic numbers
---------------------------------+------------------------------------------
Reporter: rbeezer | Owner: davidloeffler
Type: defect | Status: new
Priority: major | Milestone: sage-5.1
Component: number fields | Resolution:
Keywords: sd40.5 | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
---------------------------------+------------------------------------------
Comment (by rbeezer):
This is not the same bug as above, but could be another symptom. Still
hunting for an example to grab onto.
Run
{{{
n = 3
A = matrix(QQbar, n, n, [i/(i+1) + (i^2+3)*I for i in range(n^2)])
ev = A.eigenvalues()
(A-ev[0]).rref()
}}}
which should hangup and when killed with Ctrl-C yields
{{{
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "_sage_input_5.py", line 10, in <module>
exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8
-*-\\n" +
_support_.preparse_worksheet_cell(base64.b64decode("biA9IDMKQSA9IG1hdHJpeChRUWJhciwgbiwgbiwgW2kvKGkrMSkgKyAoaV4yKzMpKkkgZm9yIGkgaW4gcmFuZ2Uobl4yKV0pCmV2ID0gQS5laWdlbnZhbHVlcygpCihBLWV2WzBdKS5ycmVmKCk="),globals())+"\\n");
execfile(os.path.abspath("___code___.py"))
File "", line 1, in <module>
File "/tmp/tmpiYGsKk/___code___.py", line 6, in <module>
exec compile(u'(A-ev[_sage_const_0 ]).rref()
File "", line 1, in <module>
File "matrix2.pyx", line 5526, in sage.matrix.matrix2.Matrix.rref
(sage/matrix/matrix2.c:28497)
File "matrix2.pyx", line 5839, in
sage.matrix.matrix2.Matrix.echelon_form (sage/matrix/matrix2.c:29882)
File "matrix2.pyx", line 5745, in sage.matrix.matrix2.Matrix.echelonize
(sage/matrix/matrix2.c:29429)
File "matrix2.pyx", line 5897, in
sage.matrix.matrix2.Matrix._echelon_in_place_classical
(sage/matrix/matrix2.c:30399)
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 3831, in __nonzero__
return self.__nonzero__()
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 3834, in __nonzero__
self.exactify()
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 3466, in exactify
self._set_descr(self._descr.exactify())
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 7591, in exactify
left.exactify()
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 3466, in exactify
self._set_descr(self._descr.exactify())
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 7591, in exactify
left.exactify()
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 3466, in exactify
self._set_descr(self._descr.exactify())
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 7593, in exactify
gen = left._exact_field().union(right._exact_field())
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 2251, in union
pari_nf = self.pari_field()
File "/sage/sage-5.1.beta0/local/lib/python2.7/site-
packages/sage/rings/qqbar.py", line 2145, in pari_field
self._pari_field = pari_pol.nfinit(1)
File "gen.pyx", line 7470, in sage.libs.pari.gen.gen.nfinit
(sage/libs/pari/gen.c:31383)
File "gen.pyx", line 7484, in sage.libs.pari.gen.gen._nfinit_with_prec
(sage/libs/pari/gen.c:31587)
KeyboardInterrupt
__SAGE__
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13054#comment:1>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
