Skip to Main content Skip to Navigation
Journal articles

Time and space efficient generators for quasiseparable matrices

Clement Pernet 1, 2 Arne Storjohann 3
1 ARIC - Arithmetic and Computing
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : The class of quasiseparable matrices is defined by the property that any submatrix entirely below or above the main diagonal has small rank, namely below a bound called the order of quasiseparability. These matrices arise naturally in solving PDE's for particle interaction with the Fast Multi-pole Method (FMM), or computing generalized eigenvalues. From these application fields, structured representations and algorithms have been designed in numerical linear algebra to compute with these matrices in time linear in the matrix dimension and either quadratic or cubic in the quasiseparability order. Motivated by the design of the general purpose exact linear algebra library LinBox, and by algorithmic applications in algebraic computing, we adapt existing techniques introduce novel ones to use quasiseparable matrices in exact linear algebra, where sub-cubic matrix arithmetic is available. In particular, we will show, the connection between the notion of quasiseparability and the rank profile matrix invariant, that we have introduced in 2015. It results in two new structured representations, one being a simpler variation on the hierarchically semiseparable storage, and the second one exploiting the generalized Bruhat decomposition. As a consequence, most basic operations, such as computing the quasiseparability orders, applying a vector, a block vector, multiplying two quasiseparable matrices together, inverting a quasiseparable matrix, can be at least as fast and often faster than previous existing algorithms.
Document type :
Journal articles
Complete list of metadata

Cited literature [33 references]  Display  Hide  Download
Contributor : Clément Pernet <>
Submitted on : Sunday, October 20, 2019 - 4:36:48 PM
Last modification on : Friday, June 25, 2021 - 3:40:05 PM
Long-term archiving on: : Tuesday, January 21, 2020 - 1:26:57 PM




Clement Pernet, Arne Storjohann. Time and space efficient generators for quasiseparable matrices. Journal of Symbolic Computation, Elsevier, 2018, Special issue on the 41th International Symposium on Symbolic and Alge-braic Computation (ISSAC’16), 85, pp.224-246. ⟨10.1016/j.jsc.2017.07.010⟩. ⟨hal-01424252v2⟩



Record views


Files downloads