E. Kowalski's blog

Comments on mathematics, mostly.

Archive for the ‘Computers’ Category

Update on a bijective challenge

leave a comment

A long time ago, I presented in a post a fun property of the family of curves given by the equations
x+1/x+y+1/y=t,
where t is the parameter: if we consider the number (say N_p(t)) of solutions over a finite field with p\geq 3 elements, then we have
N_p(t)=N_p(16/t),
provided t is not in the set \{0, 4, -4\}. This was a fact that Fouvry, Michel and myself found out when working on our paper on algebraic twists of modular forms.

At the time, I knew one proof, based on computations with Magma: the curves above “are” elliptic curves, and Magma found an isogeny between the curves with parameters t and 16/t, which implies that they have the same number of points by elementary properties of elliptic curves over finite fields.

By now, however, I am aware of three other proofs:

  • The most elementary, which motivates this post, was just recently put on arXiv by T. Budzinski and G. Lucchini Arteche; it is based on the methods of Chevalley’s proof of Warning’s theorem: it computes N_p(t) modulo p, proving the desired identity modulo p, and then concludes using upper bounds for the number of solutions, and for its parity, to show that this is sufficient to have the equality as integers.  Interestingly, this proof was found with the help of high-school students participating in the Tournoi Français des Jeunes Mathématiciennes et Mathématiciens. This is a French mathematical contest for high-school students, created by D. Zmiaikou in 2009, which is devised to look much more like a “real” research experience than the Olympiads: groups of students work together with mentors on quite “open-ended” questions, where sometimes the answer is not clear (see for instance the 2016 problems here).
  • A proof using modular functions was found by D. Zywina, who sent it to me shortly after the first post (I have to look for it in my archives…)
  • Maybe the most elegant argument comes by applying a more general result of Deligne and Flicker on local systems of rank 2 on the projective line minus four points with unipotent tame ramification at the four missing points (the first cohomology of the curves provide such a local system, the missing points being 0, 4, -4, \infty). Deligne and Flicker prove (Section 7 of their paper, esp. Prop. 7.6), using a very cute game with products of matrices and invariant theory, that if \mathcal{V} is such a local system and \sigma is any automorphism of the projective line that permutes the four points, then \sigma^*\mathcal{V} is isomorphic to \mathcal{V}.

Not too bad a track record for such a simple-looking question… Whether there is a bijective proof remains open, however!

Written by Kowalski

September 24th, 2016 at 10:47 am

L-functions database!

leave a comment

People studying automorphic forms, automorphic representations, number fields, diophantine equations, function fields, algebraic curves, equidistribution and many other arithmetic objects (j’en passe, et des meilleurs), often end up with some “L-function” to deal with — indeed, probably equally often, with a whole family of them, sometimes not so well-behaved… These objects are fascinating, mystifying, exhilarating, random and possibly spooky. Where they really come from is still a mystery, even with buzzwords aplenty ringing around our ears. But one remarkable thing was already known to Euler and to Riemann: one can compute with L-functions. One impressive research project has been building, for quite a few years, a very sophisticated website presenting enormous amounts of data about L-functions of many kinds. The L-functions Database is now out of its beta status: go see it, and have a look at the list of editors to see who should be thanked for this amazing work!

Written by Kowalski

May 9th, 2016 at 9:46 pm

New versions, new bugs

leave a comment

Ph. Michel found the first bug in the new version of Kloostermania: when entering a modulus for a sum from the menu, there was no check that it is prime (or adjustment to make it so). This is now corrected in version 1.01, to be found in the usual place.

As another bonus feature, a double tap on the screen will cycle between the types of sums presented: Kloosterman sum to Birch sum to both to Kloosterman…

Written by Kowalski

May 7th, 2015 at 9:22 am

Kloostermania turns 1.0

one comment

Since the Pocket Kloostermania has remained unchanged for almost five years, an update is probably welcome! This is now available, in time to honor N. Katz on the occasion of this 71st birthday (which I will help celebrate in China next week).

Important changes in the new version:

(1) It was recompiled — hence a new modern look instead of a style reminiscent of the dark ages –, and the resulting binary should work on any Android system with version 4.0.3 or later;

(2) To keep up with recent progress, the program now displays Kloosterman sums and/or Birch sums, instead of Kloosterman and/or Salié sums. On the other hand, the moduli used are now restricted to primes, to simplify things a bit, and only one parameter is used: the sums displayed are
\sum_{x\in\mathbf{F}_p^{\times}}e\Bigl(\frac{ax+\bar{x}}{p}\Bigr)
and
\sum_{x\in\mathbf{F}_p}e\Bigl(\frac{ax+x^3}{p}\Bigr)

