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'');']);
|
daniele@221
|
40 % danieleb: added call to omp functions with fast implementation.
|
daniele@221
|
41 elseif (strcmpi(solver.toolbox, 'ompbox_fast'))
|
daniele@221
|
42 DtX=A'*b;
|
daniele@221
|
43 XtX = sum(b.*b);
|
daniele@221
|
44 G=A'*A;
|
daniele@221
|
45 epsilon=solver.param.epsilon;
|
daniele@221
|
46 maxatoms=solver.param.maxatoms;
|
daniele@221
|
47 y = eval([solver.name,'(DtX, XtX, G,epsilon,''maxatoms'',maxatoms,''checkdict'',''off'');']);
|
luis@199
|
48 elseif (strcmpi(solver.toolbox, 'ompsbox'))
|
luis@199
|
49 basedict = Problem.basedict;
|
luis@199
|
50 if issparse(Problem.A)
|
luis@199
|
51 A = Problem.A;
|
luis@199
|
52 else
|
luis@199
|
53 A = sparse(Problem.A);
|
luis@199
|
54 end
|
luis@199
|
55 G = dicttsep(basedict,A,dictsep(basedict,A,speye(size(A,2))));
|
luis@199
|
56 epsilon=solver.param.epsilon;
|
luis@199
|
57 maxatoms=solver.param.maxatoms;
|
luis@199
|
58 y = eval([solver.name,'(basedict, A, b, G,epsilon,''maxatoms'',maxatoms,''checkdict'',''off'');']);
|
luis@199
|
59 Problem.sparse=1;
|
luis@199
|
60 elseif (strcmpi(solver.toolbox, 'ALPS'))
|
luis@199
|
61 if ~isa(Problem.A,'float')
|
luis@199
|
62 % ALPS does not accept implicit dictionary definition
|
luis@199
|
63 A = opToMatrix(Problem.A, 1);
|
luis@199
|
64 end
|
luis@199
|
65 [y, numiter, time, y_path] = wrapper_ALPS_toolbox(b, A, solver.param);
|
luis@199
|
66 elseif (strcmpi(solver.toolbox, 'MMbox'))
|
luis@199
|
67 if ~isa(Problem.A,'float')
|
luis@199
|
68 % MMbox does not accept implicit dictionary definition
|
luis@199
|
69 A = opToMatrix(Problem.A, 1);
|
luis@199
|
70 end
|
luis@199
|
71
|
luis@199
|
72 [y, cost] = wrapper_mm_solver(b, A, solver.param);
|
luis@199
|
73
|
luis@199
|
74 % To introduce new sparse representation algorithm put the files in
|
luis@199
|
75 % your Matlab path. Next, unique name <TolboxID> for your toolbox and
|
luis@199
|
76 % prefferd API <Preffered_API> needs to be defined.
|
luis@199
|
77 %
|
luis@199
|
78 % elseif strcmpi(solver.toolbox,'<ToolboxID>')
|
luis@199
|
79 %
|
luis@199
|
80 % % - Evaluate the function (solver.name - defined in the main) with
|
luis@199
|
81 % % parameters given above
|
luis@199
|
82 %
|
luis@199
|
83 % y = eval([solver.name,'(<Preffered_API>);']);
|
luis@199
|
84
|
luis@199
|
85 else
|
luis@199
|
86 printf('\nToolbox has not been registered. Please change SMALL_solver file.\n');
|
luis@199
|
87 return
|
luis@199
|
88 end
|