Matlabは以下を最適化していませんか？

Question

私は非常に長いベクトル1xr v、非常に長いベクトルw 1xs、および行列A rxsを持っています。これはスパースですが（寸法は非常に大きいです）。

私はMATLABが最適化されるように、次のように私はメモリとのトラブルに実行されません期待していた：私はので、MATLABは、実際にフルv'*w行列を生成しようとしているよう

A./(v'*w) 

が、それはそうですメモリ不足の問題が発生しています。これを克服する方法はありますか？ Aの多くの値が0であるため、すべてv'*wを計算する必要はありません。

EDIT：それができなかった場合は、それを行うための一つの方法は、A(find(A))./(v'*w)(find(A));

を行うことであろうが、あなたが最初にそれを計算し、それをかけることなく、マトリックス（この場合はv'*w）のサブセットを選択することはできません変数。

Answer 1

• bsxfunを使用できます。

  bsxfun(@rdivide, bsxfun(@rdivide, A, v'), w) 
  

• 別の可能性：これは行列v.'*wを発生させることなくA./(v'*w)と同じ結果を与えるだけゼロ以外の値、使用したい場合：

  [ii jj Anz] = find(A); 
  Anz./v(ii)'./w(jj).' 
  

  これはあなたのA(find(A))./(v'*w)(find(A))に対応する列ベクトルを与えます再びv.'*wを生成せずに。スパース行列A./(v'*w)（ゼロ以外の値の列ベクトルの場合）が必要な場合は、sparse(ii,jj,Anz./v(ii)'./w(jj).')を使用します。