---
title: Stepwise Relaxation Algorithm
url: https://www.emergentmind.com/topics/stepwise-relaxation-algorithm-sr
type: topic
---

# Stepwise Relaxation Algorithm

Stepwise Relaxation Algorithm (SR) denotes, in the cited literature, a class of iterative procedures in which relaxation is applied step by step rather than through a single stationary update. In numerical linear algebra, SR most directly refers to Successive (or Stepwise) Relaxation for solving \(Ax=b\), either with a fixed relaxation factor \(\omega\), with self-adapted \(\omega_i\) carried by an evolutionary population, or with a scheduled cycle \((\omega_1,\dots,\omega_M)\) of Jacobi weights. In adjacent literatures, related stepwise procedures include forward–backward graph construction in Gaussian graphical models and coarse-grained switching between trapped and high-dissipation regimes in dispersively coupled oscillator–particle systems [1304.2097][2112.06854][1808.06016][1706.02992].

## 1. Terminological scope and core abstraction

In iterative numerical linear algebra, “SR” stands for *Successive (or Stepwise) Relaxation*. The basic idea is to replace a full Jacobi update by a relaxed step that moves only a fraction of the way from the current iterate toward the Jacobi iterate. For \(A=D+L+U\), with \(D\) diagonal, \(L\) strictly lower triangular, and \(U\) strictly upper triangular, the Jacobi fixed-point form is
\[
x = Hx + V,\qquad H = D^{-1}(-L-U),\qquad V = D^{-1}b.
\]
The plain Jacobi iteration is
\[
x^{(k+1)} = Hx^{(k)} + V,
\]
and the relaxed step is
\[
x^{(k+1)} = x^{(k)} + \omega\big(x^{J,(k+1)}-x^{(k)}\big).
\]
Componentwise, the Jacobi-SR update is
\[
x_i^{(k+1)} = x_i^{(k)} + \omega\frac{1}{a_{ii}}\left(b_i-\sum_{j=1}^n a_{ij}x_j^{(k)}\right),\qquad i=1,\dots,n.
\]
The cases \(\omega=1\), \(0<\omega<1\), and \(1<\omega<2\) correspond respectively to Jacobi, under-relaxation, and over-relaxation [1304.2097].

The Scheduled Relaxation Jacobi (SRJ) method is a non-stationary extension of the same principle. Instead of using one fixed \(\omega\), SRJ applies a prescribed sequence of relaxation factors \(\{\omega_i\}_{i=1}^M\) within a cycle and then repeats that cycle. The paper itself uses the name SRJ rather than “SR”, but conceptually SRJ is a stepwise relaxation algorithm: each iteration uses a possibly different relaxation factor \(\omega_k\), according to a precomputed schedule [2112.06854].

A broader use of “stepwise relaxation” appears outside classical solvers. In high-dimensional Gaussian graphical models, the graphical stepwise algorithm starts from an empty graph and iteratively adds or removes edges based on residual correlations; the paper does not explicitly use the label “Stepwise Relaxation (SR)”, but states that the method fits naturally the idea of a stepwise relaxation of the structural zeros in \(\Omega\) [1808.06016]. In nonlinear oscillator–particle systems, stepwise relaxation refers instead to a physical relaxation pattern with alternating plateaus and rapid drops in energy dissipation [1706.02992].

## 2. Classical SR for linear systems and its spectral properties

For the stationary Jacobi-SR method, the relaxed iteration can be written in matrix form as
\[
x^{(k+1)} = H_\omega x^{(k)} + V_\omega,
\]
with
\[
H_\omega = D^{-1}\big((1-\omega)I-\omega(L+U)\big),\qquad V_\omega = \omega D^{-1}b.
\]
Here \(H_\omega\) is the Jacobi-SR iteration matrix. The paper constrains \(\omega\) to a domain
\[
\omega\in(\omega_L,\omega_U),\qquad 0<\omega_L<\omega_U<2,
\]
because a necessary condition for convergence of SR is \(0<\omega<2\). More generally, convergence of the fixed-point iteration requires
\[
\rho(H_\omega)<1,
\]
so the asymptotic rate is governed by the spectral radius of \(H_\omega\) [1304.2097].

