三元组稀疏矩阵加法不能直接for循环硬扫，因为会退化为o(nnz_a × nnz_b)时间复杂度；正确做法是双指针归并两个按(row, col)字典序升序排列的三元组数组，实现o(nnz_a + nnz_b)高效相加，并需过滤浮点零值以保证稀疏性。
三元组稀疏矩阵加法为什么不能直接 for 循环硬扫？
因为两个三元组数组 triplets_A 和 triplets_B 都是按行主序（或列主序）排序的，但非零元位置完全不重合。硬扫会退化成 O(nnz_A × nnz_B)，而实际只需单趟合并——类似归并排序的双指针遍历，时间复杂度压到 O(nnz_A + nnz_B)。
常见错误是忽略坐标比较顺序：必须先比 row，再比 col，否则相同位置元素可能被漏合并；还有人用 std::map<:pair>, double></:pair> 累加，看似简洁，但插入和查找都是 O(log N)，整体变慢且内存开销大。
- 确保输入三元组已按
(row, col)字典序升序排列（可用std::sort预处理，但生产环境建议构造时就维护） - 用两个索引
i、j分别遍历 A 和 B，比较triplets_A[i].row与triplets_B[j].row，相等再比col - 结果写入预分配的 vector（如
std::vector<triplet></triplet>），避免频繁push_back触发 realloc
如何避免重复坐标导致的格式化输出错乱？
加法后可能出现多个三元组指向同一 (row, col)（比如 A 和 B 在同一位置都有非零元），必须合并为单个项；但若未清除零值（如 1.0 + (-1.0) = 0.0），后续转 CSR 或打印时就会混入冗余零元，破坏稀疏性假设。
典型错误是在合并阶段只累加值，却不检查是否为零；或者用 fabs(val) 判零，但 <code>eps 设得过大（如 1e-6）会误删小量，过小（如 1e-16）又无法过滤浮点误差累积。
- 合并逻辑中，只要
row和col相同，就累加value；之后立即判断std::abs(sum) > 1e-12（对 double 类型足够安全） - 格式化输出（如打印为 “row col value” 表格）前，务必过滤掉所有零值项——这不是可选项，是稀疏矩阵语义要求
- 若需导出为 Matrix Market 格式，header 中的非零元数量
nnz必须是过滤后的实际个数，否则解析器会读错
CSR 转换时最容易被忽略的内存布局陷阱
从三元组转 CSR（Compressed Sparse Row）不是简单排序+计数。关键在 row_ptr 数组长度是 n_rows + 1，且 row_ptr[i] 表示第 i 行第一个非零元在 col_idx/values 中的起始下标——这个定义容易和 0-based/1-based 混淆，尤其当矩阵行数为 0 或某行全零时。
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
常见翻车点：循环写成 for (int i = 0; i 填 <code>row_ptr，却忘了最后一项 row_ptr[n_rows] 必须等于总非零元数；或者用 std::lower_bound 查每行起始位置，但没处理空行边界。
- 先用
std::vector<int> row_count(n_rows, 0)</int>统计每行非零元个数，再做前缀和生成row_ptr -
row_ptr[0] = 0，然后row_ptr[i] = row_ptr[i-1] + row_count[i-1]（i 从 1 到 n_rows），最后row_ptr[n_rows]自动等于总 nnz - 填充
col_idx和values时，按三元组当前顺序（已按 row 主序）逐个写入，不要重新排序
性能瓶颈往往卡在 std::vector 的移动和内存分配上
高频调用场景（如迭代求解器内多次矩阵加）下，反复构造/析构 std::vector<triplet></triplet> 会触发大量堆分配。即使用了 reserve()，小对象（如 3×double 的 Triplet）仍受 malloc/glibc 分配器碎片影响。
有人改用 std::array 或栈内存，但尺寸固定不灵活；也有人套用内存池，却因线程安全或生命周期管理出错。最稳妥的优化是复用缓冲区 + move 语义。
- 把三元组容器封装进类（如
SparseTriplet），内部持有一个std::vector<triplet></triplet>成员，并提供clear()和reserve()接口 - 加法函数接收左值引用并 move 出结果：
auto result = add_triplets(std::move(A), std::move(B))，避免拷贝 - 若确定最大 nnz 上限，可配合
std::unique_ptr<triplet></triplet>+ 定长分配器（如std::pmr::monotonic_buffer_resource），但需 C++17 及以上
真正难的不是算法逻辑，而是让三元组在加、转、输之间不产生隐式拷贝，也不留零值垃圾——这两点漏掉任意一个，高效就只剩字面意思。
C++免费学习笔记（深入）：立即使用
在学习笔记中，你将探索 C++ 的入门与实战技巧！
