Skip to Main content Skip to Navigation
New interface
Journal articles

On the Complexity of the Block Low-Rank Multifrontal Factorization

Abstract : Matrices coming from elliptic Partial Differential Equations have been shown to have a low-rank property: well defined off-diagonal blocks of their Schur complements can be approximated by low-rank products and this property can be efficiently exploited in multifrontal solvers to provide a substantial reduction of their complexity. Among the possible low-rank formats, the Block Low-Rank format (BLR) is easy to use in a general purpose multifrontal solver and has been shown to provide significant gains compared to full-rank on practical applications. However, unlike hierarchical formats, such as H and HSS, its theoretical complexity was unknown. In this paper, we extend the theoretical work done on hierarchical matrices in order to compute the theoretical complexity of the BLR multifrontal factorization. We then present several variants of the BLR multifrontal factorization, depending on the strategies used to perform the updates in the frontal matrices and on the constraints on how numerical pivoting is handled. We show how these variants can further reduce the complexity of the factorization. In the best case (3D, constant ranks), we obtain a complexity of the order of O(n^4/3). We provide an experimental study with numerical results to support our complexity bounds.
Document type :
Journal articles
Complete list of metadata

Cited literature [41 references]  Display  Hide  Download
Contributor : Theo Mary Connect in order to contact the contributor
Submitted on : Thursday, October 10, 2019 - 4:49:23 PM
Last modification on : Tuesday, October 25, 2022 - 4:24:34 PM


Files produced by the author(s)



Patrick Amestoy, Alfredo Buttari, Jean-Yves L'Excellent, Théo Mary. On the Complexity of the Block Low-Rank Multifrontal Factorization. SIAM Journal on Scientific Computing, 2017, 39 (4), pp.34. ⟨10.1137/16M1077192⟩. ⟨hal-01322230v3⟩



Record views


Files downloads