Beyond Polyhedral Homotopies
Anton Leykin, Josephine Yu

TL;DR
This paper introduces a novel algorithmic framework combining tropical geometry and homotopy continuation to solve polynomial systems, extending the capabilities of existing polyhedral homotopies.
Contribution
It generalizes polyhedral homotopies by incorporating tropical geometry and linear subspace structures for solving polynomial equations.
Findings
Successfully extends polyhedral homotopies
Demonstrates effectiveness on complex polynomial systems
Provides a new computational approach for algebraic geometry
Abstract
We present a new algorithmic framework which utilizes tropical geometry and homotopy continuation for solving systems of polynomial equations where some of the polynomials are generic elements in linear subspaces of the polynomial ring. This approach generalizes the polyhedral homotopies by Huber and Sturmfels.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsPolynomial and algebraic computation · Commutative Algebra and Its Applications · Advanced Numerical Analysis Techniques
Beyond Polyhedral Homotopies
Anton Leykin
School of Mathematics, Georgia Tech, Atlanta GA, USA
and
Josephine Yu
Abstract.
We present a new algorithmic framework which utilizes tropical geometry and homotopy continuation for solving systems of polynomial equations where some of the polynomials are generic elements in linear subspaces of the polynomial ring. This approach generalizes the polyhedral homotopies by Huber and Sturmfels.
AL is supported by NSF-DMS grant #1151297
JY is supported by NSF-DMS grant #1600569.
1. Introduction
The polyhedral homotopy continuation method of Huber and Sturmfels [HS95], which is implemented in PHCpack [Ver99], HOM4PS [LLT, CLL14], and PHoM [GKK*+*04], is used for computing numerical approximations of all the roots of polynomial equations in variables, where it is assumed that each equation has generic coefficients with respect to its monomial support. We generalize this to
**Setting A: **
Instead of solving for all solutions in , we would like to find solutions lying on a variety defined by the set of polynomials , away from the base locus, where the number of additional given polynomials (other than ) is equal to and each of the additional polynomials is generic with respect to its monomial support.
This can be further generalized to
**Setting B: **
Instead of monomial supports, we can consider arbitrary supports, i.e. we would like to find roots on of polynomials that are generic linear combinations of arbitrary sets of polynomials.
We propose a 3-stage framework to solve this problem in the general setting.
**Initialization: **
Reformulate the problem in Setting B into Setting A as explained in §2. Then pick—in practice, randomly—a one-parameter system on many polynomials with the specified monomial supports and coefficients that are generic rational powers of the parameter multiplied by a generic complex number.
**Stage 1 (mostly symbolic): **
Compute the tropical variety . This usually involves both polyhedral computations and Gröbner basis computations.
**Stage 2 (polyhedral): **
Compute the (transverse) intersection of with tropical hypersurfaces of polynomials in .
**Stage 3 (mostly numerical): **
Find the initial terms—in general, distinct truncations—of Puiseux series solutions of the system corresponding to the tropical points found in Stage 2. Track the homotopy paths for , which are approximated by these truncated Puiseux series in the beginning ( close to [math]) and lead to solutions of a system with generic complex coefficients at the end ().
When in Setting A, the Stage 1 is not needed, and the framework specializes to the polyhedral homotopy approach. The precise statements are in the pseudocode of Algorithm 7, which is followed by remarks on currently available tools for implementation.
2. The problem setup
The more general setup (Setting B) is as follows. Let . Let be vector subspaces of spanned by finite sets respectively. Let . Let be the base locus of the linear spaces . Let .
Our main goal is to compute all the points in for some generic element . Our enumerative goal is to compute the number of such points, which we assume to be finite.
A discussion of subtleties surrounding the base locus and genericity appears in the Appendix. We can deal with rational functions in ’s by clearing the denominators and removing the zero locus of the the denominators from .
We will now reformulate the Setting B into Setting A. Let be the set of non-monomials in . Consider
[TABLE]
which are sets of polynomials in . Solving the system on is equivalent to solving on the variety where is obtained from by replacing each with . The new polynomials are generic with respect to their monomial support. For the rest of the paper we assume Setting A, that is, forms a monomial basis of for each . We drop the primes ′ for simpler notation.
Example 1**.**
Two generic circles
[TABLE]
intersect in two points in although the mixed volume of their Newton polytopes is . We rewrite the system as:
[TABLE]
The original equations are transformed into polynomials that are generic with respect to their monomial supports, but we acquire a new equation whose coefficients may be special. ∎
To achieve the main goal, we will construct a homotopy with exactly paths to track. The number is called the intersection index in [KK12] and equals the mixed volume of Newton-Okounkov bodies associated to on .
3. Algorithmic framework
Let be the field of Puiseux series that are convergent on a punctured neighborhood of [math] in . (See [Ghy] for a proof of this fact and a historical excursion.) There is a valuation from to given by the leading (lowest) degree.
For an ideal in and a weight vector , the -initial ideal is obtained by first taking the usual initial ideal with the min-convention (leading terms are lightest), where the weight of is and the weights of ’s are given by , and then setting . The -initial ideal is an ideal in ; it does not involve . For example, .
The tropical variety of is defined as
[TABLE]
We often write to denote when , and we write to denote . The tropical variety is a polyhedral complex, and we can define multiplicities on its maximal faces. See [MS15, Chapter 3] and [JMM08] for details.
Theorem 2** (Fundamental Theorem of Tropical Algebraic Geometry).**
[MS15*, Theorem 3.2.3]**
The points in are exactly the coordinatewise valuations of the Puiseux series points in the variety of . The multiplicity of each tropical point is equal to the number of Puiseux series point with valuation , counted with multiplicities.*
We will now formulate the the polyhedral homotopy continuation method of Huber and Sturmfels using tropical geometry. Given a polynomial system with generic coefficients with respect to their monomial supports, we perturb the coefficients by throwing in extra factors of the form where are arbitrary rational numbers, to obtain a family of systems . See (2) with . The homotopy continuation approach looks to “connect” solutions of the original system to the Puiseux series solutions of convergent in some neighborhood of .
If the exponents of ’s are sufficiently generic, then the intersection of tropical hypersurfaces is finite and transverse, i.e. locally at each intersection point it is a transverse intersection of affine spaces. In this case we have
[TABLE]
The points of are most commonly computed by enumerating mixed cells of the mixed subdivision, which is the projection of lower convex hull the Newton polyhedron of the product onto the -coordinates.
The solutions of the binomial initial system give us the leading terms of the convergent Puiseux series with valuation satisfying . Take as a numerical approximation of a point satisfying for a small . We can numerically track a segment of a real curve , starting at that point and finishing with an approximation of a solution of to the original system . This is called the polyhedral homotopy because mixed subdivisions of Newton polyhedra play a crucial role.
Getting back to our Setting A, polynomials in the set , e.g. in Example 1 above, can have special coefficients, while the others have generic coefficients with respect to fixed monomial supports. As above, we wish to compute the tropical variety of the system, which should consist of finitely many points.
For each , let be complex polynomials with generic coefficients. Let
[TABLE]
for some generic . We refer to the system as a homotopy. We recover the original system simply by specializing .
We can look at from two different angles:
- (1)
Consider as a family of systems in , parameterized by . Given a path , we get a path . Suppose the path in does not go through the branch locus of the projection where
[TABLE]
is the incidence variety. Then the homotopy induces smooth homotopy paths , , that give a one-to-one correspondence between starting roots and target roots , where and .
Once we find an appropriate path from some , where is known, to , we achieve our main goal by homotopy continuation along . 2. (2)
Consider as a polynomial system over the Puiseux series with many roots over .
The second point of view relates to the first as follows. The Puiseux series roots in (2) converge in some punctured neighborhood of in . Thus, if one can approximate these roots close to , one can find a starting point for (1) along with approximations of the starting roots. This requires approximating the Puiseux series roots, which is explained in Remark 6.
Remark 3**.**
One can show that the (really) “unlucky” vectors of coefficients in the construction of — when for some intersects the ramification locus of — is contained in a Zariski closed set of real codimension 1 in the ambient (real) coefficient space (with each contributing two coordinates and ). This implies that with generic choices of , the real line segment path running over for a small is “lucky”. ∎
Lemma 4**.**
With the notation as above, if the coefficients have sufficiently generic valuations, then
[TABLE]
Proof.
The inclusion is always true by the definition of tropical varieties, but the containment may be strict. However, when the coefficients of have sufficiently generic valuations, the intersections are transverse, so the result follows from [BJS*+*07, Lemma 3.2]. ∎
The multiplicities of the transverse intersection points can be computed using integer linear algebra [JY16, formula (3)].
For ’s that are not monomial, that is, for polynomials that are generic with respect to non-monomial supports, the intersection of tropicalizations is not necessarily transverse. They become transverse to after the reformulation (1). See Figure 1.
The valuations of the Puiseux series roots are provided by tropical computations (see Theorem 2), while the leading coefficients at a tropical point are given by the root(s) of the -initial ideal
[TABLE]
The next statement follows from [JY16, Lemma 3.1].
Lemma 5**.**
Let , then
[TABLE]
The roots give us the leading coefficients of Puiseux series roots. The lemma shows that the may be generated by simpler polynomials than the original polynomials in . In the original polyhedral homotopies the ideal is binomial and all roots are regular and easy to obtain.
Remark 6**.**
In general, distinct Puiseux series roots may share the same leading terms. This happens exactly when a root is multiple.
For a multiple root , one needs to find more terms in the truncated Puiseux series
[TABLE]
so that these are distinct for distinct Puiseux series roots. Note that all and ’s are polynomials in for some positive integer .
The most comprehensive algorithmic treatment of this is can be found in [JMM08]. ∎
Algorithm 7** (Main algorithm).**
0:* , a collection of polynomials;, , sets of monomials. *
0:* generic for ;approximations to the points of . *
1:* Compute . *
2:* Construct , of the form , . *
3:* Compute . *
4:* *
5:* for every point do
6:* Construct truncations of Puiseux series roots of . *
7:* for do
8:* Pick close to [math], let . *
9:* , where is the output of a homotopy continuation algorithm tracking roots of , , starting with . *
10:* end *for
11:* return , , and . *
12:* end ***for **
The following points remained unsaid in the pseudocode:
- •
The software Gfan [Jen] can compute in Line 1 when has rational coefficients.
- •
The computation of may involve Gröbner bases, while the transverse intersection in Line 3 does not.
- •
Jensen’s tropical homotopy continuation [Jen16] may be useful for Line 3.
- •
Solutions to provide the leading coefficients of . They give distinct in Line 6, unless some solutions are multiple. See Remark 6.
- •
The ideal is a byproduct of the computation of in Line 1.
- •
A practical way to pin down in Line 8 is out of the scope of this article; we envision obtaining with heuristics that depend on and . Such exists according to Remark 3 and the discussion preceding it.
- •
We also purposefully omit the discussion of how one tracks a homotopy path in practice. The machinery of numerical homotopy continuation is well established with several books (e.g., [Mor87, AG03] and several more modern ones) devoted to its details.
Algorithm 7 achieves our main goal. The enumerative goal is achieved by executing it until Line 3 and then computing the degrees of . The sum of these degrees is .
4. Conclusion
In Setting B we construct an optimal homotopy, optimal in the sense that the number of homotopy paths is equal to the number of solutions generically. Our method combines symbolic, polyhedral, and numerical parts. One potential strength is that, for a concrete polynomial system, one can distribute the load between these parts to avoid bottlenecks or to decrease the generic solution count.
Indeed, polynomials in any given set of equations can be divided into two groups, and , fixing the variety , on which the roots of are sought. For each polynomial one can decide on which ingredients this polynomial is “made of”. The list of ingredients may be either inherent to the given problem or be a subject of experimentation. This gives a lot of flexibility.
Appendix A Appendix: Base locus
When we say a statement is true “for generic in ” or “for general in ”, we mean that it is true for all elements in some Zariski open dense subset of , not to be confused with the “generic point” of a scheme. We say that a polynomial is generic or general with respect to support if where the coefficients avoid a Zariski closed proper subset of which depends on the context. One should assume that an explicit description of this exceptional set is hard to acquire algorithmically.
In Bernstein–Khovanskii–Kushnirenko theorems and in the original polyhedral homotopies, the goal is to compute (the number of) solutions in the algebraic torus , of polynomials that are generic with respect to their monomial supports. In other words, we choose a monomial basis for each linear system and remove the union of their hypersurfaces from , obtaining if all variables appear in the monomial basis. Here we remove the base locus instead, which is the intersection of the hypersurfaces. The following argument shows that these two settings are equivalent. For generic polynomials the solution set does not depend on the choice of set removed as long as it has smaller dimension and contains the base locus.
Proposition 8**.**
Suppose that is finite for generic . Let be a variety such that and . Then
[TABLE]
for generic .
Proof.
The idea is that for generic choices of , each cuts down the dimension of by one, and we assumed that , so we should have
[TABLE]
which implies (3). We will make this precise.
If , then there is nothing to prove. Suppose . Let be a Zariski open dense subset of such that is finite for each . Let
[TABLE]
We wish to show that contains a Zariski open dense subset of .
Note that is constructible, since for the incidence variety
[TABLE]
with the projection onto we have .
For a point , the set of satisfying is the complement of a hyperplane in . Since the degrees of the ’s in are bounded, the condition that does not contain any component of can be translated as not vanishing on a finite set of points; hence the set of such ’s is a Zariski open dense subset of . By Krull’s Hauptidealsatz, if does not contain any component of , then is either empty or has dimension .
For a fixed , by a similar argument, there is a Zariski open dense subset of such that has dimension , and so on. This shows that the set is dense in . Since it is constructible, it contains a Zariski open dense subset. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[AG 03] E.L. Allgower and K. Georg, Introduction to numerical continuation methods , Classics in Applied Mathematics, vol. 45, SIAM, 2003.
- 2[BJS + 07] T. Bogart, A. N. Jensen, D. Speyer, B. Sturmfels, and R. R. Thomas, Computing tropical varieties , J. Symbolic Comput. 42 (2007), no. 1-2, 54–73. MR 2284285 (2007 j:14103)
- 3[CLL 14] Tianran Chen, Tsung-Lin Lee, and Tien-Yien Li, Hom 4PS-3: a parallel numerical solver for systems of polynomial equations based on polyhedral homotopy continuation methods , International Congress on Mathematical Software, Springer, 2014, pp. 183–190.
- 4[Ghy] Étienne Ghys, A singular mathematical promenade , ar Xiv preprint ar Xiv:1612.06373.
- 5[GKK + 04] Takayuki Gunji, Sunyoung Kim, Masakazu Kojima, Akiko Takeda, Katsuki Fujisawa, and Tomohiko Mizutani, P Ho M: a polyhedral homotopy continuation method for polynomial systems , Computing 73 (2004), no. 1, 57–77.
- 6[HS 95] Birkett Huber and Bernd Sturmfels, A polyhedral method for solving sparse polynomial systems , Math. Comp. 64 (1995), no. 212, 1541–1555. MR 1297471 (95m:65100)
- 7[Jen] Anders N. Jensen, Gfan, a software system for Gröbner fans and tropical varieties , Available at http://home.imf.au.dk/jensen/software/gfan/gfan.html .
- 8[Jen 16] Anders Nedergaard Jensen, Tropical homotopy continuation , 2016.
