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;
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.