The following code is fine if 10^17 is replaced by (say) 10^16, so the
underlying cause is probably due to being close to pari's limitations.
But the error message is rather unacceptable! Here gp is yielding an
error message (see end)
[EMAIL PROTECTED]
----------------------------------------------------------------------
| SAGE Version 2.10.2, Release Date: 2008-02-22 |
| Type notebook() for the GUI, and license() for information. |
----------------------------------------------------------------------
sage: p=next_prime(10^17); p
100000000000000003
sage: E=EllipticCurve(GF(p),[123,324])
sage: P=E(53989543867592054, 16386953280180785)
sage: P.order()
---------------------------------------------------------------------------
<type 'exceptions.TypeError'> Traceback (most recent call last)
/home/jec/sage-2.10.2/<ipython console> in <module>()
/home/jec/sage-2.10.2/local/lib/python2.5/site-packages/sage/schemes/elliptic_curves/ell_point.py
in order(self, disable_warning)
633 lb,ub=ell_generic.Hasse_bounds(q)
634 if K.is_prime_field():
--> 635 e = E._gp()
636 N = rings.Integer(e.ellzppointorder(list(self.xy())))
637 else:
/home/jec/sage-2.10.2/local/lib/python2.5/site-packages/sage/schemes/elliptic_curves/ell_finite_field.py
in _gp(self)
126 if not F.is_prime_field():
127 raise NotImplementedError
--> 128 self.__gp = gp_cremona.ellinit(self.a_invariants(),
F.characteristic())
129 return self.__gp
130
/home/jec/sage-2.10.2/local/lib/python2.5/site-packages/sage/schemes/elliptic_curves/gp_cremona.py
in ellinit(e, p)
74 """
75 init()
---> 76 return gp("e=ellzpinit(%s,%s);"%(e,p))
77
78
/home/jec/sage-2.10.2/local/lib/python2.5/site-packages/sage/interfaces/expect.py
in __call__(self, x)
736 return x
737 if isinstance(x, basestring):
--> 738 return cls(self, x)
739 try:
740 return self._coerce_from_special_method(x)
/home/jec/sage-2.10.2/local/lib/python2.5/site-packages/sage/interfaces/expect.py
in __init__(self, parent, value, is_name)
1005 except (TypeError, KeyboardInterrupt,
RuntimeError, ValueError), x:
1006 self._session_number = -1
-> 1007 raise TypeError, x
1008 self._session_number = parent._session_number
1009
<type 'exceptions.TypeError'>: Error executing code in GP/PARI:
CODE:
sage[1]=e=ellzpinit([0, 0, 0, 123, 324],100000000000000003);;
GP/PARI ERROR:
*** error opening input file: /home/jec/.sage//tem
^--------------------
------------
I think gp is showing an error message like this:
*** vector: the PARI stack overflows !
current stack size: 10000000 (9.537 Mbytes)
[hint] you can increase GP stack with allocatemem()
I will be writing around this anyway since the new code now in 2.10.2
for elliptic curves over finite fields means that we will no longer
need to use my old gp scripts for any of this (ellzpinit() is not part
of standard gp, but is part of the external scripts included in Sage).
I guess this should be another trac ticket?
John
--
John Cremona
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-devel@googlegroups.com
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at http://groups.google.com/group/sage-devel
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
