使用 UMFPACK 求解 Ax =B

solving Ax =B using UMFPACK

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

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

This routine does not solve for multiple right-hand-sides, nor does it allow b to be stored in a sparse-column form.

此外,在 "Future work" 下:

The solve routines could handle multiple right-hand sides, and sparse right-hand sides.

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