Direct Schur complement method by domain decomposition based on H-matrix approximation

摘要

The goal of this paper is the construction of a data-sparse approximation to the Schur complement on the interface corresponding to FEM and BEM approximations of an elliptic equation by domain decomposition. Using the hierarchical (H-matrix) formats we elaborate the approximate Schur complement inverse in an explicit form. The required cost O (N Γlog q N Γ) is almost linear in N Γ – the number of degrees of freedom on the interface. As input, we require the Schur complement matrices corresponding to subdomains and represented in the H-matrix format. In the case of piecewise constant coefficients these matrices can be computed via the BEM representation with the cost O(N Γlog q N Γ), while in the general case the FEM discretisation leads to the complexity O(N Ωlog q N Ω), where N Ω is the number of degrees of freedom in the domain.