在没有内存开销的情况下设置特征::稀疏矩阵的稀疏模式

我已经知道了 (我有唯一的排序列索引和行偏移量)。显然,通过setFromTriplets是可能的,但不幸的是,setFromTriplets需要很多额外的内存(至少在我的例子中是这样)

const long nRows = 5000000;
const long nCols = 100000;
const long nCols2Skip = 1000;
//It's quite big!
const long nTriplets2Reserve = nRows * (nCols / nCols2Skip) * 1.1;
Eigen::SparseMatrix<double, Eigen::RowMajor, long> mat(nRows, nCols);

std::vector<Eigen::Triplet<double, long>> triplets;

triplets.reserve(nTriplets2Reserve);
for(long row = 0; row < nRows; ++row){
    for(long col = 0; col < nCols; col += nCols2Skip){
        triplets.push_back(Eigen::Triplet<double, long>(row, col, 1));
    }
}
std::cout << "filling mat" << std::endl << std::flush;
mat.setFromTriplets(triplets.begin(), triplets.end());

std::cout << "Finished! nnz " << mat.nonZeros() << std::endl;
//Stupid way to check memory consumption
std::cin.get();

htop

所以我有两个问题:

1. 如何以尽可能少的开销填充特征::稀疏矩阵的稀疏模式

如果我的例子错了,请告诉我。

我的Eigen版本是3.3.2

编辑: 看起来像是 inserting (带 preallocation