GPU-accelerated Monte Carlo convolution/superposition implementation for dose calculation.

PURPOSE Dose calculation is a key component in radiation treatment planning systems. Its performance and accuracy are crucial to the quality of treatment plans as emerging advanced radiation therapy technologies are exerting ever tighter constraints on dose calculation. A common practice is to choose either a deterministic method such as the convolution/superposition (CS) method for speed or a Monte Carlo (MC) method for accuracy. The goal of this work is to boost the performance of a hybrid Monte Carlo convolution/superposition (MCCS) method by devising a graphics processing unit (GPU) implementation so as to make the method practical for day-to-day usage. METHODS Although the MCCS algorithm combines the merits of MC fluence generation and CS fluence transport, it is still not fast enough to be used as a day-to-day planning tool. To alleviate the speed issue of MC algorithms, the authors adopted MCCS as their target method and implemented a GPU-based version. In order to fully utilize the GPU computing power, the MCCS algorithm is modified to match the GPU hardware architecture. The performance of the authors' GPU-based implementation on an Nvidia GTX260 card is compared to a multithreaded software implementation on a quad-core system. RESULTS A speedup in the range of 6.7-11.4x is observed for the clinical cases used. The less than 2% statistical fluctuation also indicates that the accuracy of the authors' GPU-based implementation is in good agreement with the results from the quad-core CPU implementation. CONCLUSIONS This work shows that GPU is a feasible and cost-efficient solution compared to other alternatives such as using cluster machines or field-programmable gate arrays for satisfying the increasing demands on computation speed and accuracy of dose calculation. But there are also inherent limitations of using GPU for accelerating MC-type applications, which are also analyzed in detail in this article.

[1]  I. Rosen,et al.  Beam‐commissioning methodology for a three‐dimensional convolution/superposition photon dose algorithm , 2000, Journal of applied clinical medical physics.

[2]  J J Battista,et al.  Generation of photon energy deposition kernels using the EGS Monte Carlo code. , 1988, Physics in medicine and biology.

[3]  Shuang Luan,et al.  Stochastic versus deterministic kernel-based superposition approaches for dose calculation of intensity-modulated arcs. , 2008, Physics in medicine and biology.

[4]  Cedric X. Yu,et al.  Arc-modulated radiation therapy (AMRT): a single-arc form of intensity-modulated arc therapy , 2008, Physics in medicine and biology.

[5]  He Wang,et al.  Use of deformed intensity distributions for on-line modification of image-guided IMRT to account for interfractional anatomic changes. , 2005, International journal of radiation oncology, biology, physics.

[6]  Ergun Ahunbay,et al.  An On-line Adaptive Planning Strategy for Inter-fraction Image Guidance , 2007 .

[7]  A L Boyer.  Relationship between attenuation coefficients and dose-spread kernels. , 1988, Radiation research.

[8]  C Cameron,et al.  Sweeping-window arc therapy: an implementation of rotational IMRT with automatic beam-weight calculation , 2005, Physics in medicine and biology.

[9]  T. Mackie,et al.  A dual source photon beam model used in convolution/superposition dose calculations for clinical megavoltage x-ray beams. , 1997, Medical physics.

[10]  Wen-mei W. Hwu,et al.  Optimization principles and application performance evaluation of a multithreaded GPU using CUDA , 2008, PPoPP.

[11]  Benoît Ozell,et al.  Fast convolution-superposition dose calculation on graphics hardware. , 2009, Medical physics.

[12]  A Bel,et al.  Accelerated ray tracing for radiotherapy dose calculations on a GPU. , 2009, Medical physics.

[13]  J. Battista,et al.  A convolution method of calculating dose for 15-MV x rays. , 1985, Medical physics.

[14]  M A Earl,et al.  An arc-sequencing algorithm for intensity modulated arc therapy. , 2007, Medical physics.

[15]  George L. Yuan,et al.  GPU compute memory systems , 2009 .

[16]  John Amanatides,et al.  A Fast Voxel Traversal Algorithm for Ray Tracing , 1987, Eurographics.

[17]  Naga K. Govindaraju,et al.  A Survey of General‐Purpose Computation on Graphics Hardware , 2007 .

[18]  Yonghui Chen.  Using FPGA to Accelerate Monte Carlo Superposition Based Radiation Dose Calculation , 2009 .

[19]  Steve B. Jiang,et al.  Development of a GPU-based Monte Carlo dose calculation code for coupled electron–photon transport , 2009, Physics in medicine and biology.

[20]  A Fenster,et al.  X-ray sources of medical linear accelerators: focal and extra-focal radiation. , 1993, Medical physics.

[21]  Russell H. Taylor,et al.  Towards real-time radiation therapy: GPU accelerated superposition/convolution , 2010, Comput. Methods Programs Biomed..

[22]  N Papanikolaou,et al.  Investigation of the convolution method for polyenergetic spectra. , 1993, Medical physics.

[23]  Justin P. Haldar,et al.  Accelerating advanced mri reconstructions on gpus , 2008, CF '08.

[24]  R. Mohan,et al.  Differential pencil beam dose computation model for photons. , 1986, Medical physics.

[25]  A. Ahnesjö.  Collapsed cone convolution of radiant energy for photon dose calculation in heterogeneous media. , 1989, Medical physics.

[26]  Raj Shekhar,et al.  Direct aperture deformation: an interfraction image guidance strategy. , 2006, Medical physics.