{"title":"Punctual Hilbert scheme and certified approximate singularities","authors":"Angelos Mantzaflaris, B. Mourrain, Á. Szántó","doi":"10.1145/3373207.3404024","DOIUrl":"https://doi.org/10.1145/3373207.3404024","url":null,"abstract":"In this paper we provide a new method to certify that a nearby polynomial system has a singular isolated root and we compute its multiplicity structure. More precisely, given a polynomial system f = (f1, ..., fN) ∈ C[x1, ..., xn]N, we present a Newton iteration on an extended deflated system that locally converges, under regularity conditions, to a small deformation of f such that this deformed system has an exact singular root. The iteration simultaneously converges to the coordinates of the singular root and the coefficients of the so-called inverse system that describes the multiplicity structure at the root. We use α-theory test to certify the quadratic convergence, and to give bounds on the size of the deformation and on the approximation error. The approach relies on an analysis of the punctual Hilbert scheme, for which we provide a new description. We show in particular that some of its strata can be rationally parametrized and exploit these parametrizations in the certification. We show in numerical experimentation how the approximate inverse system can be computed as a starting point of the Newton iterations and the fast numerical convergence to the singular root with its multiplicity structure, certified by our criteria.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"5 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"123690899","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Sub-quadratic time for riemann-roch spaces: case of smooth divisors over nodal plane projective curves","authors":"S. Abelard, Alain Couvreur, Grégoire Lecerf","doi":"10.1145/3373207.3404053","DOIUrl":"https://doi.org/10.1145/3373207.3404053","url":null,"abstract":"We revisit the seminal Brill-Noether algorithm in the rather generic situation of smooth divisors over a nodal plane projective curve. Our approach takes advantage of fast algorithms for polynomials and structured matrices. We reach sub-quadratic time for computing a basis of a Riemann-Roch space. This improves upon previously known complexity bounds.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"3 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"126899067","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On the apolar algebra of a product of linear forms","authors":"Michael DiPasquale, Zachary Flores, C. Peterson","doi":"10.1145/3373207.3404014","DOIUrl":"https://doi.org/10.1145/3373207.3404014","url":null,"abstract":"Apolarity is a important tool in commutative algebra and algebraic geometry which studies a form, f, by the action of polynomial differential operators on f. The quotient of all polynomial differential operators by those which annihilate f is called the apolar algebra of f. In general, the apolar algebra of a form is useful for determining its Waring decomposition, which consists of writing the form as a sum of powers of linear forms with as few summands as possible. In this article we study the apolar algebra of a product of linear forms, which generalizes the case of monomials and connects to the geometry of hyperplane arrangements. In the first part of the article we provide a bound on the Waring rank of a product of linear forms under certain genericity assumptions; for this we use the defining equations of so-called star configurations due to Geramita, Harbourne, and Migliore. In the second part of the article we use the computer algebra system Bertini, which operates by homotopy continuation methods, to solve certain rank equations for catalecticant matrices. Our computations suggest that, up to a change of variables, there are exactly six homogeneous polynomials of degree six in three variables which factor completely as a product of linear forms defining an irreducible multi-arrangement and whose apolar algebras have dimension six in degree three. As a consequence of these calculations, we find six cases of such forms with cactus rank six, five of which also have Waring rank six. Among these are products defining subarrangements of the braid and Hessian arrangements.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"19 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"129508659","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Signature-based algorithms for Gröbner bases over tate algebras","authors":"X. Caruso, Tristan Vaccon, Thibaut Verron","doi":"10.1145/3373207.3404035","DOIUrl":"https://doi.org/10.1145/3373207.3404035","url":null,"abstract":"Introduced by Tate in [Ta71], Tate algebras play a major role in the context of analytic geometry over the p-adics, where they act as a counterpart to the use of polynomial algebras in classical algebraic geometry. In [CVV19] the formalism of Gröbner bases over Tate algebras has been introduced and effectively implemented. One of the bottlenecks in the algorithms was the time spent on reduction, which are significantly costlier than over polynomials. In the present article, we introduce two signature-based Gröbner bases algorithms for Tate algebras, in order to avoid many reductions. They have been implemented in SageMath. We discuss their superiority based on numerical evidence.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"20 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125694756","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Cyrille Chenavier, Clemens Hofstadler, C. Raab, G. Regensburger
{"title":"Compatible rewriting of noncommutative polynomials for proving operator identities","authors":"Cyrille Chenavier, Clemens Hofstadler, C. Raab, G. Regensburger","doi":"10.1145/3373207.3404047","DOIUrl":"https://doi.org/10.1145/3373207.3404047","url":null,"abstract":"The goal of this paper is to prove operator identities using equalities between noncommutative polynomials. In general, a polynomial expression is not valid in terms of operators, since it may not be compatible with domains and codomains of the corresponding operators. Recently, some of the authors introduced a framework based on labelled quivers to rigorously translate polynomial identities to operator identities. In the present paper, we extend and adapt the framework to the context of rewriting and polynomial reduction. We give a sufficient condition on the polynomials used for rewriting to ensure that standard polynomial reduction automatically respects domains and codomains of operators. Finally, we adapt the noncommutative Buchberger procedure to compute additional compatible polynomials for rewriting. In the package OperatorGB, we also provide an implementation of the concepts developed.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"42 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122075172","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Efficient ECM factorization in parallel with the lyness map","authors":"A. Hone","doi":"10.1145/3373207.3404044","DOIUrl":"https://doi.org/10.1145/3373207.3404044","url":null,"abstract":"The Lyness map is a birational map in the plane which provides one of the simplest discrete analogues of a Hamiltonian system with one degree of freedom, having a conserved quantity and an invariant symplectic form. As an example of a symmetric Quispel-Roberts-Thompson (QRT) map, each generic orbit of the Lyness map lies on a curve of genus one, and corresponds to a sequence of points on an elliptic curve which is one of the fibres in a pencil of biquadratic curves in the plane. Here we present a version of the elliptic curve method (ECM) for integer factorization, which is based on iteration of the Lyness map with a particular choice of initial data. More precisely, we give an algorithm for scalar multiplication of a point on an arbitrary elliptic curve over Q, which is represented by one of the curves in the Lyness pencil. In order to avoid field inversion (I), and require only field multiplication (M), squaring (S) and addition, projective coordinates in P1 ×P1 are used. Neglecting multiplication by curve constants (assumed small), each addition of the chosen point uses 2M, while each doubling step requires 15M. We further show that the doubling step can be implemented efficiently in parallel with four processors, dropping the effective cost to 4M. In contrast, the fastest algorithms in the literature use twisted Edwards curves (equivalent to Montgomery curves), which correspond to a subset of all elliptic curves. Scalar muliplication on twisted Edwards curves with suitable small curve constants uses 8M for point addition and 4M+4S for point doubling, both of which can be run in parallel with four processors to yield effective costs of 2M and 1M + 1S, respectively. Thus our scalar multiplication algorithm should require, on average, roughly twice as many multiplications per bit as state of the art methods using twisted Edwards curves. In our conclusions, we discuss applications where the use of Lyness curves may provide potential advantages.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130479086","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Separating variables in bivariate polynomial ideals","authors":"Manfred Buchacher, Manuel Kauers, G. Pogudin","doi":"10.1145/3373207.3404028","DOIUrl":"https://doi.org/10.1145/3373207.3404028","url":null,"abstract":"We present an algorithm which for any given ideal I ⊆ K[x, y] finds all elements of I that have the form f(x) - g(y), i.e., all elements in which no monomial is a multiple of xy.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"11 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"114887253","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"General witness sets for numerical algebraic geometry","authors":"F. Sottile","doi":"10.1145/3373207.3403995","DOIUrl":"https://doi.org/10.1145/3373207.3403995","url":null,"abstract":"Numerical algebraic geometry has a close relationship to intersection theory from algebraic geometry. We deepen this relationship, explaining how rational or algebraic equivalence gives a homotopy. We present a general notion of witness set for subvarieties of a smooth complete complex algebraic variety using ideas from intersection theory. Under appropriate assumptions, general witness sets enable numerical algorithms such as sampling and membership. These assumptions hold for products of flag manifolds. We introduce Schubert witness sets, which provide general witness sets for Grassmannians and flag manifolds.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"519 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"116262188","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pascal Giorgi, Bruno Grenet, Armelle Perret du Cray
{"title":"Essentially optimal sparse polynomial multiplication","authors":"Pascal Giorgi, Bruno Grenet, Armelle Perret du Cray","doi":"10.1145/3373207.3404026","DOIUrl":"https://doi.org/10.1145/3373207.3404026","url":null,"abstract":"We present a probabilistic algorithm to compute the product of two univariate sparse polynomials over a field with a number of bit operations that is quasi-linear in the size of the input and the output. Our algorithm works for any field of characteristic zero or larger than the degree. We mainly rely on sparse interpolation and on a new algorithm for verifying a sparse product that has also a quasi-linear time complexity. Using Kronecker substitution techniques we extend our result to the multivariate case.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"121131115","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On fast multiplication of a matrix by its transpose","authors":"J. Dumas, Clément Pernet, A. Sedoglavic","doi":"10.1145/3373207.3404021","DOIUrl":"https://doi.org/10.1145/3373207.3404021","url":null,"abstract":"We present a non-commutative algorithm for the multiplication of a 2 × 2-block-matrix by its transpose using 5 block products (3 recursive calls and 2 general products) over C or any field of prime characteristic. We use geometric considerations on the space of bilinear forms describing 2 × 2 matrix products to obtain this algorithm and we show how to reduce the number of involved additions. The resulting algorithm for arbitrary dimensions is a reduction of multiplication of a matrix by its transpose to general matrix product, improving by a constant factor previously known reductions. Finally we propose schedules with low memory footprint that support a fast and memory efficient practical implementation over a prime field. To conclude, we show how to use our result in L · D · LT factorization.","PeriodicalId":186699,"journal":{"name":"Proceedings of the 45th International Symposium on Symbolic and Algebraic Computation","volume":"23 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2020-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"123000259","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}