idamnjanovic@1: function SMALL_solver_test
ivan@128: %%  Example test of solvers from different toolboxes on Sparco problem 6
idamnjanovic@1: %
idamnjanovic@1: %   The main purpose of this example is to show how to use SMALL structure
idamnjanovic@1: %   to solve SPARCO compressed sensing problems (1-11) and compare results
idamnjanovic@1: %   from different solvers.
idamnjanovic@1: %   To generate SMALL.Problem part of structure you can use generateProblem
idamnjanovic@1: %   function from Sparco toolbox giving the problem number and any
idamnjanovic@1: %   additional parameters you might want to change. Alternatively, you can
idamnjanovic@1: %   might want to consult sparco documentation to write a problem by
idamnjanovic@1: %   yourself. There are four fields the must be specified in SMALL.Problem 
idamnjanovic@1: %   - A, b, sizeA and reconstruct.
idamnjanovic@1: %   
idamnjanovic@1: %   To generate SMALL.solver part of the structure you must specify three
idamnjanovic@1: %   fields:
idamnjanovic@1: %   
idamnjanovic@1: %       SMALL.solver.toolbox - string with toolbox name is needed because
idamnjanovic@1: %                              different toolboxes are calling solver 
idamnjanovic@1: %                              functions in different ways.
idamnjanovic@1: %       SMALL.solver.name - its string representing solver name (e.g.
idamnjanovic@1: %                           SolveOMP)
idamnjanovic@1: %       SMALL.solver.param - string that contains optional parameters for
idamnjanovic@1: %                            particular solver (all parameters you want to
idamnjanovic@1: %                            specify except A, b and size of solution)
idamnjanovic@1: %                            
idamnjanovic@1: %   Every call to SMALL_solve function will generate following output:
idamnjanovic@1: %
idamnjanovic@1: %       SMALL.solver.solution - contains solution vector x
idamnjanovic@1: %       SMALL.solver.reconstructed - vector containing signal reconstructed
idamnjanovic@1: %                                    from the solution
idamnjanovic@1: %       SMALL.solver.time - time that solver spent to find the solution
idamnjanovic@1: %           
idamnjanovic@1: %   SMALL_plot function plots the SMALL.solver.solution and reconstructed
idamnjanovic@1: %   against original signal.
idamnjanovic@1: %   
idamnjanovic@1: %   In this particular example we are testing SMALL_cgp, SMALL_chol, 
idamnjanovic@1: %   SolveOMP form SparseLab and greed_pcgp form Sparsify against "PROB006  
idamnjanovic@1: %   Daubechies basis, Gaussian ensemble measurement basis, piecewise cubic 
idamnjanovic@1: %   polynomial signal" from Sparco. 
idamnjanovic@1: %   
idamnjanovic@1: %
idamnjanovic@25: 
idamnjanovic@25: 
maria@83: %   Centre for Digital Music, Queen Mary, University of London.
maria@83: %   This file copyright 2009 Ivan Damnjanovic.
maria@83: %
maria@83: %   This program is free software; you can redistribute it and/or
maria@83: %   modify it under the terms of the GNU General Public License as
maria@83: %   published by the Free Software Foundation; either version 2 of the
maria@83: %   License, or (at your option) any later version.  See the file
maria@83: %   COPYING included with this distribution for more information.
maria@83: %%   
idamnjanovic@1: 
idamnjanovic@20: fprintf('\n\nExample test of SMALL solvers against their counterparts on Sparco problems.\n\n');
idamnjanovic@1: 
idamnjanovic@1: %%
idamnjanovic@1: % Generate SPARCO problem 
idamnjanovic@41: clear  
idamnjanovic@4: 
idamnjanovic@1: SMALL.Problem = generateProblem(6, 'P', 6, 'm', 270,'n',1024, 'show');
idamnjanovic@1: %%
idamnjanovic@4: i=1;
idamnjanovic@20: % %%
idamnjanovic@20: % % SMALL Conjugate Gradient test 
idamnjanovic@20: % SMALL.solver(i)=SMALL_init_solver;
idamnjanovic@20: % SMALL.solver(i).toolbox='SMALL';    
idamnjanovic@20: % SMALL.solver(i).name='SMALL_cgp';
idamnjanovic@20: % 
idamnjanovic@20: % % In the following string all parameters except matrix, measurement vector
idamnjanovic@20: % % and size of solution need to be specified. If you are not sure which
idamnjanovic@20: % % parameters are needed for particular solver type "help <Solver name>" in
idamnjanovic@20: % % MATLAB command line
idamnjanovic@20: % 
idamnjanovic@20: % SMALL.solver(i).param='200, 1e-14';
idamnjanovic@20: % 
idamnjanovic@20: % SMALL.solver(i)=SMALL_solve(SMALL.Problem,SMALL.solver(i));
idamnjanovic@20: % 
idamnjanovic@20: % 
idamnjanovic@20: % i=i+1;
idamnjanovic@1: %%
idamnjanovic@41: % SMALL Conjugate Gradient test 
idamnjanovic@4: SMALL.solver(i)=SMALL_init_solver;
idamnjanovic@4: SMALL.solver(i).toolbox='SMALL';    
ivan@163: SMALL.solver(i).name='SMALL_pcgp';
idamnjanovic@1: 
idamnjanovic@1: % In the following string all parameters except matrix, measurement vector
idamnjanovic@1: % and size of solution need to be specified. If you are not sure which
idamnjanovic@1: % parameters are needed for particular solver type "help <Solver name>" in
idamnjanovic@1: % MATLAB command line
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i).param='200, 1e-14';
idamnjanovic@1: 
idamnjanovic@41: SMALL.solver(i)=SMALL_solve(SMALL.Problem,SMALL.solver(i));
idamnjanovic@41: 
idamnjanovic@1: 
idamnjanovic@4: i=i+1;
idamnjanovic@41: 
idamnjanovic@1: %%
idamnjanovic@1: % SolveOMP from SparseLab test 
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i)=SMALL_init_solver;
idamnjanovic@4: SMALL.solver(i).toolbox='SparseLab';  
idamnjanovic@4: SMALL.solver(i).name='SolveOMP';
idamnjanovic@1: 
idamnjanovic@1: % In the following string all parameters except matrix, measurement vector
idamnjanovic@1: % and size of solution need to be specified. If you are not sure which
idamnjanovic@1: % parameters are needed for particular solver type "help <Solver name>" in
idamnjanovic@1: % MATLAB command line
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i).param='200, 0, 0, 0, 1e-14';
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i)=SMALL_solve(SMALL.Problem, SMALL.solver(i));
idamnjanovic@1: 
idamnjanovic@4: i=i+1;
idamnjanovic@1:   
idamnjanovic@1: %%
idamnjanovic@41: % SMALL OMP with Cholesky update test 
idamnjanovic@41: SMALL.solver(i)=SMALL_init_solver;
idamnjanovic@41: SMALL.solver(i).toolbox='SMALL';    
idamnjanovic@41: SMALL.solver(i).name='SMALL_chol';
idamnjanovic@41: 
idamnjanovic@41: % In the following string all parameters except matrix, measurement vector
idamnjanovic@41: % and size of solution need to be specified. If you are not sure which
idamnjanovic@41: % parameters are needed for particular solver type "help <Solver name>" in
idamnjanovic@41: % MATLAB command line
idamnjanovic@41: 
idamnjanovic@41: SMALL.solver(i).param='200, 1e-14';
idamnjanovic@41: 
idamnjanovic@41: SMALL.solver(i)=SMALL_solve(SMALL.Problem, SMALL.solver(i));
idamnjanovic@41: 
idamnjanovic@41: i=i+1;
idamnjanovic@41: 
idamnjanovic@41: %%
idamnjanovic@1: % greed_pcgp from Sparsify test 
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i)=SMALL_init_solver;
idamnjanovic@4: SMALL.solver(i).toolbox='Sparsify';  
idamnjanovic@4: SMALL.solver(i).name='greed_pcgp';
idamnjanovic@1: 
idamnjanovic@1: % In the following string all parameters except matrix, measurement vector
idamnjanovic@1: % and size of solution need to be specified. If you are not sure which
idamnjanovic@1: % parameters are needed for particular solver type "help <Solver name>" in
idamnjanovic@1: % MATLAB command line
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i).param='''stopCrit'', ''M'', ''stopTol'', 200';
idamnjanovic@1: 
idamnjanovic@4: SMALL.solver(i)=SMALL_solve(SMALL.Problem, SMALL.solver(i));
idamnjanovic@4: 
idamnjanovic@4: %%
idamnjanovic@1: 
idamnjanovic@1: SMALL_plot(SMALL);
idamnjanovic@1:   
idamnjanovic@4: 
idamnjanovic@1:   
idamnjanovic@1:  
idamnjanovic@1: end % function SMALL_solver_test