English

On matrices with displacement structure: generalized operators and faster algorithms

Symbolic Computation 2017-03-13 v1 Computational Complexity

Abstract

For matrices with displacement structure, basic operations like multiplication, inversion, and linear system solving can all be expressed in terms of the following task: evaluate the product AB\mathsf{A}\mathsf{B}, where A\mathsf{A} is a structured n×nn \times n matrix of displacement rank α\alpha, and B\mathsf{B} is an arbitrary n×αn\times\alpha matrix. Given B\mathsf{B} and a so-called "generator" of A\mathsf{A}, this product is classically computed with a cost ranging from O(α2M(n))O(\alpha^2 \mathscr{M}(n)) to O(α2M(n)log(n))O(\alpha^2 \mathscr{M}(n)\log(n)) arithmetic operations, depending on the type of structure of A\mathsf{A}; here, M\mathscr{M} is a cost function for polynomial multiplication. In this paper, we first generalize classical displacement operators, based on block diagonal matrices with companion diagonal blocks, and then design fast algorithms to perform the task above for this extended class of structured matrices. The cost of these algorithms ranges from O(αω1M(n))O(\alpha^{\omega-1} \mathscr{M}(n)) to O(αω1M(n)log(n))O(\alpha^{\omega-1} \mathscr{M}(n)\log(n)), with ω\omega such that two n×nn \times n matrices over a field can be multiplied using O(nω)O(n^\omega) field operations. By combining this result with classical randomized regularization techniques, we obtain faster Las Vegas algorithms for structured inversion and linear system solving.

Keywords

Cite

@article{arxiv.1703.03734,
  title  = {On matrices with displacement structure: generalized operators and faster algorithms},
  author = {Alin Bostan and Claude-Pierre Jeannerod and Christophe Mouilleron and Éric Schost},
  journal= {arXiv preprint arXiv:1703.03734},
  year   = {2017}
}

Comments

To appear in SIAM Journal on Matrix Analysis and Applications

R2 v1 2026-06-22T18:42:27.221Z