稀疏对称矩阵的LDL<sup>T</sup>分解在GPU上的高效实现

稀疏对称矩阵的LDL^T分解在GPU上的高效实现

陈鑫峰,王武

An Effective Implementation of LDLT Decomposition of Sparse Symmetric Matrix on GPU

Chen Xinfeng,Wang Wu

伪代码5. 清空当前列多线程分解所用内存

1: function cleartmpMem(Lp, Li, tmpMem,
tmpMem1, n,k,col)
2: id=blockIdx.x*blockDim.x+threadIdx.x;
3: if id < Lp[k+1]-Lp[k]-1 then
4: tmpMem1[n*col+Li[Lp[k]+1+id]] = 0;
5: tmpMem[n*col+Li[Lp[k]+1+id]] = 0;
6: end if
7: end function