当矩阵大小超过特征矩阵类型的某个极限时，c++分配不好

在c++动态库中，我使用特征库解决了一个最小二乘问题。这个Dll是在python软件中调用的，在该软件中可以解决问题配置。在小型问题中，代码工作正常，并返回正确的解决方案。如果点数增加，则库抛出std::bad_alloc。

更准确地说，产生简化到最大错误的代码是

try {
    matrixA = new Eigen::MatrixXd(sizeX,NvalidBtuple); // initialize A
    for (int i=0;i<sizeX;++i) {
        int secondIndex = 0;
        for (int k=0;k<btermSize;++k) {
            if (bterm[k] == 1) {                        // select btuple that are validated by density exclusion
                // product of terms
                (*matrixA)(i,secondIndex) = 1.0;
                secondIndex += 1;
            }
        }
    }
} catch (std::bad_alloc& e) {
    errorString = "Error 3: bad allocation in computation of coefficients!";
    std::cout<<errorString<<" "<<e.what()<<std::endl;
    return;
} catch (...) {
    errorString = "Error 4: construction of matrix A failed! Unknown error.";
    std::cout<<errorString<<std::endl;
    return;
}

其中CCD_ 2被定义在具有CCD_ 3的头文件中。

如果sizeX和NvalidBtuple小于约20'000x3'000，则矩阵定义有效。如果大小更大，它就会崩溃。

我做测试的电脑有足够的可用内存，大约有15G的可用内存。

这是堆/堆栈问题吗？如何使库接受更大的矩阵？

欢迎发表评论。谢谢

编辑：正如下面的回答所说，我不清楚NvalidBtuple的定义：

NvalidBtuple = 0;
for (int i=0;i<btermSize;++i) {NvalidBtuple += bterm[i];}

其中CCD_ 7是布尔向量。因此，由于在循环中我们检查if (bterm[k] == 1)，所以secondIndex总是小于NvalidBtuple。