Alphabeta Math
Pipeline-generated
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.

✓ 13 results · all verified · 11 also independently AI-judged
Every result on this page is machine-checked by a proof checker and read in full and owner-audited; the judge is an additional, independent cross-model AI review of the proofs. The 2 not AI-judged were verified by owner audit (typically over a confirmed judge false positive), not failures.

Coordinate and block descent methods

1 · Prerequisites

2 · Summary

This page defines randomized coordinate-gradient and proximal updates, cyclic exact block minimization, and the Lasso coordinate formula. It proves finite separable block penalties are proper, closed, and convex by a local neighborhood argument that does not invoke sequential closure. Its smooth-gradient branch proves a finite-horizon expected O(nLmax⁡R2/k) bound and importance-sampling constants; the proximal branch proves well-defined updates and error recurrences both under strong convexity of the composite and, with a sharper factor, of the smooth term. It gives no penalty-independent R-only rate. Cyclic convergence is established only for convex separable composites with compact sublevels and unique continuous block-minimizer maps. A finite closed nonconvex domain shows why arbitrary tie selections need a cycling caveat. The canonical Lasso objective, minimizer existence theorem, and orthogonal-design soft-threshold example remain on the sparse-estimation page.

3 · Logical flowchart

4 · Definitions, theorems and proofs

AlgorithmAlgorithm: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Randomized coordinate gradient descent

Definition

Whenever a subscript i indexes a block, xi is shorthand for xBi.

Let d,n≥1, let B1,…,Bn be a partition of {1,…,d} into nonempty blocks, and let Ui:R∣Bi∣→Rd insert a block vector in coordinates Bi and put zero in the other coordinates. Coordinate label j∈{1,…,d} here denotes standard coordinate j−1∈{0,…,d−1}. Let f:Rd→R be convex and continuously differentiable in the sense of Continuously differentiable maps, local inverses, and local diffeomorphisms. Its gradient ∇f is the one defined by The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case; total differentiability supplies its partial derivatives by A total derivative computes every directional derivative, and its matrix is the Jacobian. Define ∇if(x) to be the entries of ∇f(x) with labels in Bi, in increasing order. For each i, an Li block-gradient Lipschitz bound is a number Li>0 such that ∥∇if(x+Uih)−∇if(x)∥2≤Li∥h∥2(x∈Rd, h∈R∣Bi∣), Choose majorizers Mi≥Li, probabilities pi>0 with ∑ipi=1, an initial point x0, and a finite update budget K∈N.

At each update k<K, draw Ik∈{1,…,n} with probabilities (pi), independently of earlier draws, and set xBIkk+1=xBIkk−1MIk∇Ikf(xk),xBjk+1=xBjk(j≠Ik). Uniform sampling means pi=1/n. The exact coordinate update is the minimizer of its quadratic model. If the caller supplies a stopping predicate, it is checked after each accepted update; on the first true result the procedure returns the current iterate with status tolerance-met. Otherwise it returns xK with status budget-limit. The latter status does not certify a specified finite objective gap. Every finite horizon uses only a finite product distribution on index histories.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Coordinate-block smoothness descent model

Statement

Let f:Rd→R be continuously differentiable, let Bi be one block of a finite coordinate partition, and let Ui be its coordinate insertion map. Suppose Li>0 satisfies ∥∇if(x+Uih)−∇if(x)∥2≤Li∥h∥2 for every x,h. Then f(x+Uih)≤f(x)+⟨∇if(x),h⟩+Li2∥h∥22.

Facts & Assumptions

Given: The function, block and Lipschitz bound in the Statement.

[F1]

The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case and Laws of finite sums and finite products: identify block label j∈{1,…,d} with standard coordinate j−1 and order each block increasingly; then ∇if(x) is the corresponding restriction of ∇f(x) and ⟨∇f(x),Uih⟩=⟨∇if(x),h⟩ by the finite-coordinate definition of the Euclidean inner product and finite-sum splitting.

[F2]

Continuously differentiable maps, local inverses, and local diffeomorphisms: continuously differentiable means totally differentiable at every point, with continuous derivative entries.

[F3]

The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a): along a differentiable curve, the derivative of the composition is the derivative of the outer function applied to the curve derivative.

[F4]

For a differentiable scalar field, Dvf(a)=⟨∇f(a),v⟩ and the unit direction of steepest ascent is the normalized gradient: for a totally differentiable scalar function, the directional derivative in direction v is ⟨∇f,v⟩.

[F6]

The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a): if a real function is continuous on a closed interval and differentiable in its interior, its endpoint difference is its derivative at an interior point times the interval length.

[F8]

The principle of mathematical induction: an identity that holds at zero and is preserved at successors holds for every natural number.

[F9]

Finite sums and finite products, by recursion and Laws of finite sums and finite products: finite sums satisfy the zero base case, successor recursion, additivity, telescoping, scaling, and the constant-sum law.

[F10]

For every ε>0 in a complete ordered field there is a natural n≥1 with 1/n<ε: for every η>0 there is an integer N≥1 with 1/N<η.

[F11]

Every natural-number-indexed list of nonempty sets has a choice function on its family of values: a finite indexed list of nonempty witness sets has a choice function in ZF.

Proof

1.1givenF1F2F3F4F5F7

If h=0 the claim is equality. Otherwise set ϕ(t)=f(x+tUih) on [0,1]. By [F2] and [F3], ϕ is totally differentiable at every point of [0,1], so [F5] makes it continuous there; [F3] and [F4] give ϕ′(t)=⟨∇if(x+tUih),h⟩. Hence the Lipschitz hypothesis and [F7] give ∣ϕ′(t)−ϕ′(0)∣≤Lit∥h∥22.

2.1F6F8F9F11step 1.1algebra

Fix N≥1 and apply [F6] on each interval [j/N,(j+1)/N]. Each mean-value witness set is nonempty, so [F11] selects cj∈(j/N,(j+1)/N) simultaneously for the finite list j<N, with ϕ((j+1)/N)−ϕ(j/N)=ϕ′(cj)/N. By [F9], summing and telescoping gives ϕ(1)−ϕ(0)−ϕ′(0)=1N∑j<N(ϕ′(cj)−ϕ′(0))≤Li∥h∥22N2∑j<N(j+1). To evaluate the last sum, set TN=∑j<N(j+1); [F9] gives T0=0 and TN+1=TN+(N+1), and [F8] proves TN=N(N+1)/2 by the base case and the algebraic successor identity N(N+1)/2+(N+1)=(N+1)(N+2)/2. Thus the upper bound is Li2(1+1N)∥h∥22.

3.1F10step 2.1assume-contraalgebradischarge-contradiction∎

The left side of step 2.1 is independent of N. If it exceeded Li∥h∥22/2 by δ>0 and Li∥h∥22>0, choose N with 1/N<2δ/(Li∥h∥22) by [F10]; step 2.1 would then give a strict smaller upper bound, a contradiction. If Li∥h∥22=0, step 2.1 already gives the result. Substituting ϕ(1)−ϕ(0)=f(x+Uih)−f(x) and ϕ′(0)=⟨∇if(x),h⟩ proves the inequality.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

The one-coordinate absolute-value quadratic minimizer

Statement

Let a>0, b∈R, and τ>0. The scalar convex function q(u)=a2u2−bu+τ∣u∣(u∈R) has the unique global minimizer u∗=Sτ(b)/a, where Sτ is the soft-threshold map of Soft-thresholding map. At b=τ and b=−τ, the unique minimizer is 0.

Facts & Assumptions

Given: The parameters in the Statement.

[F1]

Soft-thresholding map: Sτ(b)=b−τ for b>τ, 0 for −τ≤b≤τ, and b+τ for b<−τ.

