如何在Eigen中将块矩阵添加到稀疏矩阵

Question

例如，我有一个10x10的SparseMatrix A，我想在A的左上角添加一个3x3的单位矩阵。

在这3个条目中，已知A已经不为零。

如果必须一一添加值，也可以，但是我没有找到在Eigen中对稀疏矩阵的元素进行操作的方法。

我错过了什么？

Answer 1

如果只想一次将操作应用于特定元素，则可以像下面这样使用coeffRef ：

typedef Eigen::Triplet<double> T;
std::vector<T> coefficients;
for (int i = 0; i < 9; i++) coefficients.push_back(T(i, i, (i % 3) + 1));

Eigen::SparseMatrix<double> A(10, 10);
A.setFromTriplets(coefficients.begin(), coefficients.end());

std::cout << A << "\n\n";

for (int i = 0; i < 3; i++) A.coeffRef(i,i) += 1;
std::cout << A << "\n\n";