English

An exact solver for simple ${\mathcal H}$-matrix systems

Numerical Analysis 2014-02-24 v1

Abstract

Hierarchical matrices (usually abbreviated H{\mathcal H}-matrices) are frequently used to construct preconditioners for systems of linear equations. Since it is possible to compute approximate inverses or LULU factorizations in H{\mathcal H}-matrix representation using only O(nlog2n){\mathcal O}(n \log^2 n) operations, these preconditioners can be very efficient. Here we consider an algorithm that allows us to solve a linear system of equations given in a simple H{\mathcal H}-matrix format \emph{exactly} using O(nlog2n){\mathcal O}(n \log^2 n) operations. The central idea of our approach is to avoid computing the inverse and instead use an efficient representation of the LULU factorization based on low-rank updates performed with the well-known Sherman-Morrison-Woodbury equation.

Keywords

Cite

@article{arxiv.1402.5398,
  title  = {An exact solver for simple ${\mathcal H}$-matrix systems},
  author = {Steffen Börm and Jessica Gördes},
  journal= {arXiv preprint arXiv:1402.5398},
  year   = {2014}
}
R2 v1 2026-06-22T03:13:23.634Z