稀疏矩陣乘積的特征求解器

Question

我想使用 Eigen 來計算 L_1^{-1}L_2，其中 L_1 和 L_2 都是下三角矩陣，並在 Eigen 中存儲為面向列的稀疏矩陣。 我嘗試了本征三角求解器。 但是，這需要 L_2 是密集的。

Answer 1

實際上，對於稀疏rhs， solve方法並沒有重載，但是您可以像這樣使用solveInPlace方法（我實際上沒有嘗試過）：

Eigen::SparseMatrix<double> foo(Eigen::SparseMatrix<double> const& L1, Eigen::SparseMatrix<double> const& L2)
{
    Eigen::SparseMatrix<double> res = L2;
    L1.triangularView<Eigen::Lower>().solveInPlace(res);
    return res;
}

您仍然應該考慮是否確實需要完整的矩陣。