期刊
SIAM JOURNAL ON SCIENTIFIC COMPUTING
卷 39, 期 4, 页码 A1710-A1740出版社
SIAM PUBLICATIONS
DOI: 10.1137/16M1077192
关键词
sparse linear algebra; multifrontal factorization; block low-rank
资金
- [ANR-11-LABX-0040-CIMI]
- [ANR-11-IDEX-0002-02]
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 (BLR) format 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.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据