An exact solver for simple ${\mathcal H}$-matrix systems
Numerical Analysis
2014-02-24 v1
Abstract
Hierarchical matrices (usually abbreviated -matrices) are frequently used to construct preconditioners for systems of linear equations. Since it is possible to compute approximate inverses or factorizations in -matrix representation using only 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 -matrix format \emph{exactly} using operations. The central idea of our approach is to avoid computing the inverse and instead use an efficient representation of the factorization based on low-rank updates performed with the well-known Sherman-Morrison-Woodbury equation.
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}
}