[F2]

The triangle inequality and Basic properties of the absolute value: absolute value is subadditive and positively homogeneous.

Proof

1.1givenF2algebra

For 0≤s≤1, the square identity ((1−s)u+sv)2=(1−s)u2+sv2−s(1−s)(u−v)2 and a>0 make the quadratic term convex; [F2] makes the absolute-value term convex, and the linear term is affine. Thus q is convex. On u≥0, completing the square gives q(u)=a2(u−(b−τ)/a)2−(b−τ)22a, so the unique minimizer on that half-line is u+=max⁡{(b−τ)/a,0}.

1.2givenalgebra

On u≤0, q(u)=a2(u−(b+τ)/a)2−(b+τ)22a, so the unique minimizer on that half-line is u−=min⁡{(b+τ)/a,0}.

2.1F1step 1.1step 1.2algebra∎

If b>τ, then u+=(b−τ)/a>0 and q(u+)=−(b−τ)22a<0=q(0), while the negative-half-line minimum is 0; hence this is the unique global minimizer. If b<−τ, then u−=(b+τ)/a<0 and q(u−)=−(b+τ)22a<0=q(0), while the positive-half-line minimum is 0; hence this is the unique global minimizer. If −τ≤b≤τ, then u+=u−=0, the unique minimizer on each half-line by steps 1.1–1.2. These cases equal Sτ(b)/a by [F1], including the two kink values.

Remarks

For the finite sanity case b=3 and τ=1, the formula gives S1(3)=2 and the minimizer of 12u2−3u+∣u∣ is u=2.

AlgorithmAlgorithm: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Cyclic exact block minimization

Definition

Let d,n≥1, let B1,…,Bn partition the coordinates of Rd into nonempty blocks, let F:Rd→(−∞,+∞], and let x0 satisfy F(x0)<+∞. A cyclic exact block-minimization run with sweep budget K∈N processes blocks in the fixed order 1,…,n for at most K sweeps. At the update of block i, hold all other coordinates at their current values and form Mi(x−i)=arg⁡min⁡u∈R∣Bi∣F(x−i,u). If Mi(x−i)=∅, return the last iterate and status block-minimum-failure. Otherwise replace the current block by any member of Mi(x−i). If a minimizer is nonunique, the procedure is relation-valued and every member is an admissible choice; the convergence theorem on this page assumes unique continuous block-minimizer maps. After K complete sweeps, return the current iterate and status budget-limit. The finite budget does not certify a stationary point.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Rates for randomized coordinate gradient descent and importance sampling

Statement

Let f:Rd→R be convex and continuously differentiable, with a partition into n≥1 nonempty blocks and block-gradient Lipschitz constants Li>0 as in Coordinate-block smoothness descent model. Let Mi≥Li, choose probabilities pi>0 with ∑ipi=1, and set ρ=min⁡ipi/Mi. Assume f∗=min⁡xf(x) is attained at x∗ and the initial sublevel set has finite radius R=sup⁡{∥x−x∗∥2:f(x)≤f(x0)}<+∞. Run Randomized coordinate gradient descent from the initial point x0 with exact updates, no stopping predicate, and finite budget K∈N. For every integer k with 1≤k≤K, E[f(xk)−f∗]≤2R2ρk. If, in addition, f is globally μ-strongly convex in the sense that f(y)≥f(x)+⟨∇f(x),y−x⟩+μ2∥y−x∥22(x,y∈Rd) for some μ>0, then E[f(xk)−f∗]≤(1−μρ)k(f(x0)−f∗). The strong-convexity factor is in [0,1]. Uniform sampling with Mi=Lmax⁡:=max⁡iLi gives the factor 1−μ/(nLmax⁡) and the bound 2nLmax⁡R2/k. For nonuniform sampling, choosing Mi=Li and pi=Li/∑jLj maximizes min⁡ipi/Mi, giving ρ=1/∑iLi; its constants are not the uniform constants.

Facts & Assumptions

Given: The objective, block partition, smoothness constants, sampling probabilities, and finite-history expectation in the Statement.

[F1]

Randomized coordinate gradient descent: at a selected block i, the update is xi+=xi−Mi−1∇if(x) and the block draw has probability pi.

[F2]

Coordinate-block smoothness descent model: f(x+Uih)≤f(x)+⟨∇if(x),h⟩+(Li/2)∥h∥22.

[F3]

The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn: the product Euclidean norm satisfies ∥x∥22=∑i∥xBi∥22.

[F5]

Differentiable convex functions are characterized by the gradient inequality: convex differentiability gives f(y)≥f(x)+⟨∇f(x),y−x⟩.

[F6]

Fermat's theorem: an interior differentiable local extremum has zero gradient: an unconstrained differentiable local minimizer has zero gradient.

[F7]

Laws of finite sums and finite products: finite weighted sums preserve nonnegativity and obey the ordinary algebraic sum laws.

[F8]

The principle of mathematical induction: a base case and successor implication prove a claim for every natural number.

[A1]

f is convex, continuously differentiable, and satisfies the block-gradient Lipschitz hypotheses stated above.

[A2]

The initial sublevel radius R is finite, the algorithm has no stopping predicate and k≤K, and every expectation at a fixed k is the finite weighted sum over the k-step index histories, so no infinite product probability space is needed.

[A3]

For the linear rate only, f satisfies the displayed global μ-strong-convexity inequality and has a minimizer.

Proof

1.1A1F1F2algebra

Fix x and a block i. With h=−∇if(x)/Mi, [F2] gives f(x+Uih)≤f(x)−1Mi∥∇if(x)∥22+Li2Mi2∥∇if(x)∥22≤f(x)−12Mi∥∇if(x)∥22, because Mi≥Li.

1.2A3F2algebra

Under [A3], strong convexity at x and x∗ gives, for e=x−x∗, f(x)−f∗≤⟨∇f(x),e⟩−(μ/2)∥e∥22≤∥∇f(x)∥22/(2μ); the last inequality follows by expanding 0≤∥∇f(x)−μe∥22. Hence ∥∇f(x)∥22≥2μδ(x). Also μ≤Li≤Mi for each block, since the strong-convexity lower model and [F2]'s upper model applied to any nonzero displacement in block i bound the same function difference.

2.1A2F1F7step 1.1algebra

Conditional on the first k draws, the next draw has probabilities pi. By the finite-sum laws [F7], averaging step 1.1 over those finitely many outcomes yields E[f(xk+1)∣xk]≤f(xk)−12∑ipiMi∥∇if(xk)∥22≤f(xk)−ρ2∥∇f(xk)∥22. Taking the finite-history expectation gives the same inequality with all terms averaged.

3.1A1A2F3F4F5F7F8step 2.1algebra

Put δ(x)=f(x)−f∗. If R=0, then x0=x∗ because x0 belongs to the initial sublevel set, and [F6] gives ∇f(x∗)=0; all later updates leave x∗ fixed and all gaps vanish. Now suppose R>0. By [F5], δ(x)≤⟨∇f(x),x−x∗⟩. Since every exact step decreases f, each iterate lies in the initial sublevel set; thus [F4] and the definition of R give δ(xk)≤R∥∇f(xk)∥2. For each finite k-step history ω, write wω for its probability and zω=∥∇f(xωk)∥2; the weights are nonnegative and sum to 1. With zˉ=∑ωwωzω, [F7] gives 0≤∑ωwω(zω−zˉ)2=∑ωwωzω2−zˉ2. Hence E∥∇f(xk)∥22≥(E∥∇f(xk)∥2)2≥ϕk2/R2. Step 2.1 therefore implies ϕk+1≤ϕk−ρϕk2/(2R2).

