eigen values of 3D matrix
9 次查看(过去 30 天)
显示 更早的评论
Hello,
I need to compute eigen values of 3D matrix, something like 3x3x30000, means 30000 different 3x3 matrices. Do you have any suggestions how to do it without the "for" loop?
I've tried convert the large matrix to the cell format and do it through cellfun(@eig...) and it takes for 17 seconds but problem is that I need to do it for 150 times. So, it's still 45 minutes. Any improvement will be appreciated.
Thanks
回答(2 个)
James Tursa
2012-4-17
You can try Bruno Luong's small size multi-eigenvalue solver on the FEX:
0 个评论
Matt Tearle
2012-4-17
I have suggestions how to do it without the for loop, but not in a way that will speed things up! :) Making a block diagonal or such-like just slows things down with the memory it chews up.
So: sorry, it looks like looping is your best option. If you have Parallel Computing Toolbox, you can at least distribute this task, which might give you some improvement.
I wondered about writing your own eig function for the special case of 3-by-3, but it's hard to beat the built-in eig.
Here's my experiment:
n = 30000;
x = randi(10,[3,3,n]);
tic
ev1 = zeros(3*n,1);
for k = 1:n
j = 3*k;
ev1(j-2:j) = eig(x(:,:,k));
end
toc
tic
ev2 = zeros(3*n,1);
for k = 1:n
j = 3*k;
ev2(j-2:j) = eig3by3(x(:,:,k)); % my own eig for 3x3 matrices
end
toc
tic
y = num2cell(x,[1,2]);
z = cellfun(@eig,y,'unif',false);
ev3 = cat(1,z{:});
toc
Output:
Elapsed time is 0.315167 seconds.
Elapsed time is 0.814756 seconds.
Elapsed time is 0.281028 seconds.
Out of curiosity, what version of MATLAB are you running? (Note the significant difference in run times I get from what you were getting.)
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!