使用Eigen setFromTriplets填充大稀疏矩阵的内存优化问题

你提到的内存计算存在明显错误：100k个非零元素的三元组，根本不需要400GB内存。Eigen的Triplet对象每个约占24字节（两个64位索引+一个64位数值），100k个Triplet总内存仅约2.4MB，远低于你给出的数值。内存崩溃的原因大概率是计算逻辑误解，或是生成Triplet时的冗余/错误操作，而非Triplet本身占用过大。

以下是降低内存占用的具体方法：

• 避免冗余Triplet生成
  若代码中先生成大量重复或无效的Triplet再去重，会额外消耗内存。应在生成Triplet阶段就确保每个非零元素仅生成一次，跳过重复的(i,j)索引对。

• 改用单精度数值类型
  若你的场景允许使用单精度（float）而非双精度（double），可将Triplet<double>替换为Triplet<float>，每个Triplet占用16字节（两个64位索引+32位数值），进一步压缩内存占用。

• 分批次插入元素
  无需一次性将所有Triplet传入setFromTriplets，可通过insert方法边生成边插入，跳过完整Triplet集合的存储：

  SparseMatrix<double> mat(300000, 300000);
  mat.reserve(100000); // 提前预留非零元素空间，避免动态扩容
  // 假设这里是你的Triplet生成逻辑，直接插入而非存储
  for (int i = 0; i < 100000; ++i) {
      int row = ...; // 合法行索引
      int col = ...; // 合法列索引
      double val = ...; // 元素值
      mat.insert(row, col) = val;
  }
  mat.makeCompressed(); // 压缩矩阵，释放冗余内存
  

• 提前预留矩阵内存
  无论使用setFromTriplets还是insert，都要先调用reserve(non_zero_count)，让矩阵提前分配足够空间，避免动态扩容带来的内存碎片和额外占用。

• 校验索引合法性
  若Triplet中存在超出[0, 299999]范围的行/列索引，Eigen处理时可能会分配远超预期的内存，导致崩溃。务必确保所有索引都在矩阵尺寸范围内。

内容的提问来源于stack exchange，提问作者Cong Pan