|
nikcleju@10
|
1 # -*- coding: utf-8 -*-
|
|
nikcleju@10
|
2 """
|
|
nikcleju@10
|
3 Created on Sat Nov 05 18:08:40 2011
|
|
nikcleju@10
|
4
|
|
nikcleju@10
|
5 @author: Nic
|
|
nikcleju@10
|
6 """
|
|
nikcleju@10
|
7
|
|
nikcleju@10
|
8 import numpy
|
|
nikcleju@10
|
9 import pyCSalgos
|
|
nikcleju@10
|
10
|
|
nikcleju@10
|
11 def gap_paramsetup(y,M,Omega,epsilon,lbd):
|
|
nikcleju@10
|
12 gapparams = dict(num_iteration = 1000,
|
|
nikcleju@10
|
13 greedy_level = 0.9,
|
|
nikcleju@10
|
14 stopping_coefficientstopping_coefficient_size = 1e-4,
|
|
nikcleju@10
|
15 l2solver = 'pseudoinverse',
|
|
nikcleju@10
|
16 noise_level = epsilon)
|
|
nikcleju@10
|
17 return y,M,M.T,Omega,Omega.T,gapparams,numpy.zeros(Omega.shape[1])
|
|
nikcleju@10
|
18
|
|
nikcleju@10
|
19 def omp_paramsetup(y,M,Omega,epsilon,lbd):
|
|
nikcleju@10
|
20 gapparams = dict(num_iteration = 1000,
|
|
nikcleju@10
|
21 greedy_level = 0.9,
|
|
nikcleju@10
|
22 stopping_coefficientstopping_coefficient_size = 1e-4,
|
|
nikcleju@10
|
23 l2solver = 'pseudoinverse',
|
|
nikcleju@10
|
24 noise_level = epsilon)
|
|
nikcleju@10
|
25 return y,M,M.T,Omega,Omega.T,gapparams,numpy.zeros(Omega.shape[1])
|
|
nikcleju@10
|
26
|
|
nikcleju@10
|
27 gap = (pyCSalgos.GAP, gap_paramsetup)
|
|
nikcleju@10
|
28
|
|
nikcleju@10
|
29
|
|
nikcleju@10
|
30
|
|
nikcleju@10
|
31 gap = (pyCSalgos.GAP, gap_paramsetup)
|
|
nikcleju@10
|
32
|
|
nikcleju@10
|
33
|
|
nikcleju@10
|
34
|
|
nikcleju@10
|
35
|
|
nikcleju@10
|
36 def mainrun():
|
|
nikcleju@10
|
37
|
|
nikcleju@10
|
38 algos = (gap, sl0)
|
|
nikcleju@10
|
39
|
|
nikcleju@10
|
40 for algofunc,paramsetup in algos:
|
|
nikcleju@10
|
41 xrec = algofunc(algosetup(y, Omega, epsilon, lbd))
|
