Не улучшает ли Matlab следующее?

Question

У меня очень длинный вектор 1xr v, и очень длинный вектор w 1xs и матрица A rxs, которая является редкой (но очень большой по размеру).

Я ожидал следующее быть оптимизированы Matlab, так что я не столкнуться с проблемами с памятью:

A./(v'*w) 

но похоже, Matlab на самом деле пытается создать полный v'*w матрицу, потому что я проблемы с памятью. Есть ли способ преодолеть это? Обратите внимание, что нет необходимости вычислять все v'*w, потому что многие значения A составляют 0.

EDIT: Если бы это было возможно, один из способов сделать это будет делать A(find(A))./(v'*w)(find(A));

, но вы не можете выбрать подмножество матрицы (v'*w в данном случае) без первого расчета его и положить его в Переменная.

Answer 1

• Вы можете использовать bsxfun. Это дает тот же результат, как и A./(v'*w) без генерации матрицы v.'*w:

  bsxfun(@rdivide, bsxfun(@rdivide, 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).').