4.1F6F7step 3.1algebra

If some ϕj=0, then f(xj)=f∗ for every history of positive probability, [F6] gives zero gradient at each such unconstrained minimizer, and all later gaps are zero. Otherwise all ϕj>0, and ϕj+1≤ϕj gives 1ϕj+1−1ϕj=ϕj−ϕj+1ϕjϕj+1≥ϕj−ϕj+1ϕj2≥ρ2R2. Summing for j=0,…,k−1 telescopes by [F7] to 1/ϕk≥kρ/(2R2), which is the stated O(1/k) bound.

5.1A2A3F7F8step 2.1step 1.2algebra∎

Substituting step 1.2 in step 2.1 gives E[δ(xk+1)∣xk]≤(1−μρ)δ(xk). Because 0<μ≤Mi and pi≤1, 0≤μρ≤1. Taking expectations and using induction [F8] proves the linear bound. For any probabilities, [F7] gives 1=∑ipi≥ρ∑iMi, so ρ≤1/∑iMi; equality holds for pi=Mi/∑jMj. Since Mi≥Li, this is at most 1/∑iLi, attained by Mi=Li and pi=Li/∑jLj. The uniform and importance-sampling specializations follow by substitution.

LemmaStatement: AI-adaptedProof: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

Continuity of the cyclic block-minimum value

Statement

Fix a block i and define mi(a)=min⁡u{f(a,u)+gi(u)}, where a denotes all other blocks, f is finite and continuous, and gi is proper closed convex. Suppose the minimizer map Ti(a) exists, is unique, and is continuous on the parameter set in question. Then mi is continuous there.

Facts & Assumptions

Given: The composite block problem and minimizer map in the Statement.

[F1]

Cyclic exact block minimization: Ti(a) is the exact global minimizer of the ith block slice.

[F2]

Proper closed convex extended-real functions on Rn: closed epigraph in finite dimensions is the lower-semicontinuity convention for gi.

[A1]

f is continuous and Ti is single-valued and continuous on the parameter set.

Proof

1.1A1F2givenalgebra

Fix a parameter a, let u=Ti(a), and put mi(a)=f(a,u)+gi(u). For ε>0, set c=gi(u)−ε/3. The weak sublevel Cc={v:gi(v)≤c} is closed because v↦(v,c) is a continuous slice of the closed epigraph in [F2]; since u∉Cc, its open complement contains a neighborhood of u on which gi>c. Continuity of Ti and f gives a relative neighborhood V of a on which Ti(a′) lies in that neighborhood and both f(a′,Ti(a′)) and f(a′,u) differ from f(a,u) by less than ε/3. Hence mi(a′)>mi(a)−2ε/3.

2.1A1F1step 1.1algebra∎

For each a′∈V, the point u is a competitor, so mi(a′)≤f(a′,u)+gi(u)<mi(a)+ε/3. Together with step 1.1, this gives ∣mi(a′)−mi(a)∣<ε, which is continuity of mi at the arbitrary parameter a. The argument uses neighborhoods and closed epigraphs directly; it does not infer topological continuity from a sequential criterion.

CounterexampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

Cyclic exact block minimization can cycle outside the convex separable hypotheses

Statement refuted

The claim that every infinite continuation of the cyclic exact block-minimization update rule for an arbitrary nonconvex, nonseparable extended-real objective converges to one point is false when block minimizers are nonunique and the allowed tie choices vary by update.

Witness

Let C={(0,0),(2,0),(2,2),(1,2),(1,1),(0,1)}⊂R2 and let F be its indicator: F=0 on C and F=+∞ off C. This is a finite, closed but nonconvex, nonseparable domain. Use singleton blocks in the order x,y and start at (0,0). Any finite sweep budget follows the corresponding prefix of the periodic tie-selection rule; three sweeps complete one six-update period.

Facts & Assumptions

Given: The objective, blocks and starting point in the Counterexample.

[F1]

Cyclic exact block minimization: at each update the selected coordinate is replaced by any global minimizer of the current coordinate slice.

[F2]

The extended real line R‾=R∪{−∞,+∞}, its order, and the arithmetic that is left undefined: the indicator values 0 and +∞ give the extended-real objective and its finite domain.

Proof

1.1F1F2givenchoose

At y=0, the feasible x values are 0 and 2, and both minimize the slice with value 0; choose x=2. At x=2, the feasible y values are 0 and 2, and both minimize; choose y=2.

2.1F1F2step 1.1choose

At y=2, choose x=1 from the two minimizing values 1,2. At x=1, choose y=1 from the two minimizing values 1,2. At y=1, choose x=0 from the two minimizing values 0,1. At x=0, choose y=0 from the two minimizing values 0,1. Thus one six-update cycle returns to (0,0), and repeating the same admissible tie choices gives a period-six trajectory.

3.1F1F2F3step 1.1step 2.1algebra

Every point in the trajectory lies in C and has objective 0, so every coordinate update described in steps 1.1 and 2.1 is an exact global block minimization. The set C is closed: for any z∉C, the finitely many distances {∥z−c∥2:c∈C} are positive, so [F3] gives a positive minimum δ; if ∥w−z∥2<δ/2, the triangle inequality gives ∥w−c∥2≥∥z−c∥2−∥w−z∥2>δ/2 for every c∈C. Thus the open ball of radius δ/2 about z misses C. The set is nonconvex because (0,0),(2,0)∈C but their midpoint (1,0)∉C. It is not a Cartesian product of its coordinate projections: both projections are {0,1,2}, whose product has nine points, whereas C has six. If F(x,y)=h(x)+k(y) for functions h,k taking values in R∪{+∞}, its finite domain consists exactly of pairs where h and k are each finite, hence is a Cartesian product. Thus F is nonseparable.

4.1F1step 1.1step 2.1step 3.1algebra∎

In the infinite continuation of the finite-budget prefixes, the trajectory does not converge: the subsequences at updates 6m and 6m+1 are constantly (0,0) and (2,0), respectively, and those limits differ. Repeating the six admissible choices therefore gives a nonconvergent cyclic exact-minimization trajectory.

Remarks

The example has a nonconvex, nonseparable objective and nonunique block minimizers. It lies outside Cyclic exact block descent has globally optimal accumulation points; the convergence theorem does not assert convergence for this case.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Exact Lasso coordinate update, including a zero column

Statement

Use the Lasso objective Qy and notation of Lasso estimator for fixed design, with n,d≥1 and τ>0. Fix a coordinate j and all other coefficients, and set rj=y−∑ℓ≠jXℓθℓ. If Xj≠0, the unique minimizing coordinate is θj+=Snτ(XjTrj)∥Xj∥22. If Xj=0, the unique minimizing coordinate is θj+=0.

Facts & Assumptions

Given: The fixed-design Lasso objective, a coordinate j, and the other coordinates.

[F1]

Lasso estimator for fixed design: Qy(θ)=n−1∥y−Xθ∥22+2τ∥θ∥1, with τ>0.

[F2]

The one-coordinate absolute-value quadratic minimizer: for a>0 and λ>0, a2u2−bu+λ∣u∣ has unique minimizer Sλ(b)/a.

[F3]

Soft-thresholding map: Sλ(b)=b−λ for b>λ, equals 0 for ∣b∣≤λ, and equals b+λ for b<−λ.

[A1]

The residual excluding coordinate j is rj=y−∑ℓ≠jXℓθℓ.

Proof

1.1F1F4A1algebra

For a candidate value u in coordinate j, y−Xθ=rj−Xju. By bilinearity in [F4], expanding the square in [F1] shows that the terms depending on u are n−1∥Xj∥22u2−2n−1(XjTrj)u+2τ∣u∣.

2.1F2F3F4step 1.1algebra

If Xj≠0, [F4] gives ∥Xj∥22>0. Multiply these terms by the positive constant n/2. The result is 12∥Xj∥22u2−(XjTrj)u+nτ∣u∣. Apply [F2] with a=∥Xj∥22, b=XjTrj, and λ=nτ>0 to obtain the stated unique minimizer; [F3] fixes its threshold convention at both kink values.

3.1F1A1givenalgebra∎

If Xj=0, the squared-residual term in [F1] is independent of u, while the only varying term is 2τ∣u∣. Since τ>0, this term has the unique minimum at u=0.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Finite separable block penalties are proper, closed, and convex

Statement

Let n≥1 be finite, let di≥1, and let each gi:Rdi→(−∞,+∞] be proper, closed, and convex in the epigraph convention of Proper closed convex extended-real functions on Rn. On the product Rd1×⋯×Rdn define g(x)=∑i=1ngi(xi). Then g is proper, closed, and convex. The proof below uses only finite construction and the closed epigraphs of the separate terms.

Facts & Assumptions

Given: A finite nonempty block list and the functions in the Statement.

[F1]

Proper closed convex extended-real functions on Rn: each gi has nonempty effective domain, satisfies the extended-real convexity inequality, and has closed epigraph.

[F2]

The principle of mathematical induction: a property true at the base index and preserved at each successor holds for every natural index in the finite block list.

[A1]

Each block dimension di is finite and positive, and the index set has size n≥1.

Proof

1.1A1F1F2construct

For each i, [F1] gives a point in dom⁡gi. The empty tuple is a choice for the first zero blocks; if a tuple has been constructed for the first k<n blocks, append one point from the nonempty dom⁡gk+1. Induction [F2] gives a tuple x in all n domains, so g(x)<+∞ and g is proper. Only finitely many choices are made.

1.2F1givenalgebra

If g(x),g(y)<+∞ and 0≤θ≤1, apply the convexity inequality in [F1] to each pair (xi,yi) and add the finitely many inequalities. If either endpoint has infinite g value, the extended-real inequality is automatic under the endpoint convention. Thus g is convex.

1.3F1givenalgebra

Fix a block i and a real c. Its weak sublevel Ci,c={ui:gi(ui)≤c} is closed: the map ui↦(ui,c) preserves Euclidean distances and is continuous, so Ci,c is its preimage of the closed set epi⁡gi from [F1]. Therefore, whenever gi(xi)>c, the open complement of Ci,c is a neighborhood of xi throughout which gi>c.

2.1A1step 1.3givenalgebra

Take any (x,t)∉epi⁡g. If all gi(xi) are finite, put s=∑igi(xi)>t and ci=gi(xi)−(s−t)/(2n). By step 1.3, each xi has an open neighborhood on which gi>ci, so their product is a neighborhood of x on which g(u)>∑ici=(s+t)/2. With η=(s−t)/4>0, the product of these neighborhoods and (t−η,t+η) misses epi⁡g, because its height coordinate is below t+η<(s+t)/2.

2.2A1step 1.3givenalgebra

If some gi(xi)=+∞, choose such an index i. For each j≠i, put cj=gj(xj)−1 when that value is finite, and cj=0 when it is infinite. Choose B=t+1−∑j≠icj. Since gi(xi)>B and gj(xj)>cj, step 1.3 supplies an open neighborhood of each block on which the corresponding strict lower bound holds. Throughout their product, g(u)>B+∑j≠icj=t+1. Restrict the epigraph-height coordinate to (t−1/2,t+1/2); this neighborhood also misses epi⁡g.

3.1step 1.1step 1.2step 2.1step 2.2algebra∎

Every point outside epi⁡g is covered by one of steps 2.1 or 2.2, so its complement is open and epi⁡g is closed. Together with properness from step 1.1 and convexity from step 1.2, this proves the Statement without converting an arbitrary sequential-closure claim into topological closedness.

AlgorithmAlgorithm: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Randomized proximal coordinate descent with certified inner solves

Definition

Whenever a subscript i indexes a block, xi is shorthand for xBi.

Let F(x)=f(x)+∑i=1ngi(xBi) on Rd, where B1,…,Bn is a disjoint nonempty partition of coordinate labels {1,…,d}, with label j denoting standard coordinate j−1. Let f:Rd→R be convex and continuously differentiable in the sense of Continuously differentiable maps, local inverses, and local diffeomorphisms, and let each gi:R∣Bi∣→(−∞,+∞] be proper, closed and convex. Define ∇if(x) from The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case as the entries of ∇f(x) on the labels in Bi, in increasing order; their existence follows from total differentiability and A total derivative computes every directional derivative, and its matrix is the Jacobian. Let Li>0 satisfy the block-gradient Lipschitz bound of Coordinate-block smoothness descent model, choose majorizers Mi≥Li, probabilities pi>0 with ∑ipi=1, a feasible initial point x0, an update budget K∈N, and nonnegative inner tolerances ϵ0,…,ϵK−1.

At update k, draw Ik with probabilities (pi) independently of previous draws. For the selected block i=Ik, form qik(u)=⟨∇if(xk),u−xBik⟩+Mi2∥u−xBik∥22+gi(u). In exact mode, set xBik+1 to the global minimizer of qik; leave every other block unchanged. In certified-inexact mode, an inner oracle returns (uk,δk) with a certificate 0≤qik(uk)−inf⁡uqik(u)≤δk≤ϵk; accept uk only when that certificate is available and valid. Otherwise return the last accepted point with status inner-solve-failure. The output is the current point, its objective value, the last certified inner gap when available, and status tolerance-met, budget-limit, or inner-solve-failure.

If a proximal-gradient stopping tolerance η≥0 is requested, after each accepted update compute the block mapping GM(x)=(Mi[xBi−prox⁡gi/Mi(xBi−Mi−1∇if(x))])i=1n. Here prox⁡gi/Mi abbreviates prox⁡(1/Mi)gi, a singleton by Proximal mapping and Existence, resolvent characterization, and firm nonexpansiveness of the proximal map; Mi>0 by the input conditions. Stop with tolerance-met when ∥GM(x)∥2≤η; otherwise stop after at most K accepted updates with budget-limit. A finite budget or tolerance status is not an uncomputed objective-gap guarantee. Uniform sampling is pi=1/n; nonuniform choices require the separate probability/step relation proved for the smooth gradient branch on this page.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Coordinatewise minima of a separable convex composite are global minima

Statement

Let f:Rd→R be differentiable and convex, and let B1,…,Bn be a partition of the coordinate labels {1,…,d} into nonempty blocks, with n≥1. Let g(x)=∑i=1ngi(xBi) be proper, closed and convex, with each gi proper, closed and convex on its block. Set F=f+g. If F(x)<+∞ and x globally minimizes every block slice u↦F(x−i,u), where (x−i,u) means replace block Bi of x by u and keep all other blocks fixed, then −∇if(x)∈∂gi(xBi)(1≤i≤n),−∇f(x)∈∂g(x), and x is a global minimizer of F. For an extended-real convex h, p∈∂h(a) means h(z)≥h(a)+⟨p,z−a⟩ for every z.

Facts & Assumptions

Given: The finite block partition, the functions and the coordinatewise minimizing point in the Statement.

[F1]

Subdifferentials of proper extended-real convex functions: p∈∂h(a) exactly when h(z)≥h(a)+⟨p,z−a⟩ for every z in the full Euclidean space.

[F2]

Proper closed convex extended-real functions on Rn: each gi is convex and has nonempty effective domain; closed epigraph is lower semicontinuity in this convention.

[F3]

