optimize inefficient piece of code

4 次查看（过去 30 天）
Wesley Ooms
Wesley Ooms 2014-1-6
评论： Simon 2014-1-7
Hi,
I have a piece of code as part of a bigger piece that is inefficient. The code is called many times and is a bottleneck in the proces also according to the matlab profiler. I feel it can be done faster. The code looks as follows
clear all
% initialize variables of interest
l1=rand(1,10);
l2=rand(1,10);
p1=rand(3,10);
p2=rand(3,10);
% start inefficient code
[c,ind]=min([l1;l2]);
r(:,:,1)=p1;
r(:,:,2)=p2;
for it=1:10
p3(:,it)=r(:,it,ind(it));
end
I would like to do the for loop in one go (so remove the for loop). l1 and l2 are distances (l1 is norm of p1 and l2 is norm of p2), and p1, p2 are 3d coordinates. The goal of the code is to pick the coordinates of p1, p2 from the one with the smallest distance.
  1 个评论
Wesley Ooms
Wesley Ooms 2014-1-6
Ican do something like this:
p31=p2.*([ind;ind;ind]-1)+p1.*~([ind;ind;ind]-1)
which results in a correct answer but it feels like a waste of memory and a lot of unnecessary operations.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Simon
Simon 2014-1-6
Hi!
You can of course avoid the loop with
p3(1:3, ind==1) = p1(:, ind==1);
p3(1:3, ind==2) = p2(:, ind==2);
This avoids the r-lines before the loop as well. Use '~' instead of 'c' if you are not interessted in the value of the 'min' call.
  4 个评论
显示 2更早的评论
Wesley Ooms
Wesley Ooms 2014-1-6
编辑：Wesley Ooms 2014-1-6
I run the (complete) code in profile viewer. I also use tic toc in a for loop. In both cases your code is indeed the fastest. Strangely the ~ instead of c makes the code slower.
I optimized allready quite a bit. Matlab functions i substituted to decrease cputime are for example circshift, repmat, cross, dot, mean. Why are these handy and very common functions slower than just straightforward implementing raw code?
Simon
Simon 2014-1-7
I think because of the overhead, like checking the input for correct type and dimension and things like that. But that implies that you know what you do and what kind of input to your substitutes you have.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Code Performance 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by