c++ - solving Ax =B using UMFPACK

Question

I'm using UMFPACK to solve Ax = B where both A and B are matrices. Currently in my implementation, I'm solving x(:,i) = A\b(:,i) one by one to get x. Is there a way to directly get x without having to loop through each columns of B like LAPACK does?

Answer 1

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

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

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

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

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