CSR与CSC的核心区别在于indptr的语义:CSR中indptr[i]表示第i行首个非零元在values中的下标,适合按行访问;CSC中indptr[j]表示第j列首个非零元下标,适合按列访问。

CSR 和 CSC 的核心区别在哪?别只看存储顺序
CSR(Compressed Sparse Row)和 CSC(Compressed Sparse Column)本质不是“行优先 vs 列优先”的简单类比,而是**访问模式决定结构**。CSR 适合按行遍历(比如矩阵-向量乘 y = A * x),CSC 适合按列遍历(比如 A.transpose() * x 或求解器中 LU 的列主元选主)。误用会导致缓存失效、循环嵌套反直觉、甚至逻辑错误——比如用 CSR 写一个按列求和的函数,性能可能差 5 倍以上。
二者共享三数组结构:values(非零值)、indices(对应行号或列号)、indptr(行/列起始偏移)。关键差异在 indptr 含义:
– CSR 中 indptr[i] 是第 i 行第一个非零元在 values 中的下标;
– CSC 中 indptr[j] 是第 j 列第一个非零元在 values 中的下标。
手写 CSR 类时,indptr 容易越界或漏填 1 个元素
indptr 长度必须是 rows + 1(CSR)或 cols + 1(CSC),最后一个元素恒为非零元总数 nnz。新手常犯两个错:
– 把 indptr 初始化成长度 rows,导致访问 indptr[rows] 时越界;
– 构造时只填了前 rows 个位置,忘了设 indptr[rows] = nnz,后续遍历某行会读到错误范围。
正确做法:先统计每行非零个数,做前缀和:
std::vector<int> indptr(rows + 1, 0);
for (int i = 0; i < rows; ++i) {
indptr[i + 1] = indptr[i] + row_nnz[i]; // row_nnz[i] 是第 i 行非零元数
}注意:所有索引默认从 0 开始,别混用 1-based 的文献伪代码。
插入新元素?别直接 push_back 到 values 和 indices
C++ 稀疏矩阵一旦用 CSR/CSC 格式构建完成,就**不支持高效随机插入**。因为插入会破坏 indptr 的连续性,还可能触发整行数据搬移。常见错误场景:
– 读取 COO 格式(三元组)时边读边插,结果时间复杂度退化成 O(nnz²);
– 在已构造好的 CSR 上调用 insert(i, j, val),内部反复 resize 向量。
正确路径只有两条:
– 批量构建:先收集全部三元组,排序(按行→列对 CSR),再一次性构造三数组;
– 用中间格式过渡:COO → 排序 → CSR。标准库没提供现成排序,得自己写:std::sort(coo_triplets.begin(), coo_triplets.end(), [](const auto& a, const auto& b) { return a.row < b.row || (a.row == b.row && a.col < b.col); });
用 Eigen 或 Intel MKL?先确认你真需要它们
Eigen 的 SparseMatrix<double, RowMajor> 默认就是 CSR,ColMajor 就是 CSC。但它封装太深:比如想直接访问 indptr,得调 innerIndexPtr()、outerIndexPtr()、valuePtr(),且这些指针只读;MKL 的 sparse_matrix_t 更是要求显式 convert,调试时看不到原始数组。
组合式C++代码评审方案,融合静态分析、AI推理、多轮迭代评审和C++专项检查,适用于PR审查、增量代码审查、全项目评审和代码质量评分,触发词包括review cpp、cpp代码评审、C++review、代码审查。
立即学习“C++免费学习笔记(深入)”;
如果你只是做教学、小规模实验、或需深度定制(比如带符号压缩、混合精度索引),手写轻量 CSR 类更可控。示例骨架:
struct CSRMatrix {
std::vector<double> values;
std::vector<int> indices; // 行索引(CSR)或列索引(CSC)
std::vector<int> indptr; // 长度 rows+1(CSR)或 cols+1(CSC)
int rows, cols, nnz;
};真正难的不是存数据,而是写对 matvec 和 transpose —— 这俩函数里 indptr 和 indices 的配合稍有偏差,结果就全错。建议先用小矩阵(3×4)手算一遍每步下标,再写代码验证。


















