view SL0_pt_noise.py @ 21:d395461b92ae tip

Lots and lots of modifications. Approximate recovery script working.

# -*- coding: utf-8 -*-
"""
Test SL0 behavior in the presence of strong noise

Author: Nicolae Cleju
"""

from algos import *
import ABSapprox
import stdparams


def run_mp():
    
  SNRs = numpy.array([4., 7., 10., 13., 16., 19.])
  
  # Default params
  algosN,algosL,d,sigma,deltas,rhos,lambdas,numvects,SNRdb,dosavedata,\
    savedataname,doshowplot,dosaveplot,saveplotbase,saveplotexts = stdparams.std1()
    
  # Set only SL0 algorithms
  algosN = (sl0analysis,)
  algosL = (sl0,)
  
  # Overwrite SNRdb:
  for SNRdb in SNRs:
    savedataname = 'approx_pt_sl0_SNR'+ str(int(SNRdb)) +'.mat'
    saveplotbase = 'approx_pt_sl0_SNR'+ str(int(SNRdb)) +'_'
        
    # Run
    ABSapprox.run_multi(algosN,algosL,d,sigma,deltas,rhos,lambdas,numvects,SNRdb,\
      dosavedata=dosavedata,savedataname=savedataname,doparallel=True, ncpus=None,\
      doshowplot=doshowplot,dosaveplot=dosaveplot,saveplotbase=saveplotbase,saveplotexts=saveplotexts)
  
# Script main
if __name__ == "__main__":
  #import cProfile
  #cProfile.run('mainrun()', 'profile')    
  run_mp()
  #runsingleexampledebug()