nikcleju@4: % Run BP and save parameters and solutions as reference test data
nikcleju@4: % to check if other algorithms are correct
nikcleju@4: 
nikcleju@5: % numA = 10;
nikcleju@5: % numY = 100;
nikcleju@5: % 
nikcleju@5: % sizesA{1} = [50 100];
nikcleju@5: % sizesA{2} = [20 25];
nikcleju@5: % sizesA{3} = [10 120];
nikcleju@5: % sizesA{4} = [15 100];
nikcleju@5: % sizesA{5} = [70 100];
nikcleju@5: % sizesA{6} = [80 100];
nikcleju@5: % sizesA{7} = [90 100];
nikcleju@5: % sizesA{8} = [99 100];
nikcleju@5: % sizesA{9} = [100 100];
nikcleju@5: % sizesA{10} = [250 400];
nikcleju@5: % 
nikcleju@5: % for i = 1:numA
nikcleju@5: %     sz = sizesA{i};
nikcleju@5: %     cellA{i} = randn(sz);
nikcleju@5: %     cellY{i} = randn(sz(1), numY);
nikcleju@5: %     for j = 1:numY
nikcleju@5: %         cellTol{i}(j) = rand / 5; % restrict from 0 to 20% if measurements
nikcleju@5: %         %cellX0{i}(:,j) = cellA{i} \ cellY{i}(:,j);
nikcleju@5: %         cellX0{i}(:,j) = zeros(size(cellA{i},2),1);
nikcleju@5: %     end
nikcleju@5: % end
nikcleju@5: % %load BPtestdata
nikcleju@5: tic
nikcleju@4: for i = 1:numA
nikcleju@4:     for j = 1:numY
nikcleju@4:         cellXr{i}(:,j) = RecommendedTST(cellA{i}, cellY{i}(:,j), 300, cellTol{i}(j), cellX0{i}(:,j));
nikcleju@4:     end
nikcleju@4: end
nikcleju@5: toc
nikcleju@4: 
nikcleju@5: % save RecomTSTtestdata