7122
Comment:
|
12762
minor edits
|
Deletions are marked like this. | Additions are marked like this. |
Line 9: | Line 9: |
* Compute regulators of elliptic curves over function fields | * Compute regulators of elliptic curves over function fields: * This worksheet does it using Tate's suggestion (i.e., use the definition): http://sagenb.org/home/pub/1198/ |
Line 11: | Line 12: |
* Is there an algorithm to enumerate all elliptic curves over a function field of a given conductor? | * Is there an algorithm to enumerate all elliptic curves over a function field of a given conductor? * Victor Miller reports that there are some papers that give efficient algorithms for solving "S-unit equations" over function fields, which seems relevant. It's well known that finding all Elliptic Curves with good reduction outside a finite set, $S$ of places is equivalent to solving S-unit equations. There are a series of four papers by Gaal and Pohst that give efficient algorithms for this when dealing with function fields over finite fields, especially $\mathbb{F}_q(t)$. [[http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6WKD-4HVF15H-2&_user=10&_rdoc=1&_fmt=&_orig=search&_sort=d&_docanchor=&view=c&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=ec0b143bb1c45f3bdcd4bea8841e8735|Diophantine Equations over Global Function Fields I: The Thue Equation]] |
Line 13: | Line 16: |
* Implement Tate's algorithm for elliptic curves over the function field F_p(t) | [[http://akpeters.metapress.com/content/458402lu66634164/|Diophantine Equations over Global Function Fields II: R-Integral Solutions of Thue Equations]] [[http://projecteuclid.org/DPubS?verb=Display&version=1.0&service=UI&handle=euclid.facm/1229696556&page=record|Diophantine equations over global function fields III: An application to resultant form equations]] [[http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6WKD-4TP7HDM-3&_user=10&_rdoc=1&_fmt=&_orig=search&_sort=d&_docanchor=&view=c&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=c50fd011b673de5511d49b10b3fe9935|Diophantine equations over global function fields IV: S-unit equations in several variables with an application to norm form equations]] * Implement Tate's algorithm for elliptic curves over the function field $\mathbb{F}_p(t)$. |
Line 17: | Line 26: |
* Robert Bradshaw implemented faster arithmetic for $F_p(t)$: [[http://trac.sagemath.org/sage_trac/ticket/7585|#7585]] | |
Line 18: | Line 28: |
People: ''Sal Baig'', William Stein, David Roe, Ken Ribet, Kiran Kedlaya, Robert Bradshaw, Victor Miller (S-unit equations) | * Let $K=\mathbb{F}_5(t)$ and $E/K: y^2 = x(x+t)(x+t^2)$ be the Legendre curve and consider its twist by $f=1+t$, $E_f$. $E_f$ has Weierstrass equation $$ y^2 = x^3 + \left(3 t^{6} + 3 t^{5} + 3 t^{3} + 3 t^{2}\right)x +\left(4 t^{9} + t^{8} + 3 t^{7} + 2 t^{6} + 3 t^{5} + t^{4} + 4 t^{3}\right). $$ The L-funtion of E_f is: $$ L(E_f/K,T) = 1-5T-25T^2+125T^3.$$ The analytic rank is 2 with $L^{(2)}(E_f/K,1/5) = 100$. The points $P=[3*t^3 + t^2 + 4*t: t^3 + 3*t^2:1]$, $Q=[3*t^3 + 2*t^2 + 3*t: t^3:1]$ are independent on $E_F(K)$, so its algebraic rank is 2. Thus the refined BSD conjecture is true. We do not if $P,Q$ generate $E_f(K)$ mod torsion, so the determinant of their height matrix (=5) gives an upper bound on the regulator of $E_f/K$. The Tamagawa number was computed by hand to be $2^7/5$. The torsion subgroup has size 4 (all 2-torsion). Thus $$ \frac{100}{2} = \frac{L^{(r)}(E_f/K,1/q)}{r!} = \frac{|Sha|R\tau}{|E(K)_{tors}|^2}\leq\frac{|Sha|\cdot 5\cdot 2^7\cdot 5^{-1}}{2^4}.$$ Note that there is probably a factor of 2 going unaccounted for in the computation of the regulator. From this, we see that $|Sha|\geq \frac{25}{4}$ (again, this may be off by a power of 2). People: ''Sal Baig'', William Stein, David Roe, Ken Ribet, Kiran Kedlaya, Robert Bradshaw, Victor Miller (S-unit equations), Thomas Barnet-Lamb |
Line 27: | Line 51: |
* Implement the algorithm in the reducible case due to Jeechul Woo. | * Implement the algorithm in the reducible case due to Jeechul Woo: [[attachment:ThreeDescent.gp|Jeechul Woo's GP script]] |
Line 31: | Line 55: |
People: ''Robert Miller'', William Stein, Victor Miller, Jeechul Woo (Noam's student; around only Thu, Fri) | People: ''Robert Miller'', William Stein, Victor Miller, Jeechul Woo (Noam's student; around only Thu, Fri) |
Line 38: | Line 62: |
People: ''William Stein'', Dimitar Jetchev, Drew Sutherland, Mirela Ciperiani, Ken Ribet, Victor Miller | * Stein: I posted a bundle based against Sage-4.2.1 here (called heegner-4.2.1.hg), which has highly relevant code: http://trac.sagemath.org/sage_trac/ticket/6616 People: ''William Stein'', Dimitar Jetchev, Drew Sutherland, Mirela Ciperiani, Ken Ribet, Victor Miller |
Line 58: | Line 84: |
People: ''Drew Sutherland'', Ken Ribet, William Stein, Kiran Kedlaya, David Roe | People: ''Drew Sutherland'', Ken Ribet, William Stein, Kiran Kedlaya, David Roe |
Line 65: | Line 91: |
People: ''William Stein'', Robert Bradshaw, Jen Balakrishnan | People: ''William Stein'', Robert Bradshaw, Jen Balakrishnan |
Line 70: | Line 96: |
This will be based on Jared Weinstein's talk. | This will be based on Jared Weinstein's talk. See http://wstein.org/misc/sagedays18_papers/weinstein-kolyvagin_classes_for_higher_rank_elliptic_curves.pdf |
Line 72: | Line 98: |
People: ''Jared Weinstein'', Mirela Ciperiani, William Stein | People: ''Jared Weinstein'', Mirela Ciperiani, William Stein |
Line 78: | Line 104: |
* Figure out exactly what needs to be computed and what might be an optimal curve and quadratic imaginary field to work with. | * Figure out exactly what needs to be computed and what might be an optimal curve and quadratic imaginary field to work with: [[/5077aestimates|Some details for 5077a]] |
Line 84: | Line 110: |
People: ''William Stein'', Dimitar Jetchev, Victor Miller (sparse linear algebra), Jen Balakrishnan | People: ''William Stein'', Dimitar Jetchev, Victor Miller (sparse linear algebra), Jen Balakrishnan, Robert Bradshaw |
Line 93: | Line 119: |
People: ''Matthew Greenberg'', Cameron Frank, Kiran Kedlaya, Robert Pollack, Avner Ash, David Roe, Jay Pottharst | People: ''Matthew Greenberg'', Cameron Frank, Kiran Kedlaya, Robert Pollack, Avner Ash, David Roe, Jay Pottharst, Thomas Barnet-Lamb |
Line 112: | Line 138: |
On the one hand, one of Kolyvagin's suite of conjectures is that the classes he constructs out of Heegner points ought to generate the entirety of Sha(E/K). On the other hand each element of Sha(E/K) is "visible" in some other abelian variety, in the sense of Mazur (see for instance http://www.williamstein.org/home/was/www/home/wstein/www/papers/visibility_of_sha/). I wonder if the Kolyvagin classes d(n) contributing to Sha become visible in abelian varieties in some *systematic* way; ie, in a modular Jacobian of possibly higher level. There are already some examples out there of computation with visibility and computation with Kolyvagin classes, so maybe we can gather some data. | * Jared describes the problem thus: On the one hand, one of Kolyvagin's suite of conjectures is that the classes he constructs out of Heegner points ought to generate the entirety of Sha(E/K). On the other hand each element of Sha(E/K) is "visible" in some other abelian variety, in the sense of Mazur (see for instance http://www.williamstein.org/home/was/www/home/wstein/www/papers/visibility_of_sha/). I wonder if the Kolyvagin classes d(n) contributing to Sha become visible in abelian varieties in some *systematic* way; ie, in a modular Jacobian of possibly higher level. There are already some examples out there of computation with visibility and computation with Kolyvagin classes, so maybe we can gather some data. |
Line 114: | Line 140: |
* I wrote a short paper with Dimitar Jetchev in 2005 (which I forgot about until just now, and never published), which I think was motivated by this question: [[attachment::jetchev-stein-congruences_and_unramified_cohomology.pdf]] |
|
Line 119: | Line 146: |
Consider the Heegner $y_{p^n}$ over the anticyclotomic tower for a rank >= 2 curve, Sha trivial, etc. | * Consider the Heegner $y_{p^n}$ over the anticyclotomic tower for a rank >= 2 curve, Sha trivial, etc. * Barry Mazur remarks that this may be connected to his notion (with Rubin) of "Shadow lines" in Mordell-Weil. |
Line 127: | Line 156: |
* might use David Harvey's super fast Kedlaya zeta function implementation in Sage (his C++ code) | * might use David Harvey's super fast Kedlaya zeta function implementation in Sage (his C++ code), but unfortunately the curves that typically arise aren't hyperelliptic, and Harvey's code only applies to hyperelliptic curves. |
Line 129: | Line 158: |
* could lead to questions of an "anabelian sort"... e.g., take elliptic curve over finite field, take n-torsion points, extract roots, get coverings, etc. Compute Frobenius eigenvalues of all these covers. | * could lead to questions of an "anabelian sort"... e.g., take elliptic curve over finite field, take n-torsion points, extract roots, get coverings, etc. Compute Frobenius eigenvalues of all these covers. * To do a useful computation, one needs to compute Frobenius, say for the curve $y^{\ell} = f_N(x)$, where $f_N(x)$ is the $N$-th division polynomial on an elliptic curve. This following paper is about how to count points on superelliptic curves: [[attachment:gaudry-gurel-an_extension_of_Kedlayas_point-counting_algorithm_to_superelliptic_curves.pdf|Gaudry-Gurel]] More from Barry: " I 've thought about my off-the-cuff suggestion for computing Frobenius eigenvalues, and my worry is that one won't really see anything terribly interesting if one works only with N=2; but maybe when one works with N=3. For example, take an elliptic curve E over F_q and let f_3 be (``3-division") function on E: meaning that it vanishes to order one at all nontrivial 3-torsion points and has the appropriate order pole at the origin (and no other poles or zeroes). Then (for small n) extract n-th powers of that f_3 to get curves C_n mapping to E (of unfortunately high genus). But it could be that the collection of their Frobenius eigenvalues tells us more than just the isogeny class of E? From what I learn by talking with William, this might be a very difficult computational problem though..." * There is a better paper maybe about this super-elliptic algorithm, with a Magma implementation here: http://www.math.tu-berlin.de/~minzlaff/. And here is a Sage worksheet that *wraps* that Magma code, so if you have Magma you can compute quickly the zeta function of $y^\ell = f(x)$ for f of odd degree. And indeed Minzlaff's code seems to work fine. * Kiran has a short Sage program to compute the zeta function of y^4 = (cubic), to be published soon. |
Sage Days 18 Coding Sprint Projects
Contents
-
Sage Days 18 Coding Sprint Projects
- Elliptic curves over function fields
- Implement computation of the 3-Selmer rank of an elliptic curve over QQ
- Compute statistics about distribution of Heegner divisors and Kolyvagin divisors modulo primes p
- Create a table of images of Galois representations, for elliptic curves
- Fast computation of Heegner points
- Implement code to compute the asymptotic distribution of Kolyvagin classes
- Verify Kolyvagin's conjecture for a specific rank 3 curve
- Implement an algorithm in Sage to compute Stark-Heegner points
- Compute the higher Heegner point y_5 on the curve 389a provably correctly
- Compute a Heegner point on the Jacobian of a genus 2 curve
- Visibility of Kolyvagin Classes
- Find an algorithm to decide if y_{p^n} is divisible by (g-1) and run it for a curve of rank >= 2
- Compute Frobenius eigenvalues for a bunch of curves to illustrate Katz-Sarnak
- Sage Tutorials
Elliptic curves over function fields
This project will include the following topics:
- Compute regulators of elliptic curves over function fields:
This worksheet does it using Tate's suggestion (i.e., use the definition): http://sagenb.org/home/pub/1198/
- Is there an algorithm to enumerate all elliptic curves over a function field of a given conductor?
Victor Miller reports that there are some papers that give efficient algorithms for solving "S-unit equations" over function fields, which seems relevant. It's well known that finding all Elliptic Curves with good reduction outside a finite set, S of places is equivalent to solving S-unit equations. There are a series of four papers by Gaal and Pohst that give efficient algorithms for this when dealing with function fields over finite fields, especially \mathbb{F}_q(t).
Implement Tate's algorithm for elliptic curves over the function field \mathbb{F}_p(t).
- Verify BSD for elliptic curves over function fields of analytic rank 3 or higher
Robert Bradshaw implemented faster arithmetic for F_p(t): #7585
Let K=\mathbb{F}_5(t) and E/K: y^2 = x(x+t)(x+t^2) be the Legendre curve and consider its twist by f=1+t, E_f. E_f has Weierstrass equation
The L-funtion of E_f is:
The analytic rank is 2 with L^{(2)}(E_f/K,1/5) = 100. The points P=[3*t^3 + t^2 + 4*t: t^3 + 3*t^2:1], Q=[3*t^3 + 2*t^2 + 3*t: t^3:1] are independent on E_F(K), so its algebraic rank is 2. Thus the refined BSD conjecture is true. We do not if P,Q generate E_f(K) mod torsion, so the determinant of their height matrix (=5) gives an upper bound on the regulator of E_f/K. The Tamagawa number was computed by hand to be 2^7/5. The torsion subgroup has size 4 (all 2-torsion). Thus
Note that there is probably a factor of 2 going unaccounted for in the computation of the regulator. From this, we see that |Sha|\geq \frac{25}{4} (again, this may be off by a power of 2).
People: Sal Baig, William Stein, David Roe, Ken Ribet, Kiran Kedlaya, Robert Bradshaw, Victor Miller (S-unit equations), Thomas Barnet-Lamb
Implement computation of the 3-Selmer rank of an elliptic curve over QQ
Some projects:
- Implement the algorithm from Schaefer-Stoll which does the general case
Implement the algorithm in the reducible case due to Jeechul Woo: Jeechul Woo's GP script
- Compute the 3-Selmer ranks of all curves of conductor up to 1000
People: Robert Miller, William Stein, Victor Miller, Jeechul Woo (Noam's student; around only Thu, Fri)
Compute statistics about distribution of Heegner divisors and Kolyvagin divisors modulo primes p
Compute the reduction of x_1 using ternary quadratic forms, then use distribution relations and hit by Hecke operators to get reduction of all x_n. There is a theorem of Jetchev-Kane about the asymptotic distribution of x_n; compare our new data with that.
Stein: I posted a bundle based against Sage-4.2.1 here (called heegner-4.2.1.hg), which has highly relevant code: http://trac.sagemath.org/sage_trac/ticket/6616
People: William Stein, Dimitar Jetchev, Drew Sutherland, Mirela Ciperiani, Ken Ribet, Victor Miller
Create a table of images of Galois representations, for elliptic curves
The goals of this project are:
Compute and record in a nice table the exact image of Galois in GL_2(F_p) for all p<60 and all curves in Cremona's tables, using Drew's new code/algorithm.
Extend the above to all p by using the explicit bound coded in Sage.
Extend the above to all p^k.
Compute the exact image for all curves of conductor up to 10^8 from the Stein-Watkins database. Add this data with some nice key to that database (i.e., change all the files to include a new field).
Think about images in GL(Z/mZ).
- Think about statistics resulting from the above computation.
People: Drew Sutherland, Ken Ribet, William Stein, Kiran Kedlaya, David Roe
Fast computation of Heegner points
Implement the algorithm of Delauny/Watkins's algorithm for fast computation of Heegner points y_K \in E(K).
People: William Stein, Robert Bradshaw, Jen Balakrishnan
Implement code to compute the asymptotic distribution of Kolyvagin classes
This will be based on Jared Weinstein's talk. See http://wstein.org/misc/sagedays18_papers/weinstein-kolyvagin_classes_for_higher_rank_elliptic_curves.pdf
People: Jared Weinstein, Mirela Ciperiani, William Stein
Verify Kolyvagin's conjecture for a specific rank 3 curve
This is done for examples of rank 2 curves. Nobody has ever checked that Kolyvagin's conjecture holds for a rank 3 curve.
Figure out exactly what needs to be computed and what might be an optimal curve and quadratic imaginary field to work with: Some details for 5077a
- Verify that one Kolyvagin class for that curve is nonzero.
Possibly verify the conjecture for the first (known) rank 4 curve, which has conductor 234446. This would be computationally hard, but not impossible!
People: William Stein, Dimitar Jetchev, Victor Miller (sparse linear algebra), Jen Balakrishnan, Robert Bradshaw
Implement an algorithm in Sage to compute Stark-Heegner points
There is a new algorithm due to Darmon and Pollack for computing Stark Heegner point using overconvergent modular symbols. So this project would involve:
- Implementing computation of overconvergent modular symbols.
- Using an implementation of overconvergent modular symbols to implement the Stark-Heegner point algorithm.
People: Matthew Greenberg, Cameron Frank, Kiran Kedlaya, Robert Pollack, Avner Ash, David Roe, Jay Pottharst, Thomas Barnet-Lamb
Compute the higher Heegner point y_5 on the curve 389a provably correctly
Implement an algorithm to compute the Gross-Zagier Rankin-Selberg convolution L-functions L(f,\chi,s), and use Zhang's formula to deduce heights of Heegner points.
- Apply in the particular curves 389a for n=5.
- Come up with an algorithm that is definitely right for provably computing Heegner points given the height.
- Implement algorithm and run for 389a and n=5.
- Make a table of heights of higher Heegner points. (Search to find any of height 0!)
- Make a table of heights of derived Kolyvagin points.
People: Robert Bradshaw, William Stein, Jen Balakrishnan
Compute a Heegner point on the Jacobian of a genus 2 curve
People: Noam Elkies, virtually via his comments in this thread.
Visibility of Kolyvagin Classes
Jared describes the problem thus: On the one hand, one of Kolyvagin's suite of conjectures is that the classes he constructs out of Heegner points ought to generate the entirety of Sha(E/K). On the other hand each element of Sha(E/K) is "visible" in some other abelian variety, in the sense of Mazur (see for instance http://www.williamstein.org/home/was/www/home/wstein/www/papers/visibility_of_sha/). I wonder if the Kolyvagin classes d(n) contributing to Sha become visible in abelian varieties in some *systematic* way; ie, in a modular Jacobian of possibly higher level. There are already some examples out there of computation with visibility and computation with Kolyvagin classes, so maybe we can gather some data.
I wrote a short paper with Dimitar Jetchev in 2005 (which I forgot about until just now, and never published), which I think was motivated by this question: :jetchev-stein-congruences_and_unramified_cohomology.pdf
People: Jared Weinstein, Mirela Ciperiani, William Stein, Dimitar Jetchev, Ken Ribet, Barry Mazur
Find an algorithm to decide if y_{p^n} is divisible by (g-1) and run it for a curve of rank >= 2
Consider the Heegner y_{p^n} over the anticyclotomic tower for a rank >= 2 curve, Sha trivial, etc.
- Barry Mazur remarks that this may be connected to his notion (with Rubin) of "Shadow lines" in Mordell-Weil.
People: Mirela Ciperiani, William Stein, Barry Mazur, Jay Pottharst
Compute Frobenius eigenvalues for a bunch of curves to illustrate Katz-Sarnak
- curves that vary in various ways, e.g., coverings of one curve...
- might use David Harvey's super fast Kedlaya zeta function implementation in Sage (his C++ code), but unfortunately the curves that typically arise aren't hyperelliptic, and Harvey's code only applies to hyperelliptic curves.
- could lead to questions of an "anabelian sort"... e.g., take elliptic curve over finite field, take n-torsion points, extract roots, get coverings, etc. Compute Frobenius eigenvalues of all these covers.
To do a useful computation, one needs to compute Frobenius, say for the curve y^{\ell} = f_N(x), where f_N(x) is the N-th division polynomial on an elliptic curve. This following paper is about how to count points on superelliptic curves: Gaudry-Gurel More from Barry: " I 've thought about my off-the-cuff suggestion for computing
Frobenius eigenvalues, and my worry is that one won't really see anything terribly interesting if one works only with N=2; but maybe when one works with N=3. For example, take an elliptic curve E over F_q and let f_3 be (3-division") function on E: meaning that it vanishes to order one at all nontrivial 3-torsion points and has the appropriate order pole at the origin (and no other poles or zeroes). Then (for small n) extract n-th powers of that f_3 to get curves C_n mapping to E (of unfortunately high genus). But it could be that the collection of their Frobenius eigenvalues tells us more than just the isogeny class of E? From what I learn by talking with William, this might be a very difficult computational problem though..."
There is a better paper maybe about this super-elliptic algorithm, with a Magma implementation here: http://www.math.tu-berlin.de/~minzlaff/. And here is a Sage worksheet that *wraps* that Magma code, so if you have Magma you can compute quickly the zeta function of y^\ell = f(x) for f of odd degree. And indeed Minzlaff's code seems to work fine.
- Kiran has a short Sage program to compute the zeta function of y^4 = (cubic), to be published soon.
People: Barry Mazur, Kiran Kedlaya, Thomas Barnet-Lamb, David Harvey, Mirela Ciperiani, Sal Baig (lots of possibly relevant data over function fields)
Sage Tutorials
We would like to have a sequence of informal Sage tutorials on the following topics:
- Introduction to Python/Sage (Kiran Kedlaya, 10am on Wednesday)
- Linear algebra modulo p (Robert Bradshaw, 11am on Wednesday)
- Tate's algorithm over number fields (David Roe, 10am on Thursday)
- 2-descent in Sage (Robert Miller, 11am on Thursday)
- Computing images of Galois representations (Drew Sutherland and William Stein, 11am on Friday)