Ticket #1234 (closed defect: fixed)
[with patch; with positive review] analytic_rank crashes
|Reported by:||zimmerma||Owned by:||was|
Maybe the following is not feasible, but analytic_rank could crash smoothly.
sage: d=100032426715415089/251987961355200625 sage: E = EllipticCurve([0, -d^3+5*d^2, 0, -8*d^5+8*d^4, 4*d^8-8*d^7+4*d^6]) sage: F = E.minimal_model() sage: F.analytic_rank(algorithm='cremona') <type 'exceptions.RuntimeError'>: Error: ' *** elltors: precision too low in torsell. sage: F.analytic_rank(algorithm='rubinstein') <type 'exceptions.TypeError'>: unable to convert x (= 6.90579e+20 and is too large) to an integer sage: F.analytic_rank(algorithm='sympow') sympow 1.018 RELEASE (c) Mark Watkins --- see README and COPYING for details **ERROR** c4 invariant is too large
From John Cremona: the "cremona" version just wraps a gp script, which needs to have sufficient precision even for ellinit() to work ok here. Unfortunately this call will run in its own gp session, and it is not possible for the user to set the precision. (I found the same while running a lot of examples through Denis Simon's gp scripts). The solution is to change the wrapper to have a precision parameter, with some reasonable default for backwards compatibility, which gets passes through to gp.
- Summary changed from analytic_rank crashes to [with patch; needs review] analytic_rank crashes
- Summary changed from [with patch; needs review] analytic_rank crashes to [with patch; with positive review] analytic_rank crashes
comment:12 Changed 4 years ago by mabshoff
- Summary changed from [with patch; with positive review] analytic_rank crashes to [with patch; with positive review] analytic_rank crashes