Visual Studio下实现Eigen+Intel MKL的矩阵乘法

最新推荐文章于 2024-02-28 17:34:32 发布

zhujiahui622

最新推荐文章于 2024-02-28 17:34:32 发布

阅读量1w

点赞数 4

分类专栏：项目配置文章标签： Eigen Intel-MKL 矩阵乘法 visual studio

本文链接：https://blog.youkuaiyun.com/zhujiahui622/article/details/50504678

版权

项目配置专栏收录该内容

1 篇文章

订阅专栏

科学计算过程中经常会用到线性代数运算，尤其是矩阵乘法，自己手写的C++代码想要赶上MATLAB几乎是不可能的。在C++中，谈到高性能计算就往往要用Intel MKL，但是用过的人都知道，这东西烦人的很，函数风格也怪异，不好用，一个简单的矩阵乘法函数dgemm的参数列表就老长老长。Eigen是一个相对来说十分优雅和好用的线性代数运算库（尤其是在Windows下，Blas，Lapack太麻烦，基本不考虑），缺点就是密集型计算“不那么快”。

所以，如果既能有Intel MKL的速度，又能有Eigen的优雅，还是在Windows下，那真真是极好的。

操作系统：Windows 7 64位

安装Visual Studio 2013
下载C++下的一个线性代数运算库Eigen并解压：http://eigen.tuxfamily.org/index.php?title=Main_Page
安装Intel Parallel Studio XE 2015

在Visual Studio下配置Intel MKL：详见 http://blog.youkuaiyun.com/zhujiahui622/article/details/50504164
在Visual Studio下配置Eigen。假设Eigen解压在D:\Local\eigen下，则在项目属性中的“配置属性”下的“C/C++”面板中将“附加包含目录”设置为D:\Local\eigen
写代码

#define EIGEN_USE_MKL_ALL
#define EIGEN_VECTORIZE_SSE4_2

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

using namespace std;
using namespace Eigen;

// 使用Eigen+Intel MKL
int main(int argc, char *argv[])
{
    MatrixXd a = MatrixXd::Random(3000, 3000);  // 随机初始化矩阵
    MatrixXd b = MatrixXd::Random(3000, 3000);

    double start = clock();
    MatrixXd c = a * b;    // 乘法好简洁
    double endd = clock();
    double thisTime = (double)(endd - start) / CLOCKS_PER_SEC;

    cout << thisTime << endl;

    system("PAUSE");
    return 0;
}