---
title: 'FROSch: Fast & Robust Overlapping Schwarz'
url: https://www.emergentmind.com/topics/fast-and-robust-overlapping-schwarz-frosch
type: topic
---

# FROSch: Fast & Robust Overlapping Schwarz

Searching arXiv for FROSch and closely related overlapping Schwarz papers to ground the article in current literature.
Fast and Robust Overlapping Schwarz (FROSch) is a parallel domain decomposition framework for overlapping Schwarz preconditioning, described in the literature as a package inside Trilinos and, in particular, inside ShyLU, with support for one- and two-level methods, GDSW- and RGDSW-type coarse spaces, and later extensions to multilevel, monolithic saddle-point, GPU, ORAS, and nonlinear formulations [2304.04876, 2508.04356, 2410.24138, 2410.22871]. At its core, FROSch realizes additive overlapping Schwarz operators that combine local overlapping subdomain solves with a coarse correction built from interface-based basis functions extended by discrete harmonic or energy-minimizing solves. Closely related work on algebraic coarse spaces, high-contrast multiscale basis construction, vector-field Schwarz theory, and virtual-element coarse spaces has enlarged the conceptual scope of what counts as a FROSch-style method, even when the name itself is not used explicitly [2408.08187, 1901.00112, 2404.11986, 2512.07181].

## 1. Definition and position within domain decomposition

FROSch is presented as a framework or library for scalable overlapping Schwarz methods rather than as a single fixed preconditioner. In the GPU study, it is introduced as a package in Trilinos implementing GDSW-type two-level overlapping Schwarz preconditioners for CPU and GPU clusters [2304.04876]. In the monolithic fluid study, it is described as a framework inside ShyLU, itself part of Trilinos, providing additive overlapping Schwarz preconditioners, GDSW coarse spaces, multilevel support, and algebraic construction of local matrices, restriction and interpolation operators, coarse problems, and overlap information [2508.04356]. The deal.II coupling paper makes the same point from the operator side: FROSch is a Schwarz framework able to build one-level and two-level operators and to combine Schwarz operators additively or multiplicatively [2410.22871].

This framing places FROSch within the overlapping branch of domain decomposition rather than the substructuring branch represented by BDDC or FETI-DP. The common ingredients are an overlapping cover of the domain, local restriction and prolongation operators, local subdomain solves, and a low-dimensional coarse correction. In published FROSch work, the coarse level is usually based on extension-based GDSW or RGDSW spaces, but later papers broaden this to monolithic saddle-point coarse spaces, ORAS local operators, nonlinear coarse corrections, and multilevel recursion [2506.16179, 2508.04356, 2410.24138].

A recurring distinction in the literature is between one-level and two-level overlapping Schwarz. One-level methods retain only the sum of local overlapping subdomain solves and are repeatedly described as not scalable as the number of subdomains grows. Two-level methods add a coarse correction and thereby provide the global communication path missing from purely local overlap [2508.04356, 2410.22871, 2512.01706].

## 2. Algebraic form of the preconditioner

The standard two-level overlapping Schwarz operator used throughout the FROSch literature has the form
\[
M^{-1}=\Phi A_0^{-1}\Phi^T+\sum_{i=1}^N R_i^T A_i^{-1}R_i,
\]
or equivalently
\[
M^{-1}=R_0^T A_0^{-1}R_0+\sum_{i=1}^N R_i^T A_i^{-1}R_i,
\]
with \(A_i=R_i A R_i^T\) the local overlapping subdomain matrices and \(A_0=\Phi^T A\Phi\) the coarse operator [2408.08187, 2304.04876]. This form is used for scalar elliptic systems, elasticity, reduced-integration Stokes discretizations, and algebraic coarse-space comparisons [2304.04876, 2512.01706, 2408.08187].

A defining feature of GDSW-type FROSch is that the coarse basis is obtained from interface traces and then extended into subdomain interiors by an energy-minimizing or discrete harmonic solve. After ordering unknowns into interior and interface blocks,
\[
A = \begin{bmatrix} A_{II} & A_{I\Gamma}\\ A_{\Gamma I} & A_{\Gamma\Gamma} \end{bmatrix},
\]
the coarse basis satisfies
\[
\Phi = \begin{bmatrix} -A_{II}^{-1}A_{I\Gamma}\\ I \end{bmatrix}\Phi_\Gamma,
\]
so the coarse space is determined entirely by the choice of interface traces \(\Phi_\Gamma\) [2408.08187, 2304.04876, 2310.12659]. For scalar diffusion, the traces are restrictions of the null space of the global Neumann problem; for vector problems such as elasticity they are multiplied by translational, and in some settings rotational, null-space modes [2304.04876, 2502.19975].

