Speed up a program

12 views (last 30 days)
Vlad
Vlad on 4 Dec 2014
Commented: Vlad on 5 Dec 2014
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.

Answers (2)

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

Sign in to comment.


Edric Ellis
Edric Ellis on 5 Dec 2014
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 Comment
Vlad
Vlad on 5 Dec 2014
OK. I will try and write about results later.

Sign in to comment.

Products

Community Treasure Hunt

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

Start Hunting!