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.