Related Experiment Video
Updated: Jun 3, 2026

Novel 3D/VR Interactive Environment for MD Simulations, Visualization and Analysis
Published on: December 18, 2014
GPUMCD: A new GPU-oriented Monte Carlo dose calculation platform
Sami Hissoiny1, Benoît Ozell, Hugo Bouchard
1Département de Génie Informatique et Génie Logiciel, Ecole Polytechnique de Montréal, 2500 Chemin de Polytechnique, Montréal, Québec H3T 1J4, Canada. sami.hissoiny@polymtl.ca
Purpose:
Monte Carlo methods are considered as the gold standard for dosimetric computations in radiotherapy. Their execution time is, however, still an obstacle to the routine use of Monte Carlo packages in a clinical setting. To address this problem, a completely new, and designed from the ground up for the GPU, Monte Carlo dose calculation package for voxelized geometries is proposed: GPUMCD.
Method:
GPUMCD implements a coupled photon-electron Monte Carlo simulation for energies in the range of 0.01-20 MeV. An analog simulation of photon interactions is used and a class II condensed history method has been implemented for the simulation of electrons. A new GPU random number generator, some divergence reduction methods, as well as other optimization strategies are also described. GPUMCD was run on a NVIDIA GTX480, while single threaded implementations of EGSnrc and DPM were run on an Intel Core i7 860.
Results:
Dosimetric results obtained with GPUMCD were compared to EGSnrc. In all but one test case, 98% or more of all significant voxels passed the gamma criteria of 2%-2 mm. In terms of execution speed and efficiency, GPUMCD is more than 900 times faster than EGSnrc and more than 200 times faster than DPM, a Monte Carlo package aiming fast executions. Absolute execution times of less than 0.3 s are found for the simulation of 1M electrons and 4M photons in water for monoenergetic beams of 15 MeV, including GPU-CPU memory transfers.
Conclusion:
GPUMCD, a new GPU-oriented Monte Carlo dose calculation platform, has been compared to EGSnrc and DPM in terms of dosimetric results and execution speed. Its accuracy and speed make it an interesting solution for full Monte Carlo dose calculation in radiation oncology.
Related Concept Videos
Maxwell-Boltzmann Distribution: Problem Solving
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
Numerical Calculations
The solution to a problem is obtained using different methods. While manually solving algebraic symbols is one of the most common methods, the graphical method is often preferred. Computers...
Parallel Processing
Computed Tomography
The technique was invented in the 1970s and is based on the principle that as X-rays pass through the body, they are absorbed or reflected at different levels. In the technique, a patient lies on a motorized platform while a computerized axial tomography (CAT) scanner rotates...
Multimachine Stability
In analyzing the system, the nodal equations represent the relationship between bus voltages, machine voltages, and machine currents. The nodal equation is given by:
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...