C++ (特征库)中的分割故障

Question

我正在尝试使用特征库来确定超定线性系统的最小二乘解。

当矩阵很大时，如下面的程序所示，如6000行，代码将进入分段错误。当rows = 5000时，相同的程序也在运行。我读了很多关于修复分段错误的帖子，但我找不到任何有用的东西。

#include <cstdlib>
#include <iostream>
#include <Eigen/Dense>
#include <complex>
#define cols 45
#define rows 6000

using namespace Eigen;

int main() {
    int i,j;
    Matrix<std::complex<double>, rows,cols> zMat;

    for(i=0;i<rows;i++) {
        for(j=0;j<cols;j++) {
            zMat(i,j)=std::complex<double>(rand(),rand());
        }
    }

    VectorXcd ampIn(rows);
    for(i=0;i<rows;i++)
        ampIn(i)=std::complex<double>(rand(),rand());

    //Solving Matrix equation using QR Decomposition
    VectorXcd solution = zMat.colPivHouseholderQr().solve(ampIn);

    std::cout << "The least-squares solution is:\n" << solution << std::endl;

    return 0;
}