Four types of maps between elliptic curves may be constructed: isogenies, isomorphisms, translations, and rational maps. Isogenies and isomorphisms are by far the most important and have the most functions associated to them. Isogenies are always surjective as scheme maps, even though the Magma parlance of a map from E(Q) -> F(Q) may seem to indicate that (for instance) the multiplication-by-two map is not surjective. There is an internal limit that the degrees of the polynomials defining an isogeny cannot be more than 107.
> FF := FiniteField(167); > a := FF!2; b := FF!3; r := a^2 - 4*b; > E1 := EllipticCurve([0, a, 0, b, 0]); > E2 := EllipticCurve([0, -2*a, 0, r, 0]); > _<x> := PolynomialRing(BaseRing(E1)); > Ff, f := IsogenyFromKernel(E1, x); > Fg, g := IsogenyFromKernel(E2, x); > b, m1 := IsIsomorphic(Ff, E2); assert b; > b, m2 := IsIsomorphic(Fg, E1); assert b; > // Verify that f and g are dual isogenies of degree 2 > &and[ m2(g(m1(f(P)))) eq 2*P : P in RationalPoints(E1) ]; true > &and[ m1(f(m2(g(Q)))) eq 2*Q : Q in RationalPoints(E2) ]; true
Given elliptic curves E and F defined over the same field K, and four elements r, s, t, u of K with u ≠0, this function returns the isomorphism E -> F mapping OE |-> OF and mapping (x, y) |-> (u2x + r, u3y + su2x + t). This function returns an error if the values passed do not define such an isomorphism.
Given elliptic curves E and F defined over the same field K, this function computes and returns an isomorphism from E to F where such an isomorphism exists. The map returned will be the same as the second return value of IsIsomorphic.
Given an elliptic curve E defined over a field K, and four elements r, s, t, u of K with u ≠0, this function returns the automorphism E -> E mapping OE |-> OE and mapping (x, y) |-> (u2x + r, u3y + su2x + t). This function returns an error if the values passed do not define such an automorphism.
The sequence [r, s, t, u] of elements defining the isomorphism I.
> K := GF(73); > E1 := EllipticCurve([K | 3, 4, 2, 5, 1]); > E2 := EllipticCurve([K | 8, 2, 29, 45, 28]); > IsIsomorphic(E1, E2); true > m := Isomorphism(E1, E2, [3, 2, 1, 4]); > m; Elliptic curve isomorphism from: CrvEll: E1 to CrvEll: E2 Taking (x : y : 1) to (16*x + 3 : 64*y + 32*x + 1 : 1) > P1 := Random(E1); > P2 := Random(E1); > m(P1 + P2) eq m(P1) + m(P2); trueFrom the isomorphism data we can apply the map by hand if desired:
> r, s, t, u := Explode(IsomorphismData(Inverse(m))); > P3 := E2![ 69, 64 ]; > x, y := Explode(Eltseq(P3)); > E1 ! [ u^2*x + r, u^3*y + s*u^2*x + t ]; (68 : 32 : 1) > m($1) eq P3; true
This function returns true if and only if the isogeny I has the same action as some isomorphism; if so, the isomorphism is also returned.
This function takes a map I of type isomorphism and returns an equivalent map with type isogeny.
> FF := FiniteField(23); > E0 := EllipticCurve([FF | 1, 1]); > E1 := EllipticCurve([FF | 3, 2]); > b, iso := IsIsomorphic(E0, E1); > b; true > iso; Elliptic curve isomorphism from: CrvEll: E0 to CrvEll: E1 Taking (x : y : 1) to (16*x : 18*y : 1) > isog := IsomorphismToIsogeny(iso); > isog; Elliptic curve isogeny from: CrvEll: E0 to CrvEll: E1 taking (x : y : 1) to (16*x : 18*y : 1) > b, new_iso := IsIsomorphism(isog); > b; true > inv := Inverse(new_iso); > P := Random(E0); > inv(isog(P)) eq P; true
Given a rational point P on the elliptic curve E, this function returns the morphism tP : E -> E defined by tP(Q) = P + Q for all rational points Q of E.
Let i be an isogeny and t be a translation map tP: E -> E where tP(Q) = P + Q for some rational point P ∈E. This function returns the rational map φ : E -> F obtained by composing i and t (applying t first). Any rational map E -> F can be represented in this form.
Given a 2-torsion point P of the elliptic curve E, this function returns a 2-isogeny on E with P as its kernel.
> E1 := EllipticCurve([ GF(23) | 1, 1 ]); > E2 := EllipticCurve([ GF(23, 2) | 1, 1 ]);The doubling map on E1 lifted to E2:
> f := map<E1 -> E2 | P :-> 2*P>;Two slightly different ways to define the negation map on E1 lifted to E2:
> f := map<E1 -> E2 | P :-> E2![ P[1], -P[2], P[3] ]>; > f := map<E1 -> E2 | P :-> (P eq E1!0) select E2!0 > else E2![ P[1], -P[2], 1 ]>;
Let G be a subgroup scheme of an elliptic curve E. There is a separable isogeny f: E -> Ef to some other elliptic curve Ef, which has kernel G. This function returns the curve Ef and the map f using Velu's formulae.
Returns a sequence of isogenies whose product is the isogeny returned by the invocation IsogenyFromKernel(G). These isogenies have either degree 2 or odd degree. This function was introduced because composing isogenies can be computationally expensive. The generic Expand function on a composition of maps can then be used if desired.
Given an elliptic curve E and a univariate polynomial psi which defines a subgroup scheme of E, compute an isogeny using Velu's formulae as above.
Returns a sequence of isogenies whose product is the isogeny returned by the invocation IsogenyFromKernel(E, psi). These isogenies have either degree 2 or odd degree. This function was introduced because composing isogenies can be computationally expensive. The generic Expand function on a composition of maps can then be used if desired.
Given an isogeny I and a subgroup G which contains the kernel of I, find the image of G under the action of I. The subgroup G may be replaced by its defining polynomial v.
Given an isogeny φ : E1 -> E2, the function returns another isogeny φ * : E2 -> E1 such that φ * φ is multiplication by the degree of φ. The result is remembered and DualIsogeny(DualIsogeny(phi)) returns phi.
> E := EllipticCurve([GF(97) | 2, 3]); > E1, f := IsogenyFromKernel(E, DivisionPolynomial(E, 5));The image curve E1 is isomorphic, but not equal, to the original curve E. We proceed to find the dual of f.
> deg := Degree(f); > psi := DivisionPolynomial(E, deg); > f1 := PushThroughIsogeny(f, psi); > E2, g := IsogenyFromKernel(E1, f1); > // Velu's formulae give an isomorphic curve, not the curve itself. > IsIsomorphic(E2, E); true > h := Isomorphism(E2, E); > f_dual := g*IsomorphismToIsogeny(h);The latter isogeny is the dual of f, as we verify:
> &and [ f_dual(f(P)) eq deg*P : P in RationalPoints(E) ]; trueA simpler way to verify this is just to check equality (note that both constructions of the dual involved an arbitrary choice of isomorphism):
> f_dual eq DualIsogeny(f); true > f_dual eq DualIsogeny(f) * NegationMap(E); false;
Returns true if and only if the image of the isogeny I is the zero element of its codomain.
Returns true if and only if the isogenies I and J are equal.
Returns the univariate polynomial ψ used in defining the isogeny I. The roots of ψ determine the kernel of I.
Returns the polynomial ψ used in defining the isogeny I as a multivariate polynomial.
Returns the denominator of the fraction giving the image of the x-coordinate of a point under the isogeny I (the numerator is returned by IsogenyMapPhi(I)).
Returns the univariate polynomial φ used in defining the isogeny I.
Returns the polynomial φ used in defining the isogeny I as a multivariate polynomial.
Returns the multivariate polynomial ω used in defining the isogeny I.
Returns the subgroup of the domain consisting of the elements that map to zero under the isogeny I.
Returns the degree of the morphism I.
Let E be an elliptic curve defined over the field K. The set of endomorphisms (End)(E) of E to itself forms a ring under composition and addition of maps.
The endomorphism ring of E always contains a subring isomorphic to Z. A curve E whose endomorphism ring contains additional isogenies is said to have complex multiplication. If E is defined over a finite field then E always has complex multiplication. The endomorphism ring of E is isomorphic to an order in either a quadratic number field or a quaternion algebra.
Returns the multiplication-by-m endomorphism [m] : E -> E of the elliptic curve E, such that [m](P) = m * P.
Returns the identity map E -> E of the elliptic curve E as an isogeny.
Returns the identity map E -> E of the elliptic curve E as an isomorphism. The defining coefficients are [r, s, t, u] = [0, 0, 0, 1].
Returns the isomorphism of the elliptic curve E which maps P to -P, for all P ∈E.
Returns the Frobenius isogeny (x : y : 1) |-> (xpi : ypi : 1) of an elliptic curve E defined over a finite field of characteristic p.
Equivalent to FrobeniusMap(E, 1).
> p := 23;
> FF1 := FiniteField(p);
> FF2 := FiniteField(p, 2);
> E1 := EllipticCurve([FF1 | 1, 3]);
> E2 := BaseExtend(E1, FF2);
> frob := FrobeniusMap(E2, 1);
> #{ E1!P : P in RationalPoints(E2) | P eq frob(P) } eq #E1;
true
The functions in this section deal with automorphisms in the category of elliptic curves. For an elliptic curve E these are the automorphisms of the underlying curve that also preserve the group structure (equivalently, that send the identity OE to itself).
Warning: The behaviour of these functions depends on the type of the input curve. For a general curve in Magma (an object that is a Crv but not a CrvEll) over a finite field, AutomorphismGroup and Automorphisms give automorphisms in the category of curves.
For an elliptic curve E over a general field K, this function determines the group of automorphisms of E that are defined over K. It returns an abstract group A (an abelian group or a polycyclic group), together with a map which sends elements of A to actual automorphisms of E (maps of schemes with domain and codomain E).
For an elliptic curve E over a general field K, this function returns a sequence containing the elements of the AutomorphismGroup of E.