Lemma 2.1: dual sparse reconstruction property, version
ProvedCandesTao.Decoding.dual_reconstruction_l2Let be a real matrix with columns spanning . Let be such that , and let be a real vector supported on with . Then there exists a vector such that
Furthermore, there is an exceptional set , disjoint from , of size , with the properties
In addition, for some constant depending only on .
This is the first half of the dual certificate construction: a vector that interpolates prescribed values on and whose inner products with the columns outside are small in an sense, uniformly small outside an exceptional set of controlled size. Lemma 2.2 iterates it to remove the exceptional set.
Formalization Note The paper prints the second bound with ; the inequality (2.3) established in its proof, and the use of the lemma in Lemma 2.2 with , give , which is what is stated here (the two agree when ). "A constant depending only on " is formalized as a positive function of the real number , chosen before , , , and . The hypothesis is the domain on which Definition 1.1 defines .
import Definitions.Def_CandesTao_Decoding_RestrictedIsometry
namespace CandesTao.Decoding
theorem dual_reconstruction_l2 :
∃ K : ℝ → ℝ, (∀ δ : ℝ, 0 < K δ) ∧
∀ {p m : ℕ} (F : Matrix (Fin p) (Fin m) ℝ) (S S' : ℕ),
1 ≤ S → 1 ≤ S' → S + S' ≤ m → restrictedIsometryConst F S < 1 →
∀ (T : Finset (Fin m)) (c : Fin m → ℝ), T.card ≤ S → SupportedOn c T →
∃ w : Fin p → ℝ, w ∈ columnSpan F ∧
(∀ j ∈ T, dotProduct w (column F j) = c j) ∧
(∃ E : Finset (Fin m), Disjoint E T ∧ E.card ≤ S' ∧
(∀ j, j ∉ T → j ∉ E →
|dotProduct w (column F j)| ≤
restrictedOrthogonalityConst F S S' /
((1 - restrictedIsometryConst F S) * Real.sqrt S') * l2Norm c) ∧
Real.sqrt (∑ j ∈ E, dotProduct w (column F j) ^ 2) ≤
restrictedOrthogonalityConst F S S' / (1 - restrictedIsometryConst F S) * l2Norm c) ∧
l2Norm w ≤ K (restrictedIsometryConst F S) * l2Norm c := by sorry
end CandesTao.DecodingRead-back
What the Lean code literally says, in plain math · claude-fable-5-1
Read-back of dual_reconstruction_l2
Standing notation. Throughout, are natural numbers (any values, including ); is a real matrix whose rows are indexed by a set of labels and whose columns are indexed by a set of labels; is the -th column of , so ; is the matrix–vector product ; is the ordinary dot product; and
is the Euclidean norm (the square root is applied to a non-negative number, so it is the usual one; Mathlib's real square root would return on a negative input, which cannot happen here). For a finite set and , " is supported on " means for every ; nothing is required of on , so may also vanish on part or all of .
Definitions used, unfolded.
Restricted isometry constant. For ,
Restricted orthogonality constant. For ,
In both definitions the sets range over all subsets of the stated sizes, including the empty set and sets strictly smaller than or , and the first size parameter of goes with , the second with .
About the infima. Mathlib's infimum of a set of real numbers returns when the set is empty or has no lower bound. Neither situation occurs here: both sets are contained in , and both contain every sufficiently large real number (by Cauchy–Schwarz, every lies in the first set and every lies in the second). So and are genuine infima and are . Moreover each set is upward closed and defined by non-strict inequalities, hence is a closed ray ; the infimum is therefore attained, i.e. itself satisfies the two-sided inequality for every admissible , and itself satisfies the orthogonality inequality for every admissible .
Column span. , i.e. the set of vectors of the form with .
The statement.
There exists a function such that
-
for every real number (including and , values at which is never evaluated below); and
-
for all natural numbers , every matrix , and all natural numbers satisfying
and for every finite set and every vector with
there exists a vector for which all four of the following hold:
(a) , i.e. for some ;
(b) for every ;
(c) there exists a finite set such that
- ,
- ,
- for every index with and ,
- and
(d) .
Here is the -th coordinate of , and is the real square root of the natural number after casting it to a real; since , this is the ordinary positive square root and is .
Order of quantifiers and dependencies. is chosen first, before : a single function must serve every dimension, every matrix and every pair , and it enters the conclusion only through its value at the number , which under the hypotheses lies in . Nothing beyond positivity is asked of : no monotonicity, continuity, or boundedness as , and no relation to the explicit constants in (c). The vector may depend on all of ; the set may depend on all of these and on . Only existence is asserted for and for , not uniqueness.
Edge cases and conventions.
- The hypotheses , , force ; is not constrained by any explicit hypothesis. If (no rows) then for every , and with , this gives ; so the hypothesis excludes , and the statement is vacuous there.
- Because and , the denominators and are strictly positive, so Lean's convention is never triggered in the two displayed bounds. Both right-hand sides in (c) are since .
- , not , and may be empty. If , or more generally if , then every right-hand side in (c) and (d) equals ; (d) then forces , and does satisfy (a)–(c) (e.g. with ), so in that case the conclusion holds trivially.
- , not , and may be empty, in which case the pointwise bound in (c) applies to every and the sum over is .
- Indices are pinned exactly by (b). Indices are exempt from the pointwise bound and are controlled only through the aggregate bound over ; no bound is stated for any individual with . Indices outside get only the pointwise bound.
- The constant in (c) is with first parameter (matching ) and second parameter (matching ); no hypothesis is placed on .
- No assumption on appears beyond : no normalisation of columns, no rank condition, no relation between and .
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.