Search

We introduce the concept of monodromy coordinates for representing solutions to large polynomial systems. Representing solutions this way provides a time-memory trade-off in a monodromy solving algorithm. We describe an algorithm, which interpolates the usual monodromy solving algorithm, for computing such a representation and analyze its space and time complexity., Comment: 9 pages, 3 Figures

Area variables are intrinsic to connection formulations of general relativity, in contrast to the fundamental length variables prevalent in metric formulations. Within 4D discrete gravity, particularly based on triangulations, the area-length system establishes a relationship between area variables associated with triangles and the edge length variables. This system is comprised of polynomial equations derived from Heron's formula, which relates the area of a triangle to its edge lengths. Using tools from numerical algebraic geometry, we study the area-length systems. In particular, we show that given the ten triangular areas of a single 4-simplex, there could be up to 64 compatible sets of edge lengths. Moreover, we show that these 64 solutions do not, in general, admit formulae in terms of the areas by analyzing the Galois group, or monodromy group, of the problem. We show that by introducing additional symmetry constraints, it is possible to obtain such formulae for the edge lengths. We take the first steps toward applying our results within discrete quantum gravity, specifically for effective spin foam models., Comment: 15 pages, 9 figures, 4 tables: New Version matches published article on Class. Quant. Grav. Journal

We introduce the n-th Heron variety as the realization space of the (squared) volumes of faces of an n-simplex. Our primary goal is to understand the extent to which Heron's formula, which expresses the area of a triangle as a function of its three edge lengths, can be generalized. Such a formula for one face volume of an n-simplex in terms of other face volumes expresses a dependence in the algebraic matroid of the Heron variety. Whether the volume is expressible in terms of radicals is controlled by the monodromy groups of the coordinate projections of the Heron variety onto coordinates of bases. We discuss a suite of algorithms, some new, for determining these matroids and monodromy groups. We apply these algorithms toward the smaller Heron varieties, organize our findings, and interpret the results in the context of our original motivation., Comment: 25 pages, 11 figures

We express the maximum likelihood (ML) degrees of a family toric varieties in terms of Mobius invariants of matroids. The family of interest are those parametrized by monomial maps given by Lawrence lifts of totally unimodular matrices with even circuits. Specifying these matrices to be vertex-edge incidence matrices of bipartite graphs gives the ML degrees of some hierarchical models and three dimensional quasi-independence models. Included in this list are the no-three-way interaction models with one binary random variable, for which, we give closed formulae., Comment: 23 pages, 5 figures. Comments welcome!

It is a classical result that there are $12$ (irreducible) rational cubic curves through $8$ generic points in $\mathbb{P}_{\mathbb{C}}^2$, but little is known about the non-generic cases. The space of $8$-point configurations is partitioned into strata depending on combinatorial objects we call quatroids, a higher-order version of representable matroids. We compute all $779777$ quatroids on eight distinct points in the plane, which produces a full description of the stratification. For each stratum, we generate several invariants, including the number of rational cubics through a generic configuration. As a byproduct of our investigation, we obtain a collection of results regarding the base loci of pencils of cubics and positive certificates for non-rationality., Comment: 34 pages, 11 figures, 5 tables. Comments are welcome!

OSCAR is an innovative new computer algebra system which combines and extends the power of its four cornerstone systems - GAP (group theory), Singular (algebra and algebraic geometry), Polymake (polyhedral geometry), and Antic (number theory). Assuming little familiarity with the subject, we give an introduction to computations in polyhedral geometry using OSCAR, as a chapter of the upcoming OSCAR book. In particular, we define polytopes, polyhedra, and polyhedral fans, and we give a brief overview about computing convex hulls and solving linear programs. Three detailed case studies are left for experts in polyhedral geometry. These are concerned with face numbers of random polytopes, constructions and properties of Gelfand-Tsetlin polytopes, and secondary polytopes., Comment: 20 pages, 8 figures

The maximum function, on vectors of real numbers, is not differentiable. Consequently, several differentiable approximations of this function are popular substitutes. We survey three smooth functions which approximate the maximum function and analyze their convergence rates. We interpret these functions through the lens of tropical geometry, where their performance differences are geometrically salient. As an application, we provide an algorithm which computes the max-convolution of two integer vectors in quasi-linear time. We show this algorithm's power in computing adjacent sums within a vector as well as computing service curves in a network analysis application., Comment: 24 pages, 21 Figures, 2 Tables

We establish how the coefficients of a sparse polynomial system influence the sum (or the trace) of its zeros. As an application, we develop numerical tests for verifying whether a set of solutions to a sparse system is complete. These algorithms extend the classical trace test in numerical algebraic geometry. Our results rely on both the analysis of the structure of sparse resultants as well as an extension of Esterov's results on monodromy groups of sparse systems.