The same algebraic template extends to coupled saddle-point systems. For the monolithic Stokes and Navier–Stokes operators,
\[
\mathcal A=\begin{bmatrix} A & B^T\\ B & 0 \end{bmatrix},
\]
the preconditioner becomes
\[
\mathcal{B}_{\mathrm GDSW}=\Phi \mathcal A_0^{-1}\Phi^T+\sum_{i=1}^N \mathcal R_i^T \mathcal A_i^{-1}\mathcal R_i,
\]
with block restrictions in velocity and pressure and a coarse basis matrix \(\Phi\) containing both components [2508.04356, 2506.16179]. This monolithic formulation is central to later fluid papers because the local and coarse problems preserve the full saddle-point coupling instead of approximating only a velocity block or Schur complement factor.

The ORAS extension modifies the first level rather than the coarse space. In that case the local matrices \(A_i\) are replaced by optimized local operators \(B_i\) assembled with Robin or impedance-type transmission conditions:
\[
M_{\rm ORAS}^{-1}=\sum_{i=1}^N P_i D_i B_i^{-1}R_i.
\]
The deal.II coupling paper emphasizes that this step is no longer purely algebraic, because the optimized local matrices must be assembled separately from geometric interface data [2410.22871].

## 3. Coarse-space families

The literature around FROSch is largely a literature on coarse-space design. Standard FROSch papers treat GDSW and RGDSW as the reference family, while later work introduces intermediate, monolithic, algebraic multiscale, and energy-minimizing multiscale variants [2304.04876, 2408.08187, 2506.16179, 1901.00112].

| Coarse-space family | Interface construction | Role reported in the literature |
|---|---|---|
| GDSW | Vertices, edges, and faces with null-space restrictions | Standard FROSch coarse space [2304.04876] |
| RGDSW | Reduced interface components or nodal equivalence classes with partition-of-unity weighting | Smaller coarse space, often weaker robustness [2408.08187, 2512.01706] |
| GDSW\(^\star\) | Vertices with adjacent edges grouped, faces kept separate | Intermediate space between GDSW and RGDSW [2506.16179] |
| AMS / algebraic MsFEM | Recursive edge/face traces, then harmonic extension | Algebraic multiscale alternative in the same two-level architecture [2408.08187] |
| Energy-minimizing multiscale basis | Local generalized eigenproblems plus constrained energy minimization | FROSch-related high-contrast coarse-space design [1901.00112] |

In classical GDSW, interface components are faces, edges, and vertices, and the interface basis functions form a partition of unity. RGDSW reduces coarse dimension by merging interface components according to nodal equivalence classes and using inverse-multiplicity-type scaling. The computational study of algebraic coarse spaces interprets GDSW as essentially the standard FROSch coarse space and RGDSW as a reduced-dimension FROSch-type alternative [2408.08187]. The GPU paper and the direct-solver comparison both use reduced GDSW spaces in practice, reflecting the persistent tradeoff between coarse-space richness and coarse-problem cost [2304.04876, 2310.12659].

A more specialized development is GDSW\(^\star\), introduced for monolithic Navier–Stokes preconditioning. It keeps faces separate, but groups each vertex together with its adjacent edges, yielding a space larger than RGDSW and smaller than full GDSW. The reported motivation is to preserve much of the iteration-count benefit of GDSW while reducing setup and coarse-solve cost. In the cited Navier–Stokes experiments, the combination “velocity: GDSW\(^\star\), pressure: RGDSW” is repeatedly identified as particularly effective [2506.16179].

A different line of development replaces classical null-space traces by coefficient-aware algebraic or multiscale traces. The algebraic coarse-space study places AMS, GDSW, and RGDSW inside the same two-level additive Schwarz architecture, differing only in the choice of interface traces \(\Phi_\Gamma\) before harmonic extension [2408.08187]. The multiscale coarse-basis paper goes further: local generalized eigenproblems identify low-energy modes, and the final basis functions are obtained by constrained energy minimization rather than by direct partition-of-unity injection. Although that paper does not use the term FROSch, its overlap, partition-of-unity, spectral selection, and robustness philosophy is explicitly described as closely aligned with modern fast and robust overlapping Schwarz methods [1901.00112].

For coupled multiphysics, several papers argue that different fields need different coarse spaces. Thermo-elastoplastic laser beam welding uses different GDSW variants for temperature and displacement and reports that richer mechanics coarse spaces with reduced thermal coarse spaces are advantageous [2502.19975]. Monolithic Navier–Stokes preconditioners likewise use different GDSW-type spaces for velocity and pressure [2506.16179]. This suggests that, in FROSch practice, the decisive question is often not whether a coarse level exists, but how its interface objects are chosen for each field.

