2

英特尔MKL库提供了优化的线程函数集,但对于迭代稀疏求解器(ISS),预条件共轭梯度方法似乎并不简单。更精确地说,使用预处理技术如不完全Cholesky分解或ILU,在某些点需要稀疏三角解算器,但是相应的MKL函数执行三角解法不是线程化的。我的问题是,是否有任何稀疏的求解器库在多核处理器上击败当前版本的MKL(它不是完全线程的)?CPU上最快的多线程迭代稀疏求解器?

+0

MKL是否改进了最新版本?你有没有尝试将它链接到MATLAB? – Royi

回答

1

我不知道它是否击败MKL--毕竟它是一个研究项目 - 但有一个混合模式OpemMP+MPI branch of PETSc可用。在我的小组对SpMV实现的实验中,它对于简单的矩阵很好,但对于我们经常处理的更复杂的矩阵并不适用。根据你的问题,你可能会有更好的运气。