---
title: Non-Iterative Domain Decomposition Integrator
url: https://www.emergentmind.com/topics/non-iterative-domain-decomposition-integrator
type: topic
---

# Non-Iterative Domain Decomposition Integrator

In the cited literature, a non-iterative domain decomposition integrator denotes a domain decomposition scheme for an evolutionary problem in which the solution is advanced from one time level to the next by a fixed number of algebraic or subdomain-local substeps, rather than by inner Schwarz, fixed-point, or Krylov iterations on every time slab. The common construction begins from an evolution equation such as $\frac{du}{dt}+Au=f(t)$ or its finite element form $B\frac{dy}{dt}+Ay=\varphi(t)$, decomposes the spatial operator into subdomain contributions, and embeds the coupling into an operator-splitting, factorized, projection, or Robin-transmission update. This viewpoint is explicit in the regionally additive schemes for parabolic and first-order evolutionary equations, in additive Stokes solvers in primitive variables, in partitioned fluid–structure interaction algorithms, and in recent acoustic-wave integrators based on local prediction plus subdomain Crank–Nicolson solves [1401.0798], [1101.2395], [1102.0642], [2007.00781], [2507.19379].

## 1. Definition and distinguishing criteria

A central distinction in this literature is between domain decomposition as an **iterative solver** and domain decomposition as a **time integrator**. For time-dependent problems, the iteration-free approach constructs a splitting so that “a transition to a new time level is performed via solving problems in particular subdomains,” with no convergence loop at fixed time level [1102.0642]. In the parabolic setting, this is described as an “iteration-free domain decomposition algorithm” in which the advance from $t^n$ to $t^{n+1}$ is realized by a fixed number of substeps after an additive decomposition of the operator and the right-hand side [1401.0798]. The same formulation appears in the two-component overlapping scheme for parabolic equations, where each time step requires only a fixed sequence of local solves and “there is no inner Schwarz iteration or fixed-point iteration to convergence” [1705.03434].

The classification of such methods is tied to three choices: the method of domain decomposition, the choice of decomposition operators, and the splitting scheme employed [1401.0798]. Within this framework, overlapping subdomain methods are preferred when one seeks homogeneous numerical algorithms, because the overlap and the associated weighting functions encode interface exchange directly inside the decomposed operators. Substructuring variants reorganize the same idea by separating interior subdomain nodes from common boundary or interface nodes, still within a fixed-step, non-iterative update [1105.3448].

This definition excludes several nearby but distinct classes. BDD and FETI, for example, are iterative Krylov substructuring solvers for interface problems, even when their iterates admit useful local error decompositions [1312.4505]. Likewise, some nonlocal substructuring frameworks establish exact multi-domain equivalence without yet providing an iteration-free solver [2008.11780]. The term therefore has a stricter meaning than “domain decomposition without a global monolithic matrix”; it refers to the absence of inner domain-decomposition iterations in the update itself.

## 2. Abstract operator formulation

The basic analytical template is an evolution equation with a nonnegative operator,
\[
\frac{du}{dt}+Au=f(t), \qquad u(0)=u^0,
\]
or, after finite element semidiscretization of a parabolic problem,
\[
B\frac{dy}{dt}+Ay=\varphi(t), \qquad y(0)=y^0,
\]
with $B=B^*>0$ and $A=A^*>0$ [1401.0798]. A standard two-level weighted scheme is then
\[
\frac{y^{n+1}-y^n}{\tau}+A y^{n+\sigma}=\varphi^{n+\sigma},\qquad
y^{n+\sigma}=\sigma y^{n+1}+(1-\sigma)y^n,
\]
and, in the non-self-adjoint first-order setting, unconditional stability holds for $\sigma\ge \frac12$ with transition operator
\[
S=(E+\sigma\tau A)^{-1}(E-(1-\sigma)\tau A),
\qquad \|S\|\le 1
\]
[1101.2395]. The self-adjoint parabolic theory uses the analogous weighted estimate and then embeds domain decomposition into the operator itself [1401.0798].

