GPU加速剂量计算中微分卷积积分算法的实现.PDFVIP

第36 卷 第5 期 核 技 术 Vol. 36, No.5 2013 年5 月 NUCLEAR TECHNIQUES May 2013 GPU加速剂量计算中微分卷积/积分算法的实现 1,2 1 1 1 1 王先良 刘乐乐 吴章文 勾成俊 侯 氢 1 （四川大学原子核科学技术研究所辐射物理及技术教育部重点实验室 成都 610064 ） 2 （四川省肿瘤医院放疗科 成都 610041 ） 摘要 微分卷积/积分算法是计算精度较高的光子线剂量计算算法，较长的计算时间限制了该算法在临床上的 使用。本文对微分卷积/积分算法中最耗时的部分实现了基于GPU 的并行化计算，与基于CPU 的计算相比， 在Tesla C1060 上计算速度提高可达30X−60X 。利用γ 因子对计算结果的准确性进行了分析，结果显示，无论 是均匀水模还是非均匀头模，在单照射野还是多照射野情况下，加速后的结果都与CPU 的计算结果有相同的 准确性。通过GPU 并行加速，微分卷积/积分算法能成为日常的剂量计算算法。 关键词 CUDA ，微分卷积/积分算法，GPU，剂量计算 中图分类号 TL72 DOI: 10.11889/j.0253-3219.2013.hjs.36.050207 剂量计算是放疗计划系统的核心内容之一，快 的计算有较高的精度，是计算精度仅次于蒙特卡洛 速准确地提供感兴趣区内所受照射剂量的数据，对 方法的光子线剂量计算算法[4,8] 。另外该算法还能很 放疗计划的制定至关重要，如何在保证剂量计算精 好地进行非规则野和调强野的剂量计算，但较长计 度的前提下，减少剂量计算时间是放疗领域研究的 算时间限制了该算法在临床上的使用，为了使微分 [1] 热点之一 。 卷积/积分算法能成为日常的剂量计算算法，本文采 提高剂量计算速度的方法主要有：采用不同的 用计算能力更强的硬件(GPU)对该算法进行加速。 剂量计算算法和借助于计算能力更强的硬件。目前， GPU 最初是用于计算机图像显示和渲染的元 放疗中研究和使用的剂量计算算法[2]基本分为蒙特 件。现在的GPU 不仅是功能强大的图像处理引擎， 卡洛算法、卷积/积分算法和有限笔形束算法。蒙特 还可用于科学计算的高度并行化的可编程处理器。 卡洛算法通常被作为剂量计算的标准，模拟了粒子 与传统的并行运算方式(如集群计算)相比，GPU 并 与物质相互作用的全过程，能计算各种复杂条件下 行计算有很多优点：首先GPU 有很多计算核心，硬 的剂量分布，但计算复杂，耗时长，临床实际应用 件上适于做并行计算；其次，GPU 有较高的存储器 较少，理论研究较多[3] 。卷积/积分法又可分为三维 读写带宽和浮点数运算能力，具备很高的计算性能； [4] [5] 卷积/积分 和二维卷积/积分 ，采用空间不变的积 再次，基于GPU 的并行运算所需投入少，能在个人 分核计算的二维卷积/积分可以通过快速傅里叶变 电脑上实现，在低成本情况下就可获得不错的计算 换加速剂量计算，所以也称为快速傅里叶法。有限 效率。近几年来，GPU 在可编程和硬件架构等方面 笔束算法的核心思想是将照射野离散成许多相互独 的快速发展在各领域引发了大量关于 GPU 并行计 立的有限大小的射束元，介质中沉积在某