编译期已知尺寸时，如何动态分配Eigen矩阵的std::vector？

我正在使用Eigen线性代数库编写计算密集型程序，需要在循环中对大量约300×300的矩阵进行对角化。我希望在循环前分配一个可重复使用的STL容器存储若干Eigen矩阵，现有可运行代码如下：

std::vector<Eigen::Matrix<complex<double>, Eigen::Dynamic, Eigen::Dynamic>> A; 
for (int i = 0; i < Imax; i++){
// some calculation using A 
}

但我在编译期已知该vector的大小（含4个元素）以及每个矩阵的尺寸（300×300），能否在循环前动态分配该vector，使其拥有固定的非零大小？

当然可以，这么做能避免循环中频繁的内存分配，显著提升计算密集型程序的性能，针对你的需求有几种高效实现方式：

方法1：预分配vector大小并初始化矩阵 #

直接在vector构造时指定元素数量，再逐个将矩阵resize到目标尺寸：

// 构造时指定vector包含4个默认构造的空矩阵
std::vector<Eigen::Matrix<std::complex<double>, Eigen::Dynamic, Eigen::Dynamic>> A(4);
// 遍历每个矩阵，设置为300×300尺寸
for (auto& mat : A) {
    mat.resize(300, 300);
    // 可选：若需要初始化为零矩阵，添加此行
    // mat.setZero();
}

方法2：直接构造指定尺寸的矩阵 #

跳过默认构造的空矩阵，直接在vector中构造已指定尺寸的矩阵，同时预分配vector容量避免扩容：

std::vector<Eigen::Matrix<std::complex<double>, Eigen::Dynamic, Eigen::Dynamic>> A;
// 预分配足够容量，避免vector自动扩容
A.reserve(4);
// 构造4个300×300的矩阵并加入vector
for (int i = 0; i < 4; ++i) {
    A.emplace_back(300, 300);
    // 可选：初始化矩阵值
    // A.back().setZero();
}

方法3：使用Eigen固定尺寸矩阵（性能最优） #

既然编译期已知矩阵是300×300，推荐使用Eigen的固定尺寸矩阵类型，彻底消除动态尺寸的额外开销：

// 固定尺寸矩阵，编译期确定大小，内存布局更紧凑
std::vector<Eigen::Matrix<std::complex<double>, 300, 300>> A(4);
// 每个矩阵默认已初始化，无需resize，可直接在循环中使用

这种方式下Eigen会优化矩阵的内存管理，对角化这类密集计算的性能会比动态尺寸矩阵更好。

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