how to vectorize this for loops?

2 次查看（过去 30 天）
Siddhesh Karbhari
Siddhesh Karbhari 2018-1-10
编辑： per isakson 2018-1-13
  4 个评论
显示 2更早的评论
Siddhesh Karbhari
Siddhesh Karbhari 2018-1-10
编辑：per isakson 2018-1-13
clc,
octave=3;
lavel=3;
sigma0=sqrt(2);
for i=1:octave
temp_D=D{i};
for j=1:level
scale=sigma0*sqrt(2)^(1/level)^((i-1)*level+j)
p=(level)*(i-1);
figure(1);
subplot(octave,level,p+j);
f=fspecial('gaussian',[1,floor(6*scale)],scale);
L1=temp_img;
if(i==1&&j==1)
L2=conv2(temp_img,f,'same');
L2=conv2(L2,f','same');
temp_D(:,:,j)=L2-L1;
imshow(uint8(255 * mat2gray(temp_D(:,:,j))));
L1=L2;
else
L2=conv2(temp_img,f,'same');
L2=conv2(L2,f','same');
temp_D(:,:,j)=L2-L1;
L1=L2;
if(j==level)
temp_img=L1(2:end-1,2:end-1);
end
imshow(uint8(255 * mat2gray(temp_D(:,:,j))));
end
end
D{i}=temp_D;
temp_img=temp_img(1:2:end,1:2:end);
temp_img=padarray(temp_img,[1,1],'both','replicate');
end
Jan
Jan 2018-1-10
By the way:
% Simpler:
% sigma0 * sqrt(2) ^ (1/level) ^ ((i-1)*level+j)
sigma0 * sqrt(2) ^ (i - 1 + j / level)
With a proper indentation the code looks cleaner: Ctrl-A Ctrl-I

请先登录，再进行评论。

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

回答（1 个）

Jan
Jan 2018-1-10
Neither the inner nor the outer loop can be vectorized. The two imshow commands are impossible to vectorize. Do you really need them? Do you want a vectorization to increase the speed? Then start with using the profiler to find the most time consuming part of the code.
  3 个评论
显示 1更早的评论
Jan
Jan 2018-1-10
What do you want to implement on GPU exactly? imshow, subplot, fspecial are not suitable for this, but conv2 will be fine.
Siddhesh Karbhari
Siddhesh Karbhari 2018-1-10
编辑：Siddhesh Karbhari 2018-1-10
I want to implement SIFT descriptor on GPU.

请先登录，再进行评论。

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by