## 4. Theory, robustness, and known failure modes

The classical theoretical estimate recalled in the GPU paper is
\[
\kappa \left( M_{\rm GDSW}^{-1} A \right) \leq C \left(1+\frac{H}{\delta}\right) \left(1+\log\left(\frac{H}{h}\right)\right)^2,
\]
with \(H\) the maximum subdomain diameter, \(\delta\) the overlap width, and \(h\) the fine mesh size [2304.04876]. This is the canonical explanation for why a two-level overlapping Schwarz method can be robust and scalable while a one-level method cannot. Closely related three-dimensional vector-field analysis for Nédélec and Raviart–Thomas discretizations establishes bounds of the form
\[
\kappa(M_\xi^{-1}A_\xi)\le C\max_i\left(1+\frac{H_i}{\delta_i}\right),
\]
with robustness to topology and nonconvexity rather than to coefficient contrast [2404.11986].

For high-contrast scalar diffusion, the standard piecewise-polynomial coarse problem is explicitly described as non-robust. The multiscale energy-minimization paper addresses exactly this failure mode by selecting local generalized eigenvectors below a threshold \(\Lambda\), constructing energy-minimizing basis functions constrained by these auxiliary modes, proving exponential decay on oversampled patches, and obtaining stable-decomposition constants independent of the coefficient contrast and overlap width, provided oversampling is large enough [1901.00112]. This paper is particularly significant because it isolates a design principle that reappears in later FROSch-adjacent work: detect problematic local modes spectrally, then realize the actual coarse correction by an energy-minimizing extension.

At the same time, published FROSch work makes clear that robustness is not automatic. In the algebraic coarse-space comparison, GDSW and RGDSW are not robust for one high-conductivity inclusion-at-node configuration, while AMS is robust and scalable there; for parallel channel configurations all three coarse spaces yield condition numbers of the order of the contrast [2408.08187]. The chemo-mechanics study reaches a different but related conclusion: fully algebraic black-box FROSch can be used successfully on a monolithic coupled problem, but numerical scalability should not be expected and is not observed in that mode [2212.00801]. The direct-solver comparison emphasizes a further limitation: fully algebraic use with an incomplete nullspace sacrifices the theoretical GDSW guarantee, and iteration counts rise with core count accordingly [2310.12659].

These papers collectively correct two recurrent misconceptions. The first is that any GDSW-type coarse space is automatically coefficient-robust; the literature instead shows strong dependence on coefficient geometry, interface traces, and possible adaptive enrichment [2408.08187, 1901.00112]. The second is that “fully algebraic” is synonymous with “numerically scalable”; the chemo-mechanics and direct-solver studies explicitly state the opposite for their settings [2212.00801, 2310.12659].

## 5. Software realization and parallel performance

FROSch is embedded in the Trilinos ecosystem. The deal.II coupling paper describes it as accessible through Stratimikos and compatible with both Epetra and Tpetra through Xpetra, with Tpetra identified as the forward-looking path because it enables Kokkos and KokkosKernels [2410.22871]. Several application papers emphasize that the implementation is algebraic in the sense that only the global matrix and the parallel distribution of the degrees of freedom are needed, after which FROSch builds overlap information, local subdomain matrices, restriction and prolongation operators, and coarse matrices [2508.04356, 2512.01706].

The implementation footprint has broadened substantially. FROSch has been coupled to FEATFLOW for monolithic fluid solvers [2508.04356], to deal.II through a Tpetra-based interface supporting both Lagrange and Nédélec elements [2410.22871], to NGSolve through an interface from NGSolve’s Python frontend to Trilinos’ C++ implementation [2512.01706], and to FEDDLib for nonlinear two-level Schwarz methods [2410.24138, 2603.24542]. The chemo-mechanics paper uses deal.II and FROSch in fully algebraic mode for monolithic nonlinear deformation–diffusion problems [2212.00801].

Parallel performance studies show both the strength and the cost structure of the framework. On Summit with NVIDIA V100 GPUs, the GPU study reports that solve time was reduced by factors of about \(2\times\) using GPUs, and that GPU-friendly iterative ILU variants yielded solve-time speedups around \(2.5\times\) to \(3.8\times\) in weak scaling and up to about \(4.8\times\) in a single-node study [2304.04876]. A central implementation result there is that multiple MPI processes per GPU, enabled through NVIDIA MPS, improved both numerical setup and solve times by shrinking local problems and reducing fill-heavy factorization costs [2304.04876].

