使用 UMFPACK 求解 Ax =B

Question

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

Answer 1

恐

怕不是。 UMFPACK 用户指南 p59，关于 umfpack_*_solve 例程：

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

此外，在"未来工作"下：

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

这反而意味着它目前不受支持。