Speed up a program

3 次查看（过去 30 天）
Vlad
Vlad 2014-12-4
评论： Vlad 2014-12-5
It is required to speed up a very simple program:
for i=1:1000
for j=1:1000
A=fun(i, j) % matrix 3x3
D=eig(A);
end
end
Execution time for fun(i, j) is much less then time for eig(A). A mean execution time of whole program is about 40 minutes. I wrote a mex-file for calculating eigenvalues and improved the speed by 4 times. So now the mean time is about 10 minutes. I tried parfor-loop, it did not work. CUDA did not help as well =(
I will be very appreciated for any idea.

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

回答（2 个）

Zoltán Csáti
Zoltán Csáti 2014-12-4
Why didn't parfor and CUDA work? Do you meet the requirements?
  1 个评论
Vlad
Vlad 2014-12-5
I suppose, that matrix 3x3 is very small for gpu-calculations, that's why there is no advantage of using it.

请先登录，再进行评论。

Edric Ellis
Edric Ellis 2014-12-5
Here's what I tried with PARFOR, and this worked:
function D = pfeg
N = 100;
D = zeros(N, N, 3);
parfor i = 1:N
for j = 1:N
A = fun(i, j);
D(i, j, :) = eig(A);
end
end
end
function A = fun(i, j)
A = rand(3) .* i + j;
end
I'm not surprised that the GPU didn't help much - 3x3 matrices are simply too small to take good advantage of the GPU capabilities (unless you can use pagefun - which unfortunately you cannot in this case).
  1 个评论
Vlad
Vlad 2014-12-5
OK. I will try and write about results later.

请先登录，再进行评论。

类别

在 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