Fast Monte-Carlo Photon Transport Employing GPU-Based Parallel Computation

M. Mirzapour, K. Hadad, R. Faghihi, R. J. Hamilton, Christopher J Watchman

Research output: Contribution to journalArticlepeer-review

1 Scopus citations


Monte Carlo (MC) is known to be the most accurate dose calculation method. However, MC suffers from high computational cost as a large number of particles have to be simulated to achieve the desired statistical uncertainty. Enhancing computational power by parallelizing the simulation with multiple GPU threads reduces the time required to reach the desired uncertainty in MC simulation. In this article, we present DOSXYZgpu, a GPU implementation of EGSnrc code which is written in CUDA Fortran as an algorithm. This article relies on a well validated and popular code among medical physicists, EGSnrc/DOSXYZnrc. In order to transport particles between two consecutive interactions, we developed an algorithm to handle several thousands of histories per warp. DOSXYZgpu implementation is evaluated with the original sequential EGSnrc/DOSXYZnrc. Maximum speedup of 205 times is achieved while the statistical uncertainty of the simulation is preserved. The t-test statistical analysis indicates that for more than 95% of the voxels there is no significant difference between the results obtained from the GPU and the CPU.

Original languageEnglish (US)
Article number8985184
Pages (from-to)450-460
Number of pages11
JournalIEEE Transactions on Radiation and Plasma Medical Sciences
Issue number4
StatePublished - Jul 2020


  • Dosimetry
  • GPU-based computation
  • Monte Carlo (MC) methods
  • parallel programming

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics
  • Radiology Nuclear Medicine and imaging
  • Instrumentation


Dive into the research topics of 'Fast Monte-Carlo Photon Transport Employing GPU-Based Parallel Computation'. Together they form a unique fingerprint.

Cite this