The central practical issue is sensitivity to \(\omega\). The paper assumes an optimal relaxation factor \(\omega^*\) such that \(\rho(\omega)\) is strictly decreasing for \(\omega<\omega^*\) and strictly increasing for \(\omega>\omega^*\). Under that assumption, SR has a single best relaxation factor but is highly problem-dependent. The experimental examples make this explicit. For the system
\[
a_{ii}=2n,\qquad a_{ij}=j\ (i\neq j),\qquad b_i=i,\qquad n=100,
\]
the authors note that the optimal relaxation is near to \(1.20\), while classical Jacobi-SR with \(\omega=0.5\) or \(\omega=1.5\) failed to achieve error below \(10^{-12}\) in 1000 iterations. For a second problem with
\[
a_{ii}=n^2,\qquad a_{ij}=j,\qquad b_i=i,
\]
Jacobi-SR converges in 18 iterations for \(\omega\approx 0.81\), whereas many nearby choices diverge. The numerical message is that slight changes around \(\omega^*\) can shift the method from rapid convergence to divergence [1304.2097].

This sensitivity is the main reason SR becomes a design problem rather than a mere scalar tuning problem. A fixed-\(\omega\) scheme is simple and fully parallelizable in the Jacobi case, but robust performance depends on spectral information that is often unavailable or expensive to estimate directly.

## 3. Adaptive and scheduled Jacobi variants

Two distinct stepwise generalizations of Jacobi-SR have been studied. The first is the Jacobi-SR Based Uniform Adaptive hybrid evolutionary algorithm (JBUA). Each individual \(x_i^{(k)}\in\mathbb{R}^n\) carries its own relaxation factor \(\omega_i\), and the population is
\[
X^{(k)}=\{x_1^{(k)},x_2^{(k)},\dots,x_N^{(k)}\}.
\]
Recombination uses a stochastic mixing matrix \(R\) to form
\[
X^{(k,c)} = R\,X^{(k)}.
\]
Mutation is one Jacobi-SR step,
\[
x_i^{(k,m)} = H_{\omega_i^{(k)}}x_i^{(k,c)} + V_{\omega_i^{(k)}},\qquad i=1,\dots,N.
\]
Fitness is measured by \(\|e(x)\|=\|Ax-b\|\). The uniform adaptation rule compares adjacent ranked individuals: the worse \(\omega_x\) is moved toward the better \(\omega_y\) through
\[
\omega_x'=(0.5+p_x)(\omega_x+\omega_y),\qquad p_x\in(-0.01,0.01),
\]
while the better factor is moved away from the worse one within \((\omega_L,\omega_U)\) using \(p_y\in(0.008,0.012)\). Selection keeps the best \(N/2\) individuals and duplicates them to form the next generation. A convergence theorem states that if there exists \(0<\varepsilon<1\) such that \(\|H_{\omega_i}\|\le \varepsilon<1\) for all relevant \(\omega_i\), then every individual converges to the exact solution \(x^*\). An adaptation theorem further shows, under monotonicity assumptions on \(\rho(\omega)\), that the update rules move relaxation factors toward lower spectral radius [1304.2097].

