Alphabeta Math
TheoremStatement: Literature-sourcedProof: AI-adaptedPipeline-generatedprecheck pass
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

SOR converges exactly for positive-definite A and relaxation parameter in (0, 2)

Statement

Let n≥1, A∈Mn(C) be Hermitian with strictly positive diagonal, and let ω∈R∖{0}. Write A=D−L−L∗ by Diagonal, lower, and upper parts in the stationary-iteration convention, and define

Mω=ω−1D−L,Nω=(ω−1−1)D+L∗,Gω=Mω−1Nω.

The SOR iteration is x(k+1)=Gωx(k)+Mω−1b. It converges to the unique solution of Ax=b for every b,x(0)∈Cn if and only if A is positive definite (meaning x∗Ax>0 for every nonzero x) and 0<ω<2.

Facts & Assumptions

Given: n≥1, a Hermitian matrix A with positive diagonal, a real ω≠0, and arbitrary b,x(0)∈Cn.

[A1]

A=A∗, aii>0, and ω∈R∖{0}.

[F1]

SOR uses Mω=ω−1D−L and Nω=(ω−1−1)D+L∗ (Finite-budget successive over-relaxation stationary iteration).

[L1]

The fixed diagonal/lower/upper convention gives A=D−L−L∗ for Hermitian A (Diagonal, lower, and upper parts in the stationary-iteration convention).

[L2]

For an invertible matrix M and A=M−N, the stationary update is x(k+1)=M−1Nx(k)+M−1b (Stationary iteration from a matrix splitting A=M−N).

[L3]

This iteration converges to the unique solution for every right-hand side and every start exactly when ρ(M−1N)<1 (A stationary splitting converges for every start if and only if its iteration matrix has spectral radius below 1).

[L4]

Hermitian positive definite means A=A∗ and x∗Ax>0 for every nonzero x (Hermitian positive-definite matrices and Cholesky factorisation A = LL* with positive diagonal).

[L5]

A lower-triangular system with nonzero diagonal has a unique forward substitution solution (Forward and backward substitution are correct, unique, and quadratic in scalar operations).

[L6]

Determinants are multiplicative for same-sized square matrices over a commutative ring (For same-sized finite square matrices over a commutative ring, det⁡(AB)=det⁡(A)det⁡(B)).

[L7]

The determinant of a triangular matrix is the product of its diagonal entries (The determinant of a triangular matrix is the product of its diagonal entries).

[L9]

Every complex polynomial splits over C (The complex numbers are algebraically closed).

Proof

1.1A1F1L1L2L5L6L7algebra

By [A1] and [F1], Mω is lower triangular with diagonal aii/ω≠0, hence invertible by [L5]. Also A=Mω−Nω by [L1], so [L2] applies with iteration matrix Gω. The matrix Nω is upper triangular with diagonal entries (ω−1−1)aii. Thus [L6]–[L7] give det⁡(Mω)=ω−n∏iaii and det⁡(Nω)=(ω−1−1)n∏iaii. From MωGω=Nω and multiplicativity, det⁡(Gω)=det⁡(Nω)/det⁡(Mω)=(1−ω)n.

1.2A1F1L1L2L4algebra

For sufficiency, suppose A is positive definite and 0<ω<2. Then Q:=Mω+Mω∗−A=(2/ω−1)D is positive definite. For an eigenpair Gωv=λv with v≠0, one has λ≠1 since Av=(1−λ)Mωv and A is invertible. Put a=v∗Av>0. Then v∗Mωv=a/(1−λ), so v∗Qv=a(2Re⁡11−λ−1)=a1−∣λ∣2∣1−λ∣2>0. Therefore ∣λ∣<1 for every eigenvalue and ρ(Gω)<1.

2.1L3L8L9step 1.1algebra

If SOR converges for every right-hand side and start, [L3] gives ρ(Gω)<1. By [L8]–[L9], ∣det⁡(Gω)∣=∏j=1n∣λj∣≤ρ(Gω)n<1. Step 1.1 therefore gives ∣1−ω∣<1, which for real ω is equivalent to 0<ω<2.

2.2A1F1L1L2L3step 1.1algebra

Assume now that SOR converges for every right-hand side and start. Then ρ(Gω)<1 by [L3], so 1 is not an eigenvalue of Gω and I−Gω is invertible. By [L2], A=Mω(I−Gω) is invertible. For b=0, its unique solution is x∗=0. For arbitrary nonzero e(0):=x(0), set e(k):=x(k)−x∗=x(k) and δk=e(k+1)−e(k). The homogeneous SOR recurrence gives Mωe(k+1)=Nωe(k), hence Ae(k)+Mωδk=0. Expanding the Hermitian quadratic form gives (e(k))∗Ae(k)−(e(k+1))∗Ae(k+1)=δk∗(Mω+Mω∗−A)δk=(2/ω−1)δk∗Dδk. The last equality follows from [L1] and [F1].

3.1A1L3L4step 2.1step 2.2algebra

By step 2.1, c:=2/ω−1>0. Summing the identity in step 2.2 from k=0 to N−1 and taking N→∞ (the error tends to 0 by the assumed convergence) gives (e(0))∗Ae(0)=c∑k=0∞δk∗Dδk. Moreover δ0≠0: otherwise step 2.2 implies Ae(0)=0, contrary to invertibility of A. Since D has positive diagonal, the first summand is strictly positive and all others are nonnegative. Thus (e(0))∗Ae(0)>0 for every nonzero e(0), so [L4] gives that A is positive definite.

4.1L3step 2.1step 3.1step 1.2∎

By [L3], ρ(Gω)<1 implies convergence to the unique solution for every right-hand side and every initial vector. Steps 2.1 and 3.1 prove the reverse implication.

Remarks

The conclusion is exact-arithmetic asymptotic convergence; it is not a one-step error decrease or a finite-precision accuracy guarantee. No optimal-ω formula is asserted.

Depends on

Used by

Dependency tree · two levels

41 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources