A massively parallel nonoverlapping additive Schwarz method for discontinuous Galerkin discretization of elliptic problems

A second order elliptic problem with discontinuous coefficient in 2-D or 3-D is considered. The problem is discretized by a symmetric weighted interior penalty discontinuous Galerkin finite element method with nonmatching simplicial elements and piecewise linear functions. The resulting discrete problem is solved by a two-level additive Schwarz method with a relatively coarse grid and with local solves restricted to subdomains which can be as small as single element. In this way the method has a potential for a very high level of fine grained parallelism. Condition number estimate depending on the relative sizes of the underlying grids is provided. The rate of convergence of the method is independent of the jumps of the coefficient if its variation is moderate inside coarse grid substructures or on local solvers’ subdomain boundaries. Numerical experiments are reported which confirm theoretical results.