Fastest way to differentiate only the sides of large matricies
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I believe I have found one of the fastest ways to differentiate different portions of a large matrix using built in Matlab functions. However, there is a problem to my method, there is a remaining middle value between the side values. Allow me to demonstate with some simple code which makes use of gpuArray. For my case, I need the data on a gpu.
N = 300; % Number of points in x,y,z dimensions
S = 10; % Number of points from side of a dimension
for i = 1:20
B = gpuArray(rand(N,N,N)); % Load some random data
An = diff(B(1:S+1,:,:),1,1); % Differentation from left side
Ap = diff(B(end-S:end,:,:),1,1); % Differentation from right side
A1 = B([2:S+1,N-S+1:N],:,:)-B([1:S,N-S:N-1],:,:);
% The above is the finite difference which is applied to
% both sides of the matrix at the same time. This method
% however, is slower than using built in diff
A2 = diff(B([1:S+1,N-S:N],:,:),1,1);
% The above method is the fastest method, however there is one
% additional data point between S+1 and N - S. Deleting this
% point would add extra computational time
end
However diff is coded on the backend of Matlab, and would need to be modified for my application. If this is true, how would I go about learning how to code up my own diff function for this special case?
Thank you
0 commentaires
Réponse acceptée
Matt J
le 3 Oct 2021
Modifié(e) : Matt J
le 3 Oct 2021
This seems to be about as fast, at least on the GTX 1080 Ti.
T=reshape(B([1:S+1,N-S:N],:) ,S+1,2*N^2 );
A3 = reshape( diff(T,1,1), 2*S,N,N);
5 commentaires
Matt J
le 4 Oct 2021
Modifié(e) : Matt J
le 4 Oct 2021
No, you don't need to permute.
% along 1st axis
T = reshape(B([1:S+1,N-S:N],:) ,S+1,2,N,N );
A1 = reshape(diff(T,1,1), 2*S,N,N);
% along 2nd axis
T=reshape(B(:,[1:S+1,N-S:N],:),N,S+1,2,N);
A2=reshape( diff(T,1,2) , N,2*S,N);
% along 3rd axis
T=reshape(B(:,:,[1:S+1,N-S:N]),N,N,S+1,2);
A3=reshape( diff(T,1,3) , N,N,2*S);
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Get Started with Phased Array System Toolbox dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!