What is the fastest, circshift or indices ?

Question

Stéphane 2014-11-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/161447-what-is-the-fastest-circshift-or-indices

编辑： Matt J 2014-11-5

采纳的回答： Matt J

Hi

I'm trying to improve the efficiency of my algorithm. I need to compute something like :

Y(2:end-1,:) = X(2:end-1,:) + X(1:end-2,:) - X(3:end,:)

For Y and X n by m matrices. I was wondering if it was faster to do it as it is above, or using circshift :

Y = X + circshift(X,[-1,0]) - circshift(X,[+1,0]).

I've done a few tests but I can't get a clear idea... The aim is to do this on pretty big matrices on a cluser with multiple cores...

Any idea ? Thanks.

Stéphane

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Adam 2014-11-5

在 MATLAB Online 中打开

I regularly create small test scripts with as many alternative implementations as I can think of, each wrapped in a

timeit(f)

instruction for function handle, f. It is extremely useful for improving understanding as well as getting the (hopefully conclusive) final answer as to which method is fastest. Sounds like you did kind of do that though...

请先登录，再进行评论。

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

Answer 1

Matt J 2014-11-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/161447-what-is-the-fastest-circshift-or-indices#answer_157817

编辑：Matt J 2014-11-5

在 MATLAB Online 中打开

Just use conv2

Y=conv2(X,[-1;1;1],'valid');

As for circshift, it is not a builtin function and, if you look inside it

>>type circshift

you will see that it does pretty much the same indexing as you are doing. So, I wouldn't expect it to be a competitive approach.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Stéphane 2014-11-5

Second thought...

I've tried in 1D to understand what conv (conv2) does and it seems conv(X,[-1;1]) and the same with circshift doesn't produce the same result. Is it normal or am I doing something wrong ?

(What I want is X(2:end) - X(1:end-1) )