SU‐E‐T‐476

GPU‐Based Monte Carlo Radiotherapy Dose Calculation Using Phase‐ Space Sources

R. Townson, Xun Jia, S. Zavgorodni, S. Jiang

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

Purpose: To design an efficient method for utilizing phase‐space source models in the GPU‐based Monte Carlo (MC) dose calculation engine gDPM. Methods: In GPU‐based MC algorithms, particles are transported in parallel on different threads. Particles of different types and energies can require significantly different execution times. This can cause “thread divergence” and lower efficiency when source particles are read sequentially from a phase‐space file. We have developed a strategy for utilizing phase‐ space files in a GPU compatible manner whereby the particles are grouped into phase‐space‐lets (PSLs) by type, energy, and location in the phase‐ space plane. This allows for dose calculations using only particles inside the field opening defined by the secondary collimators. For validation, the gDPM PSL implementation is compared with DOSXYZnrc using a BEAMnrc phase‐space source model as input. Results: Two phase‐spaces were generated using a BEAMnrc head model of a 6MV Varian Clinac 21EX, one above the upper jaws used to generate PSLs for gDPM and the other below the lower jaws used for DOSXYZnrc dose calculation. Profiles and depth dose curves for a variety of field sizes were generated in a water phantom. The agreement between gDPM and DOSXYZnrc is within 2% for all field sizes. For the 10 cm × 10 cm field, the calculation times of 650 million histories were 147 CPU hours and 54 GPU seconds for DOSXYZnrc and gDPM, respectively. In addition, we have tested the gDPM PSL implementation for dose calculation in a realistic 7‐field IMRT tongue treatment plan. The calculation times were 59 CPU‐hours and 66 GPU‐ seconds for DOSXYZnrc and gDPM for 485 million histories, respectively. Gamma pass rate for the two dose distributions was 99.54% for 3 mm/3% criteria within the 10% isodose. Conclusions: Methods for the efficient use of phase‐space sources for GPU‐based MC dose calculations have been developed.

Original languageEnglish (US)
Pages (from-to)3814
Number of pages1
JournalMedical Physics
Volume39
Issue number6
DOIs
StatePublished - 2012

Fingerprint

Radiotherapy
Maxilla
Jaw
Tongue
Head
Water

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

SU‐E‐T‐476 : GPU‐Based Monte Carlo Radiotherapy Dose Calculation Using Phase‐ Space Sources. / Townson, R.; Jia, Xun; Zavgorodni, S.; Jiang, S.

In: Medical Physics, Vol. 39, No. 6, 2012, p. 3814.

Research output: Contribution to journalArticle

@article{d23e663ae1cc40c7848134e4862632a6,
title = "SU‐E‐T‐476: GPU‐Based Monte Carlo Radiotherapy Dose Calculation Using Phase‐ Space Sources",
abstract = "Purpose: To design an efficient method for utilizing phase‐space source models in the GPU‐based Monte Carlo (MC) dose calculation engine gDPM. Methods: In GPU‐based MC algorithms, particles are transported in parallel on different threads. Particles of different types and energies can require significantly different execution times. This can cause “thread divergence” and lower efficiency when source particles are read sequentially from a phase‐space file. We have developed a strategy for utilizing phase‐ space files in a GPU compatible manner whereby the particles are grouped into phase‐space‐lets (PSLs) by type, energy, and location in the phase‐ space plane. This allows for dose calculations using only particles inside the field opening defined by the secondary collimators. For validation, the gDPM PSL implementation is compared with DOSXYZnrc using a BEAMnrc phase‐space source model as input. Results: Two phase‐spaces were generated using a BEAMnrc head model of a 6MV Varian Clinac 21EX, one above the upper jaws used to generate PSLs for gDPM and the other below the lower jaws used for DOSXYZnrc dose calculation. Profiles and depth dose curves for a variety of field sizes were generated in a water phantom. The agreement between gDPM and DOSXYZnrc is within 2{\%} for all field sizes. For the 10 cm × 10 cm field, the calculation times of 650 million histories were 147 CPU hours and 54 GPU seconds for DOSXYZnrc and gDPM, respectively. In addition, we have tested the gDPM PSL implementation for dose calculation in a realistic 7‐field IMRT tongue treatment plan. The calculation times were 59 CPU‐hours and 66 GPU‐ seconds for DOSXYZnrc and gDPM for 485 million histories, respectively. Gamma pass rate for the two dose distributions was 99.54{\%} for 3 mm/3{\%} criteria within the 10{\%} isodose. Conclusions: Methods for the efficient use of phase‐space sources for GPU‐based MC dose calculations have been developed.",
author = "R. Townson and Xun Jia and S. Zavgorodni and S. Jiang",
year = "2012",
doi = "10.1118/1.4735565",
language = "English (US)",
volume = "39",
pages = "3814",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "6",

}