The decomposed form is additive:
\[
A=\sum_{\alpha=1}^p A_\alpha,\qquad \varphi=\sum_{\alpha=1}^p \varphi_\alpha,
\]
with $A_\alpha=A_\alpha^*\ge 0$ in the standard parabolic formulation [1401.0798]. For non-self-adjoint problems, the decomposition is applied after splitting
\[
A=C+D,\qquad C=\frac12(A-A^*),\qquad D=\frac12(A+A^*),
\]
so that $D=D^*\ge 0$ and $C=-C^*$, and then constructing $A_\alpha=C_\alpha+D_\alpha$ with $A_\alpha\ge 0$ preserved by design [1101.2395]. This preservation of nonnegativity is the key algebraic property behind the unconditional stability of the regularized additive schemes.

A further technical device is symmetrization. In the parabolic finite element setting one introduces
\[
w=B^{1/2}v,\qquad \widetilde A=B^{-1/2}AB^{-1/2},\qquad \widetilde\varphi=B^{-1/2}\varphi,
\]
which yields
\[
\frac{dw}{dt}+\widetilde A\,w=\widetilde\varphi(t), \qquad
\widetilde A=\sum_{\alpha=1}^p \widetilde A_\alpha,\quad
\widetilde A_\alpha=B^{-1/2}A_\alpha B^{-1/2}\ge 0
\]
[1401.0798]. This form places the DD integrator inside standard operator-splitting stability theory.

## 3. Overlap, partitions of unity, and interface corrections

The most characteristic spatial construction is an overlapping covering
\[
\Omega=\Omega_1\cup\Omega_2\cup\cdots\cup\Omega_p,
\]
together with partition-of-unity functions $\eta_\alpha(\mathbf x)$ such that
\[
\eta_\alpha(\mathbf x)>0 \text{ in }\Omega_\alpha,\qquad
\eta_\alpha(\mathbf x)=0 \text{ outside }\Omega_\alpha,\qquad
\sum_{\alpha=1}^p \eta_\alpha(\mathbf x)=1
\]
[1401.0798]. These weights localize both the operator and the forcing:
\[
a(u,v)=\sum_{\alpha=1}^p a_\alpha(u,v), \qquad (f,v)=\sum_{\alpha=1}^p (f_\alpha,v),
\]
with the standard choice
\[
(f_\alpha,v)=\int_\Omega \eta_\alpha(\mathbf x)\,f(\mathbf x,t)\,v\,d\mathbf x,\qquad
a_\alpha(u,v)=\int_\Omega \eta_\alpha(\mathbf x)\,k(\mathbf x)\,\nabla u\cdot\nabla v\,d\mathbf x
\]
[1401.0798]. Equivalent forms,
\[
a_\alpha(u,v)=\int_\Omega k(\mathbf x)\,\nabla u\cdot\nabla(\eta_\alpha(\mathbf x)v)\,d\mathbf x,
\qquad
a_\alpha(u,v)=\int_\Omega k(\mathbf x)\,\nabla(\eta_\alpha(\mathbf x)u)\cdot\nabla v\,d\mathbf x,
\]
show that the overlap is not merely geometric padding; it is part of the interface-transfer mechanism.

A different overlapping construction uses indicator functions rather than smooth weights. For two subdomains with overlap $\Omega_{12}=\Omega_1\cap\Omega_2\neq\varnothing$, define
\[
\chi_\alpha(\mathbf x)=
\begin{cases}
1,& \mathbf x\in \Omega_\alpha,\\
0,& \mathbf x\notin \Omega_\alpha,
\end{cases}
\qquad
\chi_{12}(\mathbf x)=
\begin{cases}
1,& \mathbf x\in \Omega_{12},\\
0,& \mathbf x\notin \Omega_{12},
\end{cases}
\]
so that
\[
\chi_1(\mathbf x)+\chi_2(\mathbf x)-\chi_{12}(\mathbf x)=1.
\]
This yields the corrected decomposition
\[
A=A_1+A_2-A_{12},
\]
where the overlap contribution is subtracted once so that it is not double counted [1705.03434]. Rewriting with
\[
\bar A_\alpha=A_\alpha-\frac12A_{12},\qquad
G=I+\frac{\tau}{2}A_{12},\qquad G=G^*\ge I,
\]
reduces the three-component overlap-corrected formulation to a standard positive two-component factorization [1705.03434].

