我对求解形式为AX=B的线性系统感兴趣，其中A和B是稀疏矩阵。
我通常使用scipy.sparse
X = scipy.sparse.linalg.spsolve(A, B)
但是，在我的应用程序中，我想针对不同的B多次解决该系统。
典型方法包括存储A的因式分解并使用关联的solve方法。例如，使用scipy.sparse.linalg.splu
A_inv = splu(A)
for B_i in Bs:
X_i = A_inv.solve(B_i)
但是似乎scipy.sparse.linalg.splu和scipy.sparse.linalg.factorize都不支持稀疏矩阵的RHS。
除了将B转换为密集矩阵外，python中是否还有其他稀疏的线性代数软件包可以使用此功能？