Mercurial > hg > pycsalgos
annotate tests/l1eq_gentest.m @ 68:cab8a215f9a1 tip
Minor
author | Nic Cleju <nikcleju@gmail.com> |
---|---|
date | Tue, 09 Jul 2013 14:50:09 +0300 |
parents | e684f76c1969 |
children |
rev | line source |
---|---|
nikcleju@62 | 1 % Run l1eq_pd() and save parameters and solutions as reference test data |
nikcleju@62 | 2 % to check if other algorithms are correct |
nikcleju@62 | 3 |
nikcleju@62 | 4 numA = 10; |
nikcleju@62 | 5 numY = 100; |
nikcleju@62 | 6 |
nikcleju@62 | 7 sizesA{1} = [50 100]; |
nikcleju@62 | 8 sizesA{2} = [20 25]; |
nikcleju@62 | 9 sizesA{3} = [10 120]; |
nikcleju@62 | 10 sizesA{4} = [15 100]; |
nikcleju@62 | 11 sizesA{5} = [70 100]; |
nikcleju@62 | 12 sizesA{6} = [80 100]; |
nikcleju@62 | 13 sizesA{7} = [90 100]; |
nikcleju@62 | 14 sizesA{8} = [99 100]; |
nikcleju@62 | 15 sizesA{9} = [100 100]; |
nikcleju@62 | 16 sizesA{10} = [250 400]; |
nikcleju@62 | 17 |
nikcleju@62 | 18 for i = 1:numA |
nikcleju@62 | 19 sz = sizesA{i}; |
nikcleju@62 | 20 cellA{i} = randn(sz); |
nikcleju@62 | 21 cellY{i} = randn(sz(1), numY); |
nikcleju@62 | 22 for j = 1:numY |
nikcleju@62 | 23 X0{i}(:,j) = cellA{i} \ cellY{i}(:,j); |
nikcleju@62 | 24 end |
nikcleju@62 | 25 |
nikcleju@62 | 26 end |
nikcleju@62 | 27 |
nikcleju@62 | 28 for i = 1:numA |
nikcleju@62 | 29 for j = 1:numY |
nikcleju@62 | 30 Xr{i}(:,j) = l1eq_pd(X0{i}(:,j), cellA{i}, [], cellY{i}(:,j)); |
nikcleju@62 | 31 end |
nikcleju@62 | 32 end |
nikcleju@62 | 33 |
nikcleju@62 | 34 save l1eq_testdata |