Weeks 6,7 (June 28 – July 10)

We’re two weeks into the second half of GSoC 2016. The current goal is implementing resolution of singularities for algebraic curves via blow-ups. My mentors and I have discussed ways to implement blow-ups for curves in Sage, and settled on describing the blow-ups with collections of affine charts.

Progress in weeks 6, 7 has been to finish tickets #20811 and #20848, and to work on ticket #20895 which will implement functionality for computing ordinary models of plane curves. I also worked on #20930 which implements multiplicity for arbitrary affine/projective subschemes, and multiplicity/intersection multiplicity for subschemes of products of projective spaces. Helper functions in the corresponding point classes to access this functionality are being implemented as well.

Implementing the ordinary model functionality has proven to be somewhat tricky. Theoretically, there is a straightforward method for transforming a plane curve into one with only ordinary singularities that works in a finite number of steps. Before describing it, here’s some preliminaries:

We use the birational automorphism of \mathbb{P}^{2} known as the standard Cremona transformation, and will denote it as \phi. This is the map defined (x : y : z)\mapsto (yz : xz : xy) for (x : y : z)\in\mathbb{P}^{2}, and is undefined at the points (1 : 0 : 0), (0 : 1 : 0), (0 : 0 : 1) which we will call the fundamental points of \mathbb{P}^{2} for convenience. We also call the lines V(x), V(y), V(z) the exceptional lines of \mathbb{P}^{2}. Next consider the open subset U := \mathbb{P}^{2}\setminus V(xyz)\subseteq\mathbb{P}^{2}, the complement of the union of the exceptional lines. Note that \phi is defined everywhere in U, and that the exceptional lines are exactly the fibers of \phi above the fundamental points, so the image of \phi is contained in U. We then also see that, when restricted to U, \phi is its own inverse: for (x : y : z)\in U, \phi(\phi((x : y : z))) = \phi((yz : xz : xy)) = (xzxy : yzxy : yzxz) = (x : y : z), since x,y,z \neq 0. So \phi is an isomorphism of U to itself, and thus is indeed a birational automorphism of \mathbb{P}^{2}.

Next, given a plane curve C := V(F)\subseteq \mathbb{P}^{2} for some homogeneous polynomial F\in k[x,y,z], k algebraically closed (our base field), a defining polynomial for the image of C by \phi is F(yz,xz,xy), which we can write as F(yz,xz,xy) = G(x,y,z)H(x,y,z) for a homogeneous H\in k[x,y,z] and G(x,y,z) = x^{r}y^{s}z^{t} for some r, s, t\in\mathbb{Z}_{\geq 0} such that no positive power of x,y or z divides H. We call C^{\prime} := V(H) the quadratic transformation of C.

We also define the apparent genus of a plane curve C := V(F) as g^{*}(C) = \frac{(d - 1)(d - 2)}{2} - \sum_{p\in\mathrm{Sing}(C)}\frac{m_{C}(p)(m_{C}(p) - 1)}{2}, where \mathrm{Sing} is the set of singular points of C, and m_{C}(p) is the multiplicity of p. Finally, if F has degree n, we say C is in excellent position if (0:0:1) is a point of multiplicity r > 1 of C, no exceptional line is tangent to C at a fundamental point, the exceptional lines V(x), V(y) intersect C each in n - r distinct points aside from (0 : 0 : 1), transversally, and the exceptional line V(z) intersects C transversally in n distinct, non-fundamental points.

The main result is that if (0 : 0 : 1) is a non-ordinary singularity of C, and if C is in excellent position, then its quadratic transformation C^{\prime} either has one less non-ordinary singularity, or has a smaller apparent genus. The apparent genus can be shown to always be a nonnegative integer, thus we have the following algorithm to resolve the non-ordinary singularities of a plane curve:

Input: a plane curve C

  • while C has non-ordinary singularities,
    • pick a non-ordinary singularity p of C
    • find a change of coordinates sending p to (0 : 0 : 1), and putting C into excellent position
    • compute the quadratic transformation C^{\prime}
    • replace C with C^{\prime}
  • return the plane curve C

The output plane curve will be birational to the original curve, and will have only ordinary singularities. By the above note, this algorithm will always terminate. It isn’t difficult to prove that such change of coordinates maps exist, but this is where the difficulties of actual implementation begin.

Currently, the approach in #20895 to constructing a change of coordinates map to put a curve C with non-ordinary singular point p into excellent position is to incrementally generate two points q, d such that p,q,d are not collinear until the lines connecting the points satisfy the excellent position conditions. That is, ideally we would map p,q,d to (0:0:1), (0:1:0), (1:0:0), and the curve would then be in excellent position.

With some care taken to incrementally choose the points q,d the selection process will eventually terminate after finitely many choices. However, checking whether the transformation constructed from a choice of points will put the curve into excellent position can be an extremely expensive operation. Most of the cost is due to finding and analyzing the intersections of the lines connecting the points p,q,d with C. We currently work over the implementation of the algebraic closure of the rationals \overline{\mathbb{Q}} in Sage so that we can find all non-ordinary singular points and all intersections of lines with C without needing to extend the field we start with, but computations over this field can be quite slow.

To get around this, our current approach is to bypass the explicit checking process; we just apply the transformation and repeat. To verify that the transformations will work, we instead check that the procedure terminates after t + g^{*}(C) steps, where C is our starting curve, and t is the number of non-ordinary singularities of C. If the change of coordinate maps all put the curve into excellent position, then the algorithm should terminate in at most t + g^{*}(C) steps, otherwise something went wrong, and we choose points differently in the next attempt.

Most of the time, the first choices of points works, but the resulting transformations can be quite messy. Here is a simple example with only one non-ordinary singularity:

sage: set_verbose(-1)
sage: P.<x,y,z> = ProjectiveSpace(QQbar, 2)
sage: C = Curve([y^2*z - x^3], P)
sage: C.singular_points()
[(0 : 0 : 1)]
sage: C.tangents(P([0,0,1])) # higher multiplicity tangent
sage: D = C.ordinary_model()
sage: D
(Scheme morphism:
   From: Projective Plane Curve over Algebraic Field defined by -x^3 +
   To:   Projective Space of dimension 2 over Algebraic Field
   Defn: Defined on coordinates by sending (x : y : z) to
         ((-1/16)*x^2 + (-1/4)*x*y + (-1/4)*y^2 + 1/4*x*z + 1/2*y*z :
1/16*x^2 + (-1/4)*y^2 + (-1/4)*x*z + 1/2*y*z : (-1/16)*x^2 + 1/4*y^2),
 Projective Plane Curve over Algebraic Field defined by x^3*y +
2*x^2*y^2 + x*y^3 + (-7)*x^3*z + 26*x^2*y*z + (-23)*x*y^2*z + 8*y^3*z)
[(0 : 0 : 1)]
D[1].tangents(P([0,0,1])) # singularity is now ordinary
[x + (-2.630294284738199?)*y,
 x + (-0.5419957147737579? + 0.3751512615427921?*I)*y,
 x + (-0.5419957147737579? - 0.3751512615427921?*I)*y]

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s