Numpy Array: Encuentra eficientemente los índices coincidentes

Tengo dos listas, una de las cuales es masiva (millones de elementos), la otra varios miles. Quiero hacer lo siguiente

<code>bigArray=[0,1,0,2,3,2,,.....]

smallArray=[0,1,2,3,4]

for i in len(smallArray):
  pts=np.where(bigArray==smallArray[i])
  #Do stuff with pts...
</code>

Lo anterior funciona, pero es lento. ¿Hay alguna forma de hacerlo de manera más eficiente sin tener que recurrir a escribir algo en C?

Respuestas a la pregunta(3)

Su respuesta a la pregunta