Pub Date : 2025-11-17DOI: 10.1016/j.jsc.2025.102528
Laura Casabella , Michael Joswig , Rafael Mohr
We study certain generic systems of real polynomial equations associated with triangulations of convex polytopes and investigate their number of real solutions. Our main focus is set on pairs of plane algebraic curves which form a so-called Wronski system. The computational tasks arising in the analysis of such Wronski pairs lead us to the frontiers of current computer algebra algorithms and their implementations, both via Gröbner bases and numerical algebraic geometry.
{"title":"Wronski pairs of honeycomb curves","authors":"Laura Casabella , Michael Joswig , Rafael Mohr","doi":"10.1016/j.jsc.2025.102528","DOIUrl":"10.1016/j.jsc.2025.102528","url":null,"abstract":"<div><div>We study certain generic systems of real polynomial equations associated with triangulations of convex polytopes and investigate their number of real solutions. Our main focus is set on pairs of plane algebraic curves which form a so-called Wronski system. The computational tasks arising in the analysis of such Wronski pairs lead us to the frontiers of current computer algebra algorithms and their implementations, both via Gröbner bases and numerical algebraic geometry.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"135 ","pages":"Article 102528"},"PeriodicalIF":1.1,"publicationDate":"2025-11-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145580224","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-17DOI: 10.1016/j.jsc.2025.102530
Hoon Hong, Irina A. Kogan
We develop a theory and an algorithm for constructing minimal-degree polynomial moving frames for polynomial curves in an affine space. The algorithm is equivariant under volume-preserving affine transformations of the ambient space and parameter shifts. We show that any matrix-completion algorithm can be turned into an equivariant moving frame algorithm via an equivariantization procedure that we develop. We prove that if a matrix-completion algorithm is of minimal degree, so is the resulting equivariant moving frame algorithm. We propose a novel minimal-degree matrix-completion algorithm, complementing the existing body of literature on this topic.
{"title":"Equi-affine minimal-degree moving frames for polynomial curves","authors":"Hoon Hong, Irina A. Kogan","doi":"10.1016/j.jsc.2025.102530","DOIUrl":"10.1016/j.jsc.2025.102530","url":null,"abstract":"<div><div>We develop a theory and an algorithm for constructing minimal-degree polynomial moving frames for polynomial curves in an affine space. The algorithm is equivariant under volume-preserving affine transformations of the ambient space and parameter shifts. We show that any matrix-completion algorithm can be turned into an equivariant moving frame algorithm via an equivariantization procedure that we develop. We prove that if a matrix-completion algorithm is of minimal degree, so is the resulting equivariant moving frame algorithm. We propose a novel minimal-degree matrix-completion algorithm, complementing the existing body of literature on this topic.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"135 ","pages":"Article 102530"},"PeriodicalIF":1.1,"publicationDate":"2025-11-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145625260","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-17DOI: 10.1016/j.jsc.2025.102525
Marcin Dumnicki , Mikołaj Le Van , Grzegorz Malara , Tomasz Szemberg , Justyna Szpond , Halszka Tutaj-Gasińska
The purpose of the present note is to provide a new proof of the well-known result due to Hartshorne and Hirschowitz to the effect that general lines in projective spaces have good postulation. Our approach uses specialization to a hyperplane and thus opens a door to study postulation of general codimension 2 linear subspaces in projective spaces.
{"title":"Postulation of lines in P3 revisited","authors":"Marcin Dumnicki , Mikołaj Le Van , Grzegorz Malara , Tomasz Szemberg , Justyna Szpond , Halszka Tutaj-Gasińska","doi":"10.1016/j.jsc.2025.102525","DOIUrl":"10.1016/j.jsc.2025.102525","url":null,"abstract":"<div><div>The purpose of the present note is to provide a new proof of the well-known result due to Hartshorne and Hirschowitz to the effect that general lines in projective spaces have <em>good postulation</em>. Our approach uses specialization to a hyperplane and thus opens a door to study postulation of general codimension 2 linear subspaces in projective spaces.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"135 ","pages":"Article 102525"},"PeriodicalIF":1.1,"publicationDate":"2025-11-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145555331","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-14DOI: 10.1016/j.jsc.2025.102526
Henri Lombardi , Stefan Neuwirth , Ihsen Yengui
Let M be a finitely generated submodule of a free module over a multivariate polynomial ring with coefficients in a discrete coherent ring. We prove that its module of leading terms is countably generated and give an algorithm for computing explicitly a generating set. This result is also useful when is not finitely generated. Suppose that the base ring is strongly discrete coherent. We provide a Buchberger-like algorithm and prove that it converges if and only if the module of leading terms is finitely generated. We also state a constructive version of Hilbert's syzygy theorem by following Schreyer's method.
{"title":"Generalised Buchberger and Schreyer algorithms for strongly discrete coherent rings","authors":"Henri Lombardi , Stefan Neuwirth , Ihsen Yengui","doi":"10.1016/j.jsc.2025.102526","DOIUrl":"10.1016/j.jsc.2025.102526","url":null,"abstract":"<div><div>Let <em>M</em> be a finitely generated submodule of a free module over a multivariate polynomial ring with coefficients in a discrete coherent ring. We prove that its module <span><math><mrow><mi>MLT</mi></mrow><mo>(</mo><mi>M</mi><mo>)</mo></math></span> of leading terms is countably generated and give an algorithm for computing explicitly a generating set. This result is also useful when <span><math><mrow><mi>MLT</mi></mrow><mo>(</mo><mi>M</mi><mo>)</mo></math></span> is not finitely generated. Suppose that the base ring is strongly discrete coherent. We provide a Buchberger-like algorithm and prove that it converges if and only if the module of leading terms is finitely generated. We also state a constructive version of Hilbert's syzygy theorem by following Schreyer's method.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"135 ","pages":"Article 102526"},"PeriodicalIF":1.1,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145625259","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-14DOI: 10.1016/j.jsc.2025.102527
Rui-Juan Jing , Marc Moreno Maza , Chirantan Mukherjee , Yan-Feng Xie , Chun-Ming Yuan
Fourier-Motzkin elimination is a fundamental operation in polyhedral geometry. It can be performed by several equivalent procedures, and can be regarded as an adaptation of Gaussian elimination to systems of linear inequalities. These procedures tend to generate large numbers of redundant inequalities. Efficiently detecting these redundancies is essential for obtaining software implementation of practical interest. In this paper, we propose a novel detection technique. We demonstrate its benefits over alternative approaches. A detailed experimentation is reported.
{"title":"Efficient detection of redundancies in systems of linear inequalities","authors":"Rui-Juan Jing , Marc Moreno Maza , Chirantan Mukherjee , Yan-Feng Xie , Chun-Ming Yuan","doi":"10.1016/j.jsc.2025.102527","DOIUrl":"10.1016/j.jsc.2025.102527","url":null,"abstract":"<div><div>Fourier-Motzkin elimination is a fundamental operation in polyhedral geometry. It can be performed by several equivalent procedures, and can be regarded as an adaptation of Gaussian elimination to systems of linear inequalities. These procedures tend to generate large numbers of redundant inequalities. Efficiently detecting these redundancies is essential for obtaining software implementation of practical interest. In this paper, we propose a novel detection technique. We demonstrate its benefits over alternative approaches. A detailed experimentation is reported.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"135 ","pages":"Article 102527"},"PeriodicalIF":1.1,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145737433","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-14DOI: 10.1016/j.jsc.2025.102531
John M. Campbell
The first known q-analogues for any of the 17 formulas for due to Ramanujan were introduced in 2018 by Guo and Liu (2018), via the q-Wilf–Zeilberger method. Through a “normalization” method, which we refer to as EKHAD-normalization, based on the q-polynomial coefficients involved in first-order difference equations obtained from the q-version of Zeilberger's algorithm, we introduce q-WZ pairs that extend WZ pairs introduced by Guillera, 2002, Guillera, 2006. We apply our EKHAD-normalization method to prove four new q-analogues for three of Ramanujan's formulas for along with q-analogues of Guillera's first two series for . Our normalization method does not seem to have been previously considered in any equivalent way in relation to q-series, and this is substantiated through our survey on previously known q-analogues of Ramanujan-type series for and of Guillera's series for . We conclude by showing how our method can be adapted to further extend Guillera's WZ pairs by introducing hypergeometric expansions for .
{"title":"q-analogues of π-formulas due to Ramanujan and Guillera","authors":"John M. Campbell","doi":"10.1016/j.jsc.2025.102531","DOIUrl":"10.1016/j.jsc.2025.102531","url":null,"abstract":"<div><div>The first known <em>q</em>-analogues for any of the 17 formulas for <span><math><mfrac><mrow><mn>1</mn></mrow><mrow><mi>π</mi></mrow></mfrac></math></span> due to Ramanujan were introduced in 2018 by <span><span>Guo and Liu (2018)</span></span>, via the <em>q</em>-Wilf–Zeilberger method. Through a “normalization” method, which we refer to as <span>EKHAD</span>-normalization, based on the <em>q</em>-polynomial coefficients involved in first-order difference equations obtained from the <em>q</em>-version of Zeilberger's algorithm, we introduce <em>q</em>-WZ pairs that extend WZ pairs introduced by <span><span>Guillera, 2002</span></span>, <span><span>Guillera, 2006</span></span>. We apply our <span>EKHAD</span>-normalization method to prove four new <em>q</em>-analogues for three of Ramanujan's formulas for <span><math><mfrac><mrow><mn>1</mn></mrow><mrow><mi>π</mi></mrow></mfrac></math></span> along with <em>q</em>-analogues of Guillera's first two series for <span><math><mfrac><mrow><mn>1</mn></mrow><mrow><msup><mrow><mi>π</mi></mrow><mrow><mn>2</mn></mrow></msup></mrow></mfrac></math></span>. Our normalization method does not seem to have been previously considered in any equivalent way in relation to <em>q</em>-series, and this is substantiated through our survey on previously known <em>q</em>-analogues of Ramanujan-type series for <span><math><mfrac><mrow><mn>1</mn></mrow><mrow><mi>π</mi></mrow></mfrac></math></span> and of Guillera's series for <span><math><mfrac><mrow><mn>1</mn></mrow><mrow><msup><mrow><mi>π</mi></mrow><mrow><mn>2</mn></mrow></msup></mrow></mfrac></math></span>. We conclude by showing how our method can be adapted to further extend Guillera's WZ pairs by introducing hypergeometric expansions for <span><math><mfrac><mrow><mn>1</mn></mrow><mrow><msup><mrow><mi>π</mi></mrow><mrow><mn>2</mn></mrow></msup></mrow></mfrac></math></span>.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"134 ","pages":"Article 102531"},"PeriodicalIF":1.1,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145568784","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-15DOI: 10.1016/j.jsc.2025.102523
Jean-Guillaume Dumas, Bruno Grenet
<div><div>This paper deals with simultaneously fast and in-place algorithms for formulae where the result has to be linearly accumulated: some output variables are also input variables, linked by a linear dependency. Fundamental examples include the in-place accumulated multiplication of polynomials or matrices, <span><math><mi>C</mi><mo>+</mo><mo>=</mo><mi>A</mi><mi>B</mi></math></span> (that is with only <span><math><mi>O</mi><mo>(</mo><mn>1</mn><mo>)</mo></math></span> extra space). The difficulty is to combine in-place computations with fast algorithms: those usually come at the expense of (potentially large) extra temporary space, but with accumulation the output variables are not even available to store intermediate values. We first propose a novel automatic design of fast and in-place accumulating algorithms for any bilinear formulae (and thus for polynomial and matrix multiplication) and then extend it to any linear accumulation of a collection of functions. For this, we relax the in-place model to any algorithm allowed to modify its inputs, provided that those are restored to their initial state afterwards. This allows us to ultimately derive unprecedented in-place accumulating algorithms for fast polynomial multiplications and for Strassen-like matrix multiplications.</div><div>We then consider the simultaneously fast and in-place computation of the Euclidean polynomial modular remainder <span><math><mi>R</mi><mo>(</mo><mi>X</mi><mo>)</mo><mo>≡</mo><mi>A</mi><mo>(</mo><mi>X</mi><mo>)</mo><mspace></mspace><mrow><mi>mod</mi></mrow><mspace></mspace><mi>B</mi><mo>(</mo><mi>X</mi><mo>)</mo></math></span>. Fast algorithms for this usually also come at the expense of a linear amount of extra temporary space. In particular, they require one to first compute and store the whole quotient <span><math><mi>Q</mi><mo>(</mo><mi>X</mi><mo>)</mo></math></span> such that <span><math><mi>A</mi><mo>=</mo><mi>B</mi><mi>Q</mi><mo>+</mo><mi>R</mi></math></span>. We here propose an <em>in-place</em> algorithm to compute the remainder only. If <em>A</em> and <em>B</em> have respective degree <span><math><mi>m</mi><mo>+</mo><mi>n</mi></math></span> and <em>n</em>, and <span><math><mi>M</mi><mo>(</mo><mi>k</mi><mo>)</mo></math></span> denotes the complexity of a (not-in-place) algorithm to multiply two degree-<em>k</em> polynomials, our algorithm uses at most <span><math><mi>O</mi><mrow><mo>(</mo><mfrac><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow></mfrac><mi>M</mi><mo>(</mo><mi>n</mi><mo>)</mo><mi>log</mi><mo></mo><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow></math></span> arithmetic operations. In this particular case this is a factor <span><math><mi>log</mi><mo></mo><mo>(</mo><mi>n</mi><mo>)</mo></math></span> more than the not-in-place algorithm. But if <span><math><mi>M</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>=</mo><mi>Θ</mi><mo>(</mo><msup><mrow><mi>n</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>ϵ</mi></mrow></msup><mo>)</mo></math></span> for some <span><ma
本文同时处理公式的快速和就地算法,其中结果必须是线性累积的:一些输出变量也是输入变量,通过线性依赖关系连接。基本的例子包括多项式或矩阵的就地累积乘法,C+=AB(即只有O(1)个额外空间)。困难之处在于将就地计算与快速算法结合起来:这些计算通常以(可能很大的)额外临时空间为代价,但是随着累积,输出变量甚至无法用于存储中间值。我们首先提出了一种新的快速和就地积累算法的自动设计,适用于任何双线性公式(从而适用于多项式和矩阵乘法),然后将其扩展到函数集合的任何线性积累。为此,我们将就地模型放宽为允许修改其输入的任何算法,前提是这些输入随后恢复到初始状态。这使我们最终能够为快速多项式乘法和Strassen-like矩阵乘法推导出前所未有的就地累积算法。然后我们考虑欧几里得多项式模余数R(X)≡A(X)modB(X)的同时快速就地计算。快速算法通常也以额外的线性临时空间为代价。特别是,它们要求首先计算并存储整个商Q(X),使A=BQ+R。我们在此提出一种只计算余数的就地算法。如果A和B的阶分别为m+n和n, m (k)表示两个k阶多项式相乘的(非原位)算法的复杂度,则我们的算法最多使用O(mnM(n)log (n))个算术运算。在这种特殊情况下,它比非原地算法大log (n)个因子。但是,如果M(n)=Θ(n1+ λ)对于某些ϵ>;0,那么我们的算法确实符合O(mnM(n))的非原位复杂性界限。我们还提出了计算变体-仍然在原地并且具有相同类型的复杂性界限-过位余数A(X)≡A(X)modB(X),累积余数R(X)+=A(X)modB(X)和累积模乘法R(X)+=A(X)C(X)modB(X),即在有限域的多项式扩展中的乘法。为了实现这一点,我们开发了Toeplitz矩阵运算,广义卷积,短乘积和幂级数除法和余数的技术,其输出也是输入的一部分。
{"title":"Fast in-place accumulation","authors":"Jean-Guillaume Dumas, Bruno Grenet","doi":"10.1016/j.jsc.2025.102523","DOIUrl":"10.1016/j.jsc.2025.102523","url":null,"abstract":"<div><div>This paper deals with simultaneously fast and in-place algorithms for formulae where the result has to be linearly accumulated: some output variables are also input variables, linked by a linear dependency. Fundamental examples include the in-place accumulated multiplication of polynomials or matrices, <span><math><mi>C</mi><mo>+</mo><mo>=</mo><mi>A</mi><mi>B</mi></math></span> (that is with only <span><math><mi>O</mi><mo>(</mo><mn>1</mn><mo>)</mo></math></span> extra space). The difficulty is to combine in-place computations with fast algorithms: those usually come at the expense of (potentially large) extra temporary space, but with accumulation the output variables are not even available to store intermediate values. We first propose a novel automatic design of fast and in-place accumulating algorithms for any bilinear formulae (and thus for polynomial and matrix multiplication) and then extend it to any linear accumulation of a collection of functions. For this, we relax the in-place model to any algorithm allowed to modify its inputs, provided that those are restored to their initial state afterwards. This allows us to ultimately derive unprecedented in-place accumulating algorithms for fast polynomial multiplications and for Strassen-like matrix multiplications.</div><div>We then consider the simultaneously fast and in-place computation of the Euclidean polynomial modular remainder <span><math><mi>R</mi><mo>(</mo><mi>X</mi><mo>)</mo><mo>≡</mo><mi>A</mi><mo>(</mo><mi>X</mi><mo>)</mo><mspace></mspace><mrow><mi>mod</mi></mrow><mspace></mspace><mi>B</mi><mo>(</mo><mi>X</mi><mo>)</mo></math></span>. Fast algorithms for this usually also come at the expense of a linear amount of extra temporary space. In particular, they require one to first compute and store the whole quotient <span><math><mi>Q</mi><mo>(</mo><mi>X</mi><mo>)</mo></math></span> such that <span><math><mi>A</mi><mo>=</mo><mi>B</mi><mi>Q</mi><mo>+</mo><mi>R</mi></math></span>. We here propose an <em>in-place</em> algorithm to compute the remainder only. If <em>A</em> and <em>B</em> have respective degree <span><math><mi>m</mi><mo>+</mo><mi>n</mi></math></span> and <em>n</em>, and <span><math><mi>M</mi><mo>(</mo><mi>k</mi><mo>)</mo></math></span> denotes the complexity of a (not-in-place) algorithm to multiply two degree-<em>k</em> polynomials, our algorithm uses at most <span><math><mi>O</mi><mrow><mo>(</mo><mfrac><mrow><mi>m</mi></mrow><mrow><mi>n</mi></mrow></mfrac><mi>M</mi><mo>(</mo><mi>n</mi><mo>)</mo><mi>log</mi><mo></mo><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow></math></span> arithmetic operations. In this particular case this is a factor <span><math><mi>log</mi><mo></mo><mo>(</mo><mi>n</mi><mo>)</mo></math></span> more than the not-in-place algorithm. But if <span><math><mi>M</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>=</mo><mi>Θ</mi><mo>(</mo><msup><mrow><mi>n</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>ϵ</mi></mrow></msup><mo>)</mo></math></span> for some <span><ma","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"134 ","pages":"Article 102523"},"PeriodicalIF":1.1,"publicationDate":"2025-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145362885","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We propose a strategy for the generation of fast and accurate versions of non-commutative recursive matrix multiplication algorithms. To generate these algorithms, we consider matrix and tensor norm bounds governing the stability and accuracy of numerical matrix multiplication. We start by a unification on known max-norm bounds on matrix multiplication stability and then extend them to further norms and more generally to recursive bilinear algorithms and the alternative basis matrix multiplication algorithms. Then our strategy has three phases. First, we reduce those bounds by minimizing a growth factor along the orbits of the associated matrix multiplication tensor decomposition. Second, we develop heuristics that minimize the number of operations required to realize a bilinear formula, while further improving its accuracy. Third, we perform an alternative basis sparsification that improves on the time complexity constant and mostly preserves the overall accuracy. Let denote a non-commutative algorithm multiplying an matrix with an matrix, using r coefficient products. Our strategy allows us to propose a algorithm that reaches simultaneously a better accuracy in practice, compared to previously known such fast Strassen-like algorithms, and a time complexity bound with the best currently known leading term (obtained via alternative basis sparsification). We also present detailed results of our technique on Smirnov's family of algorithms.
{"title":"Towards automated generation of fast and accurate algorithms for recursive matrix multiplication","authors":"Jean-Guillaume Dumas , Clément Pernet , Alexandre Sedoglavic","doi":"10.1016/j.jsc.2025.102524","DOIUrl":"10.1016/j.jsc.2025.102524","url":null,"abstract":"<div><div>We propose a strategy for the generation of fast and accurate versions of non-commutative recursive matrix multiplication algorithms. To generate these algorithms, we consider matrix and tensor norm bounds governing the stability and accuracy of numerical matrix multiplication. We start by a unification on known max-norm bounds on matrix multiplication stability and then extend them to further norms and more generally to recursive bilinear algorithms and the alternative basis matrix multiplication algorithms. Then our strategy has three phases. First, we reduce those bounds by minimizing a growth factor along the orbits of the associated matrix multiplication tensor decomposition. Second, we develop heuristics that minimize the number of operations required to realize a bilinear formula, while further improving its accuracy. Third, we perform an alternative basis sparsification that improves on the time complexity constant and mostly preserves the overall accuracy. Let <span><math><mo>〈</mo><mi>m</mi><mo>×</mo><mi>k</mi><mo>×</mo><mi>n</mi><mo>:</mo><mi>r</mi><mo>〉</mo></math></span> denote a non-commutative algorithm multiplying an <span><math><mi>m</mi><mo>×</mo><mi>k</mi></math></span> matrix with an <span><math><mi>k</mi><mo>×</mo><mi>n</mi></math></span> matrix, using <em>r</em> coefficient products. Our strategy allows us to propose a <span><math><mo>〈</mo><mn>2</mn><mo>×</mo><mn>2</mn><mo>×</mo><mn>2</mn><mo>:</mo><mn>7</mn><mo>〉</mo></math></span> algorithm that reaches simultaneously a better accuracy in practice, compared to previously known such fast <span><math><mo>〈</mo><mn>2</mn><mo>×</mo><mn>2</mn><mo>×</mo><mn>2</mn><mo>:</mo><mn>7</mn><mo>〉</mo></math></span> Strassen-like algorithms, and a time complexity bound with the best currently known leading term (obtained via alternative basis sparsification). We also present detailed results of our technique on Smirnov's <span><math><mo>〈</mo><mn>3</mn><mo>×</mo><mn>3</mn><mo>×</mo><mn>6</mn><mo>:</mo><mn>40</mn><mo>〉</mo></math></span> family of algorithms.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"134 ","pages":"Article 102524"},"PeriodicalIF":1.1,"publicationDate":"2025-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145362884","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-14DOI: 10.1016/j.jsc.2025.102521
Alexander Levin
In this paper we introduce new types of characteristic sets in rings of inversive difference polynomials over an inversive difference field K and in inversive difference K-modules. These characteristic sets are defined via reductions associated with partitions of the basic set of automorphisms σ and with a generalization of the concept of the effective order of an ordinary difference polynomial. Using these characteristic sets, we provide an effective construction of multivariate difference dimension polynomials of a new type. They describe the transcendence degrees of intermediate fields of inversive difference field extensions and dimensions of submodules of inversive difference modules that are obtained by adjoining transforms of the generators whose orders with respect to the components of the partition of σ are bounded by two sequences of natural numbers. We show that such dimension polynomials carry substantially more invariants (that is, characteristics that do not depend on the set of generators) than standard (univariate) difference dimension polynomials. We also show how the obtained results can be applied to the equivalence problem for systems of algebraic difference equations and to the isomorphism problem for inversive difference modules.
{"title":"New dimension polynomials of inversive difference field extensions and inversive difference modules","authors":"Alexander Levin","doi":"10.1016/j.jsc.2025.102521","DOIUrl":"10.1016/j.jsc.2025.102521","url":null,"abstract":"<div><div>In this paper we introduce new types of characteristic sets in rings of inversive difference polynomials over an inversive difference field <em>K</em> and in inversive difference <em>K</em>-modules. These characteristic sets are defined via reductions associated with partitions of the basic set of automorphisms <em>σ</em> and with a generalization of the concept of the effective order of an ordinary difference polynomial. Using these characteristic sets, we provide an effective construction of multivariate difference dimension polynomials of a new type. They describe the transcendence degrees of intermediate fields of inversive difference field extensions and dimensions of submodules of inversive difference modules that are obtained by adjoining transforms of the generators whose orders with respect to the components of the partition of <em>σ</em> are bounded by two sequences of natural numbers. We show that such dimension polynomials carry substantially more invariants (that is, characteristics that do not depend on the set of generators) than standard (univariate) difference dimension polynomials. We also show how the obtained results can be applied to the equivalence problem for systems of algebraic difference equations and to the isomorphism problem for inversive difference modules.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"134 ","pages":"Article 102521"},"PeriodicalIF":1.1,"publicationDate":"2025-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145321163","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-14DOI: 10.1016/j.jsc.2025.102522
Nazia Jabeen, Junaid Alam Khan
The theory of Standard bases of modules over the localization of a polynomial ring (with respect to a local ordering >) is well-developed. For a polynomial subalgebra A, the concept of “Subalgebra bases” of (localization of A with respect to >) is well-defined, known as Sasbi bases (subalgebra analogue of Standard bases for ideals). In this paper, we have developed and proposed a theory for Standard bases of modules over (admitting a finite Sasbi basis), termed as Sasbi-Standard bases. The main goal is to develop a constructive algorithm for the computation of these bases, which we implement in the computer algebra system SINGULAR. We also present an application of these bases to compute singularity invariants: the monoid of values of an algebroid curve and the monomodule of values of the Kähler's differentials.
{"title":"Sasbi-standard bases of modules","authors":"Nazia Jabeen, Junaid Alam Khan","doi":"10.1016/j.jsc.2025.102522","DOIUrl":"10.1016/j.jsc.2025.102522","url":null,"abstract":"<div><div>The theory of Standard bases of modules over the localization of a polynomial ring (with respect to a local ordering >) is well-developed. For a polynomial subalgebra <em>A</em>, the concept of “Subalgebra bases” of <span><math><msub><mrow><mi>A</mi></mrow><mrow><mo>></mo></mrow></msub></math></span> (localization of <em>A</em> with respect to >) is well-defined, known as Sasbi bases (subalgebra analogue of Standard bases for ideals). In this paper, we have developed and proposed a theory for Standard bases of modules over <span><math><msub><mrow><mi>A</mi></mrow><mrow><mo>></mo></mrow></msub></math></span> (admitting a finite Sasbi basis), termed as Sasbi-Standard bases. The main goal is to develop a constructive algorithm for the computation of these bases, which we implement in the computer algebra system SINGULAR. We also present an application of these bases to compute singularity invariants: the monoid of values of an algebroid curve and the monomodule of values of the Kähler's differentials.</div></div>","PeriodicalId":50031,"journal":{"name":"Journal of Symbolic Computation","volume":"134 ","pages":"Article 102522"},"PeriodicalIF":1.1,"publicationDate":"2025-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145362883","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}