行/列切片scipy.sparse.csc_matrix中的内存效率低下

Question

使用以下代码按行/列对scipy.sparse.csc_matrix进行切片会导致切片后的内存分配加倍。 如何在切片时避免在内存中增加内存的同时切片？

import numpy
from scipy import sparse

@profile
def slice_matrix(n=10000):
    a = sparse.csc_matrix(numpy.random.random((n,n))>0.9, dtype=float)
    mask = numpy.repeat(numpy.random.choice([True, False]), n)
    return a[mask, ][:, mask]