使用 UMFPACK 求解 Ax =B

solving Ax =B using UMFPACK

本文关键字:Ax 求解 UMFPACK 使用      更新时间:2023-10-16

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

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

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

此外,在"未来工作"下:

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

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