Papers
Topics
Authors
Recent
Search
2000 character limit reached

Greedy-Threshold Methods: Theory & Applications

Updated 4 July 2026
  • Greedy-threshold methods are techniques that replace exact maximizers with adaptive threshold rules in greedy algorithms for approximation and optimization.
  • They generalize classical approaches by incorporating parameters like τ, λ, or f(m) to control selection criteria and benchmark comparisons.
  • These methods find applications in sparse recovery, nonlinear inverse problems, and decentralized control, balancing computational cost and performance.

Greedy-threshold denotes a family of constructions in which a greedy procedure is governed by a threshold condition rather than by an exact maximizer. In the classical Thresholding Greedy Algorithm (TGA), the threshold is implicit in the requirement that the selected coefficients are the largest in modulus; in weak and generalized variants, admissibility is relaxed by parameters such as τ\tau, λ\lambda, or a comparison function f(m)f(m). In algorithmic settings outside Banach-space approximation, the same pattern appears when one accepts any atom, row, task, or candidate whose score exceeds a prescribed, adaptive, or capped threshold, and in some analyses the threshold instead marks a regime boundary below which a greedy policy is exactly optimal and above which it is not. This suggests a unifying view of greedy-threshold methods as threshold-controlled relaxations, generalizations, or phase-transition analyses of greediness (Berná et al., 2022, Berasategui et al., 2023, Yang et al., 2013, Xu et al., 2016, Wang et al., 2019).

1. Foundational definitions

In the approximation-theoretic literature, the basic object is a basis B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty with biorthogonal functionals (en)(e_n^*). For

x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,

a finite set ANA\subset\mathbb N of cardinality mm is called a greedy set when the selected coefficients dominate the unselected ones in modulus. One common formulation is

minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,

and the associated greedy operator is

Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.

Equivalent expositions use the non-strict relation λ\lambda0 together with a tie-breaking rule by natural ordering. The TGA is then the sequence λ\lambda1 obtained by repeatedly selecting the largest coefficients in modulus (Chu, 2022, Berná et al., 2022).

A first generalization introduces an explicit threshold parameter λ\lambda2. The λ\lambda3-Thresholding Greedy Algorithm chooses any set λ\lambda4 of size λ\lambda5 satisfying

λ\lambda6

and forms

λ\lambda7

When λ\lambda8, this is the classical TGA. A second generalization enlarges the greedy sum itself: instead of comparing the error after λ\lambda9 greedy terms with best f(m)f(m)0-term benchmarks, one defines

f(m)f(m)1

leading to f(m)f(m)2-almost greedy and f(m)f(m)3-partially greedy bases (Berasategui et al., 2023, Chu, 2022).

A further axis of generalization replaces the benchmark size f(m)f(m)4 by a function f(m)f(m)5. For f(m)f(m)6 in the class f(m)f(m)7 of continuous, increasing, concave maps with f(m)f(m)8, f(m)f(m)9, and B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty0 on B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty1, a basis is called B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty2-greedy if

B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty3

and B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty4-almost greedy if

B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty5

for all B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty6, B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty7, and greedy sets B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty8 of size B=(en)n=1\mathcal B=(e_n)_{n=1}^\infty9 (Chu, 2022). These variants make explicit that “threshold” may act on coefficient admissibility, on greedy-sum size, or on the benchmark against which greedy error is measured.

2. Structural theory in Banach and quasi-Banach spaces

The classical characterization of greedy bases is the Konyagin–Temlyakov theorem: a basis is greedy if and only if it is unconditional and democratic. The almost-greedy property replaces the best linear (en)(e_n^*)0-term error (en)(e_n^*)1 by the best projection error

(en)(e_n^*)2

and almost-greedy bases are exactly those that are quasi-greedy and democratic. In parallel, quasi-greediness is equivalent to convergence of greedy approximants and to uniform boundedness of greedy operators (Berná et al., 2022, Albiac et al., 2019).

Threshold enlargement does not behave uniformly across all greedy-type classes. For almost-greedy bases, the enlargement parameter (en)(e_n^*)3 does not produce a genuinely new class: a basis is almost greedy if and only if it is (en)(e_n^*)4-almost greedy for some, equivalently all, (en)(e_n^*)5. By contrast, for each (en)(e_n^*)6 there exists an unconditional basis that is (en)(e_n^*)7-partially greedy but is not (en)(e_n^*)8-partially greedy. The characterization of (en)(e_n^*)9-partially greedy bases is

x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,0

where x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,1-max-conservativity is expressed by

x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,2

for sets x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,3 with

x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,4

This distinction corrects a common oversimplification: allowing larger greedy sums preserves almost-greediness but strictly weakens partial greed when x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,5 (Chu, 2022).

The x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,6-theory yields a different phenomenon. For every non-identity x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,7,

x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,8

