luis@199
|
1 %% Configuration file used in SMALL_solve
|
luis@199
|
2 %
|
luis@199
|
3 % Use this file to change the solvers in SMALLBox
|
luis@199
|
4 % Please refer to the documentation before editing this file
|
luis@199
|
5
|
luis@199
|
6 % Centre for Digital Music, Queen Mary, University of London.
|
luis@199
|
7 % This file copyright 2009 Ivan Damnjanovic.
|
luis@199
|
8 %
|
luis@199
|
9 % This program is free software; you can redistribute it and/or
|
luis@199
|
10 % modify it under the terms of the GNU General Public License as
|
luis@199
|
11 % published by the Free Software Foundation; either version 2 of the
|
luis@199
|
12 % License, or (at your option) any later version. See the file
|
luis@199
|
13 % COPYING included with this distribution for more information.
|
luis@199
|
14 %
|
luis@199
|
15 %%
|
luis@199
|
16
|
luis@199
|
17 if strcmpi(solver.toolbox,'sparselab')
|
luis@199
|
18 y = eval([solver.name,'(SparseLab_A, b, n,',solver.param,');']);
|
luis@199
|
19 elseif strcmpi(solver.toolbox,'sparsify')
|
luis@199
|
20 if isa(Problem.A,'float')
|
luis@199
|
21 y = eval([solver.name,'(b, A, n,',solver.param,');']);
|
luis@199
|
22 else
|
luis@199
|
23 y = eval([solver.name,'(b, A, n, ''P_trans'', AT,',solver.param,');']);
|
luis@199
|
24 end
|
luis@199
|
25 elseif (strcmpi(solver.toolbox,'spgl1')||strcmpi(solver.toolbox,'gpsr'))
|
luis@199
|
26 y = eval([solver.name,'(b, A,',solver.param,');']);
|
luis@199
|
27 elseif (strcmpi(solver.toolbox,'SPAMS'))
|
luis@199
|
28 y = eval([solver.name,'(b, A, solver.param);']);
|
luis@199
|
29 elseif (strcmpi(solver.toolbox,'SMALL'))
|
luis@199
|
30 if isa(Problem.A,'float')
|
luis@199
|
31 y = eval([solver.name,'(A, b, n,',solver.param,');']);
|
luis@199
|
32 else
|
luis@199
|
33 y = eval([solver.name,'(A, b, n,',solver.param,',AT);']);
|
luis@199
|
34 end
|
luis@199
|
35 elseif (strcmpi(solver.toolbox, 'ompbox'))
|
luis@199
|
36 G=A'*A;
|
luis@199
|
37 epsilon=solver.param.epsilon;
|
luis@199
|
38 maxatoms=solver.param.maxatoms;
|
luis@199
|
39 y = eval([solver.name,'(A, b, G,epsilon,''maxatoms'',maxatoms,''checkdict'',''off'');']);
|
luis@199
|
40 elseif (strcmpi(solver.toolbox, 'ompsbox'))
|
luis@199
|
41 basedict = Problem.basedict;
|
luis@199
|
42 if issparse(Problem.A)
|
luis@199
|
43 A = Problem.A;
|
luis@199
|
44 else
|
luis@199
|
45 A = sparse(Problem.A);
|
luis@199
|
46 end
|
luis@199
|
47 G = dicttsep(basedict,A,dictsep(basedict,A,speye(size(A,2))));
|
luis@199
|
48 epsilon=solver.param.epsilon;
|
luis@199
|
49 maxatoms=solver.param.maxatoms;
|
luis@199
|
50 y = eval([solver.name,'(basedict, A, b, G,epsilon,''maxatoms'',maxatoms,''checkdict'',''off'');']);
|
luis@199
|
51 Problem.sparse=1;
|
luis@199
|
52 elseif (strcmpi(solver.toolbox, 'ALPS'))
|
luis@199
|
53 if ~isa(Problem.A,'float')
|
luis@199
|
54 % ALPS does not accept implicit dictionary definition
|
luis@199
|
55 A = opToMatrix(Problem.A, 1);
|
luis@199
|
56 end
|
luis@199
|
57 [y, numiter, time, y_path] = wrapper_ALPS_toolbox(b, A, solver.param);
|
luis@199
|
58 elseif (strcmpi(solver.toolbox, 'MMbox'))
|
luis@199
|
59 if ~isa(Problem.A,'float')
|
luis@199
|
60 % MMbox does not accept implicit dictionary definition
|
luis@199
|
61 A = opToMatrix(Problem.A, 1);
|
luis@199
|
62 end
|
luis@199
|
63
|
luis@199
|
64 [y, cost] = wrapper_mm_solver(b, A, solver.param);
|
luis@199
|
65
|
luis@199
|
66 % To introduce new sparse representation algorithm put the files in
|
luis@199
|
67 % your Matlab path. Next, unique name <TolboxID> for your toolbox and
|
luis@199
|
68 % prefferd API <Preffered_API> needs to be defined.
|
luis@199
|
69 %
|
luis@199
|
70 % elseif strcmpi(solver.toolbox,'<ToolboxID>')
|
luis@199
|
71 %
|
luis@199
|
72 % % - Evaluate the function (solver.name - defined in the main) with
|
luis@199
|
73 % % parameters given above
|
luis@199
|
74 %
|
luis@199
|
75 % y = eval([solver.name,'(<Preffered_API>);']);
|
luis@199
|
76
|
luis@199
|
77 else
|
luis@199
|
78 printf('\nToolbox has not been registered. Please change SMALL_solver file.\n');
|
luis@199
|
79 return
|
luis@199
|
80 end
|