How to compute fast?
5 ビュー (過去 30 日間)
古いコメントを表示
Mohammod Minhajur Rahman
2019 年 2 月 28 日
コメント済み: Mohammod Minhajur Rahman
2019 年 2 月 28 日
Hi, I have the following code that I need to vectorize to compute faster:
Here,
neq = 3197;
length(K1) = 152;
d1 is such d1{1,i} = 3197 by 1;
delNonZero = 152 by 1;
cd = zeros(neq,neq);
for i = 1:length(K1)
for j = 1:length(K1)
cd =cd+ d1{1,i}*d1{1,j}'*mean([delNonZero(i) delNonZero(j)]);
end
end
4 件のコメント
Bob Thompson
2019 年 2 月 28 日
Got it. I am not an expert in those things, so hopefully somebody else can offer you a better solution. I will say that I don't know if you can eliminate the loops because of your desire to look at just a few elements at a time, rather than just looking at entire rows or the like.
採用された回答
Jos (10584)
2019 年 2 月 28 日
編集済み: Jos (10584)
2019 年 2 月 28 日
Some suggestions:
- replace mean(A,B) by (A+B)/2
- you can have j run from i to length(K1), since everything seems symmetric (unless there are complex numbers involved):
for i = 1:..
cd = cd + d1{i,1}*d1{i,1}'*delNonZero(i) ;
for j = i+1: ..
cd = cd + 2*d1{i,1}*d1{j,1}'*(delNonZero(i)+delNonZero(j))/2; % you can remove the 2's !!
end
end
- do not use the function CD as avariable name ...
0 件のコメント
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Signal Generation and Preprocessing についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!