Substructuring schemes isolate interface degrees of freedom even more explicitly. In the two-component interior/interface split, one chooses
\[
\chi_2(\mathbf x)=
\begin{cases}
1,& \mathbf x\in\widehat\omega,\\
0,& \mathbf x\notin\widehat\omega,
\end{cases}
\qquad
\chi_1=1-\chi_2,
\]
where $\widehat\omega$ is the set of internal boundary nodes. The resulting operator split,
\[
A=A_1+A_2,\qquad A_\alpha=\chi_\alpha A,
\]
separates the interior subdomain problem from the boundary/interface problem [1105.3448]. This suggests a common structural interpretation across overlapping and substructured methods: interface exchange is encoded algebraically in the decomposition operators rather than enforced by repeated transmission iterations.

## 4. Factorized updates and stability theory

The classical non-iterative update for two subdomains is the factorized two-stage scheme. In its weighted variational form it reads
\[
\left (\frac{u^{n+1/2} - u^{n}}{\tau},v \right ) + a_1(\sigma u^{n+1/2} + (1-\sigma) u^{n},v) + a_2(u^{n},v) = (f^{n+\sigma},v),
\]
\[
\left (\frac{u^{n+1} - u^{n}}{\tau},v \right ) + a_1(\sigma u^{n+1/2} + (1-\sigma) u^{n},v)
+ a_2(\sigma u^{n+1} + (1-\sigma) u^{n},v) = (f^{n+\sigma},v).
\]
For $\sigma=\tfrac12$ this gives the Peaceman–Rachford scheme, and for $\sigma=1$ the Douglas–Rachford scheme [1401.0798]. In matrix form,
\[
(B+\sigma\tau A_1)\,B^{-1}\,(B+\sigma\tau A_2)\frac{y^{n+1}-y^n}{\tau}+Ay^n=\varphi^n,
\]
which is the paper’s “regionally additive” factorization [1401.0798]. The corresponding two-stage operator form in the two-component overlapping formulation is also described as the iteration-free analogue of ADI-type methods [1705.03434].

The principal stability statement is that the factorized regionally additive difference scheme is unconditionally stable for $\sigma\ge \frac12$ [1401.0798]. In symmetrized coordinates this becomes
\[
(E+\sigma\tau \widetilde A_1)(E+\sigma\tau \widetilde A_2)\frac{v^{n+1}-v^n}{\tau}
+\widetilde A v^n=\widetilde\varphi^n,
\]
with estimate
\[
\|(E+\sigma\tau \widetilde A_2)v^{n+1}\|\le
\|(E+\sigma\tau \widetilde A_2)v^n\|+\tau\|\widetilde\varphi^n\|.
\]
For the indicator-function three-component decomposition, unconditional stability again holds for $\sigma\ge 0.5$, with
\[
\|Q y^{n+1}\| \le \|Q w^0\| + \tau\sum_{j=0}^n \|\varphi^{j+\sigma}\|,
\]
where
\[
Q=\left(I+\frac{\tau}{2}A_{12}\right)^{-1}
\left(I+\sigma\tau A_2-\frac{\sigma-1}{2}\tau A_{12}\right)
\]
[1705.03434].

For general first-order evolutionary equations with non-self-adjoint $A\ge 0$, the same threshold $\sigma\ge \frac12$ underlies two regularized non-iterative families. The additive regularized transition operator is
\[
\widetilde S=\frac1p\sum_{\alpha=1}^p
(E+\sigma p\tau A_\alpha)^{-1}(E-(1-\sigma)p\tau A_\alpha),
\]
and the multiplicative analogue is
\[
\widetilde S=\prod_{\alpha=1}^p
(E+\sigma\tau A_\alpha)^{-1}(E-(1-\sigma)\tau A_\alpha),
\]
both unconditionally stable for $\sigma\ge \frac12$ [1101.2395]. The same paper also develops vector additive schemes in which $p$ components are advanced by successive inversions of operators $E+\tau A_\alpha$, again without subdomain iterations.

For the unsteady Stokes equations in primitive variables, non-iterativity appears through a different split: a viscous evolution step followed by a pressure projection step. After decomposition into $u_\alpha=\chi_\alpha u$, the viscous substep is implemented by ordered triangular sweeps,
\[
U^{n+1/4}-U^n+\tau A_1U^{n+1/4}=\tau F^{n+1/2},\qquad
U^{n+1/2}-U^{n+1/4}+\tau A_2U^{n+1/2}=0,
\]
followed by additive componentwise pressure correction [1102.0642]. Theorem 5.1 states unconditional stability of the additive domain decomposition scheme, with final estimate
\[
\|U^{n+1}\|_m\le e^\tau \|U^n\|_m+\tau \|F^{n+1/2}\|_m^2
\]
[1102.0642].

