Hi Tommy,
Well I tried using numpy.where. I want my biggest cluster to be removed. I get the id's to be removed using cl.cluster_keys[0]. Now I search for the bonds to be removed.
rembnd = np.empty(0,dtype=int)
for x in cl.cluster_keys[0]:
rembnd= np.append(rembnd,np.where(bonds==x)[0] ) # This will give me row number
new_bonds = np.delete(bonds,rembnd,axis=0)
But this method is extremely slow when I have 457632 particles in that cluster.
Is there a better way ? (I think the for loops slowing things down, if I am not mistaken )
Regards,
Jude