Eigen优化及并行性能测试

最新推荐文章于 2025-07-25 16:34:31 发布

原创

最新推荐文章于 2025-07-25 16:34:31 发布 · 4k 阅读

5 ·

CC 4.0 BY-SA版权

近期开展方程组求解过程中，发现求解速度随着节点数量的增加快速降低，还没有确定出具体原因在哪，于是就先对使用的Eigen库进行了研究，Eigen是一个开源的矩阵计算库，使用C++语言，大量应用模板，同时有BLAS、MKL、LAPACKE等高性能库的接口。

测试源程序如下：
#include
#include “Eigen/Dense”
#include <time.h>
#include <stdio.h>

using namespace Eigen;
using namespace std;
 
int main(int argc,char *argv[])
{
  //Eigen::initParallel();
  int dim;
  if(argc>1)
  {
    dim = atoi(argv[1]);
  }
  else{
    dim = 10;
  }
  MatrixXd A = MatrixXd::Random(dim,dim);
  MatrixXd B = MatrixXd::Random(dim,dim);
  //cout << "A =" << endl << A << endl;
  VectorXd b = VectorXd::Random(dim);
  //cout << "b =" << endl << b << endl;
  VectorXd x(dim);
  time_t begin_t  = clock();
  x = A.colPivHouseholderQr().solve(b);  
  //cout << "The solution is:\n" << x << endl;
  //A= A*B;
  time_t finish_t = clock();
  cout<<"it cost " << (double )(finish_t - begin_t )/CLOCKS_PER_SEC <<" s"<<