The second generalization is SRJ, where stepwise relaxation is scheduled rather than learned online. For the Jacobi iteration matrix \(B_{\mathrm J}=-D^{-1}(L+U)\), one relaxed step uses
\[
B_\omega=(1-\omega)I+\omega B_{\mathrm J}.
\]
A cycle of length \(M\) has iteration matrix
\[
B_{\mathrm{SRJ}}=\prod_{i=1}^{M} B_{\omega_i},
\]
and amplification polynomial
\[
G_M(\lambda;\omega_i)=\prod_{i=1}^{M}\big[(1-\omega_i)+\omega_i\lambda\big].
\]
The nonsymmetric extension constructs \(\{\omega_i\}\) by minimizing the maximum of \(|G_M(\lambda)|\) over an ellipse
\[
\mathcal E=\left\{z=x+iy:\ \frac{(x-x_c)^2}{a^2}+\frac{y^2}{b^2}\le 1\right\},
\]
with
\[
a=\frac{\lambda_{\max}+1}{2},\qquad b=ca,\qquad x_c=\frac{\lambda_{\max}-1}{2}.
\]
Numerically, the scheme is obtained from the constrained problem
\[
\min_{(\omega_i,\bar g)}\ \bar g^2
\]
subject to
\[
|G_M(x_j;\omega_i)|^2\le \bar g^2
\]
at test points on the ellipse boundary. The implementation described in the paper uses `scipy.optimize.minimize` with the Trust-Region Constrained algorithm. In 1D steady advection–diffusion with \(N=128\), \(\nu=1\), and \(a=300\), the \(c=0\) scheme offers almost no benefit compared to Jacobi, whereas the \(c=\tfrac12\) scheme is significantly faster. In 2D, for \(256\times 256\) grids with \(a_x=a_y=400\), the \(c=\tfrac13\) scheme is the best and robust, while \(c=0,\tfrac{1}{10},\tfrac{1}{5}\) typically diverge then stagnate [2112.06854].

| Variant | Stepwise mechanism | Design principle |
|---|---|---|
| JBUA | One Jacobi-SR step per generation with evolving \(\omega_i\) | Uniform self-adaptation guided by \(\|Ax-b\|\) |
| SRJ | One cycle of \(M\) relaxed Jacobi steps with \(\omega_1,\dots,\omega_M\) | Minimize \(\max_{\lambda\in\mathcal E}|G_M(\lambda;\omega_i)|\) |

A major structural distinction concerns parallelism. Because Jacobi-SR uses only the old iterate \(x^{(k)}\), all components \(x_i^{(k+1)}\) can be computed independently and simultaneously. The JBUA paper therefore argues that Jacobi-SR based hybrids inherently can be implemented in parallel processing environments efficiently, whereas Gauss-Seidel-SR based hybrids cannot [1304.2097].

## 4. Stepwise relaxation of graph structure in Gaussian graphical models

In Gaussian graphical models, the stepwise procedure is not a relaxation parameter method but a forward–backward algorithm on graph support. Let \(\mathbf X=(X_1,\dots,X_p)^\top\sim N_p(\mathbf 0,\Sigma)\), with precision matrix \(\Omega=\Sigma^{-1}\). Conditional independence is equivalent to \(\omega_{ij}=0\), and the partial correlation satisfies
\[
\rho_{il\cdot V\setminus\{i,l\}}=-\frac{\omega_{il}}{\sqrt{\omega_{ii}\omega_{ll}}}.
\]
The paper’s key parametrization expresses edges through Pearson correlations between prediction errors of nodewise best linear predictors. If
\[
X_i=\sum_{k\in\mathcal A_i}\beta_{ik}X_k+\varepsilon_i,
\]
then
\[
E=\{(i,l):\ \operatorname{cov}(\varepsilon_i,\varepsilon_l)\neq 0\}
=\{(i,l):\ \rho_{\varepsilon_i,\varepsilon_l}\neq 0\}.
\]
The graphical stepwise algorithm starts from empty neighborhoods \(\widehat{\mathcal A}_j^0=\emptyset\). In the forward step it computes empirical residual correlations
\[
f_{jl}^k=\operatorname{corr}(\mathbf e_j^k,\mathbf e_l^k)
\]
for non-neighbors and adds the edge with largest \(|f_{jl}^k|\) if it exceeds \(\alpha_f\). In the backward step it computes
\[
b_{jl}^k=\operatorname{corr}(\mathbf r_j^k,\mathbf r_l^k)
\]
for existing edges after leaving that edge out of the conditioning sets, and removes the edge with smallest \(|b_{jl}^k|\) if it is below \(\alpha_b\). The final precision estimate uses residual variances and covariances:
\[
\widehat\omega_{ii}=\frac{n}{\mathbf e_i^\top\mathbf e_i},
\qquad
\widehat\omega_{il}=
\frac{n\,\mathbf e_i^\top\mathbf e_l}{(\mathbf e_i^\top\mathbf e_i)(\mathbf e_l^\top\mathbf e_l)}
\]
for selected edges, and \(0\) otherwise. Thresholds \((\alpha_f,\alpha_b)\) are chosen by \(K\)-fold cross-validation minimizing
\[
CV(\alpha_f,\alpha_b)
=
\frac{1}{n}\sum_{t=1}^K \sum_{j=1}^p
\big\|\mathbf X_j^{(t)}-\widehat{\mathbf X}_j^{(t)}(\alpha_f,\alpha_b)\big\|_2^2.
\]
The paper states that this procedure fits naturally the idea of a “stepwise relaxation”: the support of \(\Omega\) starts empty and is gradually relaxed and refined, rather than being determined in one global optimization [1808.06016].

