How to get a sub matrix of a large sparse matrix / array efficiently?

I stumbled upon the same problem and it is not super easy to get what you want because of the compressed format. I think the simplest is to work at the level of the vectors I,J,K given by findnz.