For a differentiable scalar field, Dvf(a)=⟨∇f(a),v⟩ and the unit direction of steepest ascent is the normalized gradient: for a totally differentiable scalar function, the directional derivative in direction v is ⟨∇f(x),v⟩.

[F4]

Differentiable convex functions are characterized by the gradient inequality: convex differentiability gives f(z)≥f(x)+⟨∇f(x),z−x⟩.

[F5]

Identify block label j∈{1,…,d} with standard coordinate j−1; let Ui insert a block displacement into those coordinates, and define ∇if(x) as the entries of ∇f(x) on Bi, in increasing order. By Laws of finite sums and finite products and The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, finite block sums add and the full inner product is the sum of the block inner products, so ⟨∇f(x),Uih⟩=⟨∇if(x),h⟩.

[A1]

f is differentiable and convex, g=∑igi is proper closed convex, F(x)<+∞, and x minimizes each exact one-block slice.

Proof

1.1A1F1F2F3F5givenalgebra

Fix a block i, a point zi with gi(zi)<+∞, and put h=zi−xBi. For 0<t<1, block minimality and convexity of gi give 0≤f(x+Ui(th))−f(x)+gi(xBi+th)−gi(xBi)≤f(x+Ui(th))−f(x)+t[gi(zi)−gi(xBi)]. Rearranging, dividing by t>0, and taking t↓0, [F3] and [F5] give gi(zi)−gi(xBi)≥−⟨∇f(x),Uih⟩=−⟨∇if(x),zi−xBi⟩. For gi(zi)=+∞ the same subgradient inequality is automatic. Thus −∇if(x)∈∂gi(xBi) by [F1].

2.1F1F5step 1.1algebra

Fix any z=(zBi)i. If g(z)=+∞, the subgradient inequality at z is automatic. Otherwise each gi(zBi) is finite, so summing the block inequalities from step 1.1 and using [F5] gives g(z)≥g(x)−∑i⟨∇if(x),zBi−xBi⟩=g(x)−⟨∇f(x),z−x⟩. Thus the subgradient inequality holds for every z, and [F1] gives −∇f(x)∈∂g(x).

3.1A1F4step 2.1algebra∎

By [F4], f(z)≥f(x)+⟨∇f(x),z−x⟩. Adding the inequality in step 2.1 gives F(z)≥F(x) for every z, so x is globally minimizing and all the stated stationarity inclusions hold.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Coordinate proximal subproblems are attained uniquely and descend

Statement

For the model and inputs of Randomized proximal coordinate descent with certified inner solves, every selected exact block subproblem has a unique minimizer. If the update is exact, the objective does not increase. If its accepted inner point has certified model gap at most ϵk, then the new objective is at most the old objective plus ϵk.

Facts & Assumptions

Given: A feasible iterate x, a selected block i, positive majorizer Mi, and proper closed convex block penalty gi.

[F1]

Randomized proximal coordinate descent with certified inner solves: the block model is qi(u)=⟨∇if(x),u−xi⟩+(Mi/2)∥u−xi∥22+gi(u), and an accepted inexact point satisfies qi(u)−inf⁡qi≤ϵk.

[F2]

Existence, resolvent characterization, and firm nonexpansiveness of the proximal map: for proper closed convex g, positive t, and any v, the proximal set prox⁡tg(v) is a singleton.

[F3]

Proximal mapping: prox⁡tg(v) is the set of global minimizers of g(u)+∥u−v∥22/(2t).

[F4]

Proper closed convex extended-real functions on Rn: each gi is proper, closed and convex on its Euclidean block space.

[F5]

Coordinate-block smoothness descent model: changing one block by h raises f by at most ⟨∇if(x),h⟩+(Li/2)∥h∥22.

Proof

1.1F2F3F4givenalgebra

Complete the square in qi: with vi=xi−Mi−1∇if(x), its u-dependent terms are gi(u)+(Mi/2)∥u−vi∥22. Taking t=1/Mi>0 in [F2]–[F4] proves the exact block minimizer exists and is unique.

2.1F1F5step 1.1algebra∎

For an exact update, qi(u∗)≤qi(xi)=gi(xi). For an accepted inexact update, [F1] gives qi(u)≤qi(u∗)+ϵk≤qi(xi)+ϵk. Apply [F5] with h=u−xi and cancel the unchanged penalties on the other blocks; the resulting objective is at most F(x) in the exact case and at most F(x)+ϵk in the inexact case.

CounterexampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

An R-only proximal coordinate rate fails for nonsmooth penalties

Statement refuted

The claim that the uniform randomized proximal update for every proper closed convex separable penalty satisfies E[F(xk)−F∗]≤2nLmax⁡R2k+4 under only a bounded initial sublevel radius R and positive coordinate majorizers Li is false without an error term or other dependence on the penalty.

Witness

Let n=2, f(x)=0, g1(u)=g2(u)=4∣u∣, L1=L2=Lmax⁡=1, and x0=(1,1). The objective is F(x)=4(∣x1∣+∣x2∣) and F∗=0 at x∗=(0,0). The initial sublevel is the closed diamond ∣x1∣+∣x2∣≤2, whose radius about x∗ is R=2.

Facts & Assumptions

[F1]

Randomized proximal coordinate descent with certified inner solves: for selected block i the exact update minimizes ⟨∇if(x),u−xi⟩+(Mi/2)∣u−xi∣2+gi(u).

[F2]

The one-coordinate absolute-value quadratic minimizer: a2u2−bu+τ∣u∣ has unique minimizer Sτ(b)/a when a,τ>0.

[F3]

Proper closed convex extended-real functions on Rn: a block penalty is proper, closed and convex when its effective domain is nonempty, its epigraph is closed, and it satisfies the segment inequality.

[F4]

Absolute value in an ordered field, Basic properties of the absolute value and The triangle inequality: absolute value is nonnegative, positively homogeneous, and subadditive.

[F5]

The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn and Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation: the two-coordinate Euclidean norm is a norm and its triangle inequality gives ∥(x1,x2)∥2≤∣x1∣+∣x2∣ by decomposing into coordinate vectors.

[A1]

The proximal coordinate algorithm uses common model majorizers M1=M2=1, samples uniformly, and runs at least one update with no stopping tolerance.

Proof

1.1A1F1F3F4F5givenalgebra

The assumptions hold: f=0 is convex and continuously differentiable and has zero gradient, so Li=Mi=1 is a valid positive block-gradient bound and majorizer. Each gi(u)=4∣u∣ is finite everywhere and hence proper; [F4] gives convexity; and epi⁡gi={(u,t):t≥4u, t≥−4u} is an intersection of closed half-planes, so it is closed in the epigraph convention [F3]. Thus F(x)=4(∣x1∣+∣x2∣) and its minimum is F∗=0 at x∗=(0,0). The initial value is F(x0)=8 and its sublevel is ∣x1∣+∣x2∣≤2. There ∥x−x∗∥2=x12+x22≤∣x1∣+∣x2∣≤2 by [F4]--[F5], with equality at (2,0); hence R=2.

2.1F1F2step 1.1algebra

For either selected coordinate the subproblem at its current value 1 is 4∣u∣+12(u−1)2, which differs by a constant from 12u2−u+4∣u∣. By [F2] its unique minimizer is S4(1)=0. The other coordinate remains 1, so after one update F(x1)=4 for either draw and E[F(x1)−F∗]=4.

3.1step 1.1step 2.1algebra∎

The proposed upper bound at k=1 is 2⋅2⋅1⋅22/(1+4)=16/5<4, contradicting step 2.1. Therefore the claimed penalty-independent composite bound is false.

Remarks

