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.

Strong RRQR block singular values compare to the full matrix spectrum

Statement

Let A∈Fm×n have full column rank, m≥n, and let the exact-arithmetic version of Strong rank-revealing QR with finite swaps and explicit statuses return success with full QR factorization AΠ=Q[R11R120R22], where Q∈Fm×m is unitary and the displayed R is upper trapezoidal with m rows. Let the target satisfy 0≤k<n and the factor satisfy f≥1. For 1≤k<n, assume R11∈Fk×k is nonsingular and put X=R11−1R12. Put g=(1+f2k(n−k))1/2. Then σi(R11)≥σi(A)g(1≤i≤k),σj(R22)≤g σk+j(A)(1≤j≤n−k). For k=0, the second comparison is the equality R22=Q∗AΠ with g=1. These conclusions use the combined score in the algorithm; the bound on X alone is not used as a substitute for control of R22. For a computed candidate R^ with the algorithm's verified certificate (W,(bij),ε), put A~=WR^Π∗. No full-column-rank assumption on A~ is required beyond nonsingularity of R^11 when k>0. Then the corresponding residual-adjusted bounds are σi(R^11)≥σi(A)−εg(1≤i≤k),σj(R^22)≤g(σk+j(A)+ε)(1≤j≤n−k). For k=0, these reduce to σj(R^)≤σj(A)+ε. Computed scores alone do not certify either conclusion.

Facts & Assumptions

Given: A full-column-rank matrix and a successful factorization returned by the finite-swap algorithm.

[F1]

In the exact-arithmetic version of Strong rank-revealing QR with finite swaps and explicit statuses, success means ∣Xij∣2+(ωiγj)2≤f2 for every i,j, where ωi=∥ei∗R11−1∥2 and γj=∥R22(:,j)∥2. The finite-precision branch tests computed scores and is not used here.

[F2]

Singular values contract under left and right multiplication: σi(XY)≤σi(X)∥Y∥2 for compatible matrices.

[F3]

Spectral and Frobenius norms are unitarily invariant, are given by singular values, and satisfy the sharp rank comparison: the spectral and Frobenius norms are unitarily invariant, ∥M∥2=σ1(M), and ∥M∥2≤∥M∥F.

[F4]

Singular values are Lipschitz in the spectral norm bounds each singular-value change by the spectral norm of the perturbation.

Proof

1.1F1F2F3algebra

First, square unitary factors preserve every ordered singular value: [F3] gives ∥T∥2=∥T∗∥2=1 for a positive-dimensional unitary T, and applying both inequalities of [F2] to PMT gives σi(PMT)≤σi(M); applying them to M=P∗(PMT)T∗ gives the reverse inequality. Thus if k=0, R22=Q∗AΠ has the singular values of A, proving that case. Assume henceforth 1≤k<n. Write A0=R11, C=R22, and X=A0−1R12. Summing the k(n−k) inequalities in [F1] gives ∥X∥F2+∥A0−1∥F2∥C∥F2≤f2k(n−k). Indeed, ∑iωi2=∥A0−1∥F2 and ∑jγj2=∥C∥F2.

2.1step 1.1F1F3algebra

Since A has full column rank and k<n, C has full column rank: if Cz=0 for nonzero z, then R(−A0−1R12zz)=0, contradicting full column rank of R=Q∗AΠ, since Q and Π are invertible. Thus α=∥C∥2/σk(A0)>0. Since σk(A0)=1/∥A0−1∥2, one has α2=∥C∥22∥A0−1∥22≤∥C∥F2∥A0−1∥F2, so step 1.1 gives ∥X∥F2+α2≤f2k(n−k).

3.1step 1.1step 2.1F2F3algebra

By step 1.1, R=Q∗AΠ has the same singular values as A. Define D1=diag⁡(A0,C/α) and W1=[IkX0αIn−k], so R=D1W1. For vectors x,y, ∥W1(x,y)∥22=∥x+Xy∥22+α2∥y∥22 is bounded by the quadratic form of M1=[1∥X∥2∥X∥2∥X∥22+α2] at (∥x∥2,∥y∥2). For scalars u,v its quadratic form is (u+∥X∥2v)2+α2v2≥0, so it is positive semidefinite and its largest eigenvalue is at most its trace; thus ∥W1∥22≤1+∥X∥F2+α2≤g2. Since ∥C/α∥2=σk(A0), the ordered first k singular values of D1 are those of A0. By [F2], σi(R)≤σi(D1)∥W1∥2, proving σi(A0)≥σi(A)/g for i≤k.

3.2step 1.1step 2.1F2F3algebra

Define D2=diag⁡(αA0,C) and W2=[αIk−X0In−k], so D2=RW2. For vectors x,y, ∥W2(x,y)∥22=∥αx−Xy∥22+∥y∥22 is bounded by the quadratic form of M2=[α2α∥X∥2α∥X∥2∥X∥22+1] at (∥x∥2,∥y∥2). Its determinant is α2>0 and its diagonal entries are positive, so it is positive semidefinite; its largest eigenvalue is at most its trace. Thus ∥W2∥22≤1+∥X∥F2+α2≤g2. The smallest singular value of αA0 is ∥C∥2, so the first k singular values of D2 are those of αA0 and its remaining n−k singular values are those of C. Therefore σj(C)=σk+j(D2)≤σk+j(R)∥W2∥2≤gσk+j(A) by [F2] and step 1.1.

4.1step 1.1step 2.1step 3.1step 3.2F2F3F4algebra∎

The preceding block proof also applies to any upper trapezoidal R whose R11 is nonsingular and whose exact combined scores are at most f, even when R is rank deficient. If C≠0, α>0 and steps 2.1–3.2 apply without the full-column-rank argument. If C=0, write R=diag⁡(A0,0)[IX0I]. Directly R=A0[I  X] with zero lower rows and ∥[I  X]∥22=1+∥X∥22≤g2. The product bound proves the first comparison, and the second is 0≤gσk+j(R). Thus the certificate gives both comparisons with A replaced by A~ and blocks R^11,R^22. Since ∥A−A~∥2≤ε, [F4] gives σi(A~)≥σi(A)−ε and σk+j(A~)≤σk+j(A)+ε. Combining proves the computed bounds. The k=0 case follows directly from unitary invariance and [F4].

Depends on

Used by

Dependency tree · two levels

15 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