The paper also characterizes x=n=1en(x)en,x=\sum_{n=1}^\infty e_n^*(x)\,e_n,9-almost greedy as quasi-greedy plus Property ANA\subset\mathbb N0, and ANA\subset\mathbb N1-greedy as ANA\subset\mathbb N2-unconditional plus Property ANA\subset\mathbb N3. This collapses the distinction between best linear and best threshold approximants once the comparison size is strictly sublinear. A plausible implication is that sublinear benchmark growth regularizes the gap between linear and coordinate projection models of approximation (Chu, 2022).

Several later refinements restrict the admissible comparison supports. Consecutive greedy bases compare greedy error only with approximants supported on intervals of ANA\subset\mathbb N4, while CGPCC restricts further to interval-supported polynomials with constant coefficients. On Schauder bases, consecutive-greedy is equivalent to almost-greedy. The sequential theory replaces intervals by families ANA\subset\mathbb N5 generated by a fixed gap sequence; the ANA\subset\mathbb N6-almost greedy property and the ANA\subset\mathbb N7-strong partially greedy property are equivalent to their classical counterparts if and only if ANA\subset\mathbb N8 is bounded. Under additional conditions, the ANA\subset\mathbb N9-minimum partially greedy property lies strictly between almost-greedy and strong partially greedy (Berasategui et al., 2023, Berasategui et al., 2023).

Quantitative control is given by Lebesgue-type constants and related parameters. The Lebesgue constants

mm0

satisfy

mm1

where mm2 are unconditionality parameters and mm3 are the squeeze-symmetry constants. This answers Temlyakov’s 2011 question about a natural greedy-type parameter complementing mm4 in an exact linear rule for mm5. Rate-of-convergence results then bound the TGA error by a product of mm6 and mm7, and in mm8 the trigonometric system and the Haar basis admit sharp rates (Albiac et al., 2021, Temlyakov, 2023). Related work also characterizes greedy bases by approximation with polynomials of constant coefficients, including the RGPCC and URGPCC formulations, and extends these equivalences to quasi-Banach spaces (Berasategui et al., 2023).

3. Sparse recovery and thresholded pursuit

In compressive sensing, thresholding often replaces the exact “largest inner product” step of a greedy pursuit. Orthogonal Matching Pursuit with Thresholding (OMPT) takes measurements mm9, a normalized dictionary minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,0, and a threshold parameter minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,1. At iteration minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,2, instead of selecting the atom of maximum correlation with the residual minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,3, OMPT chooses any index minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,4 satisfying

minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,5

The support is updated, a least-squares projection is performed on the selected support, and the residual is recomputed. The analysis uses the global minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,6-coherence

minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,7

which bridges mutual coherence and the restricted isometry constant via

minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,8

Under

minnAen(x)>maxnAen(x),\min_{n\in A}|e_n^*(x)|>\max_{n\notin A}|e_n^*(x)|,9

and

Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.0

OMPT recovers a Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.1-sparse signal exactly in Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.2 steps, and the paper states that these conditions coincide with the best known bounds for standard OMP while reducing the cost of the selection step (Yang et al., 2013).

Thresholding Greedy Pursuit (TGP) applies the same principle to a CoSaMP-style framework. At each iteration it computes the proxy

Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.3

keeps all coordinates above a scalar threshold,

Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.4

merges them into the support estimate, solves a least-squares problem on the merged support, and updates the residual. The theory emphasizes asymptotic sparse recovery with additive noise. With

Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.5

the “No Phantom Signal” theorem yields Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.6 when Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.7 with probability at least Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.8, and the “No False Discoveries” theorem gives

Gm(x)=nAen(x)en.G_m(x)=\sum_{n\in A}e_n^*(x)\,e_n.9

under λ\lambda00. In the exact-recovery regime, if

λ\lambda01

and

λ\lambda02

then TGP recovers exactly

λ\lambda03

with probability λ\lambda04. The paper stresses that the threshold can be chosen without the knowledge of sparsity level of the signal and strength of the noise (Le et al., 2021).

Taken together, OMPT and TGP show that thresholding is not merely a heuristic weakening of greedy pursuit. In both cases, the thresholded rule replaces a more expensive exact maximization while retaining recovery guarantees under coherence- or concentration-based assumptions (Yang et al., 2013, Le et al., 2021).

4. Learning and nonlinear inverse problems

Orthogonal Greedy Learning (OGL) provides a learning-theoretic version of the same idea. Classical OGL selects

λ\lambda05

and updates the predictor by orthogonal projection onto the span of selected atoms. The λ\lambda06-greedy-threshold metric declares an atom λ\lambda07 to be active when

λ\lambda08

or, in the closely related formulation,

λ\lambda09

The λ\lambda10-Threshold OGL algorithm chooses any active atom, projects orthogonally onto the enlarged span, and stops when no active atom remains or when

λ\lambda11

Under the source condition

λ\lambda12

the output λ\lambda13 satisfies, with probability at least λ\lambda14,

λ\lambda15

and choosing

λ\lambda16