The controlling plan's O(1/k) rate is for the smooth coordinate-gradient branch gi=0. This witness does not refute that branch.

CounterexampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Strong convexity of the composite alone does not give the stated coordinate factor

Statement refuted

The claim that F=f+∑igi being μ-strongly convex by itself implies the contraction factor 1−μ/(nLmax⁡) for the randomized proximal coordinate update is false when the curvature comes only from the nonsmooth-separable term.

Witness

Take n=2, f=0, g1(u)=g2(u)=12u2, common positive majorizers L1=L2=Lmax⁡=1, and x0=(1,1). Then F(x)=12(x12+x22) is globally μ=1 strongly convex and has minimum F∗=0.

Facts & Assumptions

Given: The update rule of Randomized proximal coordinate descent with certified inner solves for the stated two-block objective.

[F1]

Randomized proximal coordinate descent with certified inner solves: in exact mode, with f=0 and Mi=1, the selected block minimizes gi(u)+12(u−xi)2.

[F2]

Proper closed convex extended-real functions on Rn: each block penalty must have nonempty finite domain, closed epigraph, and satisfy the convexity inequality.

[A1]

The run uses exact mode, the common model majorizers are M1=M2=1, the two coordinates are sampled uniformly, and the run executes at least one update with no stopping tolerance.

Proof

1.1A1F1F2givenalgebra

Each gi(u)=u2/2 is finite and its epigraph {(u,t):u2/2≤t} is closed; hence it is proper under [F2]. The identity gi((1−s)u+sv)=(1−s)gi(u)+sgi(v)−s(1−s)(u−v)2/2 proves convexity. The same identity summed over coordinates shows F(x)=12∥x∥22 is globally μ=1 strongly convex, with unique minimum at x∗=(0,0). Since f=0 has zero gradient, Li=1 are valid positive block-gradient bounds and Mi=Lmax⁡=1. For either selected coordinate the model is q(u)=12u2+12(u−1)2=(u−12)2+14, so its unique minimizer is u=1/2. The other coordinate remains 1. Thus every one-step outcome has objective gap 12((1/2)2+1)=5/8, and the expected gap is also 5/8.

2.1step 1.1givenalgebra∎

The claimed factor would give (1−μ/(nLmax⁡))(F(x0)−F∗)=(1−1/2)⋅1=1/2, which is strictly smaller than 5/8. Hence composite strong convexity alone does not imply that factor. The positive-rate theorem on this page assumes the smooth term f itself is μ-strongly convex.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Cyclic exact block descent has globally optimal accumulation points

Statement

Let F(x)=f(x)+g(x) on Rd, where f:Rd→R is convex and continuously differentiable, and g(x)=∑igi(xBi) is proper closed convex with each gi proper closed convex. Assume F(x0)<+∞ and the initial sublevel set S={x:F(x)≤F(x0)} is compact. Suppose every exact block minimizer on the visited sublevel exists uniquely and the resulting block-minimizer maps are continuous. Then every accumulation point of any indefinitely continued cyclic full-sweep run of Cyclic exact block minimization is coordinatewise minimizing, satisfies 0∈∇f(x)+∂g(x), and is a global minimizer of F. If the global minimizer is unique, the full sequence, including within-sweep iterates, converges to it. The theorem uses the unique-continuous-minimizer branch of the pair design; it makes no claim for arbitrary nonconvex or nonseparable nonsmooth updates.

Facts & Assumptions

Given: The objective, compact sublevel and block-minimizer hypotheses in the Statement.

[F1]

Cyclic exact block minimization: a sweep minimizes each block in order and leaves the other blocks unchanged.

[F2]

Continuity of the cyclic block-minimum value: with continuous unique block-minimizer maps, the value of each block minimum is continuous in the other coordinates.

[F3]

Coordinatewise minima of a separable convex composite are global minima: a coordinatewise minimum of this convex separable composite is stationary and globally minimizing.

[F4]

A finite continuous convex summand preserves proper closed convexity: adding the finite continuous convex function f to proper closed convex g preserves properness, convexity and closed epigraph.

[F5]

Semicontinuous extreme value theorem on compact Euclidean sets: a real-valued lower-semicontinuous function on a nonempty compact Euclidean set is bounded below and attains its minimum.

[F7]

Every nonempty set bounded below has an infimum and Epsilon characterisation of the infimum: a nonempty real set bounded below has an infimum, and elements lie arbitrarily close above it.

[F8]

A compact subset of a metric space is closed and bounded: a compact subset of a metric space is closed, so ambient limits of iterates in S remain in S.

[F10]

The well-ordering principle and The recursion theorem: a recursively defined subsequence can choose the least next natural index in any nonempty set of admissible indices.

[A1]

The maps for the exact block minimizers are single-valued and continuous on the states in S; every finite update is exact.

Proof

1.1A1F1F4F5F7F9givenalgebra

By [F1], F is nonincreasing at every block update, so all cycle and within-cycle iterates remain in S. By [F9], f is continuous, so [F4] makes F lower semicontinuous; it is real-valued on S because F≤F(x0)<+∞. Thus [F5] bounds it below on S. By [F7], L=inf⁡kF(xk) is real; for each ε>0, [F7] gives j with F(xj)<L+ε, and monotonicity then gives L≤F(xk)≤F(xj)<L+ε for all k≥j. Hence the full-sweep values tend to L.

2.1A1F1F8step 1.1given

Fix any accumulation point xˉ of the full-sweep iterates and a subsequence xkj→xˉ. Since S is compact, [F8] makes it closed and therefore xˉ∈S. The intermediate states yj0=xkj,yj1,…,yjn=xkj+1 satisfy F(xkj)≥F(yji)≥F(xkj+1). Both endpoint values tend to L, so every intermediate value tends to L. By [A1] and continuity of the finitely many block maps, each yji converges to a limit yˉi, with yˉ0=xˉ.

3.1A1F1F2step 2.1algebra

Fix i. Let aj be the coordinates other than block i in yji−1, and let uj− and uj+ be its block-i coordinates before and after the update. Then aj→aˉ, uj−→uˉ−, and uj+→uˉ+ by step 2.1; exactness gives uj+=Ti(aj). Write mi(a)=min⁡u{f(a,u)+gi(u)}. Since the other penalties cancel, the objective decrease is dj=f(aj,uj−)+gi(uj−)−mi(aj), and dj→0 because the two intermediate objective values both tend to L. By [F2], mi(aj)→mi(aˉ); continuity of f and lower semicontinuity of gi imply f(aˉ,uˉ−)+gi(uˉ−)≤mi(aˉ). The reverse inequality follows from the minimum definition, so u−=Ti(aˉ) by uniqueness. Continuity of Ti gives uˉ+=Ti(aˉ)=uˉ−, hence yˉi=yˉi−1.

4.1A1F3step 3.1algebra

Applying step 3.1 for i=1,…,n shows every within-sweep limit equals xˉ and that xˉ minimizes every block slice. Since xˉ∈S and F(x0)<+∞, F(xˉ)<+∞. By [F3], xˉ satisfies the composite stationarity condition and is globally minimizing. Since xˉ was arbitrary, every accumulation point is global.

5.1F6F10F11step 2.1step 3.1step 4.1algebra∎

