nikcleju@10: # -*- coding: utf-8 -*-
nikcleju@10: """
nikcleju@10: Created on Sat Nov 05 18:08:40 2011
nikcleju@10: 
nikcleju@10: @author: Nic
nikcleju@10: """
nikcleju@10: 
nikcleju@10: import numpy
nikcleju@10: import pyCSalgos
nikcleju@10: 
nikcleju@10: def gap_paramsetup(y,M,Omega,epsilon,lbd):
nikcleju@10:   gapparams = dict(num_iteration = 1000,
nikcleju@10:                    greedy_level = 0.9,
nikcleju@10:                    stopping_coefficientstopping_coefficient_size = 1e-4,
nikcleju@10:                    l2solver = 'pseudoinverse',
nikcleju@10:                    noise_level = epsilon)
nikcleju@10:   return y,M,M.T,Omega,Omega.T,gapparams,numpy.zeros(Omega.shape[1])
nikcleju@10: 
nikcleju@10: def omp_paramsetup(y,M,Omega,epsilon,lbd):
nikcleju@10:   gapparams = dict(num_iteration = 1000,
nikcleju@10:                    greedy_level = 0.9,
nikcleju@10:                    stopping_coefficientstopping_coefficient_size = 1e-4,
nikcleju@10:                    l2solver = 'pseudoinverse',
nikcleju@10:                    noise_level = epsilon)
nikcleju@10:   return y,M,M.T,Omega,Omega.T,gapparams,numpy.zeros(Omega.shape[1])
nikcleju@10: 
nikcleju@10: gap = (pyCSalgos.GAP, gap_paramsetup)
nikcleju@10: 
nikcleju@10: 
nikcleju@10: 
nikcleju@10: gap = (pyCSalgos.GAP, gap_paramsetup)
nikcleju@10:   
nikcleju@10: 
nikcleju@10: 
nikcleju@10: 
nikcleju@10: def mainrun():
nikcleju@10:   
nikcleju@10:   algos = (gap, sl0)
nikcleju@10:   
nikcleju@10:   for algofunc,paramsetup in algos:
nikcleju@10:     xrec = algofunc(algosetup(y, Omega, epsilon, lbd))