Efficient Matrix Multiplication
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
I have A(2000x5000). I need to perform the following:
P1 = A(:,1)*A(:,1)';
for i=2:5000
P1 = P1 + AA(:,i)*A(:,i)'
end
What is the most efficient way to do above? It takes so much time to do it right now due to size of the arrays.
3 Kommentare
the cyclist
am 26 Feb. 2011
From his initialization step, I would infer that "AA" is just a typo of "A."
Jan
am 27 Feb. 2011
Just an actually too obvious comment: If AA is not typo, A*A' is not a matching solution. So, Sam Da, we need your help.
Akzeptierte Antwort
the cyclist
am 26 Feb. 2011
P1 = A * A';
On my machine, that cut the execution time from 330 seconds to 1.5. :-)
3 Kommentare
Oleg Komarov
am 26 Feb. 2011
A = rand(10);
P1 = A(:,1)*A(:,1)';
for i=2:10
P1 = P1 + A(:,i)*A(:,i)';
end
P2 = A * A';
abs(P1-P2) < eps*3
Cyclist's method is essentially the same.
James Tursa
am 26 Feb. 2011
Yep. I went back & checked the code I used to double check the result & saw my mistake. Thanks.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Logical finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!