Fast Solution of 3D Elastodynamic Boundary Element Problems by Hierarchical Matrices

In this paper a fast solver for three-dimensional elastodynamic BEM problems formulated in the Laplace transform domain is presented, implemented and tested. The technique is based on the use of hierarchical matrices for the representation of the collocation matrix for each value of the Laplace parameter of interest and uses a preconditioned GMRES for the solution of the algebraic system of equations. The preconditioner is built exploiting the hierarchical arithmetic and taking full advantage of the hierarchical format. An original strategy for speeding up the overall analysis is presented and tested. The reported numerical results demonstrate the effectiveness of the technique.