Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation最新文献

筛选
英文 中文
GBLA: Gröbner Basis Linear Algebra Package GBLA: Gröbner基础线性代数包
Brice Boyer, C. Eder, J. Faugère, Sylvain Lachartre, Fayssal Martani
{"title":"GBLA: Gröbner Basis Linear Algebra Package","authors":"Brice Boyer, C. Eder, J. Faugère, Sylvain Lachartre, Fayssal Martani","doi":"10.1145/2930889.2930914","DOIUrl":"https://doi.org/10.1145/2930889.2930914","url":null,"abstract":"This is a system paper about a new GPLv2 open source C library GBLA implementing and improving the idea [8] of Faugère and Lachartre (GB reduction). We further exploit underlying structures in matrices generated during Gröbner basis computations in algorithms like F4 or F5 taking advantage of block patterns by using a special data structure called multilines. Moreover, we discuss a new order of operations for the reduction process. In various different experimental results we show that GBLA performs better than GB reduction or Magma in sequential computations (up to 40% faster) and scales much better than GB reduction for a higher number of cores: On 32 cores we reach a scaling of up to 26. GBLA is up to 7 times faster than GB reduction. Further, we compare different parallel schedulers GBLA can be used with. We also developed a new advanced storage format that exploits the fact that our matrices are coming from Gröbner basis computations, shrinking storage by a factor of up to 4. A huge database of our matrices is freely avail- able with GBLA.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"42 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"132756142","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}
引用次数: 11
Computing with Quasiseparable Matrices 拟可分矩阵的计算
Clément Pernet
{"title":"Computing with Quasiseparable Matrices","authors":"Clément Pernet","doi":"10.1145/2930889.2930915","DOIUrl":"https://doi.org/10.1145/2930889.2930915","url":null,"abstract":"The class of quasiseparable matrices is defined by a pair of bounds, called the quasiseparable orders, on the ranks of the sub-matrices entirely located in their strictly lower and upper triangular parts. These arise naturally in applications, as e.g. the inverse of band matrices, and are widely used for they admit structured representations allowing to compute with them in time linear in the dimension. We show, in this paper, the connection between the notion of quasiseparability and the rank profile matrix invariant, presented in [Dumas & al. ISSAC'15]. This allows us to propose an algorithm computing the quasiseparable orders (rL,rU) in time O{n2sω-2} where s=max(rL,rU) and ω the exponent of matrix multiplication. We then present two new structured representations, a binary tree of PLUQ decompositions, and the Bruhat, using respectively O{ns log n/s and O{ns} field elements instead of O{ns2} for the classical generator and O{ns log n} for the hierarchically semiseparable representations. We present algorithms computing these representations in time O{n2sω-2}. These representations allow a matrix-vector product in time linear in the size of their representation. Lastly we show how to multiply two such structured matrices in time O{n2sω-2}.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"20 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125199312","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}
引用次数: 6
Algorithms for Simultaneous Padé Approximations 同时逼近的算法
J. R. N. Nielsen, A. Storjohann
{"title":"Algorithms for Simultaneous Padé Approximations","authors":"J. R. N. Nielsen, A. Storjohann","doi":"10.1145/2930889.2930933","DOIUrl":"https://doi.org/10.1145/2930889.2930933","url":null,"abstract":"We describe how to solve simultaneous Padé approximations over a power series ring K[[x]] for a field K using O~(nω - 1 d) operations in K, where d is the sought precision and $n$ is the number of power series to approximate. We develop two algorithms using different approaches. Both algorithms return a reduced sub-bases that generates the complete set of solutions to the input approximations problem that satisfy the given degree constraints. Our results are made possible by recent breakthroughs in fast computations of minimal approximant bases and Hermite Padé approximations.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"23 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125384666","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}
引用次数: 13
Linear Time Interactive Certificates for the Minimal Polynomial and the Determinant of a Sparse Matrix 最小多项式和稀疏矩阵行列式的线性时间交互证书
J. Dumas, E. Kaltofen, Emmanuel Thomé, G. Villard
{"title":"Linear Time Interactive Certificates for the Minimal Polynomial and the Determinant of a Sparse Matrix","authors":"J. Dumas, E. Kaltofen, Emmanuel Thomé, G. Villard","doi":"10.1145/2930889.2930908","DOIUrl":"https://doi.org/10.1145/2930889.2930908","url":null,"abstract":"Computational problem certificates are additional data structures for each output, which can be used by a---possibly randomized---verification algorithm that proves the correctness of each output. In this paper, we give an algorithm that computes a certificate for the minimal polynomial of sparse or structured matrices over an abstract field, of sufficiently large cardinality, whose Monte Carlo verification complexity requires a single matrix-vector multiplication and a linear number of extra field operations. We also propose a novel preconditioner that ensures irreducibility of the characteristic polynomial of the generically preconditioned matrix. This preconditioner takes linear time to be applied and uses only two random entries. We then combine these two techniques to give algorithms that compute certificates for the determinant, and thus for the characteristic polynomial, whose Monte Carlo verification complexity is therefore also linear.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"18 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122361479","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}
引用次数: 15
Fast Computation of Minimal Interpolation Bases in Popov Form for Arbitrary Shifts 任意移位波波夫形式最小插值基的快速计算
C. Jeannerod, Vincent Neiger, É. Schost, G. Villard
{"title":"Fast Computation of Minimal Interpolation Bases in Popov Form for Arbitrary Shifts","authors":"C. Jeannerod, Vincent Neiger, É. Schost, G. Villard","doi":"10.1145/2930889.2930928","DOIUrl":"https://doi.org/10.1145/2930889.2930928","url":null,"abstract":"We compute minimal bases of solutions for a general interpolation problem, which encompasses Hermite-Pade approximation and constrained multivariate interpolation, and has applications in coding theory and security. This problem asks to find univariate polynomial relations between m vectors of size σ; these relations should have small degree with respect to an input degree shift. For an arbitrary shift, we propose an algorithm for the computation of an interpolation basis in shifted Popov normal form with a cost of O~(mω-1 σ) field operations, where ω is the exponent of matrix multiplication and the notation O~(·) indicates that logarithmic terms are omitted. Earlier works, in the case of Hermite-Pade approximation and in the general interpolation case, compute non-normalized bases. Since for arbitrary shifts such bases may have size Θ(m2 σ), the cost bound O~(mω-1 σ) was feasible only with restrictive assumptions on the shift that ensure small output sizes. The question of handling arbitrary shifts with the same complexity bound was left open. To obtain the target cost for any shift, we strengthen the properties of the output bases, and of those obtained during the course of the algorithm: all the bases are computed in shifted Popov form, whose size is always O(m σ). Then, we design a divide-and-conquer scheme. We recursively reduce the initial interpolation problem to sub-problems with more convenient shifts by first computing information on the degrees of the intermediate bases.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"123 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124190003","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}
引用次数: 23
Fast Computation of Shifted Popov Forms of Polynomial Matrices via Systems of Modular Polynomial Equations 用模多项式方程组快速计算多项式矩阵的移位波波夫形式
Vincent Neiger
{"title":"Fast Computation of Shifted Popov Forms of Polynomial Matrices via Systems of Modular Polynomial Equations","authors":"Vincent Neiger","doi":"10.1145/2930889.2930936","DOIUrl":"https://doi.org/10.1145/2930889.2930936","url":null,"abstract":"We give a Las Vegas algorithm which computes the shifted Popov form of an m x m nonsingular polynomial matrix of degree d in expected ~O(mω d) field operations, where ω is the exponent of matrix multiplication and ~O(·) indicates that logarithmic factors are omitted. This is the first algorithm in ~O(mω d) for shifted row reduction with arbitrary shifts. Using partial linearization, we reduce the problem to the case d ≤ ⌈ σ/m ⌉ where σ is the generic determinant bound, with σ / m bounded from above by both the average row degree and the average column degree of the matrix. The cost above becomes ~O(mω ⌈ σ/m ⌉), improving upon the cost of the fastest previously known algorithm for row reduction, which is deterministic. Our algorithm first builds a system of modular equations whose solution set is the row space of the input matrix, and then finds the basis in shifted Popov form of this set. We give a deterministic algorithm for this second step supporting arbitrary moduli in ~O(mω-1 σ) field operations, where m is the number of unknowns and σ is the sum of the degrees of the moduli. This extends previous results with the same cost bound in the specific cases of order basis computation and M-Pade approximation, in which the moduli are products of known linear factors.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"11 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130235261","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}
引用次数: 22
Fast Computation of the Nth Term of an Algebraic Series over a Finite Prime Field 有限素数域上代数级数第n项的快速计算
A. Bostan, G. Christol, P. Dumas
{"title":"Fast Computation of the Nth Term of an Algebraic Series over a Finite Prime Field","authors":"A. Bostan, G. Christol, P. Dumas","doi":"10.1145/2930889.2930904","DOIUrl":"https://doi.org/10.1145/2930889.2930904","url":null,"abstract":"We address the question of computing one selected term of an algebraic power series. In characteristic zero, the best algorithm currently known for computing the~Nth coefficient of an algebraic series uses differential equations and has arithmetic complexity quasi-linear in √N. We show that over a prime field of positive characteristic p, the complexity can be lowered to O(log N). The mathematical basis for this dramatic improvement is a classical theorem stating that a formal power series with coefficients in a finite field is algebraic if and only if the sequence of its coefficients can be generated by an automaton. We revisit and enhance two constructive proofs of this result for finite prime fields. The first proof uses Mahler equations, whose sizes appear to be prohibitively large. The second proof relies on diagonals of rational functions; we turn it into an efficient algorithm, of complexity linear in log N and quasi-linear in p.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"57 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"128315060","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}
引用次数: 7
Solving Rank-Constrained Semidefinite Programs in Exact Arithmetic 用精确算法求解秩约束半定规划
Simone Naldi
{"title":"Solving Rank-Constrained Semidefinite Programs in Exact Arithmetic","authors":"Simone Naldi","doi":"10.1145/2930889.2930925","DOIUrl":"https://doi.org/10.1145/2930889.2930925","url":null,"abstract":"We consider the problem of minimizing a linear function over an affine section of the cone of positive semidefinite matrices, with the additional constraint that the feasible matrix has prescribed rank. When the rank constraint is active, this is a non-convex optimization problem, otherwise it is a semidefinite program. Both find numerous applications especially in systems control theory and combinatorial optimization, but even in more general contexts such as polynomial optimization or real algebra. While numerical algorithms exist for solving this problem, such as interior-point or Newton-like algorithms, in this paper we propose an approach based on symbolic computation. We design an exact algorithm for solving rank-constrained semidefinite programs, whose complexity is essentially quadratic on natural degree bounds associated to the given optimization problem: for subfamilies of the problem where the size of the feasible matrix is fixed, the complexity is polynomial in the number of variables. The algorithm works under assumptions on the input data: we prove that these assumptions are generically satisfied. We also implement it in Maple and discuss practical experiments.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"255 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122153978","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}
引用次数: 11
Reduction-Based Creative Telescoping for Algebraic Functions 基于约简的代数函数创造性伸缩
Shaoshi Chen, Manuel Kauers, C. Koutschan
{"title":"Reduction-Based Creative Telescoping for Algebraic Functions","authors":"Shaoshi Chen, Manuel Kauers, C. Koutschan","doi":"10.1145/2930889.2930901","DOIUrl":"https://doi.org/10.1145/2930889.2930901","url":null,"abstract":"Continuing a series of articles in the past few years on creative telescoping using reductions, we develop a new algorithm to construct minimal telescopers for algebraic functions. This algorithm is based on Trager's Hermite reduction and on polynomial reduction, which was originally designed for hyperexponential functions and extended to the algebraic case in this paper.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"22 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125428582","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}
引用次数: 32
Validating the Completeness of the Real Solution Set of a System of Polynomial Equations 验证多项式方程组实解集的完备性
Daniel A. Brake, J. Hauenstein, Alan C. Liddell
{"title":"Validating the Completeness of the Real Solution Set of a System of Polynomial Equations","authors":"Daniel A. Brake, J. Hauenstein, Alan C. Liddell","doi":"10.1145/2930889.2930910","DOIUrl":"https://doi.org/10.1145/2930889.2930910","url":null,"abstract":"Computing the real solutions to a system of polynomial equations is a challenging problem, particularly verifying that all solutions have been computed. We describe an approach that combines numerical algebraic geometry and sums of squares programming to test whether a given set is \"complete\" with respect to the real solution set. Specifically, we test whether the Zariski closure of that set is indeed equal to the solution set of the real radical of the ideal generated by the given polynomials. Examples with finitely and infinitely many real solutions are provided, along with an example having polynomial inequalities.","PeriodicalId":169557,"journal":{"name":"Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation","volume":"7 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2016-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"127285591","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}
引用次数: 15
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
相关产品
×
本文献相关产品
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信