annotate stdparams_exact.py @ 21:d395461b92ae tip

Lots and lots of modifications. Approximate recovery script working.
author Nic Cleju <nikcleju@gmail.com>
date Mon, 23 Apr 2012 10:54:57 +0300
parents 4a967f4f18a0
children
rev   line source
nikcleju@14 1 # -*- coding: utf-8 -*-
nikcleju@14 2 """
nikcleju@17 3 Defines standard parameters for exact reconstruction simulation
nikcleju@17 4 Author: Nicolae Cleju
nikcleju@17 5 """
nikcleju@17 6 __author__ = "Nicolae Cleju"
nikcleju@17 7 __license__ = "GPL"
nikcleju@17 8 __email__ = "nikcleju@gmail.com"
nikcleju@14 9
nikcleju@14 10
nikcleju@14 11 import numpy
nikcleju@17 12
nikcleju@17 13 # Solver algorithms to run
nikcleju@14 14 from algos import *
nikcleju@14 15
nikcleju@17 16
nikcleju@17 17 # Test parameters
nikcleju@15 18 paramstest = dict()
nikcleju@15 19 paramstest['algos'] = exact_gap,exact_sl0,exact_bp,exact_ompeps,exact_tst # tuple of algorithms
nikcleju@15 20 #paramstest['algos'] = exact_bp_cvxopt, # tuple of algorithms
nikcleju@18 21 paramstest['d'] = 200.0
nikcleju@15 22 paramstest['sigma'] = 1.2
nikcleju@15 23 paramstest['deltas'] = numpy.array([0.05, 0.45, 0.95])
nikcleju@15 24 paramstest['rhos'] = numpy.array([0.05, 0.45, 0.95])
nikcleju@15 25 #deltas = numpy.array([0.6])
nikcleju@15 26 #deltas = numpy.arange(0.05,1.,0.05)
nikcleju@15 27 #rhos = numpy.array([0.05])
nikcleju@15 28 paramstest['numvects'] = 10; # Number of vectors to generate
nikcleju@15 29 paramstest['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
nikcleju@15 30 paramstest['savedataname'] = 'exact_pt_stdtest.mat'
nikcleju@15 31 paramstest['saveplotbase'] = 'exact_pt_stdtest_'
nikcleju@15 32 paramstest['saveplotexts'] = ('png','pdf','eps')
nikcleju@14 33
nikcleju@14 34 # Standard parameters 1
nikcleju@14 35 # All algorithms, 100 vectors
nikcleju@17 36 # d = 200, sigma = 1.2, delta and rho full resolution (0.05 step)
nikcleju@17 37 # Virtually no noise (100db)
nikcleju@15 38 params1 = dict()
nikcleju@15 39 params1['algos'] = exact_gap,exact_sl0,exact_bp_cvxopt,exact_ompeps,exact_tst # tuple of algorithms
nikcleju@17 40 params1['d'] = 200.0;
nikcleju@15 41 params1['sigma'] = 1.2
nikcleju@15 42 params1['deltas'] = numpy.arange(0.05,1.,0.05)
nikcleju@15 43 params1['rhos'] = numpy.arange(0.05,1.,0.05)
nikcleju@15 44 params1['numvects'] = 100; # Number of vectors to generate
nikcleju@15 45 params1['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
nikcleju@15 46 params1['savedataname'] = 'exact_pt_std1.mat'
nikcleju@15 47 params1['saveplotbase'] = 'exact_pt_std1_'
nikcleju@15 48 params1['saveplotexts'] = ('png','pdf','eps')
nikcleju@14 49
nikcleju@15 50
nikcleju@14 51 # Standard parameters 2
nikcleju@14 52 # All algorithms, 100 vectors
nikcleju@17 53 # d = 20, sigma = 10, delta and rho full resolution (0.05 step)
nikcleju@17 54 # Virtually no noise (100db)
nikcleju@15 55 params2 = dict()
nikcleju@15 56 params2['algos'] = exact_gap,exact_sl0,exact_bp_cvxopt,exact_ompeps,exact_tst # tuple of algorithms
nikcleju@15 57 params2['d'] = 20.0
nikcleju@15 58 params2['sigma'] = 10.0
nikcleju@15 59 params2['deltas'] = numpy.arange(0.05,1.,0.05)
nikcleju@15 60 params2['rhos'] = numpy.arange(0.05,1.,0.05)
nikcleju@15 61 params2['numvects'] = 100; # Number of vectors to generate
nikcleju@15 62 params2['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
nikcleju@15 63 params2['savedataname'] = 'exact_pt_std2.mat'
nikcleju@15 64 params2['saveplotbase'] = 'exact_pt_std2_'
nikcleju@15 65 params2['saveplotexts'] = ('png','pdf','eps')