Mercurial > hg > pycsalgos
annotate tests/SL0approxgentest.m @ 39:0ff08ae833be
utils.py:
Mark the areas with less than 10% error in the figures.
Joined loadsavematrices() and loadshowmatrices() together in the same function loadmatrices()
author | nikcleju |
---|---|
date | Fri, 18 Nov 2011 15:12:34 +0000 |
parents | 66d8654fb73e |
children |
rev | line source |
---|---|
nikcleju@11 | 1 % Run SL0 and save parameters and solutions as reference test data |
nikcleju@11 | 2 % to check if other algorithms are correct |
nikcleju@11 | 3 |
nikcleju@11 | 4 numA = 10; |
nikcleju@11 | 5 numY = 100; |
nikcleju@11 | 6 |
nikcleju@11 | 7 sizesA{1} = [50 100]; |
nikcleju@11 | 8 sizesA{2} = [20 25]; |
nikcleju@11 | 9 sizesA{3} = [10 120]; |
nikcleju@11 | 10 sizesA{4} = [15 100]; |
nikcleju@11 | 11 sizesA{5} = [70 100]; |
nikcleju@11 | 12 sizesA{6} = [80 100]; |
nikcleju@11 | 13 sizesA{7} = [90 100]; |
nikcleju@11 | 14 sizesA{8} = [99 100]; |
nikcleju@11 | 15 sizesA{9} = [100 100]; |
nikcleju@11 | 16 sizesA{10} = [250 400]; |
nikcleju@11 | 17 |
nikcleju@11 | 18 sigmamin = [0.00001 0.01 0.2 0.3 0.4 0.0001 0.1 0.001 0.1 0.1]; |
nikcleju@11 | 19 |
nikcleju@11 | 20 for i = 1:numA |
nikcleju@11 | 21 sz = sizesA{i}; |
nikcleju@11 | 22 cellA{i} = randn(sz); |
nikcleju@11 | 23 cellY{i} = randn(sz(1), numY); |
nikcleju@11 | 24 for j = 1:numY |
nikcleju@11 | 25 cellEps{i}(j) = rand / 100; % restrict from 0 to 1% io measurements |
nikcleju@11 | 26 end |
nikcleju@11 | 27 end |
nikcleju@11 | 28 |
nikcleju@11 | 29 %load SL0approxtestdata |
nikcleju@11 | 30 tic |
nikcleju@11 | 31 for i = 1:numA |
nikcleju@11 | 32 for j = 1:numY |
nikcleju@11 | 33 cellXr{i}(:,j) = SL0_approx(cellA{i}, cellY{i}(:,j), cellEps{i}(j), sigmamin(i)); |
nikcleju@11 | 34 end |
nikcleju@11 | 35 end |
nikcleju@11 | 36 toc |
nikcleju@11 | 37 |
nikcleju@11 | 38 save SL0approxtestdata |