Overlapping Domain Decomposition Preconditioner for Integral Equations
Abstract
The discretization of certain integral equations, e.g., the first-kind Fredholm equation of Laplace's equation, leads to symmetric positive-definite linear systems, where the coefficient matrix is dense and often ill-conditioned. We introduce a new preconditioner based on a novel overlapping domain decomposition that can be combined efficiently with fast direct solvers. Empirically, we observe that the condition number of the preconditioned system is , independent of the problem size. Our domain decomposition is designed so that we can construct approximate factorizations of the subproblems efficiently. In particular, we apply the recursive skeletonization algorithm to subproblems associated with every subdomain. We present numerical results on problem sizes up to in 2D and in 3D, which were solved in less than 16 hours and three hours, respectively, on an Intel Xeon Platinum 8280M.
Cite
@article{arxiv.2108.12574,
title = {Overlapping Domain Decomposition Preconditioner for Integral Equations},
author = {Chao Chen and George Biros},
journal= {arXiv preprint arXiv:2108.12574},
year = {2022}
}