Source author record

Anton Leykin

Anton Leykin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

$u$-generation: solving systems of polynomials equation-by-equation

We develop a new method that improves the efficiency of equation-by-equation algorithms for solving polynomial systems. Our method is based on a novel geometric construction, and reduces the total number of homotopy paths that must be numerically continued. These improvements may be applied to the basic algorithms of numerical algebraic geometry in the settings of both projective and multiprojective varieties. Our computational experiments demonstrate significant savings obtained on several benchmark systems. We also present an extended case study on maximum likelihood estimation for rank-constrained symmetric $n\times n$ matrices, in which multiprojective $u$-generation allows us to complete the list of ML degrees for $n\le 6.$

preprint2020arXiv

A numerical toolkit for multiprojective varieties

A numerical description of an algebraic subvariety of projective space is given by a general linear section, called a witness set. For a subvariety of a product of projective spaces (a multiprojective variety), the corresponding numerical description is given by a witness collection, whose structure is more involved. We build on recent work to develop a toolkit for the numerical manipulation of multiprojective varieties that operates on witness collections, and use this toolkit in an algorithm for numerical irreducible decomposition of multiprojective varieties. The toolkit and decomposition algorithm are illustrated throughout in a series of examples.

preprint2020arXiv

Noetherian operators and primary decomposition

Noetherian operators are differential operators that encode primary components of a polynomial ideal. We develop a framework, as well as algorithms, for computing Noetherian operators with local dual spaces, both symbolically and numerically. For a primary ideal, such operators provide an alternative representation to one given by a set of generators. This description fits well with numerical algebraic geometry, taking a step toward the goal of numerical primary decomposition.

preprint2020arXiv

Numerical Schubert Calculus via the Littlewood-Richardson Homotopy Algorithm

We develop the Littlewood-Richardson homotopy algorithm, which uses numerical continuation to compute solutions to Schubert problems on Grassmannians and is based on the geometric Littlewood-Richardson rule. One key ingredient of this algorithm is our new optimal formulation of Schubert problems in local Stiefel coordinates as systems of equations. Our implementation can solve problem instances with tens of thousands of solutions.

preprint2020arXiv

PL${}_{1}$P -- Point-line Minimal Problems under Partial Visibility in Three Views

We present a complete classification of minimal problems for generic arrangements of points and lines in space observed partially by three calibrated perspective cameras when each line is incident to at most one point. This is a large class of interesting minimal problems that allows missing observations in images due to occlusions and missed detections. There is an infinite number of such minimal problems; however, we show that they can be reduced to 140616 equivalence classes by removing superfluous features and relabeling the cameras. We also introduce camera-minimal problems, which are practical for designing minimal solvers, and show how to pick a simplest camera-minimal problem for each minimal problem. This simplification results in 74575 equivalence classes. Only 76 of these were known; the rest are new. In order to identify problems that have potential for practical solving of image matching and 3D reconstruction, we present several smaller natural subfamilies of camera-minimal problems as well as compute solution counts for all camera-minimal problems which have less than 300 solutions for generic data.

preprint2016arXiv

Determinantal representations of hyperbolic curves via polynomial homotopy continuation

A smooth curve in the real projective plane is hyperbolic if its ovals are maximally nested. By the Helton-Vinnikov Theorem, any such curve admits a definite symmetric determinantal representation. We use polynomial homotopy continuation to compute such representations numerically. Our method works by lifting paths from the space of hyperbolic polynomials to a branched cover in the space of pairs of symmetric matrices.

preprint2016arXiv

Hilbert series of symmetric ideals in infinite polynomial rings via formal languages

Let $R$ be the polynomial ring $K[x_{i,j}]$ where $1 \le i \le r$ and $j \in \mathbb{N}$, and let $I$ be an ideal of $R$ stable under the natural action of the infinite symmetric group $S_{\infty}$. Nagel--Römer recently defined a Hilbert series $H_I(s,t)$ of $I$ and proved that it is rational. We give a much shorter proof of this theorem using tools from the theory of formal languages and a simple algorithm that computes the series.

preprint2015arXiv

Noetherianity for infinite-dimensional toric varieties

We consider a large class of monomial maps respecting an action of the infinite symmetric group, and prove that the toric ideals arising as their kernels are finitely generated up to symmetry. Our class includes many important examples where Noetherianity was recently proved or conjectured. In particular, our results imply Hillar-Sullivant's Independent Set Theorem and settle several finiteness conjectures due to Aschenbrenner, Martin del Campo, Hillar, and Sullivant. We introduce a matching monoid and show that its monoid ring is Noetherian up to symmetry. Our approach is then to factorize a more general equivariant monomial map into two parts going through this monoid. The kernels of both parts are finitely generated up to symmetry: recent work by Yamaguchi-Ogawa-Takemura on the (generalized) Birkhoff model provides an explicit degree bound for the kernel of the first part, while for the second part the finiteness follows from the Noetherianity of the matching monoid ring.

