c++ - 如何计算 Eigen 库中稀疏矩阵的逆

Question

我对 C++ 中的 Eigen 库有疑问。实际上，我想计算稀疏矩阵的逆矩阵。当我在 Eigen 中使用密集矩阵时，我可以使用 .inverse() 操作来计算密集矩阵的逆。但是在稀疏矩阵中，我在任何地方都找不到逆运算。有谁知道计算稀疏矩阵的逆？帮我。

Answer 1

您不能直接执行此操作，但您始终可以使用一种稀疏求解器来计算它。思路是求解A*X=I，其中I是单位矩阵。如果有解，X 将是你的逆矩阵。 The eigen documentation有一个关于稀疏求解器及其使用方法的页面，但基本步骤如下:

SolverClassName<SparseMatrix<double> > solver;
solver.compute(A);
SparseMatrix<double> I(n,n);
I.setIdentity();
auto A_inv = solver.solve(I);