让我们考虑稀疏矩阵。

import numpy as np
import scipy.sparse as sps

A = sps.csc_matrix([[0, 1], [0, 0]])
L = sps.csc_matrix([[1, 1],[0, 0]]).astype(bool)

具有 A 作为任意数值矩阵和 L 稀疏逻辑矩阵，用于通过选择进行连续约简。
产量 print(A[L]) ```
matrix(0, 1)

然后在数字上看起来还可以，但是输出类型 `numpy.matrix` 这是非常不幸的，因为它既是密集矩阵类型，也是不推荐使用的类型。
有没有办法通过另一个稀疏矩阵来索引稀疏矩阵，从而获得稀疏矩阵？我在考虑最小化计算开销。