Numerical behavior in the parabolic papers is more nuanced than the unconditional stability theory. Reported experiments show $\mathcal O(\tau^2)$ behavior for larger $\tau$ and $\mathcal O(\tau)$ for smaller $\tau$, decomposition error close to $\mathcal O}(h^{-1})$ in the reported mesh-refinement studies, and strong dependence on overlap width [1401.0798]. In the indicator-function study, reducing the overlap width decreases accuracy, refining the spatial grid also decreases accuracy, and the results indicate conditional convergence of the domain decomposition approximations [1705.03434].

## 5. Representative equation classes and modern extensions

The non-iterative template has been adapted beyond self-adjoint parabolic diffusion to incompressible flow, fluid–structure interaction, and wave propagation. The common pattern remains a fixed-step coupling: a local or subdomain solve uses previously available interface data, then a complementary local solve uses the newly computed state, and no interface convergence loop is introduced.

| Problem class | Characteristic integrator | Reported property |
|---|---|---|
| Parabolic equations | Factorized regionally additive Douglas–Rachford / Peaceman–Rachford updates | Unconditionally stable for $\sigma\ge \tfrac12$ |
| First-order evolutionary equations with non-self-adjoint operators | Additive, multiplicative, and vector regularized schemes | Unconditionally stable for $\sigma\ge \tfrac12$ |
| Stokes equations | Viscous step plus pressure projection in decomposed product space | Unconditionally stable |
| Thick-structure FSI | Sequential structure solve, geometry update, fluid solve with generalized Robin conditions | Unconditionally stable; $\mathcal O(\Delta t^\frac12)$ in time and optimal convergence in space |
| Acoustic wave equation | Interface prediction plus local Crank–Nicolson solves and averaging | Second-order in time; global convergence $\mathcal O}(h+\tau^2)$ under a CFL-type condition |

In moving-domain fluid–structure interaction, the non-iterative construction is based on generalized Robin transmission conditions obtained by combining the kinematic and dynamic interface conditions with a parameter $\alpha>0$:
\[
\alpha \boldsymbol\xi+\boldsymbol\sigma_S \boldsymbol n_S
=
\alpha \boldsymbol v\circ\mathcal A
-
J\,\boldsymbol\sigma_F \boldsymbol F^{-T}\boldsymbol n_F
\quad\text{on }\hat\Gamma.
\]
At each step, one performs: a structure solve using the previous fluid data, a geometry update by harmonic extension into the fluid domain, and a fluid solve on the updated domain using the new structure state [2007.00781]. The method is non-iterative because each time step requires only one structure solve, one geometry update, and one fluid solve. Its semi-discrete stability theorem provides an energy estimate with total discrete energy
\[
\mathcal E^n=\frac{\rho_F}{2}\|\boldsymbol v^n\|_{L^2(\Omega_F(t^n))}^2
+\frac{\rho_S}{2}\|\boldsymbol\xi^n\|_{L^2(\hat\Omega_S)}^2
+\frac12\|\boldsymbol\eta^n\|_S^2,
\]
and the convergence analysis proves $\mathcal O(\Delta t^{1/2})$ convergence in time and optimal spatial convergence under the stated regularity assumptions, while the fully discrete proof uses the CFL-like condition
\[
\Delta t \le \frac{\rho_F}{\alpha C_{TI}k^2}h
\]
[2007.00781].

