N = 40; %行数
D = 32; %块数
m = 4;  %块的列数
d = 2;  %块内rank/维度
resulta = zeros(5000,1);
resultb = zeros(5000,1);
%resulta = zeros(5,20,50);    %保留结果
%resultb = zeros(5,20,50);
%按行数循环
A = zeros(N,m*D);
B = zeros(N,d*D);
%按稀疏度循环
n = N;
parfor i = 0:199
    t = floor(i/1000) + 1;
    k = floor(mod(i,1000)/50)+1;
    j = mod(i,50)+1;
    [a,b] = test(t,k,j);
    resulta(i+1) = a;
    resultb(i+1) = b;
end
resulta = reshape(resulta,5,20,50);
resultb = reshape(resultb,5,20,50);
save('gaussnoise.mat');