我会假设接口在过去 5 年内发生了变化和发展。截至目前，可以使用scipy（如评论中要求的）的标准功能来实现结果。

• 使用scipy.io.mmread将 MatrixMarket 文件中的矩阵加载到稀疏或密集的 scipy 矩阵中。我会假设你的矩阵最初是稀疏的，因此你很高兴。
• 现在，可以使用 scipy.sparse.linalg.slu （如果您选择不完全分解，也可以使用 iLU）。

相反，您可以使用

• scipy.io.mmwrite将您的矩阵写入 MatrixMarket 文件。