C++替换fftw_malloc数组为Eigen矩阵报static assertion failed错误

报错信息已经直接点明原因：THE_BRACKET_OPERATOR_IS_ONLY_FOR_VECTORS__USE_THE_PARENTHESIS_OPERATOR_INSTEAD，即Eigen对[]运算符做了编译期限制：方括号访问只允许用于一维向量（行/列向量），二维/高维矩阵访问元素必须使用圆括号()传入行列索引。
你之前用fftw_malloc申请的是裸连续内存，手动计算j + ny*i偏移访问是C语言操作二维数组的常规写法，但Eigen的Matrix类封装了存储逻辑，专门禁掉了非向量类型的[]运算符，就是为了避免开发者混淆线性索引、行列索引写出隐蔽bug，触发的是编译期静态断言，不会留到运行时才崩。
另外你写的代码里还有两个隐藏问题，原来的C版本就存在，Eigen版本里也没改：

• 内层循环j的上限写了ny+1，10x10的矩阵列索引合法范围是0~9，j跑到10的时候已经越界，属于未定义行为，之前C代码没报错只是运气好没触发内存保护。
• Eigen矩阵默认是列优先存储，和你手动计算偏移时用的行优先布局不一致，就算强行绕过[]的限制，元素存放的位置也和你预期的不一样。

最符合Eigen使用规范的写法是直接用(行号, 列号)的圆括号语法访问元素，不需要手动计算线性偏移，Eigen会自动根据存储顺序算对内存位置，也不用关心底层是行优先还是列优先：

static const int nx = 10;
static const int ny = 10; 

Eigen::Matrix<double, nx, ny> eXX;
eXX.setZero();

for(int i = 0; i < nx; i++){
    // 修正循环边界，去掉+1避免越界
    for(int j = 0; j < ny; j++){
        // 圆括号传入(i行,j列)坐标
        eXX(i, j) = i * 2 * EIGEN_PI / nx;
    }       
}

如果你需要兼容原来FFTW相关的逻辑，想继续用手动算线性偏移的裸指针操作，可以直接拿Eigen矩阵的底层数据指针，和你之前操作fftw_malloc返回的指针没有区别。如果要和原来行优先的内存布局完全对齐，记得把矩阵定义成行优先存储：

// 最后加RowMajor参数指定行优先存储，和C二维数组布局一致
Eigen::Matrix<double, nx, ny, Eigen::RowMajor> eXX;
double* XX = eXX.data(); // 直接拿底层连续内存指针，可以直接传给fftw接口
// 后面的逻辑和你原来的C代码完全一致即可，记得修正j的循环边界
memset(XX, 0, nx*ny*sizeof(double));
for(int i = 0; i< nx; i++){
    for(int j = 0; j< ny; j++){
        XX[j + ny*i] = i*2*EIGEN_PI/nx;
    }       
}

• Eigen默认分配的内存已经满足SIMD和FFTW要求的对齐规则，不需要额外用fftw_malloc申请内存再做拷贝，直接拿data()返回的指针传给FFTW接口即可。
• 不要尝试找绕过[]运算符限制的歪方法，这个限制是Eigen刻意加的防呆设计，按规范用(i,j)访问能减少很多索引写错的隐蔽bug。

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