Computing all critical points of a monomial on a very affine variety is a fundamental task in algebraic statistics, particle physics and other fields. The number of critical points is known as the maximum likelihood (ML) degree. When the variety is smooth, it coincides with the Euler characteristic. We introduce degeneration techniques that are inspired by the soft limits in CEGM theory, and we answer several questions raised in the physics literature. These pertain to bounded regions in discriminantal arrangements and to moduli spaces of point configurations. We present theory and practise, connecting complex geometry, tropical combinatorics, and numerical nonlinear algebra., Comment: 33 pages, updated to reflect reviewers' comments and added link to Zenodo to certify numerical results

We introduce a new algorithm computing the characteristic polynomials of hyperplane arrangements which exploits their underlying symmetry groups. Our algorithm counts the chambers of an arrangement as a byproduct of computing its characteristic polynomial. We showcase our julia implementation, based on OSCAR, on examples coming from hyperplane arrangements with applications to physics and computer science., Comment: 21 pages, 6 figures, 8 tables, 1 appendix. Updated version based on referee reports

We classify transversal quintic spectrahedra by the location of 20 nodes on the respective real determinantal surface of degree 5. We identify 65 classes of such surfaces and find an explicit representative in each of them.

We examine quadratic surfaces in 3-space that are tangent to nine given figures. These figures can be points, lines, planes or quadrics. The numbers of tangent quadrics were determined by Hermann Schubert in 1879. We study the associated systems of polynomial equations, also in the space of complete quadrics, and we solve them using certified numerical methods. Our aim is to show that Schubert's problems are fully real., Comment: 13 pages

The Macaulay2 package DecomposableSparseSystems implements methods for studying and numerically solving decomposable sparse polynomial systems. We describe the structure of decomposable sparse systems and explain how the methods in this package may be used to exploit this structure, with examples., Comment: 7 pages, software available at https://www.math.tamu.edu/~thomasjyahl/research/DSS/DSSsite.html

We develop a collection of numerical algorithms which connect ideas from polyhedral geometry and algebraic geometry. The first algorithm we develop functions as a numerical oracle for the Newton polytope of a hypersurface and is based on ideas of Hauenstein and Sottile. Additionally, we construct a numerical tropical membership algorithm which uses the former algorithm as a subroutine. Based on recent results of Esterov, we give an algorithm which recursively solves a sparse polynomial system when the support of that system is either lacunary or triangular. Prior to explaining these results, we give necessary background on polytopes, algebraic geometry, monodromy groups of branched covers, and numerical algebraic geometry., Comment: 150 pages, 65 figures, contains content from arXiv:1811.12279 and arXiv:2001.04228

Amendola et al. proposed a method for solving systems of polynomial equations lying in a family which exploits a recursive decomposition into smaller systems. A family of systems admits such a decomposition if and only if the corresponding Galois group is imprimitive. When the Galois group is imprimitive we consider the problem of computing an explicit decomposition. A consequence of Esterov's classification of sparse polynomial systems with imprimitive Galois groups is that this decomposition is obtained by inspection. This leads to a recursive algorithm to solve decomposable sparse systems, which we present and give evidence for its efficiency., Comment: 20 pages

We compute the degree of Stiefel manifolds, that is, the variety of orthonormal frames in a finite dimensional vector space. Our approach employs techniques from classical algebraic geometry, algebraic combinatorics, and classical invariant theory., Comment: 24 pages, Final version accepted in ECA

We present our implementation of an algorithm which functions as a numerical oracle for the Newton polytope of a hypersurface in the Macaulay2 package NumericalNP.m2. We propose a tropical membership test, relying on this algorithm, for higher codimension varieties based on ideas from Hept and Theobald. To showcase this software, we investigate the Newton polytope of both a hypersurface coming from algebraic vision and the Luroth invariant., Comment: 16 pages, 5 figures

We consider the problem of geometrically approximating a complex analytic curve in the plane by the image of a polynomial parametrization $t \mapsto (x_1(t),x_2(t))$ of bidegree $(d_1,d_2)$. We show the number of such curves is the number of primitive necklaces on $d_1$ white beads and $d_2$ black beads. We show that this number is odd when $d_1=d_2$ is squarefree and use this to give a partial solution to a conjecture by Rababah. Our results naturally extend to a generalization regarding hypersurfaces in higher dimensions. There, the number of parametrized curves of multidegree $(d_1,\ldots,d_n)$ which optimally osculate a given hypersurface are counted by the number of primitive necklaces with $d_i$ beads of color $i$.

We provide a closed formula for the degree of $\text{SO}(n)$ over an algebraically closed field of characteristic zero. In addition, we describe symbolic and numerical techniques which can also be used to compute the degree of $\text{SO}(n)$ for small values of $n$. As an application of our results, we give a formula for the number of critical points of a low-rank semidefinite programming optimization problem. Finally, we provide some evidence for a conjecture regarding the real locus of $\text{SO}(n)$., Comment: 21 pages, 3 figures