Finding absolute difference between each row of a tall array

8 visualizzazioni (ultimi 30 giorni)
I am trying to find sum of absolute differences of each row. So if the input matrix A is M x N, the output matrix B should be M x M and
B( i , j ) = sum(abs( A ( i , : ) - A ( j , : ))). For small matrices the below two codes work.
lng=size(A,1);
B=reshape(sum(abs(repelem(A,lng,1)-repmat(A,lng,1)),2),lng,lng);
OR
lng=size(A,1); wdt=size(A,2);
B=zeros(lng);
for i=1:wdt
B=B+abs((A(:,i)-A(:,i)'));
end
The problem is that A is a tall matrix with dimension M larger than N. Thus, I cannot apply neither of these two methods on it. For the first code, it doesn't allow me to give the first inputs of repelem and repmat commands something different than 1. For the second code, taking the transpose fails for the tall array (is that forbidden totally for tall arrays or is there a workaround?) . Can someone please help me on this? I would appreciate if it would be a vectorized code.
Thanks in advance.
  1 Commento
madhan ravi
madhan ravi il 24 Dic 2018
Modificato: madhan ravi il 24 Dic 2018
probably Edric Ellis can share some light , he is really an expert on this (tall arrays)..

Accedi per commentare.

Risposte (1)

Cris LaPierre
Cris LaPierre il 24 Dic 2018
I would use the diff function.
First, check that the function(s) you want to use are compatible with tall arrays. See this page.
I would probably use abs(diff()).
There is a note on the limitations for diff with tall arrays:
  • You must use the three-input syntax Y = diff(X,N,dim).
  1 Commento
Osman
Osman il 24 Dic 2018
Unfortunately, diff does not work for me in this case. It only takes the difference between consecutive rows but I need the differences of a single row to all other rows of the matrix. Basically, I am trying to implement what pdist function does but it is not compatible with tall arrays. At the end, the code will be used for a knn Imputation method so if there is an easier way to apply it to tall arrays that will be really helpful. Still thanks for the advice :)

Accedi per commentare.

Prodotti


Release

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by