comparison solvers/SMALL_chol.m @ 63:e08af264ef93

(none)

 % This Dictionary check is based on Thomas Blumensath work in sparsify 0_4 greedy solvers 
 explicitD=0;
 if     isa(Dict,'float')      
             D =@(z) Dict*z;  
             Dt =@(z) Dict'*z;
-            explicitD=1;
+            explicitD=0;
 elseif isobject(Dict) 
             D =@(z) Dict*z;  
             Dt =@(z) Dict'*z;
 elseif isa(Dict,'function_handle') 
     try
         DictT=varargin{1};
         if isa(DictT,'function_handle'); 
-                D=Dict; 
-                Dt=DictT;
+                D=@(z)Dict(z,1); 
+                Dt=@(z)Dict(z,2);
         else
                 error('If Dictionary is a function handle,Transpose Dictionary also needs to be a function handle. ');
         end
     catch
         error('If Dictionary is a function handle, Transpose Dictionary  needs to be specified. Exiting.');

     x=X(:,k);
     residual=x;
 	indx = [];
 	a = zeros(m,1);
 	currResNorm = norm(residual);
-    errorGoal=errorGoal*currResNorm;
+    errorGoal1=errorGoal*currResNorm;
 	j = 0;
     while currResNorm>errorGoal & j < maxNumCoef,
 		j = j+1;
         dir=Dt(residual);