Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 1, Page 368-396, March 2024. As a crucial first step towards finding the (approximate) common roots of a (possibly overdetermined) bivariate polynomial system of equations, the problem of determining an explicit numerical basis for the right null space of the system's Macaulay matrix is considered. If [math] denotes the total degree of the bivariate polynomials of the system, the cost of computing a null space basis containing all system roots is [math] floating point operations through standard numerical algebra techniques (e.g., a singular value decomposition, rank-revealing QR-decomposition). We show that it is actually possible to design an algorithm that reduces the complexity to [math]. The proposed algorithm exploits the Toeplitz structures of the Macaulay matrix under a nongraded lexicographic ordering of its entries and uses the low displacement rank properties to efficiently convert it into a Cauchy-like matrix with the help of fast Fourier transforms. By modifying the classical Schur algorithm with total pivoting for Cauchy-like matrices, a compact representation of the right null space is eventually obtained from a rank-revealing LU-factorization. Details of the proposed method, including numerical experiments, are fully provided for the case wherein the polynomials are expressed in the monomial basis. Furthermore, it is shown that an analogous fast algorithm can also be formulated for polynomial systems expressed in the Chebyshev basis.

Jacobi-free and Complex-Free Method for Finding Simultaneously All Zeros of Polynomials Having Only Real Zeros

A Parallel Algorithm For Finding All The Roots Of A Polynomial

New simultaneous methods for finding all zeros of a polynomial

An Iteration Method with Generally Convergent Property for Cubic Polynomials.

Recursive Algorithm Without Extra Function Evaluations for the Jacobian Matrix of Viéta's Polynomial System with Applications

Parallel computation of real solving bivariate polynomial systems by zero-matching method

Location of the Zeros of Certain Complex-Valued Harmonic Polynomials

Simple bounds for the extreme zeroes of Jacobi polynomials

Aberth's Method for the Parallel Iterative Finding of Polynomial Zeros

Zeros of univariate interval polynomials

Real Root Isolation of Polynomial Equations Based on Hybrid Computation

A New Method for Real Root Isolation of Univariate Polynomials

Polynomials with real zeros and Polya frequency sequences

Real root refinements for univariate polynomial equations

Zeros Transfer For Recursively defined Polynomials

A solution of polynomial equations

A Fast Algorithm for Computing Macaulay Null Spaces of Bivariate Polynomial Systems

The number of real zeros of polynomials with constrained coefficients

Zhuang-Zi: A New Algorithm for Solving Multivariate Polynomial Equations over a Finite Field.

Quantifier elimination for quartics

On computing zeros of a bivariate Bernstein polynomial