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
id doaj-4b218246e8a647f19fb3d062666de33e
record_format Article
spelling doaj-4b218246e8a647f19fb3d062666de33e2021-07-02T03:23:12ZengHindawi LimitedScientific Programming1058-92441875-919X2014-01-0122215717110.3233/SPR-140381Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic ProblemsPierre Jolivet0Frédéric Hecht1Frédéric Nataf2Christophe Prud'homme3Laboratoire J. Kuntzmannn, Université J. Fourier, Grenoble Cedex 9, FranceLaboratoire J.-L. Lions, Université P. et M. Curie, Paris, FranceLaboratoire J.-L. Lions, Université P. et M. Curie, Paris, FranceIRMA, Université de Strasbourg, Strasbourg Cedex, FranceDomain 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.http://dx.doi.org/10.3233/SPR-140381
collection DOAJ
language English
format Article
sources DOAJ
author Pierre Jolivet
Frédéric Hecht
Frédéric Nataf
Christophe Prud'homme
spellingShingle Pierre Jolivet
Frédéric Hecht
Frédéric Nataf
Christophe Prud'homme
Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems
Scientific Programming
author_facet Pierre Jolivet
Frédéric Hecht
Frédéric Nataf
Christophe Prud'homme
author_sort Pierre Jolivet
title Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems
title_short Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems
title_full Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems
title_fullStr Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems
title_full_unstemmed Scalable Domain Decomposition Preconditioners for Heterogeneous Elliptic Problems
title_sort scalable domain decomposition preconditioners for heterogeneous elliptic problems
publisher Hindawi Limited
series Scientific Programming
issn 1058-9244
1875-919X
publishDate 2014-01-01
description 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.
url http://dx.doi.org/10.3233/SPR-140381
work_keys_str_mv AT pierrejolivet scalabledomaindecompositionpreconditionersforheterogeneousellipticproblems
AT frederichecht scalabledomaindecompositionpreconditionersforheterogeneousellipticproblems
AT fredericnataf scalabledomaindecompositionpreconditionersforheterogeneousellipticproblems
AT christopheprudhomme scalabledomaindecompositionpreconditionersforheterogeneousellipticproblems
_version_ 1721341701581176832