c++ - Eigen 中的 vector 加和点不通过mkl加速

这是我使用的代码:

#define EIGEN_USE_MKL_ALL

#include <iostream>
#include <Eigen/Core>
#include <Eigen/Dense>
#include <time.h>

using namespace std;
using namespace Eigen;

int main(int argc, char *argv[])
{
    VectorXf a = VectorXf::Random(100000000);
    VectorXf b = VectorXf::Random(100000000);

    double start = clock();
    VectorXf c = a+b;
    float d = a.dot(b);
    double endd = clock();
    double thisTime = (double)(endd - start) / CLOCKS_PER_SEC;

    cout << thisTime << endl;

    return 0;
}

用mkl编译:

g++ mkl_test.cpp /home/tong.guo/intel/mkl/lib/intel64/libmkl_rt.so -Ieigen -Wl,--no-as-needed -lpthread -lm -ldl -m64 -I/home/tong.guo/intel/mkl/include

删除第一行代码，然后不使用mkl进行编译:

g++ mkl_test.cpp -Ieigen

时间差不多。

但是可以加快矩阵计算的速度。
将代码更改为以下代码，我可以看到速度。

    MatrixXd a = MatrixXd::Random(1000, 1000);
    MatrixXd b = MatrixXd::Random(1000, 1000);

    double start = clock();
    MatrixXd c = a * b;
    double endd = clock();
    double thisTime = (double)(endd - start) / CLOCKS_PER_SEC;

    cout << thisTime << endl;

最佳答案

从eigen page启用mkl:

由于 vector 加法和点积是1级blas例程，因此Eigen在这里将不使用外部例程。

关于c++ - Eigen 中的 vector 加和点不通过mkl加速，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/52869694/