nikcleju@3: # -*- coding: utf-8 -*-
nikcleju@3: """
nikcleju@17: Test SL0 behavior in the presence of strong noise
nikcleju@3: 
nikcleju@17: Author: Nicolae Cleju
nikcleju@3: """
nikcleju@3: 
nikcleju@3: from algos import *
nikcleju@3: import ABSapprox
nikcleju@3: import stdparams
nikcleju@3: 
nikcleju@3: 
nikcleju@3: def run_mp():
nikcleju@3:     
nikcleju@3:   SNRs = numpy.array([4., 7., 10., 13., 16., 19.])
nikcleju@3:   
nikcleju@3:   # Default params
nikcleju@3:   algosN,algosL,d,sigma,deltas,rhos,lambdas,numvects,SNRdb,dosavedata,\
nikcleju@3:     savedataname,doshowplot,dosaveplot,saveplotbase,saveplotexts = stdparams.std1()
nikcleju@3:     
nikcleju@3:   # Set only SL0 algorithms
nikcleju@3:   algosN = (sl0analysis,)
nikcleju@3:   algosL = (sl0,)
nikcleju@3:   
nikcleju@3:   # Overwrite SNRdb:
nikcleju@3:   for SNRdb in SNRs:
nikcleju@4:     savedataname = 'approx_pt_sl0_SNR'+ str(int(SNRdb)) +'.mat'
nikcleju@4:     saveplotbase = 'approx_pt_sl0_SNR'+ str(int(SNRdb)) +'_'
nikcleju@3:         
nikcleju@3:     # Run
nikcleju@3:     ABSapprox.run_multi(algosN,algosL,d,sigma,deltas,rhos,lambdas,numvects,SNRdb,\
nikcleju@3:       dosavedata=dosavedata,savedataname=savedataname,doparallel=True, ncpus=None,\
nikcleju@3:       doshowplot=doshowplot,dosaveplot=dosaveplot,saveplotbase=saveplotbase,saveplotexts=saveplotexts)
nikcleju@3:   
nikcleju@3: # Script main
nikcleju@3: if __name__ == "__main__":
nikcleju@3:   #import cProfile
nikcleju@3:   #cProfile.run('mainrun()', 'profile')    
nikcleju@3:   run_mp()
nikcleju@3:   #runsingleexampledebug()