This should just be writing a new version of the `smol_inds` function. Julia already has one. It won't be hard in python. - [x] Julia - [x] Python