使用UMFPACK解决Ax = B.

Question

我正在使用UMFPACK来解决Ax = B，其中A和B都是矩阵。目前在我的实现中，我正在逐一解决x（：，i）= A \ b（：，i）来获取x。有没有办法直接获取x，而不必像LAPACK那样遍历B的每一列？

Answer 1

我不敢。 UMFPACK用户指南p59，关于umfpack _ * _求解例程：

  

此例程不能解决多个右侧，也不允许b以稀疏列形式存储。

此外，在“未来工作”下：

  

求解程序可以处理多个右侧和稀疏的右侧。

这意味着目前不支持它。