使用 UMFPACK 求解 Ax =B

【问题标题】：solving Ax =B using UMFPACK使用 UMFPACK 求解 Ax =B
【发布时间】：2016-02-11 17:52:24
【问题描述】：

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

【问题讨论】：

【解决方案1】：

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

这个例程不能解决多个右手边，也不允许 b 以稀疏列的形式存储。

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

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

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

【讨论】：