view tests/RecomTSTgentest.m @ 18:a8ff9a881d2f

GAP test almost working. For some data the results are not the same because of representation error, so the test doesn't fully work for now. But the results seem to be accurate.

% Run BP and save parameters and solutions as reference test data
% to check if other algorithms are correct

% numA = 10;
% numY = 100;
% 
% sizesA{1} = [50 100];
% sizesA{2} = [20 25];
% sizesA{3} = [10 120];
% sizesA{4} = [15 100];
% sizesA{5} = [70 100];
% sizesA{6} = [80 100];
% sizesA{7} = [90 100];
% sizesA{8} = [99 100];
% sizesA{9} = [100 100];
% sizesA{10} = [250 400];
% 
% for i = 1:numA
%     sz = sizesA{i};
%     cellA{i} = randn(sz);
%     cellY{i} = randn(sz(1), numY);
%     for j = 1:numY
%         cellTol{i}(j) = rand / 5; % restrict from 0 to 20% if measurements
%         %cellX0{i}(:,j) = cellA{i} \ cellY{i}(:,j);
%         cellX0{i}(:,j) = zeros(size(cellA{i},2),1);
%     end
% end
% %load BPtestdata
tic
for i = 1:numA
    for j = 1:numY
        cellXr{i}(:,j) = RecommendedTST(cellA{i}, cellY{i}(:,j), 300, cellTol{i}(j), cellX0{i}(:,j));
    end
end
toc

% save RecomTSTtestdata