Python 求解具有稀疏RHS的稀疏线性系统的scipy.sparse替代方案

我感兴趣的是求解一个线性系统，它的形式是

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中是否还有其他稀疏线性代数包支持此功能