Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems

Domain decomposition methods are, alongside multigrid methods, one of the dominant paradigms in contemporary large-scale partial differential equation simulation. In this paper, a lightweight implementation of a theoretically and numerically scalable preconditioner is presented in the context of ove...

Full description

Bibliographic Details
Main Authors: Pierre Jolivet, Frédéric Hecht, Frédéric Nataf, Christophe Prud'homme
Format: Article
Language:English
Published: Hindawi Limited 2014-01-01
Series:Scientific Programming
Online Access:http://dx.doi.org/10.3233/SPR-140381
Description
Summary:Domain decomposition methods are, alongside multigrid methods, one of the dominant paradigms in contemporary large-scale partial differential equation simulation. In this paper, a lightweight implementation of a theoretically and numerically scalable preconditioner is presented in the context of overlapping methods. The performance of this work is assessed by numerical simulations executed on thousands of cores, for solving various highly heterogeneous elliptic problems in both 2D and 3D with billions of degrees of freedom. Such problems arise in computational science and engineering, in solid and fluid mechanics. While focusing on overlapping domain decomposition methods might seem too restrictive, it will be shown how this work can be applied to a variety of other methods, such as non-overlapping methods and abstract deflation based preconditioners. It is also presented how multilevel preconditioners can be used to avoid communication during an iterative process such as a Krylov method.
ISSN:1058-9244
1875-919X