如何并行化稀疏矩阵的反演

How to parallelize inversion of a sparse matrix?

本文关键字：并行化更新时间：2023-10-16

我在 c++ 中使用 eigen3 包来反转一些大型稀疏矩阵（例如 12000*12000），我需要这些矩阵进行后续操作; 但是，这真的很耗时，我无法将其扩展到更大的矩阵。是否可以并行执行此操作，例如使用 openmp？

提前谢谢。

您可以使用特征矩阵

作为稀疏矩阵容器，并使用 MKL 并行稀疏求解器函数来计算反演。

https://software.intel.com/node/521676

尽管 Eigen 不提供 MKL 备用求解器的 API。您仍然可以在低级界面中使用它们。