# HG changeset patch # User idamnjanovic # Date 1267792282 0 # Node ID 01cad25206d614dcc480cbc5771510154e6a6546 # Parent 257c289bf11daadfdf8c258477295be1433d5d3c diff -r 257c289bf11d -r 01cad25206d6 solvers/SMALL_MP.m --- a/solvers/SMALL_MP.m Fri Mar 05 12:23:57 2010 +0000 +++ b/solvers/SMALL_MP.m Fri Mar 05 12:31:22 2010 +0000 @@ -1,4 +1,4 @@ -function [A]=SMALL_MP(Dict,X, m, maxNumCoef, errorGoal1, varargin) +function [A]=SMALL_MP(Dict,X, m, maxNumCoef, errorGoal, varargin) %% %============================================= % Sparse coding of a group of signals based on a given @@ -57,12 +57,12 @@ j=0; currResNorm = norm(residual); - errorGoal=errorGoal1*currResNorm; + errGoal=errorGoal*currResNorm; a = zeros(m,1); - while currResNorm>errorGoal & j < maxNumCoef, + while currResNorm>errGoal && j < maxNumCoef, j = j+1; dir=Dt(residual); diff -r 257c289bf11d -r 01cad25206d6 solvers/SMALL_cgp.m --- a/solvers/SMALL_cgp.m Fri Mar 05 12:23:57 2010 +0000 +++ b/solvers/SMALL_cgp.m Fri Mar 05 12:31:22 2010 +0000 @@ -1,4 +1,4 @@ -function [A]=SMALL_cgp(Dict,X, m, maxNumCoef, errorGoal, varargin) +function [A, resF]=SMALL_cgp(Dict,X, m, maxNumCoef, errorGoal, varargin) %============================================= % Sparse coding of a group of signals based on a given % dictionary and specified number of atoms to use. @@ -59,7 +59,7 @@ currResNorm = residual'*residual/n; - errorGoal=errorGoal*currResNorm; + errGoal=errorGoal*currResNorm; a = zeros(m,1); p = zeros(m,1); @@ -88,12 +88,13 @@ residual=residual-alfa*Dp; currResNorm = residual'*residual/n; - if currResNorm