
Summary: Computations with Toeplitz and Toeplitz-like matrices are fundamental for many areas of algebraic and numerical computing. The list of computational problems reducible to Toeplitz and Toeplitz-like computations includes, in particular, the evaluation of the greatest common divisor (gcd), the least common multiple (lcm), and the resultant of two polynomials, computing Padé approximation and the Berlekamp-Massey recurrence coefficients, as well as numerous problems reducible to these. Transition to Toeplitz and Toeplitz-like computations is currently the basis for the design of the parallel randomized NC (RNC) algorithms for these computational problems. Our main result is in constructing nearly optimal randomized parallel algorithms for Toeplitz and Toeplitz-like computations and, consequently, for numerous related computational problems (including the computational problems listed above), where all the input values are integers and all the output values are computed exactly. This includes randomized parallel algorithms for computing the rank, the determinant, and a basis for the null-space of an \(n {\times} n\) Toeplitz or Toeplitz-like matrix \(A\) filled with integers, as well as a solution \(x\) to a linear system \(Ax=f\) if the system is consistent. Our algorithms use \(O((\log n) \log (n \log |A|))\) parallel time and \(O(n \log n)\) processors, each capable of performing (in unit time) an arithmetic operation, a comparision, or a rounding of a rational number to a closest integer. The cost bounds cover the cost of the verification of the correctness of the output. The computations by these algorithms can be performed with the precision of \(O(n \log |A|)\) bits, which matches the precision required in order to represent the output, except for the rank computation, where the precision of the computation decreases. The algorithms involve either a single random parameter or at most \(2n-1\) parameters. The cited processor bounds are less by roughly factor \(n\) than ones supported by the known algorithms that run in polylogarithmic arithmetic time and do not use rounding to the closest integers. Technically, we first devise new algorithms supporting our old nearly optimal complexity estimates for parallel computations with general matrices filled with integers. Then we decrease dramatically, by roughly factor \(n^{1.376}\), the processor bounds required in these algorithms in the case where the input matrix is Toeplitz-like. Our algorithms exploit and combine some new techniques (which may be of independent interest, e.g., in the study of parallel and sequential computation of recursive factorization of integer matrices) as well as our earlier techniques of variable diagonal (relating to each other several known algebraic and numerical methods), stream contraction, and the truncation of displacement generators in Toeplitz-like computations; our development and application of these techniques may be of independent interest.
polynomial gcd, computational complexity, Analysis of algorithms and problem complexity, randomized algorithms, displacement rank, Other matrix algorithms, parallel algorithms, block Gauss-Jordan decomposition, Complexity and performance of numerical algorithms, Grammars and rewriting systems, Toeplitz operators, Hankel operators, Wiener-Hopf operators, \(p\)-adic lifting, Parallel algorithms in computer science, Toeplitz-like matrices, Newton-Hensel's lifting, Toeplitz matrix computations
polynomial gcd, computational complexity, Analysis of algorithms and problem complexity, randomized algorithms, displacement rank, Other matrix algorithms, parallel algorithms, block Gauss-Jordan decomposition, Complexity and performance of numerical algorithms, Grammars and rewriting systems, Toeplitz operators, Hankel operators, Wiener-Hopf operators, \(p\)-adic lifting, Parallel algorithms in computer science, Toeplitz-like matrices, Newton-Hensel's lifting, Toeplitz matrix computations
| selected citations These citations are derived from selected sources. This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | 3 | |
| popularity This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network. | Average | |
| influence This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | Average | |
| impulse This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network. | Average |
