parallel-processing - 使用 RowMajor 稀疏矩阵并行 BiCGSTAB 仍然按顺序运行

Question

我正在使用BiCGSTAB来求解稀疏方程组。我已将MPFR 用于多精度和以下代码，但代码是按顺序运行的。我遵循 Make Eigen run in parallel。

using MatrixXmp = Eigen::Matrix<mpfr::mpreal, Eigen::Dynamic, Eigen::Dynamic>;
using VectorXmp = Eigen::Matrix<mpfr::mpreal, Eigen::Dynamic, 1>;
...

omp_set_num_threads(num_threads);
Eigen::setNbThreads(num_threads);
mpfr::mpreal::set_default_prec(128);

Eigen::BiCGSTAB<Eigen::SparseMatrix<mpfr::mpreal, Eigen::RowMajor>> solver;
solver.setTolerance(tol);
solver.setMaxIterations(maxiter);
solver.compute(A);
xx = solver.solve(b);

编译：

g++ -o solver main.cpp -lmpfr -lgmp -fopenmp

我检查topcpu使用情况。我的操作系统是 Ubuntu 20.04。我看过这篇几乎相同的帖子，但我仍然遇到顺序运行问题。

编辑：A我像这样填充矩阵：

Eigen::SparseMatrix<mpfr::mpreal, Eigen::RowMajor> A(num_rows, num_rows);
A.reserve(Eigen::RowVectorXi::Constant(num_elements, 3));

for (int i = 0; i < num_elements; i++)
{
    A.insert(rows[i], cols[i]) = data[i];
}
A.makeCompressed();

标签： parallel-processingeigen