#4805: S_integral points failure (possible p-adic precision problem)
---------------------------+------------------------------------------------
Reporter: cremona | Owner: was
Type: defect | Status: new
Priority: major | Milestone: sage-3.3
Component: number theory | Keywords: elliptic curve
---------------------------+------------------------------------------------
After #4741 we still have this problem:
{{{
sage: EllipticCurve("7690e1").S_integral_points([13,2])
ERROR: An unexpected error occurred while tokenizing input
The following traceback may be corrupted or invalid
The error message is: ('EOF in multi-line statement', (178, 0))
---------------------------------------------------------------------------
IndexError Traceback (most recent call
last)
/home/john/<ipython console> in <module>()
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/schemes/elliptic_curves/ell_rational_field.pyc in
S_integral_points(self, S, mw_base, both_signs, verbose, proof)
5062 while not p_prec_ok:
5063 try:
-> 5064
mw_base_p_log.append([mp_temp*(pts.padic_elliptic_logarithm(p,precision=p_prec))
for pts in mw_base])
5065 p_prec_ok=True
5066 except (PrecisionError, ZeroDivisionError,
TypeError):
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/schemes/elliptic_curves/ell_point.pyc in
padic_elliptic_logarithm(self, p, precision)
1399 t = -x/y
1400 v = t.valuation()
-> 1401 phi = Ep.formal().log(prec=1+precision//v)
1402 return phi(t)/f
1403
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/schemes/elliptic_curves/formal_group.pyc in log(self, prec)
348 -- David Harvey (2006-09-10): rewrote to use
differential
349 """
--> 350 return
self.differential(prec-1).integral().add_bigoh(prec)
351
352 def inverse(self, prec=20):
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/rings/power_series_ring_element.so in
sage.rings.power_series_ring_element.PowerSeries.add_bigoh
(sage/rings/power_series_ring_element.c:5413)()
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/rings/power_series_ring.pyc in __call__(self, f, prec,
check)
324 v = sage_eval(f.Eltseq())
325 return self(v) * (self.gen(0)**f.Valuation())
--> 326 return self.__power_series_class(self, f, prec,
check=check)
327
328 def construction(self):
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/rings/power_series_poly.so in
sage.rings.power_series_poly.PowerSeries_poly.__init__
(sage/rings/power_series_poly.c:2160)()
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/rings/polynomial/polynomial_element.so in
sage.rings.polynomial.polynomial_element.Polynomial.truncate
(sage/rings/polynomial/polynomial_element.c:27152)()
/home/john/sage-3.2.2.alpha1/local/lib/python2.5/site-
packages/sage/rings/polynomial/padics/polynomial_padic_capped_relative_dense.pyc
in __getslice__(self, i, j)
372 j = len(self._relprecs) + j
373 if j < 0:
--> 374 raise IndexError, "list index out of range"
375 if i >= j:
376 return
Polynomial_padic_capped_relative_dense(self.parent(), [])
IndexError: list index out of range
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4805>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
