次のコードがありますが、n>10^6 に対して実行する必要があります。n=10^5 の場合、記憶は終わりです。私は16GBを持っています。これを行う効率的な方法はありますか?
n = 10^6;
m = n/2;
H = sparse(m,n);
parfor jj=1:n
pos = sample(cumulative);
H(:,jj) = (rand(m,1) > 1 - degrees(pos)/m);
end