This is still a rather immature version of our toolbox. Nevertheless, you can use it to compute, for a large class of curves over the rationals, the stable reduction at primes of bad reduction.

Let Y be a smooth projective curve over a field K and let vK be a discrete valuation on K. The principal goal is to compute the semistable reduction of Y with respect to vK. This means that we want to know

• a finite Galois extension L/K,
• an extension vL of vK to L,
• the special fiber of an integral semistable model of Y over the valuation ring of vL, and
• the action of the decomposition group of vL on that special fiber.

At the moment we can do this only in certain special cases, which should nevertheless be useful.

If you have at least Sage 10 you can install the latest version of this package with

sage -pip install git+https://github.com/MCLF/mclf

The package can be loaded with

sage: from mclf import *

We create a Picard curve over the rational number field.

sage: R.<x> = QQ[]
sage: Y = SuperellipticCurve(x^4-1, 3)
sage: Y
superelliptic curve y^3 = x^4 - 1 over Rational Field

In general, the class SuperellipticCurve allows you to create a superelliptic curve of the form yⁿ = f(x), for a polynomial f over an arbitrary field K. But you can also define any smooth projective curve Y with given function field.

We define the 2-adic valuation on the rational field. Then we are able to create an object of the class SemistableModel which represents a semistable model of the curve Y with respect to the 2-adic valuation.

sage: v_2 = QQ.valuation(2)
sage: Y2 = SemistableModel(Y, v_2)
sage: Y2.is_semistable() # this may take a while
True

The stable reduction of Y at p=2 has four components, one of genus 0 and three of genus 1.

sage: [Z.genus() for Z in Y2.components()]
[0, 1, 1, 1]
sage: Y2.components_of_positive_genus()
[the smooth projective curve with Function field in y defined by y^3 + x^4 + x^2,
 the smooth projective curve with Function field in y defined by y^3 + x^2 + x,
 the smooth projective curve with Function field in y defined by y^3 + x^2 + x + 1]

We can also extract some arithmetic information on the curve Y from the stable reduction. For instance, we can compute the conductor exponent of Y at p=2:

sage: Y2.conductor_exponent()
6

Now let us compute the semistable reduction of Y at p=3:

sage: v_3 = QQ.valuation(3)
sage: Y3 = SemistableModel(Y, v_3)
sage: Y3.is_semistable()
True
sage: Y3.components_of_positive_genus()
[the smooth projective curve with Function field in y defined by y^3 + y + 2*x^4]

We see that Y has potentially good reduction at p=3. The conductor exponent is:

sage: Y3.conductor_exponent()
6

For more details on the functionality and the restrictions of the toolbox, see the Documentation. For the mathematical background see

• J. Rüth, Models of Curves and Valuations, PhD thesis, Ulm University, 2014
• I.I. Bouw, S. Wewers, Computing L-Functions and semistable reduction of superellipic curves, Glasgow Math. J., 59(1), 2017, 77-108
• I.I. Bouw, S. Wewers,Semistable reduction of curves and computation of bad Euler factors of L-functions, lecture notes for a minicourse at ICERM
• S. Wewers, Semistable reduction of superelliptic curves of degree p, preprint, 2017

See our issues list, and tell us of any bugs or omissions that are not covered there.

We recommend that you install pixi to provide all the dependencies for developing mclf. Once pixi is installed, clone this repository and use commands such as

pixi run sage                    # run SageMath with mclf installed
pixi run doctest-long            # run doctests against the latest SageMath
pixi run -e sagemath-100 doctest-long  # run doctests against SageMath 10.0
pixi run readthedocs             # validate the Read the Docs build

Most development happens on feature branches against the master branch. The test workflow runs against SageMath 10.0 through 10.7 on Ubuntu, and also tests SageMath 10.7 on Intel macOS and Apple Silicon macOS. Publishing to PyPI is triggered from a GitHub release.