スケールMATLAB

Question

のベクトルの要素によって、スパース行列のすべての行は、私は疎行列

obj.resOp = sparse(row,col,val); 

とすぐ

sums = sparse(sum(obj.resOp,2)); 

行列の各列の和を含むベクターを有します私は何をしたいです

obj.resOp = obj.resOp ./ sums; 

各行のrowsumが1になるように行列のすべての行をスケールするでしょう

は、しかし、この最後の行では、MATLABは内部obj.resOpからフル行列を構築しているようですので、私はこのエラーを取得する：十分に大きな行列のため

Error using ./ Requested 38849x231827 (17.5GB) array exceeds maximum array size preference. Creation of arrays greater than this limit may take a long time and cause MATLAB to become unresponsive. See array size limit or preference panel for more information.

。

理論的には、フルマトリクスに拡張する必要はないと私は考えています。 obj.resOpの希薄さを維持しながら、達成したいMATLABの定式化はありますか？

Answer 1

これは、described in this answerと同様の方法で実行できます。いくつかの疎行列

% Random sparse matrix: 10 rows, 4 cols, density 20% 
S = sprand(10,4, 0.2); 

sumがまばらな入力からスパース行列を返すことに注意して、行の合計を取得しますので、あなたの追加の変換（docs）を必要とせず、

スタート。

rowsums = sum(S,2); 

今要素単位分割からスパース行列を作成

out = sparse(rowidx, colidx, vals./rowsums(rowidx), size(S,1), size(S,2)); 

Answer 2

にすべての非ゼロインデックスとその値

[rowidx, colidx, vals] = find(S) 

等価計算

obj.resOp = inv(diag(sums)) * obj.resOp; 

を探します

はスムーズに動作します。