(3) In addition to being able to change the modulus by swiping horizontally, a vertical swipe on the screen scrolls among the values of the parameter a. As was the case in the previous version (and even more because phones are faster now), the scrolling is usually too fast for a single step, so tapping once close to one of the edges of the window displaying the sum will perform just one step of the corresponding move (e.g., tapping close to the right of the screen goes to the next prime modulus). The value of the sum and its parameters and then displayed quickly.

(4) The plots are presented in orthonormal configurations, to give a more faithful representation of the paths in the plane;

Kloostermania

Kloostermania

(5) In the “About” dialog, a single click will display (if a PDF viewer is installed) the paper of W. Sawin and myself that explains the limiting distribution of Kloosterman (and Birch) paths…

(6) It is possible to save (or rather to “share” in the usual Android way) a picture of the sums currently displayed, as a PNG file.

(6) And the launcher icon is better.

Icon

Icon

The installation file is available right now on the updated Kloostermania page!

Worms, brains and “graphes expanseurs”

leave a comment

Two days ago, I was in Paris and visited the École Polytechnique there (for the first time, as it turns out), in order to give the traditional introductory lecture presenting various aspects of mathematics that is given each year to the incoming class of students. (Many people would wonder how a new class can be incoming in April; this is because the students of École Polytechnique, which is a military school, begin their studies with a military/civil service from September to April).

I had selected expander graphs as the topic of my talk, as being simultaneously a very modern subject, that can be presented with very basic definitions (most students, coming from the French “Classes préparatoires” have had two very solid years of mathematical studies, but of a rather traditional kind), and that has links to many different areas, including more applied ones.

The slides of my presentation can be found here, in French, and I prepared an English translation there (up to the handwritten bits of information on certain pictures, which will remain in French…) Actually, the slides by themselves are not particularly interesting, since there are many remarks and details that I only described during the talk.

The talk begins with a discussion of graphs, continues with basic notions of expansion, and then defines expander graphs, with a bit of the history (especially the paper of Barzdin and Kolmogorov that I like very much), and concludes with two applications (among many…): Apollonian circle packings, and the Gromov-Guth knot-distorsion theorem. If this looks like a lot of ground to cover, this is because the talk was 1h30 long…

While preparing the first part of the talk, I researched some examples of graphs with more care than I had done before. Some of the things I found are extremely interesting, and since they might not be so well-known among my readers, I will present them quickly.

(1) The worm is Caenorhabditis elegans, more chummily known as C. elegans, one of the stars of neurosciences, as being the only animal whose entire nervous system has been mapped at the level of all individual neurons and connections between them. This was done in 1986 by White, Southgate, Thomson and Brenner, with minor corrections since then, and an important update and representation in 2011 by Varshney, Chen, Paniagua and Chklovskii. The resulting graph has 302 neurons (this number is apparently constant over all individuals), and about 8000 edges. (Interestingly, it is naturally a mix of directed and un-directed edges, depending on the type of biological connection). The paper of Varshney, Chen, Paniagua and Chklovskii is quite fascinating, as it investigates many mathematical invariants of the graph, and for instance compares the number of small subgraphs of various types with the expected number for random graphs with comparable parameters (certain subgraphs arise with higher frequency…)

Much more about C. elegans is found on dedicated websites, including Openworm, which seems to have as a goal to recreate the animal virtually…

(2) The brain is the humain brain; it can’t be mapped at the level of C. elegans (there are about 10^{11} neurons and 10^{15} synapses, from what I’ve seen), but I read a very interesting survey by Valiant about attempts to understand the computational model underpinning the reasoning capacities of the brain. He presents four basic tasks that must be among those that the brain performs, and explains how he succeeded in earlier work (from 1994 to 2005) in finding realistic algorithms and models for these problems. He comments:

In [5,14] it is shown that algorithms for the four random
access tasks described above can be performed on the
neuroidal model with realistic values of the numerical
parameters. The algorithms used are all of the vicinal
style. Their basic steps are all local in that they only
change synaptic strengths between pairs of neurons that
are directly connected. Yet they need to achieve the more
global objectives of random access. In order that they be
able to do this certain graph theoretic connectivity prop-
erties are required of the network. The property of
expansion [15], that any set of a certain number of
neurons have between them substantially more neighbors
than their own number, is an archetypal such property.
(This property, widely studied in computer science, was
apparently first discussed in a neuroscience setting [16].)
The vicinal algorithms for the four tasks considered here
need some such connectivity properties. In each case
random graphs with appropriate realistic parameters have
it, but pure randomness is not necessarily essential.

Here reference [16] is to the Barzdin-Kolmogorov paper.

(3) Lastly, I was wondering what is the size (number of vertices) of the largest graphs for which the first non-zero Laplace eigenvalue has been computed, approximately. The best I found is a paper of Kang, Breed, Papalexakis, and Faloutsos, which describes an algorithm that they have applied successfully to real-world graphs with about a billion vertices. This is rather impressive…

Written by Kowalski

April 11th, 2015 at 1:03 pm