Mathematical analysis of robustness of two-level domain decomposition methods with respect to inexact coarse solves
- PDF / 396,763 Bytes
- 23 Pages / 439.37 x 666.142 pts Page_size
- 88 Downloads / 211 Views
Numerische Mathematik
Mathematical analysis of robustness of two-level domain decomposition methods with respect to inexact coarse solves Frédéric Nataf1 Received: 13 November 2017 / Revised: 9 July 2019 © Springer-Verlag GmbH Germany, part of Springer Nature 2020
Abstract Convergence of domain decomposition methods rely heavily on the efficiency of the coarse space used in the second level. The GenEO coarse space has been shown to lead to a robust two-level Schwarz preconditioner which scales well over multiple cores (Spillane et al. in Numer Math 126(4):741–770, 2014. https://doi.org/10.1007/ s00211-013-0576-y; Dolean et al. in An introduction to domain decomposition methods: algorithms, theory and parallel implementation, SIAM, Philadelphia, 2015). The robustness is due to its good approximation properties for problems with highly heterogeneous material parameters. It is available in the finite element packages FreeFem++ (Hecht in J Numer Math 20(3–4):251–265, 2012), Feel++ (Prud’homme in Sci Program 14(2):81–110, 2006), Dune (Blatt et al. in Arch Numer Softw 4(100):13–29, 2016) and is implemented as a standalone library in HPDDM (Jolivet and Nataf in HPDDM: high-Performance Unified framework for Domain Decomposition methods, MPI-C++ library, 2014. https://github.com/hpddm/hpddm) and as such is available as well as a PETSc (Balay et al. in: Arge, Bruaset, Langtangen, (eds) Modern software tools in scientific computing, Birkhäuser Press, Basel, 1997) preconditioner. But the coarse component of the preconditioner can ultimately become a bottleneck if the number of subdomains is very large and exact solves are used. It is therefore interesting to consider the effect of inexact coarse solves. In this paper, robustness of GenEO methods is analyzed with respect to inexact coarse solves. Interestingly, the GenEO2 method introduced in Haferssas et al. (SIAM J Sci Comput 39(4):A1345–A1365, 2017. https://doi.org/10.1137/16M1060066) has to be modified in order to be able to prove its robustness in this context. Mathematics Subject Classification 65Y05 · 65F08
B 1
Frédéric Nataf [email protected] Laboratoire J.L. Lions, UPMC, CNRS UMR7598, Equipe LJLL-INRIA Alpines, 4 Place Jussieu, 75005 Paris, France
123
F. Nataf
1 Introduction Convergence of domain decomposition methods rely heavily on the efficiency of the coarse space used in the second level, see [16,17,23] and references therein. These methods are based on two ingredients: a coarse space (CS) and a correction formula (see e.g. [22]). The GenEO coarse space introduced in [20] has been shown to lead to a robust two-level Schwarz preconditioner which scales well over multiple cores. The robustness is due to its good approximation properties for problems with highly heterogeneous material parameters. This approach is closely related to [5]. We refer to the introduction of [20] for more details on the differences and similarities between both approaches. Here we will mainly work with a slight modification of the GenEO CS introduced in [3] for t
Data Loading...
