How to remove items from two arrays that index each other without for loop

2 visualizaciones (últimos 30 días)
This is hard to explain in words. Here's my code, and I'd like to know whether there's a way to vectorize the part with a for loop:
t = [obj.network.node(obj.posIx).type];
mPosIx = obj.posIx (t == 13 | t == 9);
mIxPos = zeros (1, max (mPosIx));
for idx = 1:length(mPosIx)
mIxPos (mPosIx (idx)) = idx;
end
This code does exactly what I need it to do. It removes all items that aren't type 9 or 13. The mPosIx array points to the correct position in the node array. The mIxPos array either contains a zero where there is no corresponding element in the position array (which does not appear in this code), or it contains the correct index in the position array. The for loop just feels clunky, but I couldn't figure out another way to do it.
  4 comentarios
madhan ravi
madhan ravi el 23 de Oct. de 2018
what do you mean by this?
mIxPos (mPosIx (idx)) = idx;
Al in St. Louis
Al in St. Louis el 23 de Oct. de 2018
It's a line of code that assigns the index of the for loop, idx, to the correct position. Here's the output of running the code on my test data. I get two vectors with the correct values in the correct positions. I just want to get rid of the for loop.
>> mIxPos
mIxPos =
1 0 0 0 2
>> mPosIx
mPosIx =
1 5

Iniciar sesión para comentar.

Respuesta aceptada

Bruno Luong
Bruno Luong el 23 de Oct. de 2018
You can vectorize LHS as with RHS, the for loop becomes one line
mIxPos(mPosIx) = 1:length(mPosIx);
That returns the inverse of permutation mPosIx
  1 comentario
Al in St. Louis
Al in St. Louis el 23 de Oct. de 2018
That's much better than calling find and saving two of its three outputs. My function has been reduced to:
t = [obj.network.node(obj.posIx).type];
mPosIx = obj.posIx (t == 13 | t == 9);
mIxPos (mPosIx) = 1:numel (mPosIx);

Iniciar sesión para comentar.

Más respuestas (1)

Al in St. Louis
Al in St. Louis el 23 de Oct. de 2018
This works:
mIxPos = zeros (1, max (mPosIx));
[~,col,v]=find(mPosIx)
mIxPos (v) = col;

Categorías

Más información sobre Logical en Help Center y File Exchange.

Productos


Versión

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by