Abstract:Purpose: Using the graphical processing units (GPU) hardware technology, an extremely fast Monte Carlo (MC) code ARCHER(RT) is developed for radiation dose calculations in radiation therapy. This paper describes the detailed software development and testing for three clinical TomoTherapy (R) cases: the prostate, lung, and head & neck.Methods: To obtain clinically relevant dose distributions, phase space files (PSFs) created from optimized radiation therapy treatment plan fluence maps were used as the input to ARCHERRT. Patient-specific phantoms were constructed from patient CT images. Batch simulations were employed to facilitate the time-consuming task of loading large PSFs, and to improve the estimation of statistical uncertainty. Furthermore, two different Woodcock tracking algorithms were implemented and their relative performance was compared. The dose curves of an Elekta accelerator PSF incident on a homogeneous water phantom were benchmarked against DOSXYZnrc. For each of the treatment cases, dose volume histograms and isodose maps were produced from ARCHER(RT) and the general-purpose code, GEANT4. The gamma index analysis was performed to evaluate the similarity of voxel doses obtained from these two codes. The hardware accelerators used in this study are one NVIDIA K20 GPU, one NVIDIA K40 GPU, and six NVIDIA M2090 GPUs. In addition, to make a fairer comparison of the CPU and GPU performance, a multithreaded CPU code was developed using OpenMP and tested on an Intel E5-2620 CPU.Results: For the water phantom, the depth dose curve and dose profiles from ARCHER(RT) agree well with DOSXYZnrc. For clinical cases, results from ARCHER(RT) are compared with those from GEANT4 and good agreement is observed. Gamma index test is performed for voxels whose dose is greater than 10% of maximum dose. For 2%/2mm criteria, the passing rates for the prostate, lung case, and head & neck cases are 99.7%, 98.5%, and 97.2%, respectively. Due to specific architecture of GPU, modified Woodcock tracking algorithm performed inferior to the original one. ARCHER(RT) achieves a fast speed for PSF-based dose calculations. With a single M2090 card, the simulations cost about 60, 50, 80 s for three cases, respectively, with the 1% statistical error in the PTV. Using the latest K40 card, the simulations are 1.7-1.8 times faster. More impressively, six M2090 cards could finish the simulations in 8.9-13.4 s. For comparison, the same simulations on Intel E5-2620 (12 hyperthreading) cost about 500-800 s.Conclusions: ARCHER(RT) was developed successfully to perform fast and accurate MC dose calculation for radiotherapy using PSFs and patient CT phantoms. c 2014 American Association of Physicists in Medicine.

Heterogeneous Concurrent Execution of Monte Carlo Photon Transport on CPU, GPU and MIC

Status of archer — a monte carlo code for the high-performance heterogeneous platforms involving gpu and mic

Archer, A New Monte Carlo Software Tool for Emerging Heterogeneous Computing Environments

Comparison of Two Accelerators for Monte Carlo Radiation Transport Calculations, NVIDIA Tesla M2090 GPU and Intel Xeon Phi 5110p Coprocessor: A Case Study for X-ray CT Imaging Dose Calculation

Scalable and massively parallel Monte Carlo photon transport simulations for heterogeneous computing platforms

GPU-Accelerated Monte Carlo Electron Transport Methods: Development and Application for Radiation Dose Calculations Using Six GPU cards

Graphics Processing Unit-Accelerated Mesh-Based Monte Carlo Photon Transport Simulations

TU-AB-BRC-10: Modeling of Radiotherapy Linac Source Terms Using ARCHER Monte Carlo Code: Performance Comparison of GPU and MIC Computing Accelerators.

Modeling of Radiotherapy Linac Source Terms Using ARCHER Monte Carlo Code: Performance Comparison for GPU and MIC Parallel Computing Devices

Efficient Execution of Microscopy Image Analysis on CPU, GPU, and MIC Equipped Cluster Systems

Testing of ARCHER-CT, A Fast Monte Carlo Code for CT Dose Calculation: Experiment Versus Simulation

Performance Study of Monte Carlo Codes on Xeon Phi Coprocessors — Testing MCNP 6.1 and Profiling ARCHER Geometry Module on the FS7ONNi Problem

SU‐F‐BRD‐02: Application of ARCHERRT‐‐ A GPU‐Based Monte Carlo Dose Engine for Radiation Therapy ‐‐ to Tomotherapy and Patient‐Independent IMRT

Archer(Rt) - A Gpu-Based And Photon-Electron Coupled Monte Carlo Dose Computing Engine For Radiation Therapy: Software Development And Application To Helical Tomotherapy

MIC acceleration of short-range molecular dynamics simulations

GPU-based Monte Carlo Simulation for Light Propagation in Complex Heterogeneous Tissues.

EVALUATION OF VECTORIZED MONTE CARLO ALGORITHMS ON GPUS FOR A NEUTRON EIGENVALUE PROBLEM

Optimizing the Monte Carlo Neutron Cross-Section Construction Code Xsbench for Mic and Gpu Platforms

Parallel 3D Deterministic Particle Transport on Intel MIC Architecture

Accelerated 3 D Full Band Self-consistent Ensemble Monte Carlo Device Simulation Utilizing Intel MIC Coprocessors on TianHe II

Performance Portable Monte Carlo Particle Transport on Intel, NVIDIA, and AMD GPUs