I need to solve a few eigenvalues of a large matrix specified by their indices. These indices are according to the whole eigenspectrum sorted in algebraic (not absolute value) ascending order.

I notice this is made available by the subset_by_index option in scipy.linalg.eigh, but not in, e.g., the sparse counterpart eigsh. Is this possible at all in the sparse context?