Doubly periodic functions

This is somewhat a continuation of How (not) to get into isogenies. The reason for the word "not" in the title is that when working with isogeny-based cryptographic protocols, you typically don’t encounter the following theorems. However, they are of great importance and form the basis for complex multiplication. They are also necessary for understanding the background of the theta model in higher-dimensional elliptic curves.

Excursion into topology

As with complex analysis last time, I will now try to refresh my knowledge of topology and intuitively understand some basic results.

Definition: A topological space is compact if every open covering has a finite sub-covering.
Theorem:

A compact subset of Cn\mathbb{C}^n with its usual metric is closed and bounded.

The topology of C\mathbb{C} is the same as that of R2,\mathbb{R}^2, so the above theorem holds for C\mathbb{C} as well.

Definition:

A function f:XYf: X \rightarrow Y is continuous if the inverse image of every open subset of YY is open in X.X.

It might not seem so at first glance, but the above definition is a generalization of the notion that continuous functions are those without abrupt changes in value.

Theorem:

If function f:XYf: X \rightarrow Y is continuous it maps compact sets to compact sets.

Doubly periodic functions and some of their properties

Let w1,w2w_1, w_2 be complex numbers that are linearly independent over R. Then

L=Zw1+Zw2={n1w1+n2w2n1,n2Z}L = \mathbb{Z}w_1 + \mathbb{Z}w_2 = \{n_1w_1 + n_2w_2 | n_1, n_2 \in \mathbb{Z}\}

is called a lattice.

The set

F={aw1+aw20ai<1,i=1,2}F = \{aw_1 + aw_2 | 0 \leq a_i < 1, i = 1, 2\}

is called a fundamental parallelogram for L.L.

Note that C/Z\mathbb{C}/Z is a torus—we can visualize it by taking a fundamental parallelogram and "gluing" the bottom to the top and the left to the right.

Definition:

A doubly periodic function is meromorphic function f:CCf: C \rightarrow C \cup \infty such that

f(z+w)=f(z)f(z + w) = f(z)

for all zCz \in \mathbb{C} and all wL.w \in L.

If ff is a meromorphic functions and wCw \in \mathbb{C}, we can write ff using a Laurent series:

f(z)=ar(zw)r+ar+1(zw)r+1+...f(z) = a_r (z - w)^r + a_{r+1} (z-w)^{r+1} + ...

The order and the residue of ff at ww are defined as follows:

r=ordw(f)r = ord_w(f)

a1=Resw(f)a_{-1} = Res_w(f)

In what follows, ff is a doubly periodic function for the lattice LL, and FF is a fundamental parallelogram for L.L.

I haven’t closely examined what a meromorphic function is, but I decided I can work with the following notion:

A meromorphic function is a function f(z)f(z) of the form f(z)=g(z)/h(z)f(z) = g(z) / h(z), where g(z)g(z) and h(z)h(z) are entire functions (holomorphic on the whole complex plane) and h(z)h(z) is not identically zero.

I will be lazy once more and simply believe that meromorphic functions can only have a finite number of zeros and poles in a compact set.

Note that the closure of FF is a compact set. Thus, ff has only a finite number of zeros and poles in F.F.

If f has no poles, then f is constant

If ff has no poles, then it’s bounded in the closure of F.F. Because ff is doubly periodic, it maps to the same values inside and outside of F,F, so it’s bounded everywhere.

It follows from Liouville’s theorem that ff is constant.

Sum of residues

It holds:

wFResw(f)=0\sum_{w\in F} Res_w(f) = 0

First, it should be noted that the sum has only finitely many nonzero terms. This is because a nonzero term occurs only when there is a zero or a pole.

Let’s recall Cauchy’s theorem:

z0FResz0(f)=12πiFf(z)dz\sum_{z_0 \in F} Res_{z_0}(f) = \frac{1}{2\pi i} \int\limits_{\partial F}f(z)dz

where F\partial F is the boundary of F.F.

Ff(z)dz=0w2f(z)dz+w2w1+w2f(z)dz+w1+w2w1f(z)dz+w10f(z)dz\int\limits_{\partial F}f(z)dz = \int_0^{w_2} f(z)dz + \int_{w_2}^{w_1+w_2} f(z)dz + \int_{w_1+w_2}^{w_1} f(z)dz + \int_{w_1}^0 f(z)dz

Note that f(z+w)=f(z).f(z + w) = f(z). Thus:

