应使用Eigen::SparseMatrix类，它专为稀疏矩阵设计，仅存储非零元素，底层采用默认列优先的压缩列存储（CSC）或可选的压缩行存储（CSR）格式，避免内存爆炸和运算低效。
稀疏矩阵该用哪个Eigen类？
别用 MatrixXd 或 MatrixXf 存稀疏矩阵——内存爆炸、运算极慢。必须用 SparseMatrix<t></t>，它只存非零元，底层用压缩列（CSC）或压缩行（CSR）格式，默认是列优先的 CSC。
常见误用：把稠密矩阵转成 SparseMatrix 后直接赋值循环填元素，性能反而更差。正确做法是预分配非零项数量，再用三元组（triplet）批量构建：
std::vector<Eigen::Triplet<double>> triplets;
triplets.reserve(nnz); // nnz 是预估非零元个数
for (auto& entry : data) {
triplets.emplace_back(entry.row, entry.col, entry.val);
}
Eigen::SparseMatrix<double> A(n, n);
A.setFromTriplets(triplets.begin(), triplets.end());注意：setFromTriplets() 会自动去重并排序，但不检查索引越界——越界会导致运行时断言失败（Debug 模式）或未定义行为（Release）。
如何高效地修改稀疏矩阵的单个元素？
直接写 A(i,j) = x 在稀疏矩阵上是 O(nnz) 操作，因为要遍历列/行查找位置。频繁单点赋值会拖慢整个流程，尤其在迭代求解器里反复更新时。
立即学习“C++免费学习笔记（深入）”；
真正高效的方案只有两个：
- 如果只是初始化阶段少量修改，先用
Triplet构建，再一次性setFromTriplets() - 如果必须动态更新（比如 ILU 预处理中填充），改用
SparseMatrix::coeffRef(i, j)—— 它会在需要时自动插入新元素，平均摊还成本较低，但首次调用仍可能触发内存重分配
别用 A.coeff(i,j) 去读——它只读，不支持写；也别在循环里反复调用 coeffRef 而不预分配空间，否则每次插入都可能引起内部缓冲区扩容。
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
稀疏矩阵乘法为什么结果变稠密了？
A * B 运算本身没问题，但结果类型默认仍是 SparseMatrix。问题出在：如果 A 和 B 的结构导致大量隐式零参与计算（例如带宽很大、pattern 重叠严重），乘积矩阵的非零元数量可能远超预期，甚至接近 n²，这时内存和时间都会失控。
关键对策：
- 用
A.prune(1e-12)清理数值上接近零的残留项（注意：这不会改变 sparsity pattern，只删掉绝对值小于阈值的现有非零元） - 乘法前用
A.isCompressed()检查是否已压缩——未压缩的矩阵（如刚用coeffRef插入过）参与运算会显著变慢 - 对称正定问题优先用
A.selfadjointView<Lower>() * B，能利用结构减少约一半计算量
另外，SparseMatrix 不支持原位乘法（如 A *= B），必须写成 A = A * B，这意味着一次额外的内存分配——大矩阵时得留意临时对象开销。
求解稀疏线性方程组选哪个模块？
Eigen 自带的 SparseLU 和 SparseQR 可直接用，但别碰 BiCGSTAB 或 ConjugateGradient 之前不设预处理器——收敛慢到无法接受是常态。
实操要点：
-
SparseLU<double>适合中小型一般矩阵（n < 50k），内存占用高但稳定；启用COLAMD重排序（默认开启）能大幅减少填充元 - 大型对称正定问题，用
ConjugateGradient<Eigen::SparseMatrix<double>, Eigen::DiagonalPreconditioner<double>>—— 对角预处理器开销小，比不加快 5–10 倍 - 非对称且病态的问题，
BiCGSTAB+ILUT预处理器效果好，但ILUT不在 Eigen 主库，需自行实现或引入第三方封装（如unsupported/Eigen/src/IterativeSolvers/IncompleteLUT.h）
最后提醒：所有稀疏求解器都要求输入矩阵已压缩（A.makeCompressed()），否则可能静默降速或触发断言；而且它们都不检查矩阵是否奇异，解出来 x 的残差得自己用 (A*x - b).norm() 验证。
