C++(特征库)中的分段错误

Question

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

当矩阵大小较大时，代码将进入分段错误，如下面的程序中所示，如6000行。当行数=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;
}

Answer 1

对于特征库，对于大于1MB 的数据，使用MatrixXcd而不是Matrix<std::complex<double>...