The sequence lies in compact S, so [F6] supplies an accumulation point. If the global minimizer is unique, step 4.1 makes it the only possible accumulation point. Fix ε>0. If infinitely many full-sweep indices lay outside the open ε-ball about that minimizer, start with the least such index and recursively choose the least later such index; infinitude makes each candidate set nonempty, while [F10] supplies the least indices and constructs the sequence. By [F6] it has a convergent subsequence with limit in S, and [F11] keeps that limit outside the ball, producing a different accumulation point. Thus all sufficiently late full-sweep iterates lie in the ball. For a fixed within-sweep position, suppose a subsequence of those states converges to y. Its cycle-start subsequence lies in compact S, so [F6] supplies a further subsequence converging to some x∈S; steps 2.1 and 3.1 show that the corresponding within-sweep states also converge to x, hence y=x and is a full-sweep accumulation point. Therefore every convergent subsequence of within-sweep states has the unique minimizer as its limit. If infinitely many states at that position lay outside the ball, the least-index recursion and [F6] would give a convergent subsequence whose limit remains outside by [F11], a contradiction. Hence each within-sweep sequence converges to the unique minimizer as well.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Strongly convex smooth-part proximal coordinate contraction

Statement

For a finite horizon k, choose an algorithm budget K≥k and do not request a proximal-gradient stopping tolerance, so all k updates are executed.

Let F(x)=f(x)+g(x) on a finite-dimensional block product Rd, where f is convex, continuously differentiable, and globally μ-strongly convex for some μ>0, and g(x)=∑i=1ngi(xBi) with each gi proper, closed and convex. Let B1,…,Bn partition the coordinates into nonempty blocks. Suppose each block gradient is Li-Lipschitz and set Lmax⁡=max⁡iLi. Start from a finite point x0 with F(x0)<+∞. Run Randomized proximal coordinate descent with certified inner solves with uniform independent block sampling and common majorizer Mi=Lmax⁡. For any finite horizon k, assume each of its first k subproblem calls returns an accepted certificate with model error at most ϵj≥0 at update j; exact updates have ϵj=0. Then F has a unique minimizer x∗, μ≤Lmax⁡, and for q=1−μ/(nLmax⁡), E[F(xk)−F(x∗)]≤qk(F(x0)−F(x∗))+∑j=0k−1qk−1−jϵj. The expectation is the finite weighted sum over the k-step block histories, with the power convention q0=1. In particular, exact updates contract the expected gap by the factor q. Curvature of F supplied only by the nonsmooth penalties is not covered.

Facts & Assumptions

Given: The finite block product, f,gi,Li,μ, a finite-valued initial point, uniform independent sampling, and certified model errors from the Statement.

[F1]

Randomized proximal coordinate descent with certified inner solves: the selected block minimizes its quadratic model, or is accepted with certified model gap at most ϵk; other blocks stay fixed.

[F2]

Coordinate proximal subproblems are attained uniquely and descend: every exact block model has a unique minimizer and an accepted error-ϵk update raises F by at most ϵk.

[F3]

Coordinate-block smoothness descent model: the smooth term obeys the block quadratic upper model with constant Li≤Lmax⁡.

[F4]

Finite separable block penalties are proper, closed, and convex: the finite separable sum g is proper, closed and convex; its proof uses only closed sublevels of the block epigraphs and finite products of neighborhoods.

[F5]

A proper closed convex function has an affine minorant: a proper closed convex function on a finite-dimensional Euclidean space has an affine minorant.

[F6]

A finite continuous convex summand preserves proper closed convexity: adding finite continuous convex f to g preserves properness, closedness and convexity.

[F8]

Semicontinuous extreme value theorem on compact Euclidean sets: a real-valued lower-semicontinuous function on a nonempty compact Euclidean set attains its minimum.

[F10]

The principle of mathematical induction: a recurrence bound holding at zero and preserved at successors holds for every finite horizon.

[F11]

Laws of finite sums and finite products: finite weighted sums obey additivity and scaling for averaging over block choices and histories.

[A1]

f(y)≥f(x)+⟨∇f(x),y−x⟩+μ2∥y−x∥22 for all x,y.

[A2]

Mi=Lmax⁡ for every block, Li≤Lmax⁡, and each block is drawn independently and uniformly at every accepted update.

[A3]

Each finite-horizon expectation is the sum over the finite product distribution on block histories.

[A4]

All update calls through the stated horizon return accepted certificates; if an inner call fails, the conclusion is asserted only for the accepted prefix before failure.

Proof

1.1A1F4F5F6F7F8F9givenalgebrachoose

By [F4], g is proper closed convex, so choose xˉ∈dom⁡g. By [F5], g(x)≥⟨a,x⟩+b for some a,b. Strong convexity at 0 yields f(x)≥f(0)+⟨∇f(0),x⟩+μ2∥x∥22, so F is coercive. By [F9], f is continuous, so [F6] gives the closed epigraph of F; the map x↦(x,F(xˉ)) is continuous, so S={x:F(x)≤F(xˉ)} is its closed preimage. Coercivity makes S bounded, and [F7] makes it compact. It is nonempty and finite-valued, so [F8] gives a global minimizer. If two distinct minimizers existed, strong convexity of f at their midpoint and convexity of g would make F there strictly less than the minimum, a contradiction; thus the minimizer is unique.

1.2A1F3givenalgebra

For any block i and nonzero block displacement h, the strong-convexity lower model at x and [F3]'s upper model for x+Uih share the same linear term. Comparing them gives μ∥h∥22/2≤Li∥h∥22/2, so μ≤Li≤Lmax⁡. Therefore α=μ/Lmax⁡ lies in (0,1] and q=1−α/n lies in [0,1).

1.3A2F1F2F3F11givenalgebra

Fix a current finite point x and let zi minimize u↦⟨∇if(x),u−xi⟩+(Lmax⁡/2)∥u−xi∥22+gi(u); [F1]--[F2] give existence. Let z=(zi)i and Hx(w)=f(x)+⟨∇f(x),w−x⟩+(Lmax⁡/2)∥w−x∥22+g(w). This separates over blocks, so z minimizes Hx. For any block candidate v, [F3] and Li≤Lmax⁡ give F(x+Ui(v−xi))≤F(x)+qik(v)−qik(xi). The accepted certificate [F1] gives qik(xik+1)≤qik(zi)+ϵk and qik(xi)=gi(xi), so F(xk+1)≤F(x)+Δi+ϵk, where Δi=qik(zi)−qik(xi). Summing Δi over i gives Hx(z)−F(x). By [F11], uniform averaging yields E[F(xk+1)∣xk=x]≤(1−1/n)F(x)+Hx(z)/n+ϵk.

2.1A1F4step 1.1step 1.2step 1.3algebra

Let x∗ be the minimizer from step 1.1 and w=(1−α)x+αx∗. Convexity of g gives g(w)≤(1−α)g(x)+αg(x∗). Strong convexity in [A1] gives f(x)+α⟨∇f(x),x∗−x⟩≤(1−α)f(x)+αf(x∗)−(αμ/2)∥x−x∗∥22. The quadratic term in Hx(w) is Lmax⁡α2∥x−x∗∥22/2; since α=μ/Lmax⁡, it cancels the negative term. Hence Hx(z)≤Hx(w)≤(1−α)F(x)+αF(x∗).

3.1A3A4F11step 1.3step 2.1algebra

By [A4], substitution of step 2.1 into step 1.3 gives, conditionally on the history, E[F(xk+1)−F(x∗)∣xk]≤q(F(xk)−F(x∗))+ϵk. By [F11] and [A3], taking the finite-history expectation yields ϕk+1≤qϕk+ϵk, where ϕk=E[F(xk)−F(x∗)].

4.1F10step 1.2step 3.1algebra∎

Iterating the scalar recurrence gives ϕ1≤qϕ0+ϵ0; if ϕk≤qkϕ0+∑j<kqk−1−jϵj, then multiplying by q≥0 and adding ϵk proves the formula at k+1. [F10] induction over k establishes the displayed bound.

TheoremStatement: AI-adaptedProof: AI-generatedOpen item page →

