scipy中的有效稀疏矩阵计算

Question

我有一个稀疏矩阵A（d乘以d），我试图从中减去密集矢量（平均值）的乘法。现在我正在做的是逐行进行减法：

for i in range(d):
    A[i,:]=A[i,:]-mean[i]*(mean.T)

但是人们可能会预料到这一点非常缓慢。 A是csr格式，我尝试将其转换为lil以加速减法，但过程仍然很慢。我想知道如何更快地完成这项操作？一个建议可能是对一组行而不是单个行进行减法。但有更好的方法吗？