求解大型稀疏矩阵的并行方法

Elliothuo · 2019 年6 月 28 日 02:49

我现在想要解一个大型稀疏矩阵线性方程组，之前矩阵小的时候直接就用 a = b \ c 这种直接解法，现在规模大了以后发现这样效率比较低。
我想通过迭代法并行（我的电脑有20核，或许可以尝试多核并行）加速求解，有关这种计算有没有哪些可以直接使用的第三方库呢,或者其他的解决办法？

（GPU加速的话也可以，但是我在这方面不太懂，如果说他的操作不需要改动太多代码的话也可以考虑）

xgdgsc · 2019 年6 月 28 日 09:15

Elliothuo · 2019 年6 月 28 日 12:47

非常感谢!我看到在MKLSparse的README文件中指出：

Solving linear systems with triangular sparse matrices is supported. These matrices should be wrapped in their corresponding type, for example LowerTriangular for lower triangular matrices.For solving general sparse linear systems using MKL we refer to Pardiso.jl.

所以我认为可能Pardiso.jl更适合我的问题。