English

Overlapping Domain Decomposition Preconditioner for Integral Equations

Numerical Analysis 2022-08-22 v2 Numerical Analysis

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 O(1)O(1), 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 16384216\,384^2 in 2D and 2563256^3 in 3D, which were solved in less than 16 hours and three hours, respectively, on an Intel Xeon Platinum 8280M.

Keywords

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}
}
R2 v1 2026-06-24T05:29:20.201Z