我正在解决一个更大的问题,当我尝试使用 OpenMP 并行化一些循环时遇到了一个错误。我用下面一些更简单的代码重现了这个问题,这些代码模仿了我自己的代码。
问题是当我运行程序时,它会随机进入某种无限循环/死锁(CPU 是 100%,但什么也不做)。从我的测试中可以看出,其中一个线程尝试计算矩阵-矩阵乘积,但由于某种原因从未完成。
我知道,如果您启用 OpenMP,Eigen 将使用 OpenMP 并行化矩阵-矩阵产品。我还在此之外添加了另一个并行循环。但是,如果我通过定义 EIGEN_DONT_PARALLELIZE 来禁用 Eigen 的并行化,仍然会出现此错误。
我在带有 Eigen 3.0.4 的 MacOS 10.6.8 上使用 gcc 版本 4.6.0 20101127。
我无法弄清楚可能出了什么问题...
#include <iostream>
#include <Eigen/Core>
using namespace std;
using namespace Eigen;
MatrixXd Test(MatrixXd const& F, MatrixXd const& G)
{
MatrixXd H(F.rows(), G.cols());
H.noalias() = F*G;
return H;
}
int main()
{
MatrixXd F = MatrixXd::Random(2,2);
MatrixXd G = MatrixXd::Random(2,2);
#pragma omp parallel for
for (unsigned int i = 0; i < 10000; ++i)
MatrixXd H = Test(F,G);
cout << "Done!" << endl;
}