The empirical comparisons emphasize graph recovery rather than convex optimality. Across AR(1), NN(2), and BG models with \(p\in\{50,100,150\}\) and \(n=100\), the method substantially outperforms graphical lasso and CLIME in MCC. Representative values include AR(1), \(p=100\): GS 0.751 versus 0.433 for Glasso and 0.464 for CLIME; NN(2), \(p=100\): GS 0.802 versus 0.382 and 0.407; BG, \(p=100\): GS 0.857 versus 0.348 and 0.461. In the breast cancer gene expression example, GS attains MCC 0.520 (0.02), CLIME 0.516 (0.02), and Glasso 0.334 (0.02), with 54 selected edges versus 4823 and 2103. The paper also notes that no formal theorems about model selection consistency or rates of convergence are provided [1808.06016].

## 5. Stepwise relaxation as a dynamical regime in dispersively coupled oscillators

A different usage appears in the study of degenerate harmonic oscillator modes dispersively coupled to particle positions. The system has Hamiltonian
\[
H = \frac{1}{2}\sum_n\left(\dot q_n^2+\Omega_n^2q_n^2\right)
+\frac{\mu}{2}\sum_k \dot x_k^2
-\frac{1}{2}\sum_{k,n,m} g_{mn}q_n q_m \phi_n(x_k)\phi_m(x_k),
\]
so the particle experiences an effective inertial potential
\[
U(\text{position},t)\propto \left(\sum_n q_n(t)\phi_n(\text{position})\right)^2.
\]
The key observation is that the oscillator energy \(E(t)\) does not decay smoothly. Instead, the dynamics shows plateaus, where the particle is trapped in a local minimum of \(U\) and dissipation is slow, alternating with rapid drops, where the particle traverses the domain in an irregular fashion and transfers energy efficiently from the oscillators to the dissipative particle bath [1706.02992].

