Opened 18 months ago
Last modified 13 months ago
#27015 new defect
bugs in integral_points for elliptic curves over QQ
Reported by: | cremona | Owned by: | |
---|---|---|---|
Priority: | major | Milestone: | |
Component: | elliptic curves | Keywords: | elliptic curves, integral points |
Cc: | Merged in: | ||
Authors: | John Cremona | Reviewers: | |
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | Stopgaps: |
Description
The code for computing all integral points on elliptic curves over QQ has bugs, some due to incorrect formulas in the literature. For example:
sage: E = EllipticCurve('2082a1') sage: E.integral_points() [(-11 : 29 : 1), (-2 : 29 : 1), (4 : 11 : 1), (13 : 29 : 1)] sage: E([507525709, 11433453531221, 1]) (507525709 : 11433453531221 : 1)
This curve has rank 1 and trivial torsion, generator P=(13,-43) and the integral points are k*P for k=1,2,3,4 and 13 (and their negatives). We miss 13*P because the upper bound computed is only 12 (as seen with verbose=True).
#10973 implements integral points on elliptic curves over arbitrary number fields, and also provides a fix for this bug since it work over QQ (the branch at #10973 has the integral_points() method over QQ call the general one). However the code there is slow over QQ (and relies on some incorrect though published assumptions). Hence there is a case for fixing the bugs over QQ independently of the general implementation, and that will be done here.
Change History (3)
comment:1 Changed 18 months ago by
comment:2 Changed 16 months ago by
- Milestone changed from sage-8.7 to sage-8.8
Ticket retargeted after milestone closed (if you don't believe this ticket is appropriate for the Sage 8.8 release please retarget manually)
comment:3 Changed 13 months ago by
- Milestone sage-8.8 deleted
As the Sage-8.8 release milestone is pending, we should delete the sage-8.8 milestone for tickets that are not actively being worked on or that still require significant work to move forward. If you feel that this ticket should be included in the next Sage release at the soonest please set its milestone to the next release milestone (sage-8.9).
The complete list of missed points with Sage 8.5, for curves in database_cremona_ellcurve (i.e. all curves of conductor<400000) is: 12 points on 11 curves:
missing 1 point in each case except 2 points on 1848f1.