基于CUDA的二维泊松方程快速直接求解  

Fast 2-dimension Poisson Direct Solver Based on CUDA

作　　者：岳小宁[1,2] 肖炳甲[2] 罗正平[2] 

机构地区：[1]中国科学技术大学核科学技术学院,合肥230027 [2]中国科学院等离子体物理研究所计算机应用研究室,合肥230031

出　　处：《计算机科学》2013年第10期21-23,38,共4页Computer Science

基　　金：国家科技部"973"项目ITER计划专项(国内配套研究)(2012GB105000);国家自然科学基金(10835009)资助

摘　　要：二维泊松方程离散化之后可以转化为一个具有特殊格式的块三对角方程的求解问题,通过对这一结构化线性方程组的研究,提出了一个适用于统一计算架构(CUDA)的泊松方程并行算法。该算法通过离散正弦变化,可以将计算任务划分为若干相互独立的部分进行求解,各部分求解完成后再通过一次离散正弦变换即可获得最终解,整个求解过程只需要两次全局通信。结合GPU的硬件特征进行优化之后,该算法相比CPU上的串行算法可以获得10倍以上的加速比。The finite difference approximation of two-dimension poisson equation would create a block-tridiagonal equa-tion system. An algorithm which is suitable with the Compute Unified Device Architecture （CUDA）was proposed. Through a discrete sine transform, the computation task could be divided into several completely independent parts,af-ter solving these parts in parallel, the final result could be obtained through another discrete sine transform. Only two global synchronizing is needed during the computation. After carefully optimized, an acceleration rate of more than 10 times is obtained.

关 键 词：泊松方程 统一计算架构 并行计算 块三对角方程 

分 类 号：TP301.6[自动化与计算机技术—计算机系统结构]