I found the following piece of code to produce the same array as r in your example.

np.take(entry,indices,axis=0) 

It seems that stacking is not necessary with the take function.