preprint2014arXiv

Equivariant lattice generators and Markov bases

It has been shown recently that monomial maps in a large class respecting the action of the infinite symmetric group have, up to symmetry, finitely generated kernels. We study the simplest nontrivial family in this class: the maps given by a single monomial. Considering the corresponding lattice map, we explicitly construct an equivariant lattice generating set, whose width (the number of variables necessary to write it down) depends linearly on the width of the map. This result is sharp and improves dramatically the previously known upper bound as it does not depend on the degree of the image monomial. In the case of of width two, we construct an explicit finite set of binomials generating the toric ideal up to symmetry. Both width and degree of this generating set are sharply bounded by linear functions in the exponents of the monomial.

preprint2013arXiv

Bertini for Macaulay2

Numerical algebraic geometry is the field of computational mathematics concerning the numerical solution of polynomial systems of equations. Bertini, a popular software package for computational applications of this field, includes implementations of a variety of algorithms based on polynomial homotopy continuation. The Macaulay2 package Bertini.m2 provides an interface to Bertini, making it possible to access the core run modes of Bertini in Macaulay2. With these run modes, users can find approximate solutions to zero-dimensional systems and positive-dimensional systems, test numerically whether a point lies on a variety, sample numerically from a variety, and perform parameter homotopy runs.

preprint2012arXiv

Robust certified numerical homotopy tracking

We describe, for the first time, a completely rigorous homotopy (path--following) algorithm (in the Turing machine model) to find approximate zeros of systems of polynomial equations. If the coordinates of the input systems and the initial zero are rational our algorithm involves only rational computations and if the homotopy is well posed an approximate zero with integer coordinates of the target system is obtained. The total bit complexity is linear in the length of the path in the condition metric, and polynomial in the logarithm of the maximum of the condition number along the path, and in the size of the input.

preprint2011arXiv

A search for an optimal start system for numerical homotopy continuation

We use our recent implementation of a certified homotopy tracking algorithm to search for start systems that minimize the average complexity of finding all roots of a regular system of polynomial equations. While finding optimal start systems is a hard problem, our experiments show that it is possible to find start systems that deliver better average complexity than the ones that are commonly used in the existing homotopy continuation software.

preprint2011arXiv

Computing localizations iteratively

Let $R=\bC[\bfx]$ be a polynomial ring with complex coefficients and $\Dx = \bC<bfx,\bfp>$ be the Weyl algebra. Describing the localization $R_f = R[f^{-1}]$ for nonzero $f\in R$ as a $\Dx$-module amounts to computing the annihilator $A = \Ann(f^a)\subset \Dx$ of the cyclic generator $f^{a}$ for a suitable negative integer $a$. We construct an iterative algorithm that uses truncated annihilators to build $A$ for planar curves.

preprint2011arXiv

Numerical Algebraic Geometry for Macaulay2

Numerical Algebraic Geometry uses numerical data to describe algebraic varieties. It is based on the methods of numerical polynomial homotopy continuation, an alternative to the classical symbolic approaches of computational algebraic geometry. We present a package, the driving idea behind which is to interlink the existing symbolic methods of Macaulay2 and the powerful engine of numerical approximate computations. The core procedures of the package exhibit performance competitive with the other homotopy continuation software.

preprint2010arXiv

Algorithms for Bernstein-Sato polynomials and multiplier ideals

The Bernstein-Sato polynomial (or global b-function) is an important invariant in singularity theory, which can be computed using symbolic methods in the theory of D-modules. After surveying algorithms for computing the global b-function, we develop a new method to compute the local b-function for a single polynomial. We then develop algorithms that compute generalized Bernstein-Sato polynomials of Budur-Mustata-Saito and Shibuta for an arbitrary polynomial ideal. These lead to computations of log canonical thresholds, jumping coefficients, and multiplier ideals. Our algorithm for multiplier ideals simplifies that of Shibuta and shares a common subroutine with our local b-function algorithm. The algorithms we present have been implemented in the D-modules package of the computer algebra system Macaulay2.

preprint2010arXiv

Certified numerical homotopy tracking

Given a homotopy connecting two polynomial systems we provide a rigorous algorithm for tracking a regular homotopy path connecting an approximate zero of the start system to an approximate zero of the target system. Our method uses recent results on the complexity of homotopy continuation rooted in the alpha theory of Smale. Experimental results obtained with the implementation in the numerical algebraic geometry package of Macaulay2 demonstrate the practicality of the algorithm. In particular, we confirm the theoretical results for random linear homotopies and illustrate the plausibility of a conjecture by Shub and Smale on a good initial pair.