For the minimal 1D prototype with a degenerate doublet and one particle, the equations of motion are
\[
\ddot q_1+\left[\Omega^2-g\phi_1^2(x)\right]q_1-g\phi_1(x)\phi_2(x)q_2=0,
\]
\[
\ddot q_2+\left[\Omega^2-g\phi_2^2(x)\right]q_2-g\phi_1(x)\phi_2(x)q_1=0,
\]
\[
\ddot x+\Gamma\dot x-\frac{g}{2\mu}\partial_x\big[(q_1\phi_1(x)+q_2\phi_2(x))^2\big]=0.
\]
During plateau phases, a rotating-wave approximation yields slow amplitude dynamics,
\[
\dot a_1=-\frac{i g}{2\Omega}\big(\phi_1^2 a_1+\phi_1\phi_2 a_2\big),\qquad
\dot a_2=-\frac{i g}{2\Omega}\big(\phi_1\phi_2 a_1+\phi_2^2 a_2\big),
\]
with slow mixing frequency
\[
\lambda\approx \frac{g}{2\Omega}\big(\phi_1^2(x_{\rm eq})+\phi_2^2(x_{\rm eq})\big).
\]
In the high-dissipation regime, the paper reports
\[
-\frac{\partial E}{\partial t}\propto E(0)\Gamma
\qquad\text{for }E(0)\Gamma<1,
\]
which leads to an effective exponential relaxation picture through particle thermalization [1706.02992].

In the membrane model, stepwise relaxation appears only in the presence of thermal noise. Noise enables escape and retrapping, and the degenerate doublet exhibits a slow-in-time stochastic precession of the mixing angle. The power spectral density of the angular velocity has a low-frequency Lorentzian whose width scales as
\[
\Delta\omega\sim \mu^2\Gamma.
\]
The paper does not explicitly propose this as an algorithm, but the logic of its analysis naturally leads to a two-regime SR scheme: micro-integration of the coupled equations, trapping-state detection, regime-dependent dissipation laws, and stochastic evolution of the mixing angle in degenerate noisy plateaus [1706.02992].

## 6. Comparative interpretation and recurrent misconceptions

Several recurrent distinctions are essential for interpreting “SR” correctly. First, stationary SR and non-stationary SR are different objects. A fixed-\(\omega\) Jacobi-SR or Gauss-Seidel-SR iteration modifies each step by the same relaxation factor, whereas SRJ uses a prescribed cycle of factors and JBUA adapts \(\omega_i\) online. The scheduled version is polynomially accelerated; the adaptive evolutionary version is population-based; the classical version is a single fixed-point iteration [1304.2097][2112.06854].

Second, SRJ is not restricted to elliptic problems. Earlier SRJ work optimized amplification on a real interval and therefore targeted symmetric positive definite elliptic systems, but the nonsymmetric extension replaces the real interval by an ellipse in the complex plane and computes \(\{\omega_i\}\) from a constrained minimization problem. The paper explicitly shows that old Chebyshev-based schemes can have \(|G_M(\lambda)|>1\) inside \(|\lambda|\le 1\), so SRJ can diverge while Jacobi still converges if the design ignores complex spectra [2112.06854].

Third, the graphical stepwise procedure is not a one-shot convex estimator of \(\Omega\). Its logic is local and structural: add one edge whose residual correlation exceeds \(\alpha_f\), then possibly remove one edge whose residual correlation falls below \(\alpha_b\). This is why the paper contrasts it with graphical lasso and CLIME, which solve global penalized or constrained optimization problems [1808.06016].

Fourth, in the oscillator–particle setting, “stepwise relaxation” names a relaxation pattern rather than a numerical solver. Plateaus and drops arise from inertial trapping, release, and retrapping of particles, with stochastic precession as a degenerate-mode effect in the noisy membrane problem [1706.02992].

Taken together, these usages suggest that “Stepwise Relaxation Algorithm” is best understood as a procedural pattern rather than a single canonical algorithm. The common structure is step-by-step control of a system that is otherwise hard to optimize or stabilize globally: by adjusting \(\omega\), by scheduling \(\omega_1,\dots,\omega_M\), by relaxing the support of \(\Omega\) through forward–backward edge updates, or by switching between coarse-grained dissipation regimes. The precise technical content of SR therefore depends on domain, but its most developed meaning remains the Jacobi-based relaxation family for linear systems and PDE discretizations [1304.2097][2112.06854].

Source: https://www.emergentmind.com/topics/stepwise-relaxation-algorithm-sr