Section: New Results

Domain decomposition preconditioning with approximate coarse solve

Convergence of domain decomposition methods relies heavily on the efficiency of the coarse space used in the second level. The GenEO coarse space has been shown to lead to a fully robust two-level Schwarz preconditioner which scales well over multiple cores  [9], [2] as has been proved rigorously in [9]. 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++ [7], Feel++  [31] and recently in Dune  [30] and is implemented as a standalone library in HPDDM [8]. 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 approximate coarse solves. In [28], robustness of GenEO methods is analyzed with respect to approximate coarse solves. Interestingly, the GenEO-2 method introduced in [6] has to be modified in order to be able to prove its robustness in this context.