Improving Multifrontal Methods by Means of Block Low-Rank Representations - Archive ouverte HAL Access content directly
Journal Articles SIAM Journal on Scientific Computing Year : 2015

Improving Multifrontal Methods by Means of Block Low-Rank Representations

(1, 2) , (3) , (4) , (1, 5) , (6) , (1, 2)
1
2
3
4
5
6

Abstract

Matrices coming from elliptic partial differential equations have been shown to have alow-rank property: well-defined off-diagonal blocks of their Schur complements can be approximatedby low-rank products. Given a suitable ordering of the matrix which gives the blocks a geometricalmeaning, such approximations can be computed using an SVD or a rank-revealing QR factorization.The resulting representation offers a substantial reduction of the memory requirement and givesefficient ways to perform many of the basic dense linear algebra operations. Several strategies,mostly based on hierarchical formats, have been proposed to exploit this property. We study asimple, nonhierarchical, low-rank format called block low-rank (BLR) and explain how it can be usedto reduce the memory footprint and the complexity of sparse direct solvers based on the multifrontalmethod. We present experimental results on matrices coming from elliptic PDEs and from variousother applications. We show that even if BLR-based factorizations are asymptotically less efficientthan hierarchical approaches, they still deliver considerable gains. The BLR format is compatiblewith numerical pivoting, and its simplicity and flexibility make it easy to use in the context of ageneral purpose, algebraic solver.

Dates and versions

hal-01237169 , version 1 (02-12-2015)

Identifiers

Cite

Patrick Amestoy, Cleve Ashcraft, Olivier Boiteau, Alfredo Buttari, Jean-Yves L'Excellent, et al.. Improving Multifrontal Methods by Means of Block Low-Rank Representations. SIAM Journal on Scientific Computing, 2015, 37 (3), pp.A1451-A1474. ⟨10.1137/120903476⟩. ⟨hal-01237169⟩
180 View
0 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More