如何用CUDA同时计算几百个实对称矩阵的eigenvalues/eigenvecot (转载) - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Joke - 肚皮舞运动

如何用CUDA同时计算几百个实对称矩阵的eigenvalues/eigenvecot (转载)

如何用CUDA同时计算几百个实对称矩阵的eigenvalues/eigenvecot (转载)# Joke - 肚皮舞运动

g*t2012-07-02 07:07

1 楼

【以下文字转载自 Computation 讨论区】
发信人: geneticdrift (不懂微积分), 信区: Computation
标题: 如何用CUDA同时计算几百个实对称矩阵的eigenvalues/eigenvecot
发信站: BBS 未名空间站 (Mon Jul 2 02:38:51 2012, 美东)
我有一个程序要反复计算几百个（约500个）64 x 64的实对称矩阵的所有的
eigenvalues/eigenvectors。自己用CUDA实现了一个Jacobi algorithm with chess
tournament ordering。具体来说，每个block(含有32个threads)处理一个矩阵，这32
个threads并行消去一个矩阵中的32个off-diagonal elements，直到算法收敛。结果无
误，计算单个矩阵所花的时间也和最近的一篇paper里的数据接近。但是这个算法和CPU
上的library比没有太大的优势。在同时处理这500个矩阵的情况下，和GSL里面高度优
化的函数比较(用单CPU)，用GPU仅仅快了一倍。我觉得主要是Jacobi algorithm对于这
个大小的矩阵效率太差，而GSL里面的函数用的好像是QR decomposition，虽然只有一
个thread但是效率很高。有没有比较适合我的问题的能在GPU上高效执行的算法？有没
有什么paper/code可以参考的？先谢谢了。

d*f2012-07-02 07:07

2 楼

http://saahpc.ncsa.illinois.edu/10/papers/paper_19.pdf
based on this paper the speed-up should be significant

32
CPU

【在 g**********t 的大作中提到】

: 【以下文字转载自 Computation 讨论区】
: 发信人: geneticdrift (不懂微积分), 信区: Computation
: 标题: 如何用CUDA同时计算几百个实对称矩阵的eigenvalues/eigenvecot
: 发信站: BBS 未名空间站 (Mon Jul 2 02:38:51 2012, 美东)
: 我有一个程序要反复计算几百个（约500个）64 x 64的实对称矩阵的所有的
: eigenvalues/eigenvectors。自己用CUDA实现了一个Jacobi algorithm with chess
: tournament ordering。具体来说，每个block(含有32个threads)处理一个矩阵，这32
: 个threads并行消去一个矩阵中的32个off-diagonal elements，直到算法收敛。结果无
: 误，计算单个矩阵所花的时间也和最近的一篇paper里的数据接近。但是这个算法和CPU
: 上的library比没有太大的优势。在同时处理这500个矩阵的情况下，和GSL里面高度优

g*t2012-07-02 07:07

3 楼

谢谢博导。我现在的程序就是基于这篇文章写的。这篇文章里的提升是对单个大矩阵来
说的，而且在这篇CPU端用的是同样的算法(Jacobi method)。其实CPU端可以用更高效
的算法，这样GPU就没有太大优势了。我想问一下有没有适合我的矩阵大小的，可以在
GPU上高效实现的算法？比如QR algorithm？

【在 d********f 的大作中提到】

: http://saahpc.ncsa.illinois.edu/10/papers/paper_19.pdf
: based on this paper the speed-up should be significant
:
: 32
: CPU