yields the rate λ\lambda17. The papers emphasize that steepest gradient descent is not the unique greedy criterion of OGL and that the threshold parameter acts simultaneously as a greediness control and an adaptive stopping rule (Xu et al., 2016, Xu et al., 2014).

A related pattern appears in nonlinear Kaczmarz methods. The greedy capped nonlinear Kaczmarz framework introduces a capped set of promising rows at each iterate. In the distance–residual capped method, one defines

λ\lambda18

and the capped set

λ\lambda19

In the residual–distance capped method, one defines

λ\lambda20

Randomization is then performed only within λ\lambda21 or λ\lambda22, with probabilities determined by residual- or distance-based scores. Under the local tangential cone condition, the resulting single-row and block methods satisfy expected contraction bounds that the paper states are strictly smaller than the NRK factor. This is an explicit exploration–exploitation design: discard low-score rows by a threshold, then sample within the retained set (Zhang et al., 2022).

5. Threshold-greedy optimization and decentralized control

For monotone submodular task allocation under a partition matroid, the Decreasing Threshold Task Allocation (DTTA) algorithm initializes a global threshold

λ\lambda23

and then sweeps

λ\lambda24

At each threshold level, each robot proposes any task with marginal gain at least λ\lambda25, conflicts are resolved by a coordination step, and many tasks can be allocated in parallel. The resulting guarantee is

λ\lambda26

with time complexity

λ\lambda27

The lazy variant LDTTA maintains cached marginals in a max-heap and uses submodularity to avoid recomputing most gains, while retaining the same λ\lambda28 approximation guarantee (Li et al., 2019).

Noisy submodular maximization uses thresholding in a statistically adaptive form. The Confident Sample subroutine decides whether an unknown mean is approximately above or below a threshold λ\lambda29 with probability at least λ\lambda30, and the monotone cardinality algorithm CTG integrates it into a threshold schedule based on

λ\lambda31

With probability at least λ\lambda32, CTG returns λ\lambda33 with

λ\lambda34

The same framework extends to unconstrained non-monotone maximization through CDG and to matroid constraints through CCTG, the latter obtaining

λ\lambda35

Here the threshold is not fixed solely by optimization geometry; it is coupled to adaptive sampling and confidence control (Chen et al., 2023).

Threshold control can also be the object being optimized. In persistent monitoring on graphs, each agent carries a threshold matrix λ\lambda36 that determines dwell times and next-hop decisions by testing whether the current node uncertainty falls below λ\lambda37 and whether some neighboring uncertainty exceeds λ\lambda38. Because local optimization over thresholds is highly initialization-dependent, the paper develops an off-line greedy initialization based on asymptotic cycle analysis and a target-cycle expansion operation. The resulting method constructs a high-performing set of initial thresholds and then refines them by on-line IPA-driven gradient descent. This use of greediness is secondary to the threshold controller itself, but it shows that threshold-greedy methodology also appears as an initialization strategy for hybrid distributed control (Welikala et al., 2019).

6. Thresholds as optimality boundaries and phase transitions

Not all greedy-threshold results are algorithmic selection rules. In some settings, the threshold is a regime boundary that exactly characterizes when a greedy policy is optimal. For a battery limited energy harvesting communication system with i.i.d. arrivals, finite battery capacity λ\lambda39, and reward λ\lambda40 increasing, concave, and continuously differentiable, the greedy policy uses all available battery energy in each slot. Its baseline throughput is

λ\lambda41

The main theorem states that there exists a critical battery size

λ\lambda42

such that

λ\lambda43

For the AWGN reward,

λ\lambda44

this becomes

λ\lambda45

The paper also provides lower and upper bounds on λ\lambda46 and asymptotic formulas for geometric, Poisson, uniform, exponential, and Rayleigh arrivals. A common misconception is that greediness should improve with larger storage; this result shows the opposite structure: greedy is exactly optimal only on the interval λ\lambda47 (Wang et al., 2019).

A second regime-threshold interpretation appears in online edge coloring. The folklore greedy algorithm uses at most

λ\lambda48

colors. The sharp-threshold results show that this guarantee is unimprovable for deterministic algorithms when λ\lambda49 and for randomized algorithms when λ\lambda50, matching the classical impossibility results. Beyond those regimes, however, the greedy barrier disappears: there is a deterministic online algorithm achieving λ\lambda51-colorings for all λ\lambda52, and a randomized algorithm achieving λ\lambda53-colorings already for λ\lambda54. The paper explicitly describes these as sharp thresholds for when greedy can be surpassed (Blikstad et al., 29 Jul 2025).

These phase-transition results use “threshold” in a different sense from TGA, OMPT, or DTTA. The threshold does not define admissible local moves; it delineates the structural regime in which greediness is exactly optimal, provably optimal only up to a factor, or asymptotically improvable. This suggests that greedy-threshold is best understood as a broader research motif rather than a single algorithmic template (Wang et al., 2019, Blikstad et al., 29 Jul 2025).

Definition Search Book Streamline Icon: https://streamlinehq.com
References (19)

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Greedy-Threshold.