Proximal coordinate contraction from strong convexity of the composite

Statement

For a finite horizon k, use an algorithm budget K≥k and do not request a proximal-gradient stopping tolerance before that horizon.

Let B1,…,Bn be a partition of Rd into nonempty coordinate blocks, with n,d≥1. Let f:Rd→R be convex and continuously differentiable, with positive block-gradient Lipschitz constants Li, and let g(x)=∑i=1ngi(xBi), where each gi is proper, closed and convex in the closed-epigraph convention. Set F=f+g and Lmax⁡=max⁡iLi>0. Assume F itself is μ-strongly convex for some μ>0: for finite-domain x,y and s∈[0,1], F((1−s)x+sy)≤(1−s)F(x)+sF(y)−μ2s(1−s)∥x−y∥22. No strong convexity of f alone is required, and no inequality μ≤Lmax⁡ is assumed.

Run Randomized proximal coordinate descent with certified inner solves from F(x0)<+∞, with common majorizer Mi=Lmax⁡ and uniform independent block draws. Fix a finite horizon k and deterministic bounds ϵ0,…,ϵk−1≥0. Assume the first k updates are executed with accepted inner-model gap certificates at most ϵj at every history in the finite block-history distribution; a stopping or failed call supplies no unexecuted later iterate. Then F has a unique minimizer x∗ and, with qc=1−μn(Lmax⁡+μ)∈(0,1), E[F(xk)−F(x∗)]≤qck(F(x0)−F(x∗))+∑j=0k−1qck−1−jϵj. The expectation is the finite weighted sum over the uniformly sampled block histories, not a distribution conditioned on successful calls. At k=0 the sum is empty and the bound is equality. Exact updates have ϵj=0. This includes curvature contributed by the separable penalties; it does not assert the sharper factor 1−μ/(nLmax⁡) without additional smooth-part curvature.

Facts & Assumptions

Given: The functions, constants, strong-convexity segment inequality, finite initial point, sampling rule and accepted per-history updates in the Statement.

[F1]

Randomized proximal coordinate descent with certified inner solves: the selected block minimizes its quadratic model, or is accepted with certified model gap at most ϵk; other blocks stay fixed.

[F2]

Coordinate proximal subproblems are attained uniquely and descend: every exact block model has a unique minimizer and an accepted error-ϵk update raises F by at most ϵk.

[F3]

Coordinate-block smoothness descent model: the smooth term obeys the block quadratic upper model with constant Li≤Lmax⁡.

[F4]

Finite separable block penalties are proper, closed, and convex: the finite separable sum g is proper, closed and convex; its proof uses only closed sublevels of the block epigraphs and finite products of neighborhoods.

[F5]

A proper closed convex function has an affine minorant: a proper closed convex function on a finite-dimensional Euclidean space has an affine minorant.

[F6]

A finite continuous convex summand preserves proper closed convexity: adding finite continuous convex f to g preserves properness, closedness and convexity.

[F8]

Semicontinuous extreme value theorem on compact Euclidean sets: a real-valued lower-semicontinuous function on a nonempty compact Euclidean set attains its minimum.

[F9]

Proper closed convex extended-real functions on Rn: properness is nonempty finite domain, convexity is the extended-real segment inequality, and closedness is closedness of the epigraph.

[F10]

Differentiable convex functions are characterized by the gradient inequality: for convex differentiable f, f(w)≥f(x)+⟨∇f(x),w−x⟩.

[F12]

The principle of mathematical induction: a recurrence bound holding at zero and preserved at successors holds for every finite horizon.

[F13]

Laws of finite sums and finite products: finite weighted sums obey additivity and scaling for averaging over block choices and histories.

[A1]

The segment strong-convexity inequality in the Statement holds for F, with μ>0; f is convex and Lmax⁡>0.

[A2]

Every executed history through the horizon has accepted certificates, with deterministic errors ϵj; block sampling remains independent and uniform, without conditioning on successful histories. Finite-sum laws [F13] justify averaging the conditional bounds over these histories.

Proof

1.1A1F4F5F6F7F8F9F11givenconstructalgebra

By [F4], g is proper closed convex, and by [F11], f is continuous. Applying [F6] makes F proper closed convex. Define H(x)=F(x)−μ∥x∥22/2, keeping +∞ off the finite domain. It is proper. The identity ∥(1−s)x+sy∥22=(1−s)∥x∥22+s∥y∥22−s(1−s)∥x−y∥22 and [A1] show H is convex, including the extended-real endpoint convention. The epigraph of H is the preimage of the closed epigraph of F under the continuous map (x,a)↦(x,a+μ∥x∥22/2); hence H is closed. By [F5] applied to H, F(x)≥μ∥x∥22/2+⟨a,x⟩+b for some a,b, so F is coercive. For a finite-domain xˉ, the nonempty sublevel S={x:F(x)≤F(xˉ)} is closed by the closed epigraph and bounded by coercivity, hence compact by [F7]. The restriction F∣S is real-valued and lower semicontinuous, so [F8] attains its minimum. Points outside S cannot improve it. The midpoint inequality [A1] forbids distinct minimizers, giving a unique x∗.

1.2A2F1F2F3F13givenalgebra

Fix a finite current point x. For block i let zi be the exact minimizer of Qi(u)=⟨∇if(x),u−xBi⟩+Lmax⁡2∥u−xBi∥22+gi(u). Existence and uniqueness follow from [F2]. Define Mx(w)=f(x)+⟨∇f(x),w−x⟩+Lmax⁡2∥w−x∥22+g(w). Separation of the finite block sum shows z=(zi)i minimizes Mx, and ∑i[Qi(zi)−Qi(xBi)]=Mx(z)−F(x). By [F3] and Li≤Lmax⁡, the actual update in a sampled block satisfies F(x+)≤F(x)+Qi(zi)−Qi(xBi)+ϵj. By [F13], averaging over the uniform block choice at this fixed history therefore gives E[F(x+)∣history]≤F(x)+Mx(z)−F(x)n+ϵj. All quantities are finite at these points because accepted model gaps and exact block minima are finite.

2.1A1F9F10step 1.1step 1.2givenalgebra

By [F10], convex differentiability of f gives f(x)+⟨∇f(x),w−x⟩≤f(w). Thus, for ws=(1−s)x+sx∗, Mx(ws)−F(x)≤F(ws)−F(x)+Lmax⁡s22∥x∗−x∥22. These segment values are finite by convexity of the effective domain. Apply [A1] to obtain Mx(z)−F(x)≤−s(F(x)−F(x∗))+Lmax⁡s2−μs(1−s)2∥x∗−x∥22. Choose s=μ/(Lmax⁡+μ)∈(0,1). The squared-distance coefficient vanishes, so Mx(z)−F(x)≤−s(F(x)−F(x∗)). This does not require μ≤Lmax⁡.

3.1A2step 1.2step 2.1algebra

Substitute step 2.1 in step 1.2 to get E[F(xj+1)−F(x∗)∣history]≤qc(F(xj)−F(x∗))+ϵj. By [A2] all histories in the finite product distribution are present; summing the inequality with their weights gives ϕj+1≤qcϕj+ϵj, where ϕj=E[F(xj)−F(x∗)]. No conditioning on the event of a successful inner solve is used.

4.1F12step 2.1step 3.1algebra∎

The proposed bound at k=0 is ϕ0=ϕ0. If ϕj≤qcjϕ0+∑ℓ<jqcj−1−ℓϵℓ, multiply by qc>0 and use step 3.1 to obtain the bound at j+1. By [F12], finite induction proves the displayed convolution. It specializes to exact-update contraction when every ϵj=0.

5 · Examples, counterexamples and false statements

None yet.

Sources