Backtrace leads into cypari2 polredbest, possibly a pari bug:
sage: a=AA(sqrt(sqrt(5)))
....: r=AA(sqrt((AA(sqrt(13))-a)^2+3))
....: c=a+r
....:
....: d= AA(sqrt(r^2-a^2))
....:
....: 2*a*c == c^2 - d^2
....:
^C---------------------------------------------------------------------------
KeyboardInterrupt Traceback (most recent call last)
<ipython-input-3-9ddbec6c7736> in <module>()
5 d= AA(sqrt(r**Integer(2)-a**Integer(2)))
6
----> 7 Integer(2)*a*c == c**Integer(2) - d**Integer(2)
/home/vbraun/Code/sage.git/src/sage/structure/element.pyx in
sage.structure.element.Element.__richcmp__
(/home/vbraun/Sage/git/src/build/cythonized/sage/structure/element.c:10360)()
1121 # an instance of Element. The explicit casts below make
1122 # Cython generate optimized code for this call.
-> 1123 return (<Element>self)._richcmp_(other, op)
1124 else:
1125 return coercion_model.richcmp(self, other, op)
/home/vbraun/Code/sage.git/src/sage/structure/element.pyx in
sage.structure.element.Element._richcmp_
(/home/vbraun/Sage/git/src/build/cythonized/sage/structure/element.c:10559)()
1156 cdef int c
1157 try:
-> 1158 c = left._cmp_(right)
1159 except NotImplementedError:
1160 # Check equality by id(), knowing that left is not right
/home/vbraun/Code/sage.git/src/sage/structure/element.pyx in
sage.structure.element.Element._cmp_
(/home/vbraun/Sage/git/src/build/cythonized/sage/structure/element.c:10948)()
1172 left_cmp = left.__cmp__
1173 if isinstance(left_cmp, MethodType):
-> 1174 return left_cmp(right)
1175 msg = LazyFormat("comparison not implemented for
%r")%type(left)
1176 raise NotImplementedError(msg)
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in __cmp__(self, other)
4408 return -other.sign()
4409 else:
-> 4410 return self._sub_(other).sign()
4411
4412 def __pow__(self, e):
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in sign(self)
4737 # OK, we'll try adding precision one more time
4738 self._more_precision()
-> 4739 return self.sign()
4740 else:
4741 # Sigh...
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in sign(self)
4740 else:
4741 # Sigh...
-> 4742 self.exactify()
4743 return self.sign()
4744
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
3270 od = self._descr
3271 if isinstance(od, (ANRational, ANExtensionElement)): return
-> 3272 self._set_descr(self._descr.exactify())
3273
3274 def _set_descr(self, new_descr):
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
7078 right = self._right
7079 left.exactify()
-> 7080 right.exactify()
7081 gen = left._exact_field().union(right._exact_field())
7082 left_value = gen(left._exact_value())
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
3270 od = self._descr
3271 if isinstance(od, (ANRational, ANExtensionElement)): return
-> 3272 self._set_descr(self._descr.exactify())
3273
3274 def _set_descr(self, new_descr):
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
7078 right = self._right
7079 left.exactify()
-> 7080 right.exactify()
7081 gen = left._exact_field().union(right._exact_field())
7082 left_value = gen(left._exact_value())
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
3270 od = self._descr
3271 if isinstance(od, (ANRational, ANExtensionElement)): return
-> 3272 self._set_descr(self._descr.exactify())
3273
3274 def _set_descr(self, new_descr):
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
7077 left = self._left
7078 right = self._right
-> 7079 left.exactify()
7080 right.exactify()
7081 gen = left._exact_field().union(right._exact_field())
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
3270 od = self._descr
3271 if isinstance(od, (ANRational, ANExtensionElement)): return
-> 3272 self._set_descr(self._descr.exactify())
3273
3274 def _set_descr(self, new_descr):
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in exactify(self)
6142 newpol_sage_y = QQy(newpol_sage)
6143
-> 6144 red_elt, red_back, red_pol = do_polred(newpol_sage_y)
6145
6146 new_nf = NumberField(red_pol, name='a', check=False)
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/rings/qqbar.pyc
in do_polred(poly)
1712 (1/4*x, 4*x, x^4 - 268435456*x^2 + 211973662764908353025)
1713 """
-> 1714 new_poly, elt_back = poly._pari_().polredbest(flag=1)
1715
1716 parent = poly.parent()
/home/vbraun/Code/sage.git/local/lib/python2.7/site-packages/sage/libs/cypari2/auto_gen.pxi
in sage.libs.cypari2.gen.Gen_auto.polredbest
(/home/vbraun/Sage/git/src/build/cythonized/sage/libs/cypari2/gen.c:82077)()
16003 '''
16004 cdef GEN _T = T.g
> 16005 sig_on()
16006 cdef GEN _ret = polredbest(_T, flag)
16007 return new_gen(_ret)
src/cysignals/signals.pyx in cysignals.signals.sig_raise_exception
(build/src/cysignals/signals.c:1303)
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-devel+unsubscr...@googlegroups.com.
To post to this group, send email to sage-devel@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