TY - JOUR

T1 - SU‐E‐T‐476

T2 - GPU‐Based Monte Carlo Radiotherapy Dose Calculation Using Phase‐ Space Sources

AU - Townson, R.

AU - Jia, Xun

AU - Zavgorodni, S.

AU - Jiang, S.

PY - 2012

Y1 - 2012

N2 - Purpose: To design an efficient method for utilizing phase‐space source models in the GPU‐based Monte Carlo (MC) dose calculation engine gDPM. Methods: In GPU‐based MC algorithms, particles are transported in parallel on different threads. Particles of different types and energies can require significantly different execution times. This can cause “thread divergence” and lower efficiency when source particles are read sequentially from a phase‐space file. We have developed a strategy for utilizing phase‐ space files in a GPU compatible manner whereby the particles are grouped into phase‐space‐lets (PSLs) by type, energy, and location in the phase‐ space plane. This allows for dose calculations using only particles inside the field opening defined by the secondary collimators. For validation, the gDPM PSL implementation is compared with DOSXYZnrc using a BEAMnrc phase‐space source model as input. Results: Two phase‐spaces were generated using a BEAMnrc head model of a 6MV Varian Clinac 21EX, one above the upper jaws used to generate PSLs for gDPM and the other below the lower jaws used for DOSXYZnrc dose calculation. Profiles and depth dose curves for a variety of field sizes were generated in a water phantom. The agreement between gDPM and DOSXYZnrc is within 2% for all field sizes. For the 10 cm × 10 cm field, the calculation times of 650 million histories were 147 CPU hours and 54 GPU seconds for DOSXYZnrc and gDPM, respectively. In addition, we have tested the gDPM PSL implementation for dose calculation in a realistic 7‐field IMRT tongue treatment plan. The calculation times were 59 CPU‐hours and 66 GPU‐ seconds for DOSXYZnrc and gDPM for 485 million histories, respectively. Gamma pass rate for the two dose distributions was 99.54% for 3 mm/3% criteria within the 10% isodose. Conclusions: Methods for the efficient use of phase‐space sources for GPU‐based MC dose calculations have been developed.

AB - Purpose: To design an efficient method for utilizing phase‐space source models in the GPU‐based Monte Carlo (MC) dose calculation engine gDPM. Methods: In GPU‐based MC algorithms, particles are transported in parallel on different threads. Particles of different types and energies can require significantly different execution times. This can cause “thread divergence” and lower efficiency when source particles are read sequentially from a phase‐space file. We have developed a strategy for utilizing phase‐ space files in a GPU compatible manner whereby the particles are grouped into phase‐space‐lets (PSLs) by type, energy, and location in the phase‐ space plane. This allows for dose calculations using only particles inside the field opening defined by the secondary collimators. For validation, the gDPM PSL implementation is compared with DOSXYZnrc using a BEAMnrc phase‐space source model as input. Results: Two phase‐spaces were generated using a BEAMnrc head model of a 6MV Varian Clinac 21EX, one above the upper jaws used to generate PSLs for gDPM and the other below the lower jaws used for DOSXYZnrc dose calculation. Profiles and depth dose curves for a variety of field sizes were generated in a water phantom. The agreement between gDPM and DOSXYZnrc is within 2% for all field sizes. For the 10 cm × 10 cm field, the calculation times of 650 million histories were 147 CPU hours and 54 GPU seconds for DOSXYZnrc and gDPM, respectively. In addition, we have tested the gDPM PSL implementation for dose calculation in a realistic 7‐field IMRT tongue treatment plan. The calculation times were 59 CPU‐hours and 66 GPU‐ seconds for DOSXYZnrc and gDPM for 485 million histories, respectively. Gamma pass rate for the two dose distributions was 99.54% for 3 mm/3% criteria within the 10% isodose. Conclusions: Methods for the efficient use of phase‐space sources for GPU‐based MC dose calculations have been developed.

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

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

U2 - 10.1118/1.4735565

DO - 10.1118/1.4735565

M3 - Article

VL - 39

SP - 3814

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 6

ER -