An OpenCL-based Monte Carlo dose calculation engine (oclMC) for coupled photon-electron transport

TL;DR

This paper introduces oclMC, a fast, cross-platform Monte Carlo dose calculation engine for radiotherapy, implemented in OpenCL, achieving comparable accuracy to CUDA-based methods across various hardware platforms.

Contribution

Developed a portable, efficient Monte Carlo dose engine in OpenCL that maintains accuracy and works across multiple hardware platforms, unlike previous CUDA-only solutions.

Findings

oCLMC achieves dose calculation accuracy within 0.5% of CUDA-based gDPM.

oCLMC runs on diverse hardware including Nvidia, AMD, and Intel devices.

The engine's efficiency is comparable across different platforms, with minor performance differences.

Abstract

Monte Carlo (MC) method has been recognized the most accurate dose calculation method for radiotherapy. However, its extremely long computation time impedes clinical applications. Recently, a lot of efforts have been made to realize fast MC dose calculation on GPUs. Nonetheless, most of the GPU-based MC dose engines were developed in NVidia CUDA environment. This limits the code portability to other platforms, hindering the introduction of GPU-based MC simulations to clinical practice. The objective of this paper is to develop a fast cross-platform MC dose engine oclMC using OpenCL environment for external beam photon and electron radiotherapy in MeV energy range. Coupled photon-electron MC simulation was implemented with analogue simulations for photon transports and a Class II condensed history scheme for electron transports. To test the accuracy and efficiency of our dose engine oclMC, we compared dose calculation results of oclMC and gDPM, our previously developed GPU-based MC code, for a 15 MeV electron beam and a 6 MV photon beam on a homogenous water phantom, one slab phantom and one half-slab phantom. Satisfactory agreement was observed in all the cases. The average dose differences within 10% isodose line of the maximum dose were 0.48-0.53% for the electron beam cases and 0.15-0.17% for the photon beam cases. In terms of efficiency, our dose engine oclMC was 6-17% slower than gDPM when running both codes on the same NVidia TITAN card due to both different physics particle transport models and different computational environments between CUDA and OpenCL. The cross-platform portability was also validated by successfully running our new dose engine on a set of different compute devices including an Nvidia GPU card, two AMD GPU cards and an Intel CPU card using one or four cores. Computational efficiency among these platforms was compared.