TU‐E‐BRB‐09: Development of a GPU‐Based Monte Carlo Dose Calculation Package for Cancer Radiotherapy

X. Jia, X. gu, Y. Graves, S. Jiang

Research output: Contribution to journalArticle

Abstract

Purpose: Monte Carlo (MC) simulation is commonly considered as the most accurate dose calculation method in radiotherapy. However, its efficiency still requires improvement for many routine clinical applications. The goal of this work is to develop a GPU‐based MC dose calculation package, gDPM. It utilizes the parallel computation ability of a GPU to achieve high efficiency, while maintaining same particle transport physics as in the original DPM code and hence simulation accuracy.Methods: The gDPM code is implemented on GPU architecture under CUDA platform. In GPU computing, divergence of execution paths between threads can considerably reduce the efficiency. Since photons and electrons undergo different physics and hence attain different execution path, we design a simulation scheme where photon transport and electron transport are separated to partially relieve the thread divergence issue. High performance random number generator and hardware linear interpolation are also utilized. We have also developed various components to account for fluence map and linac geometry, so that gDPM can be used to compute dose for real IMRT or VMAT treatment plans. The gDPM is tested for its accuracy and efficiency in both phantoms and realistic patient cases. Results: In all cases, the average relative uncertainties are less than 1%. Dose difference between CPU and GPU results is less than 3% of maximum dose in over 98% high dose region. A t‐test is also performed and the dose difference is not statistically significant in over 98% high dose region. Speed up factors of 57.3∼75.5 have been observed using an NVIDIA Tesla C2050 GPU card against a 2.27GHz Intel Xeon CPU processor. Dose calculation for a real IMRT head‐and‐neck plan and a VMAT prostate plan can be finished in 48sec and 46sec, respectively. Conclusions: The gDPM package has achieved high computation efficiency on GPU without losing accuracy. This work is supported in part by the University of California Lab Fees Research Program. GPU cards are partly provided by NVIDIA.

Original languageEnglish (US)
Pages (from-to)3768
Number of pages1
JournalMedical Physics
Volume38
Issue number6
DOIs
StatePublished - 2011

Fingerprint

Physics
Photons
Radiotherapy
Fees and Charges
Electron Transport
Uncertainty
Prostate
Neoplasms
Electrons
Research
Therapeutics

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

TU‐E‐BRB‐09 : Development of a GPU‐Based Monte Carlo Dose Calculation Package for Cancer Radiotherapy. / Jia, X.; gu, X.; Graves, Y.; Jiang, S.

In: Medical Physics, Vol. 38, No. 6, 2011, p. 3768.

Research output: Contribution to journalArticle

@article{a600f62d2bcc4cb7bb150e91391ae1b7,
title = "TU‐E‐BRB‐09: Development of a GPU‐Based Monte Carlo Dose Calculation Package for Cancer Radiotherapy",
abstract = "Purpose: Monte Carlo (MC) simulation is commonly considered as the most accurate dose calculation method in radiotherapy. However, its efficiency still requires improvement for many routine clinical applications. The goal of this work is to develop a GPU‐based MC dose calculation package, gDPM. It utilizes the parallel computation ability of a GPU to achieve high efficiency, while maintaining same particle transport physics as in the original DPM code and hence simulation accuracy.Methods: The gDPM code is implemented on GPU architecture under CUDA platform. In GPU computing, divergence of execution paths between threads can considerably reduce the efficiency. Since photons and electrons undergo different physics and hence attain different execution path, we design a simulation scheme where photon transport and electron transport are separated to partially relieve the thread divergence issue. High performance random number generator and hardware linear interpolation are also utilized. We have also developed various components to account for fluence map and linac geometry, so that gDPM can be used to compute dose for real IMRT or VMAT treatment plans. The gDPM is tested for its accuracy and efficiency in both phantoms and realistic patient cases. Results: In all cases, the average relative uncertainties are less than 1{\%}. Dose difference between CPU and GPU results is less than 3{\%} of maximum dose in over 98{\%} high dose region. A t‐test is also performed and the dose difference is not statistically significant in over 98{\%} high dose region. Speed up factors of 57.3∼75.5 have been observed using an NVIDIA Tesla C2050 GPU card against a 2.27GHz Intel Xeon CPU processor. Dose calculation for a real IMRT head‐and‐neck plan and a VMAT prostate plan can be finished in 48sec and 46sec, respectively. Conclusions: The gDPM package has achieved high computation efficiency on GPU without losing accuracy. This work is supported in part by the University of California Lab Fees Research Program. GPU cards are partly provided by NVIDIA.",
author = "X. Jia and X. gu and Y. Graves and S. Jiang",
year = "2011",
doi = "10.1118/1.3613185",
language = "English (US)",
volume = "38",
pages = "3768",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "6",

}

