用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存储在稀疏的列形式中。

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

解题例程可以处理多个右侧，以及稀疏的右侧.

这意味着目前不支持它。