nikcleju@14
|
1 # -*- coding: utf-8 -*-
|
nikcleju@14
|
2 """
|
nikcleju@14
|
3 Created on Wed Dec 07 14:04:40 2011
|
nikcleju@14
|
4
|
nikcleju@14
|
5 @author: ncleju
|
nikcleju@14
|
6 """
|
nikcleju@14
|
7
|
nikcleju@14
|
8 import numpy
|
nikcleju@14
|
9 from algos import *
|
nikcleju@14
|
10
|
nikcleju@15
|
11 paramstest = dict()
|
nikcleju@15
|
12 paramstest['algos'] = exact_gap,exact_sl0,exact_bp,exact_ompeps,exact_tst # tuple of algorithms
|
nikcleju@15
|
13 #paramstest['algos'] = exact_bp_cvxopt, # tuple of algorithms
|
nikcleju@15
|
14 paramstest['d'] = 50.0
|
nikcleju@15
|
15 paramstest['sigma'] = 1.2
|
nikcleju@15
|
16 paramstest['deltas'] = numpy.array([0.05, 0.45, 0.95])
|
nikcleju@15
|
17 paramstest['rhos'] = numpy.array([0.05, 0.45, 0.95])
|
nikcleju@15
|
18 #deltas = numpy.array([0.6])
|
nikcleju@15
|
19 #deltas = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
20 #rhos = numpy.array([0.05])
|
nikcleju@15
|
21 paramstest['numvects'] = 10; # Number of vectors to generate
|
nikcleju@15
|
22 paramstest['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
|
nikcleju@15
|
23 paramstest['savedataname'] = 'exact_pt_stdtest.mat'
|
nikcleju@15
|
24 paramstest['saveplotbase'] = 'exact_pt_stdtest_'
|
nikcleju@15
|
25 paramstest['saveplotexts'] = ('png','pdf','eps')
|
nikcleju@14
|
26
|
nikcleju@14
|
27 # Standard parameters 1
|
nikcleju@14
|
28 # All algorithms, 100 vectors
|
nikcleju@14
|
29 # d=50, sigma = 2, delta and rho full resolution (0.05 step), lambdas = 0, 1e-4, 1e-2, 1, 100, 10000
|
nikcleju@14
|
30 # Do save data, do save plots, don't show plots
|
nikcleju@15
|
31 params1 = dict()
|
nikcleju@15
|
32 params1['algos'] = exact_gap,exact_sl0,exact_bp_cvxopt,exact_ompeps,exact_tst # tuple of algorithms
|
nikcleju@15
|
33 params1['d'] = 50.0;
|
nikcleju@15
|
34 params1['sigma'] = 1.2
|
nikcleju@15
|
35 params1['deltas'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
36 params1['rhos'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
37 params1['numvects'] = 100; # Number of vectors to generate
|
nikcleju@15
|
38 params1['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
|
nikcleju@15
|
39 params1['savedataname'] = 'exact_pt_std1.mat'
|
nikcleju@15
|
40 params1['saveplotbase'] = 'exact_pt_std1_'
|
nikcleju@15
|
41 params1['saveplotexts'] = ('png','pdf','eps')
|
nikcleju@14
|
42
|
nikcleju@15
|
43
|
nikcleju@14
|
44 # Standard parameters 2
|
nikcleju@14
|
45 # All algorithms, 100 vectors
|
nikcleju@14
|
46 # d=20, sigma = 10, delta and rho full resolution (0.05 step), lambdas = 0, 1e-4, 1e-2, 1, 100, 10000
|
nikcleju@14
|
47 # Do save data, do save plots, don't show plots
|
nikcleju@15
|
48 params2 = dict()
|
nikcleju@15
|
49 params2['algos'] = exact_gap,exact_sl0,exact_bp_cvxopt,exact_ompeps,exact_tst # tuple of algorithms
|
nikcleju@15
|
50 params2['d'] = 20.0
|
nikcleju@15
|
51 params2['sigma'] = 10.0
|
nikcleju@15
|
52 params2['deltas'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
53 params2['rhos'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
54 params2['numvects'] = 100; # Number of vectors to generate
|
nikcleju@15
|
55 params2['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
|
nikcleju@15
|
56 params2['savedataname'] = 'exact_pt_std2.mat'
|
nikcleju@15
|
57 params2['saveplotbase'] = 'exact_pt_std2_'
|
nikcleju@15
|
58 params2['saveplotexts'] = ('png','pdf','eps') |