将Eigen :: BICGSTAB与不完整ILU的OpenMP配合使用不会并行运行

编程问答 2022-08-13

问题描述

我的BICGSTAB无法与不完整的LU预条件器并行运行： 我在DEBIAN下使用带有编译选项-fopenmp的g ++ 请帮我。 我给Eigen解决的电话是：

initParallel();
omp_set_num_threads(8);
setNbThreads(8);

// Fill Eigen Row Major sparse matrix            
SparseMatrix<double,RowMajor> A(n,n);
A.setFromTriplets(T.begin(),T.end());

//map std::vector B to Eigen VectorXd
Map<Eigen::VectorXd> b(B_.data(),B_.size());


// BICGSTAB
BiCGSTAB<SparseMatrix<double,RowMajor>,Eigen::IncompleteLUT<double>  >  solver(A);
solver.setTolerance(1e-15);

VectorXd X(n);
X = solver.solve(b);

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

eigenopenmpsolver