For CPU-only direct-solver studies, the picture is different. The comparison of direct solvers in FROSch on elasticity and chemo-mechanics reports good strong scaling from 64 to 1728 cores, with deterioration beyond that because the single-core coarse solve becomes dominant; among the tested solvers, MUMPS is consistently the fastest overall [2310.12659]. This coarse-bottleneck diagnosis reappears in the monolithic multilevel fluid work, where the two-level method fails at 32,768 subdomains because the coarse factorization becomes infeasible, while the three-level recursive extension remains solvable and reduces total time once the second-level coarse problem has become too large [2508.04356].

A further software development is the addition of ORAS to FROSch. The deal.II interface paper introduces new geometric classes for O(R)AS and emphasizes that optimized local operators \(B_i\) must be assembled separately from finite-element and interface data. In that sense, ORAS marks a shift from the purely algebraic local-operator model toward a geometry-aware first level, while retaining FROSch’s operator-composition viewpoint [2410.22871].

## 6. Application domains and ongoing extensions

The application range of FROSch now spans scalar elliptic problems, elasticity, incompressible flow, thermo-mechanics, chemo-mechanics, fluid–structure interaction, Maxwell equations, vector-field problems, and virtual element discretizations. In incompressible fluid mechanics, monolithic two-level and three-level overlapping Schwarz solvers implemented with FROSch and FEATFLOW are reported up to 32,768 MPI ranks for Poiseuille flow and an extrusion die geometry [2508.04356]. A broader Navier–Stokes study uses FROSch both as a monolithic preconditioner for the full saddle-point operator and as an inner solver inside block preconditioners such as PCD, SIMPLEC, and LSC, concluding that monolithic FROSch preconditioners are generally more robust, especially with respect to viscosity, Reynolds number, and CFL number [2506.16179].

For reduced-integration Stokes discretizations, FROSch is used with an RGDSW coarse space on an SPD velocity-only system, and the reported weak-scaling experiments show that the two-level RGDSW preconditioner has good numerical scalability compared with the one-level method, while setup time stays nearly constant over the tested penalty-parameter range [2512.01706]. In thermo-elastoplastic laser beam welding, monolithic two-level overlapping Schwarz with different combinations of GDSW, GDSW\(^\star\), and RGDSW coarse spaces is implemented in PETSc rather than Trilinos, but the method is explicitly described as a GDSW-family overlapping Schwarz construction for multiphysics saddle-point systems and is therefore directly relevant to FROSch-style design [2502.19975]. The earlier monolithic chemo-mechanics study, by contrast, uses FROSch in fully algebraic mode and shows that it remains applicable on the tested core counts even though numerical scalability is not observed [2212.00801].

Nonlinear two-level Schwarz is the most substantial recent conceptual extension. The 2024 FROSch paper introduces one-level, additive two-level, and hybrid two-level nonlinear Schwarz solvers on top of FROSch, with local nonlinear corrections \(T_i(u)\), a nonlinear coarse correction \(T_0(u)\), and exact tangents of the nonlinear Schwarz map [2410.24138]. The 2026 follow-up reports that the same framework, using FROSch GDSW-type coarse spaces and modifications of those coarse spaces, scales exceptionally well up to 9,000 subdomains and is more robust than standard Newton–Krylov–Schwarz for high-Reynolds-number lid-driven cavity flow and large-deformation Neo-Hookean elasticity [2603.24542]. One of the striking conclusions of that paper is that apparently small modifications of RGDSW or MsFEM interface functions can materially change nonlinear robustness, especially near Dirichlet boundaries [2603.24542].

Related research continues to expand the conceptual perimeter of FROSch. The vector-field analysis for \(H(\mathrm{curl})\) and \(H(\mathrm{div})\) problems supplies topology-robust stable-decomposition theory with linear dependence on relative overlap [2404.11986]. The VEM paper develops a two-level additive overlapping Schwarz preconditioner on irregular three-dimensional polyhedral decompositions, numerically observing weak dependence on the number of subdomains and on \(H/h\), and approximately linear dependence on \(H/\delta\), with a geometric vertex-based interface-harmonic coarse space that is close in spirit to reduced GDSW constructions [2512.07181]. This suggests that FROSch is best understood not only as a concrete Trilinos package, but also as a continuing research program centered on overlapping local solves, interface-based coarse spaces, and the search for coarse corrections that remain robust under stronger coupling, more irregular discretizations, and more demanding parallel scales.

Source: https://www.emergentmind.com/topics/fast-and-robust-overlapping-schwarz-frosch