view tests/RecomTSTgentest.m @ 25:dd0e78b5bb13

Added .squeeze() in GAP function to avoid strange error in numpy.delete(), which wasn't raised on my laptop but was raised on octave

% 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