在R中，Python相当于Matrix :: chol和Matrix :: solve？

Question

在R中我正在做以下事情：

L = ... # some sparse matrix L
chol_factor = Matrix::chol(L)

b = # some vector
z = Matrix::solve(chol_factor, b)

solve将通过Cholesky因子有效地计算事物。我一直在使用scipy.sparse代码，但似乎并不是Cholesky分解实现，也不是有效解决因素的方法。在python中有相同的方法吗？

Answer 1

这可以通过scipy的稀疏LU分解来完成。

import numpy as np
from scipy.sparse import linalg as sla

L = # some sparse matrix 
lu = sla.splu(L)

b = # some vector
z = lu.solve(b)