On the Robustness of Elliptic Resolvents Computed by Means of the Technique of Hierarchical Matrices

Hierarchical matrices provide a data-sparse way to approximate fully populated matrices. In this paper we exploit robust H-matrix techniques to approximate the resolvents of stiffness matrices as they appear in (finite element or finite difference) discretisations of elliptic partial differential equations