无法使Eigen库矩阵乘法启用多线程，求问题排查方案

我在使用Eigen库开发C++程序时，遇到矩阵乘法无法并行化的问题——无论怎么配置，线程数始终无法超过1。具体操作及验证方式如下：

已做的配置 #

• 代码中添加并行初始化逻辑：

  Eigen::initParallel();
  omp_set_num_threads(4);
  Eigen::setNbThreads(4);
  

• 编译时添加OpenMP链接参数：-fopenmp
• 运行时设置环境变量：OMP_NUM_THREADS=4 ./program

验证方式 #

1. 调用Eigen::nbThreads()接口，返回值始终为1
2. 通过htop查看进程线程数，仅显示1个工作线程

复现代码 #

void m_multiply() {
    Eigen::Matrix<double,Eigen::Dynamic, Eigen::Dynamic> m1;
    Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> m2;
    m1.resize(1000,1000);
    m2.resize(1000,1000);
    m1.setRandom();
    m2.setRandom();

    while (true) {
        m1 *= m2;
        std::cout << Eigen::nbThreads() << std::endl; // 原代码缺少函数括号，已修正
    }
}

可能的问题及解决方案 #

1. Eigen未启用OpenMP后端 #

Eigen的并行化依赖OpenMP，但默认可能未开启。解决方法：

• 在所有Eigen头文件包含前添加宏定义：

  #define EIGEN_USE_OPENMP
  #include <Eigen/Dense>
  

• 若自行编译Eigen库，编译时需添加-DEIGEN_USE_OPENMP选项。

2. 矩阵大小未达并行阈值 #

Eigen对小矩阵乘法不会触发并行（线程调度开销大于并行收益）。尝试将矩阵尺寸放大至2000x2000或4000x4000后再测试。

3. 编译链接不完整 #

确保编译和链接阶段都添加了-fopenmp参数，完整编译命令示例：

g++ -fopenmp your_program.cpp -o program

4. Eigen版本过低 #

Eigen从3.3版本开始完善并行支持，若使用版本低于3.3，建议升级至最新稳定版。

5. 代码调用错误 #

原复现代码中Eigen::nbThreads缺少括号，正确调用应为Eigen::nbThreads()，否则输出的是函数指针而非当前线程数，可能导致误判。

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