稀疏LU分解不能直接用double模拟二维数组，因非零元分布不规则，连续内存会严重浪费空间或致溢出；必须分离结构（int索引）与数值（double*），CSR格式需维护val、col_ind、row_ptr三数组，且须预留填充元空间并同步重分配。
稀疏LU分解为什么不能直接用 double* 指针模拟二维数组
因为稀疏矩阵的非零元分布不规则，用连续内存模拟满阵会浪费大量空间甚至导致内存溢出。直接用 double* 配合行主序计算索引（如 A[i * n + j]）在稀疏场景下既无法跳过零元，也无法支持高效的前向/后向代入。实际中必须分离结构与数值：结构用索引数组（int*）描述非零位置，数值用 double* 存储对应值。
CSR 格式下 LU 分解需维护哪几个指针数组
CSR（Compressed Sparse Row）是 LU 实现最常用格式，必须管理三组指针：val（非零值）、col_ind（列索引）、row_ptr（行起始偏移）。LU 过程中这些数组会被原地修改 —— val 逐步填入 L 和 U 的非零元，col_ind 和 row_ptr 可能因填充元（fill-in）而重分配。常见错误是只给 val 分配原始非零元数量的空间，没预留填充余量，导致越界写入或提前失败。
-
row_ptr[0] = 0，row_ptr[i+1] - row_ptr[i]表示第i行非零元个数（分解后会增长） -
col_ind必须与val等长，且每个col_ind[k]对应val[k]的列号 - 若预估填充量不足，
realloc时需同步更新三者，否则col_ind和val错位会导致数值错位
如何用指针安全实现 Doolittle 形式的稀疏 LU 前向扫描
Doolittle 要求 L 对角元为 1，U 保留原对角元。关键在于每轮消去时，只更新当前行（U 部分）和后续行中与当前列相关的非零元（L 部分）。不能遍历整列，必须用 col_ind 和 row_ptr 定位哪些行在该列有非零元。
例如，处理第 k 列时，先从 row_ptr[k] 开始扫描 col_ind，找出所有满足 col_ind[j] == k 的行 i > k（即 L 的待定元素），再对每个这样的 i，遍历其在 U 中已计算的行段（row_ptr[i] 到 row_ptr[i+1]-1），跳过列号 的项，对列号 <code>>= k 的项执行消去。这里极易因索引计算错误访问到未初始化的 val 位置。
立即学习“C++免费学习笔记（深入）”；
典型错误代码片段：
// 错误：假设每行都从 row_ptr[i] 开始存 U 元素，但 L 和 U 混存在 val 中
for (int j = row_ptr[i]; j < row_ptr[i+1]; ++j) {
if (col_ind[j] >= k) u_val = val[j] - l_val * u_pivot; // u_pivot 可能未定义
}
为什么 std::vector 替代裸指针反而更易出错
不是不能用，而是当接口要求传入裸指针（如绑定 C 风格求解器、对接 suitesparse 或 Intel MKL）时，用 &vec[0] 获取指针必须确保 vector 未发生 realloc —— 一旦 push_back 或 resize，原有指针失效。更隐蔽的问题是：vector 的 data() 返回 double*，但它的生命周期依赖 vector 对象；若函数内局部 vector 返回 data()，调用方拿到的是悬垂指针。
真实项目中建议：用类封装 CSR 三元组，内部用 std::unique_ptr<double[]> + std::unique_ptr<int[]> 管理内存，暴露 get_val() 等方法返回稳定指针，且禁止拷贝。裸指针操作本身不难，难的是让指针生命周期、所有权、重分配时机全部对齐数值算法逻辑。
