How can I use the Xcorr function for two matrix operation! I want to use this for gpuArray.

my code as follows

s1;     %this is a matrix for 512*403;
s2;     %this is a matrix for 512*403;
for i = 1:512
    [acor lag] = xcorr(s1(i,:),s2(i,:));
end

I want to use the gpuArray for this, and I don't want to use the for loop, how can I do this.

Answer 1

[Ridwan Alam]

The operation (xcorr) inside the for loop is not a matrix operation, rather an operation between two vectors (corresponding rows of two matrices). I don't think you can skip the for loop unfortunately.

Other than that,

S1 = gpuArray(s1);
S2 = gpuArray(s2);
acor = [];
for i = 1:512
    [acori,lag(i)] = xcorr(S1(i,:),S2(i,:));
    acor(i,:) = gather(acori);
end

[ang lee]

Do you think I can rewrite the xcorr function so that it can perform matrix operations, because I think parallel operations without the for loop will be faster.

[Ridwan Alam]

that's indeed a good idea. specially if you transpose the second matrix, it should have a syntax similar to matrix multiplication. good luck, and please keep us informed. also, if it works, please accept the response as an answer. thanks.