Skip to Main content Skip to Navigation
Journal articles

Bridging the gap between flat and hierarchical low-rank matrix formats: the multilevel BLR format

Patrick Amestoy 1 Alfredo Buttari 1 Jean-Yves l'Excellent 2 Theo Mary 3
1 IRIT-APO - Algorithmes Parallèles et Optimisation
IRIT - Institut de recherche en informatique de Toulouse
2 ROMA - Optimisation des ressources : modèles, algorithmes et ordonnancement
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : Matrices possessing a low-rank property arise in numerous scientific applications. This property can be exploited to provide a substantial reduction of the complexity of their LU or LDL T factorization. Among the possible low-rank formats, the flat Block Low-Rank (BLR) format is easy to use but achieves superlinear complexity. Alternatively, the hierarchical formats achieve linear complexity at the price of a much more complex, hierarchical matrix representation. In this paper, we propose a new format based on multilevel BLR approximations: the matrix is recursively defined as a BLR matrix whose full-rank blocks are themselves represented by BLR matrices. We call this format multilevel BLR (MBLR). Contrarily to hierarchical matrices, the number of levels in the block hierarchy is fixed to a given constant; while this format can still be represented within the H formalism, we show that applying the H theory to it leads to very pessimistic complexity bounds. We therefore extend the theory to prove better bounds, and show that the MBLR format provides a simple way to finely control the desired complexity of dense factorizations. By striking a balance between the simplicity of the BLR format and the low complexity of the hierarchical ones, the MBLR format bridges the gap between flat and hierarchical low-rank matrix formats. The MBLR format is of particular relevance in the context of sparse direct solvers, for which it is able to trade off the optimal dense complexity of the hierarchical formats to benefit from the simplicity and flexibility of the BLR format while still achieving O(n) sparse complexity. We finally compare our MBLR format with the related BLR-H (or Lattice-H) format; our theoretical analysis shows that both formats achieve the same asymptotic complexity for a given top level block size.
Document type :
Journal articles
Complete list of metadatas

Cited literature [30 references]  Display  Hide  Download
Contributor : Theo Mary <>
Submitted on : Thursday, October 10, 2019 - 4:53:12 PM
Last modification on : Tuesday, September 8, 2020 - 10:48:04 AM


Files produced by the author(s)



Patrick Amestoy, Alfredo Buttari, Jean-Yves l'Excellent, Theo Mary. Bridging the gap between flat and hierarchical low-rank matrix formats: the multilevel BLR format. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2019, 41 (3), pp.A1414-A1442. ⟨10.1137/18M1182760⟩. ⟨hal-01774642v3⟩



Record views


Files downloads