DSYEV和DSYEVD用于稀疏矩阵对角化

时间:2012-10-23 14:12:10

标签: fortran linear-algebra sparse-matrix matrix-multiplication lapack

所以我已经到了DSYEVD变得不切实际的地步,由于它的内存要求更高 - 它需要至少240GB的RAM来对齐我的矩阵,所以我考虑转移到DSYEV例程,需要较少的内存;

DSYEVD需要(1 + 6N + 3N ^ 2)字,而DSYEV只需要3N + N ^ 2.

较小矩阵的维度为~100,000,较大的维度为〜230,000。显然第二个必须使用DSYEV完成,但第一个我可以使用。

我快速查看了谷歌,但无法找到任何相关的基准 - 唯一可用的基准是N = 10-1000,并且是密集的矩阵。我的矩阵是100K维度和90%稀疏,230K维度,95%稀疏。

有谁知道我可以期待多长时间?

我已经使用DSYEVD从每个中获取了前3000个特征向量和值,但现在我需要全部使用它们。当仅采用3K向量时,使用迭代方案,100K矩阵需要~1.6天,230K矩阵需要~3天。这些工作是在Intel(R)Xeon(R)CPU E5-4640 0 @ 2.40GHz(256GB RAM)上运行的。

0 个答案:

没有答案