TY - JOUR

T1 - TU‐E‐BRB‐09

T2 - Development of a GPU‐Based Monte Carlo Dose Calculation Package for Cancer Radiotherapy

AU - Jia, X.

AU - gu, X.

AU - Graves, Y.

AU - Jiang, S.

PY - 2011

Y1 - 2011

N2 - Purpose: Monte Carlo (MC) simulation is commonly considered as the most accurate dose calculation method in radiotherapy. However, its efficiency still requires improvement for many routine clinical applications. The goal of this work is to develop a GPU‐based MC dose calculation package, gDPM. It utilizes the parallel computation ability of a GPU to achieve high efficiency, while maintaining same particle transport physics as in the original DPM code and hence simulation accuracy.Methods: The gDPM code is implemented on GPU architecture under CUDA platform. In GPU computing, divergence of execution paths between threads can considerably reduce the efficiency. Since photons and electrons undergo different physics and hence attain different execution path, we design a simulation scheme where photon transport and electron transport are separated to partially relieve the thread divergence issue. High performance random number generator and hardware linear interpolation are also utilized. We have also developed various components to account for fluence map and linac geometry, so that gDPM can be used to compute dose for real IMRT or VMAT treatment plans. The gDPM is tested for its accuracy and efficiency in both phantoms and realistic patient cases. Results: In all cases, the average relative uncertainties are less than 1%. Dose difference between CPU and GPU results is less than 3% of maximum dose in over 98% high dose region. A t‐test is also performed and the dose difference is not statistically significant in over 98% high dose region. Speed up factors of 57.3∼75.5 have been observed using an NVIDIA Tesla C2050 GPU card against a 2.27GHz Intel Xeon CPU processor. Dose calculation for a real IMRT head‐and‐neck plan and a VMAT prostate plan can be finished in 48sec and 46sec, respectively. Conclusions: The gDPM package has achieved high computation efficiency on GPU without losing accuracy. This work is supported in part by the University of California Lab Fees Research Program. GPU cards are partly provided by NVIDIA.

AB - Purpose: Monte Carlo (MC) simulation is commonly considered as the most accurate dose calculation method in radiotherapy. However, its efficiency still requires improvement for many routine clinical applications. The goal of this work is to develop a GPU‐based MC dose calculation package, gDPM. It utilizes the parallel computation ability of a GPU to achieve high efficiency, while maintaining same particle transport physics as in the original DPM code and hence simulation accuracy.Methods: The gDPM code is implemented on GPU architecture under CUDA platform. In GPU computing, divergence of execution paths between threads can considerably reduce the efficiency. Since photons and electrons undergo different physics and hence attain different execution path, we design a simulation scheme where photon transport and electron transport are separated to partially relieve the thread divergence issue. High performance random number generator and hardware linear interpolation are also utilized. We have also developed various components to account for fluence map and linac geometry, so that gDPM can be used to compute dose for real IMRT or VMAT treatment plans. The gDPM is tested for its accuracy and efficiency in both phantoms and realistic patient cases. Results: In all cases, the average relative uncertainties are less than 1%. Dose difference between CPU and GPU results is less than 3% of maximum dose in over 98% high dose region. A t‐test is also performed and the dose difference is not statistically significant in over 98% high dose region. Speed up factors of 57.3∼75.5 have been observed using an NVIDIA Tesla C2050 GPU card against a 2.27GHz Intel Xeon CPU processor. Dose calculation for a real IMRT head‐and‐neck plan and a VMAT prostate plan can be finished in 48sec and 46sec, respectively. Conclusions: The gDPM package has achieved high computation efficiency on GPU without losing accuracy. This work is supported in part by the University of California Lab Fees Research Program. GPU cards are partly provided by NVIDIA.

UR - http://www.scopus.com/inward/record.url?scp=85024796677&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85024796677&partnerID=8YFLogxK

U2 - 10.1118/1.3613185

DO - 10.1118/1.3613185

M3 - Article

AN - SCOPUS:85024796677

VL - 38

SP - 3768

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 6

ER -