For the linear acoustic wave equation, a more recent non-iterative time integrator combines a local explicit prediction step on interfaces with independent implicit Crank–Nicolson solves on overlapping subdomains, followed by nodal averaging [2507.19379]. In first-order form,
\[
\partial_t u=v,\qquad \partial_t v=Lu+f,\qquad Lu=c^2\Delta u,
\]
the interface values on artificial boundaries $\Gamma_i=\partial\Omega_i^\delta\cap\Omega$ are predicted by one explicit leapfrog step, local Crank–Nicolson problems are solved on each $\Omega_i^\delta$, and the results are averaged back to a globally continuous discrete function [2507.19379]. Under the prediction CFL bound
\[
\tau\,\|h\|_{\mathbb H_h\leftarrow \mathbb H_h}\le 4^2,
\]
the paper proves second-order accuracy in time and the final error estimate
\[
\|x(t_n)-x_n^{\text{DS}}\|\lesssim h+\tau^2.
\]
The same work reports that the allowed $\tau$ grows approximately linearly in the number of overlap layers $\ell$ and that the method permits much larger time steps than leapfrog while retaining the accuracy of Crank–Nicolson [2507.19379].

A discontinuous Galerkin variant extends this wave-equation strategy to higher-order approximations and heterogeneous material parameters. It combines a local leapfrog prediction on a narrow interface strip with a local Crank–Nicolson solve on each overlapping subdomain and then local exchange in overlap regions, avoiding global averaging and remaining non-iterative in the coupling step [2510.27314]. In the reported prism test, the wall time per step is $0.94$ s for DS with Cholesky and $1.52$ s for CN, with total simulation times $3021$ s and $4669$ s respectively, while the final-time relative $L^2$ errors are essentially the same [2510.27314]. This suggests that the non-iterative template is compatible with higher-order DG discretizations, not only with linear finite elements and mass lumping.

## 6. Limitations, adjacent methods, and terminological boundaries

A recurring misconception is that any domain decomposition method with local solves is “non-iterative.” The cited literature repeatedly rejects that equivalence. BDD and FETI remain iterative Krylov solvers for interface problems; what changes in the error-bound paper is not the solver type but the availability of a strict upper bound separating algebraic and discretization error,
\[
\|u_{ex}-u_N\|_{\Omega} \le \sqrt{r^T z} + \sqrt{\sum_s e_{\mathrm{CR},\Omega^{(s)}}^2\!\left(u_N^{(s)},\hat\sigma_N^{(s)}\right)},
\]
which yields an objective stopping criterion for the iterative solver [1312.4505]. This is a non-iterative **error decomposition**, not a non-iterative integrator.

The same boundary applies in nonlocal problems. A substructuring-based framework for nonlocal operators proves exact equivalence between the single-domain problem and a multi-domain reformulation, both continuously and after finite element discretization, but explicitly states that the subdomain equations remain coupled and that the actual parallel solution methods are deferred to later work [2008.11780]. A later nonlocal FETI-type solver performs non-iterative elimination of interior variables by Schur complements and non-iterative reconstruction by local back substitution, but the reduced interface problem is solved iteratively by projected preconditioned CG, so the overall method is hybrid rather than fully non-iterative [2306.00094].

There are also non-iterative domain decomposition methods that are not time integrators. For the Helmholtz equation, the Method of Difference Potentials reduces each subdomain problem to a Calderón boundary equation with projection, assembles a single global linear system for Dirichlet and Neumann traces, solves that system once, and reconstructs the subdomain fields independently [2103.12172]. The method is explicitly non-iterative in its interface coupling, but it is a stationary boundary solver rather than an evolutionary integrator.

Even within the true time-integration setting, the absence of inner iterations does not imply decomposition-independent accuracy. The parabolic and substructuring papers show that overlap width strongly affects the error, that increasing the number of subdomains can worsen accuracy, and that the schemes can be conditionally convergent with respect to decomposition parameters [1401.0798], [1105.3448], [1705.03434]. In thick-structure FSI, the combination parameter $\alpha$ must be tuned problem by problem; moderate values such as $10$ or $100$ behave best in the reported tests, whereas very large $\alpha$ can degrade convergence rate [2007.00781]. In wave propagation, the CFL-type restriction in the convergence proof depends on the overlap width, so larger overlap permits larger stable time steps [2507.19379].

Taken together, these results define a precise technical meaning for the subject. A non-iterative domain decomposition integrator is not merely a decomposed solver, nor merely a substructured formulation. It is a time-marching algorithm in which operator decomposition, overlap or interface correction, and factorized or partitioned updates are organized so that one advances the solution by a prescribed local sequence, with stability and convergence controlled analytically rather than by an inner interface iteration.

Source: https://www.emergentmind.com/topics/non-iterative-domain-decomposition-integrator