nikcleju@8: % Run SL0 and save parameters and solutions as reference test data
nikcleju@8: % to check if other algorithms are correct
nikcleju@8: 
nikcleju@8: numA = 10;
nikcleju@8: numY = 100;
nikcleju@8: 
nikcleju@8: sizesA{1} = [50 100];
nikcleju@8: sizesA{2} = [20 25];
nikcleju@8: sizesA{3} = [10 120];
nikcleju@8: sizesA{4} = [15 100];
nikcleju@8: sizesA{5} = [70 100];
nikcleju@8: sizesA{6} = [80 100];
nikcleju@8: sizesA{7} = [90 100];
nikcleju@8: sizesA{8} = [99 100];
nikcleju@8: sizesA{9} = [100 100];
nikcleju@8: sizesA{10} = [250 400];
nikcleju@8: 
nikcleju@8: sigmamin = [0.00001 0.01 0.2 0.3 0.4 0.0001 0.1 0.001 0.1 0.1];
nikcleju@8: 
nikcleju@8: for i = 1:numA
nikcleju@8:     sz = sizesA{i};
nikcleju@8:     cellA{i} = randn(sz);
nikcleju@8:     cellY{i} = randn(sz(1), numY);
nikcleju@8: end
nikcleju@8: 
nikcleju@8: %load SL0testdata
nikcleju@8: tic
nikcleju@8: for i = 1:numA
nikcleju@8:     for j = 1:numY
nikcleju@8:         cellXr{i}(:,j) = SL0(cellA{i}, cellY{i}(:,j), sigmamin(i));
nikcleju@8:     end
nikcleju@8: end
nikcleju@8: toc
nikcleju@8: 
nikcleju@8: save SL0testdata