我尝试使用来自 Eigen(3) 的 Intel MKL。我写了一个基本的矩阵乘法程序,如下所示。
#include <iostream>
#define EIGEN_USE_MKL_ALL
#include "../Eigen/src/Core/util/MKL_support.h"
#include "Dense"
#include "EigenTools.h"
#include <sys/time.h>
int main()
{
int rows= 10000, cols =100;
Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic> M(rows, cols);
Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic> N(cols, rows);
M.setRandom();
N.setRandom();
struct timeval tv_start;
struct timezone tz;
gettimeofday(&tv_start, &tz);
Eigen::MatrixXf P = M*N;
struct timeval tv_end;
gettimeofday(&tv_end, &tz);
float elapsedTime = (tv_end.tv_sec-tv_start.tv_sec)*1000000+tv_end.tv_usec-tv_start.tv_usec;
elapsedTime /= 1000000;
std::cout<<"rows, cols:"<<M.rows()<<","<<M.cols()<<std::endl;
std::cout<<"elpsed Time in sec is :"<<elapsedTime<<"\n";
return 0;
}
我使用“g++ EigenTest_IntelMKL.cpp -I../Eigen -DMKL_ILP64 -m64 -I /opt/intel/mkl/include”,我得到以下编译错误:
在 EigenTest_IntelMKL.cpp:3:0 中包含的文件中:../Eigen/src/Core/util/MKL_support.h:64:9: 错误:命名空间“std”中的“complex”未命名类型 typedef std::复杂的dcomplex;^ ../Eigen/src/Core/util/MKL_support.h:65:9:错误:命名空间“std”中的“complex”未命名类型typedef std::complex scomplex;