Is there a way to vectorise this vector/matrix product?

1 次查看(过去 30 天)
I have a piece of code inside a for loop, which looks like
for k = 1:size(G, 1)
G(k) = real(W(:,k)' * R * W(:,k))
end
G is preallocated (which is how I know the length of it). W and R are 2D matrices (so W(:,k) is a vector). This loop is a real bottleneck in my code, that one line taking 90% of the total time that it takes to run a script which itself calls other functions and files. It's not possible to use a parfor loop as I'm already using one in a wrapper function.
Is there a way to vectorise that loop, in order to reduce the number of times the calculation is done? I'm sure it can be done, but it's not at all obvious to me.

采纳的回答

Matt J
Matt J 2021-5-7
编辑:Matt J 2021-5-7
G=real( sum((R*W).*conj(W),1) );
  4 个评论
Matt J
Matt J 2021-5-7
Here is a test of a agreement with random data:
N=5;
W=complex(rand(N),rand(N));
R=complex(rand(N),rand(N));
%%% Original method
G=nan(N,1);
for k = 1:size(G, 1)
G(k) = real(W(:,k)' * R * W(:,k));
end
Gbase=G;
%%% Proposed alternative
G=real( sum((R*W).*conj(W),1) ).';
Discrepancy=norm(G-Gbase,'fro')
Discrepancy = 1.9860e-15
Daniel Pollard
Daniel Pollard 2021-5-7
You're right, and I just came back to your answer to say that. Thank you!

请先登录,再进行评论。

更多回答(1 个)

Bruno Luong
Bruno Luong 2021-5-7
G = real(dot(W,R*W,1)).';

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

产品


版本

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by