I have a function, call it get_block(I, J), which takes two index sets, I and J, and returns a scipy.sparse.csr_matrix wi
get_block(I, J)
I
J
scipy.sparse.csr_matrix