假设我想要计算二次形式的x^T A y,其中x和y是NumPy向量,为了提高效率,我将A存储为SciPy csc_matrix。似乎执行计算的一种方法是因为SciPy对矩阵产品使用*,而NumPy使用.dot。虽然这种方法给了我正确的答案,但它似乎与我的直觉相反,我想知道是否有更好的可读性/更有效的方法?
我想知道是否有一种方法可以简单地将密集向量添加到稀疏矩阵的所有行中,在scipy.sparse中表示为一个稀疏矩阵并返回一个稀疏矩阵,即只对稀疏矩阵的非零元素求和。如果我这样做:import scipy.sparse as sp
W = sp.csr_matrix(np.array([[0,0,1],[0,