Adaptive H-matrix computations in linear elasticity

This article deals with the efficient numerical treatment of the Lamé equations. The equations of linear elasticity are considered as boundary integral equations and solved in the setting of the boundary element method (BEM). Using BEM, one is faced with the solution of a system of equations with a fully populated system matrix, which is in general very costly. In order to overcome this difficulty, adaptive and approximate algorithms based on hierarchical matrices and the adaptive cross approximation are proposed. These new methods rely on error estimators and refinement techniques known from adaptivity but are not used here to improve the mesh. We apply these new techniques to both, the efficient solution of Lamé equations and to the multiplication with given data.