0w2f(z)dz=0w2f(z+w1)dz=w1w1+w2f(z)dz=w1+w2w1f(z)dz\int_0^{w_2} f(z)dz = \int_0^{w_2} f(z+w_1)dz = \int_{w_1}^{w_1+w_2} f(z)dz = -\int_{w_1+w_2}^{w_1} f(z)dz

We can see that the four integrals cancel each other out.

Sum of orders

wFordw(f)=0\sum_{w\in F} ord_w(f) = 0

I can’t help but be surprised again. It seems strange that you can sum orders and get 0. Summing coefficients (residues) and obtaining 0 feels slightly more intuitive to me, though I'm not sure why. The order probably seems too related to how the poles and zeros behave individually; I have difficulty grasping why you could sum orders and arrive at such a significant formula. I suppose my intuition is just lacking, as I haven’t done much of this kind of math.

But it feels like it might have something to do with derivatives (I say this after having seen the proof).

We can write f(z)f(z) as

f(z)=(zw)rg(z)f(z) = (z-w)^r g(z)

where ww is neither a zero or pole of g(z).g(z).

Now:

f(z)=r(zw)r1g(z)+(zw)rg(z)f'(z) = r(z-w)^{r-1} g(z) + (z-w)^r g'(z)

And:

f(z)/f(z)=rzw+g(z)/g(z)f'(z)/f(z) = \frac{r}{z-w} + g'(z)/g(z)

We can see:

Resw(f(z)/f(z))=rRes_w(f'(z)/f(z)) = r

This is true because g(z)/g(z)g'(z)/g(z) doesn’t have a pole at w.w.

Now, ff\frac{f'}{f} is also doubly periodic, and we know from the previous section that:

wFResw(ff)=0\sum_{w\in F} Res_w(\frac{f'}{f}) = 0

It holds Reswff=ordw(f)Res_w\frac{f'}{f} = ord_w(f), so we have proved the result.

Weighted sum of orders

If ff is not identically zero:

wFwordw(f)L\sum_{w\in F} w \cdot ord_w(f) \in L

I am not going to fully dive into this one, but I can see:

wFwordw(f)=wFwResw(ff)\sum_{w\in F} w \cdot ord_w(f) = \sum_{w\in F} w \cdot Res_w(\frac{f'}{f})

=wFResw(zff)= \sum_{w\in F} \cdot Res_w(z\frac{f'}{f})

The last step follows from observing:

Resw(zf(z)/f(z))=Resw(z(rzw+g(z)/g(z)))Res_w(zf'(z)/f(z)) = Res_w(z(\frac{r}{z-w} + g'(z)/g(z)))

=Resw((zw+w)(rzw+g(z)/g(z)))= Res_w((z-w + w)(\frac{r}{z-w} + g'(z)/g(z)))

=Resw(r+(zw)g(z)/g(z)+w(rzw+g(z)/g(z)))= Res_w(r + (z-w) g'(z)/g(z) + w(\frac{r}{z-w} + g'(z)/g(z)))

=Resw(w(rzw+g(z)/g(z)))= Res_w(w(\frac{r}{z-w} + g'(z)/g(z)))

=wResw(rzw+g(z)/g(z))= w\cdot Res_w(\frac{r}{z-w} + g'(z)/g(z))

=wordw(f)= w\cdot ord_w(f)

So:

wFwordw(f)=wFResw(zff)\sum_{w\in F} w \cdot ord_w(f) = \sum_{w\in F} \cdot Res_w(z\frac{f'}{f})

=12πiFzffdz= \frac{1}{2\pi i} \int\limits_{\partial F}z\frac{f'}{f}dz

It appears that the proof can be finished quite easily if you know what a winding number is. Intuitively, it’s easy to understand the winding number—it’s the number of times that the curve γ\gamma passes (counterclockwise) around a point. However, I wanted to understand why it’s given by the following formula:

n(γ,z0)=12πiγdzzz0n(\gamma, z_0) = \frac{1}{2\pi i} \int\limits_{\gamma}\frac{dz}{z-z_0}

But I have decided to leave this for another time. I will try to continue toward isogenies now.

Surjectivity

If ff is not constant, then f:CCf: \mathbb{C} \rightarrow \mathbb{C} \cup \infty is surjective.

Let’s observe if there is a solution for f(z)=z0.f(z) = z_0.

Remember from above that a non-constant function has poles. Now, h(z)=f(z)z0h(z) = f(z) - z_0 is doubly periodic and has the same poles as f.f.

From

wFordw(h)=0\sum_{w\in F} ord_w(h) = 0

we know if there are poles, there are zeros too.