In-crowd algorithm
The in-crowd algorithm is a numerical method for solving basis pursuit denoising quickly; faster than any other algorithm for large, sparse problems.[1] This algorithm is an active set method, which minimizes iteratively sub-problems of the global basis pursuit denoising: [math]\displaystyle{ \min_x \frac{1}{2}\|y-Ax\|^2_2+\lambda\|x\|_1. }[/math]
where [math]\displaystyle{ y }[/math] is the observed signal, [math]\displaystyle{ x }[/math] is the sparse signal to be recovered, [math]\displaystyle{ Ax }[/math] is the expected signal under [math]\displaystyle{ x }[/math], and [math]\displaystyle{ \lambda }[/math] is the regularization parameter trading off signal fidelity and simplicity. The simplicity is here measured using the sparsity of the solution [math]\displaystyle{ x }[/math], measure through its [math]\displaystyle{ \ell_1 }[/math]-norm. The active set strategies are very efficient in this context as only few coefficient are expected to be non-zero. Thus, if they can be identified, solving the problem restricted to these coefficients yield the solution. Here, the features are greedily selected based on the absolute value of their gradient at the current estimate.
Other active-set methods for the basis pursuit denoising includes BLITZ,[2] where the selection of the active set is performed using the duality gap of the problem, and The Feature Sign Search,[3] where the features are included based on the estimate of their sign.
Algorithm
It consists of the following:
- Declare [math]\displaystyle{ x }[/math] to be 0, so the unexplained residual [math]\displaystyle{ r = y }[/math]
- Declare the active set [math]\displaystyle{ I }[/math] to be the empty set, and [math]\displaystyle{ I^c }[/math] to be its complement (the inactive set)
- Calculate the usefulness [math]\displaystyle{ u_j = | \langle r A_j \rangle | }[/math] for each component in [math]\displaystyle{ I^c }[/math]
- If on [math]\displaystyle{ I^c }[/math], no [math]\displaystyle{ u_j \gt \lambda }[/math], terminate
- Otherwise, add [math]\displaystyle{ L \approx 25 }[/math] components to [math]\displaystyle{ I }[/math] based on their usefulness
- Solve basis pursuit denoising exactly on [math]\displaystyle{ I }[/math], and throw out any component of [math]\displaystyle{ I }[/math] whose value attains exactly 0. This problem is dense, so quadratic programming techniques work very well for this sub problem.
- Update [math]\displaystyle{ r = y - Ax }[/math] - n.b. can be computed in the subproblem as all elements outside of [math]\displaystyle{ I }[/math] are 0
- Go to step 3.
Since every time the in-crowd algorithm performs a global search it adds up to [math]\displaystyle{ L }[/math] components to the active set, it can be a factor of [math]\displaystyle{ L }[/math] faster than the best alternative algorithms when this search is computationally expensive. A theorem[1] guarantees that the global optimum is reached in spite of the many-at-a-time nature of the in-crowd algorithm.
Notes
- ↑ 1.0 1.1 See The In-Crowd Algorithm for Fast Basis Pursuit Denoising, IEEE Trans Sig Proc 59 (10), Oct 1 2011, pp. 4595 - 4605, [1], demo MATLAB code available [2]
- ↑ Johnson T, Guestrin C. Blitz: A principled meta-algorithm for scaling sparse optimization. In proceedings of the International Conference on Machine Learning (ICML) 2015 (pp. 1171-1179).([3])
- ↑ Lee H, Battle A, Raina R, Ng AY. Efficient sparse coding algorithms. In Advances in neural information processing systems 2007 (pp. 801-808). [4]
 |  |
