Abstract:A sparsifiable spectrally equivalent preconditioner is developed for a special finite element discretization of two-dimensional elliptic problems. We start with a spectrally equivalent preconditioner which is closely related to the discrete Dirichlet-to-Neumann (DtN) mappings in sub-domains. The fullness of discrete DtN mappings renders an essential difficulty for a direct implementation. Based on numerical evidences, we present a particular square root matrix that is spectrally equivalent to the discrete DtN mappings. Unfortunately, this square root matrix is still not sparsifiable. To remedy this problem, we resort to the Chebyshev rational approximation. Incorporated with the fast sine transforms, a sparsifiable and spectrally equivalent preconditioner is finally obtained. Numerical examples are provided to validate efficiency of the proposed preconditioning technique.

Sparsifiable Spectral Equivalence of DtN Mapping and Its Application to Elliptic Problems