Two-level preconditioning for BEM with GenEO
Abstract
Domain Decomposition Methods (DDM), such as Additive Schwarz (AS), can be used to precondition linear systems arising from Boundary Integral Equations (BIE). This approach was widely studied since then and extended in various directions. The basic idea is to adapt the classical FEM-based AS to the BIE context: this includes a two-level preconditioner relying on a coarse space, which leads to theoretical bounds on the condition number. Regarding the choice of relevant coarse spaces, important progress has been achieved in recent years for the FEM context. For the construction of coarse spaces, the Generalized Eigenproblems in the Overlaps (GenEO) has emerged as one of the most promising approach for symmetric positive definite problems. Instead of solving a coarse problem on a coarse mesh, GenEO takes eigenvectors of well chosen local eigenproblems as a basis for the coarse space. As one of its interesting features, GenEO is only based on the knowledge of the stiffness matrix elements and discretization agnostic, left apart a few reasonable assumptions. In this talk, we will present recent theoretical and numerical results in 2D and 3D aiming at adapting GenEO to the BIE context for symmetric positive definite problems on closed and open surface. Examples of applications are Laplace problems on screens or dissipative Helmholtz problems. This work is supported by the project NonlocalDD, research grant ANR-15-CE23-0017-01 from the French National Research Agency and the numerical results are obtained using HPC resources from GENCI- CINES (Grant 2017-A0020607330).
