OpenCL中的矩阵求逆

Question

我试图使用OpenCL加速一些计算,部分算法包括反转矩阵.是否有任何开源库或免费提供的代码来计算用OpenCL或CUDA编写的矩阵或一般反演的lu分解(lapack dgetrf和dgetri)？矩阵是实数和正方形,但除此之外没有任何其他特殊属性.到目前为止,我已经设法在gpu上找到了基本的blas矩阵向量运算实现.

矩阵相当小,只有大约60-100行和cols,所以它可以在cpu上更快地计算,但它在算法的中间使用,所以我必须将它转移到主机,计算逆,和然后将结果传回设备,然后在更大的计算中使用它.

Answer 1

看看ViennaCL:http://viennacl.sourceforge.net/

Answer 2

我没有在Open CL中实现,但"数字食谱"和Gil Strang的"进入应用数学"都有很好的解释,很容易编码."NR"具有您可以适应的C代码.

计算逆

这是不正确的.你没有计算LU分解的逆,你正在分解矩阵.如果你想要逆,你必须用一系列单位向量进行前向替换.这是一个小但重要的区别.