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

2 次查看（过去 30 天）
Al in St. Louis
Al in St. Louis 2018-10-23
评论： Al in St. Louis 2018-10-23
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 个评论
显示 2更早的评论
Al in St. Louis
Al in St. Louis 2018-10-23
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

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Bruno Luong
Bruno Luong 2018-10-23
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 个评论
Al in St. Louis
Al in St. Louis 2018-10-23
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);

请先登录，再进行评论。

更多回答（1 个）

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

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

产品


版本

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by