Qifan Xu, Qihui Lyu, Lu Jiang, Shusen Jing, Dan Ruan, Ke Sheng
{"title":"一种用于复杂4π放疗规划的超高性能并行(UHPP)框架。","authors":"Qifan Xu, Qihui Lyu, Lu Jiang, Shusen Jing, Dan Ruan, Ke Sheng","doi":"10.1088/1361-6560/adf2f2","DOIUrl":null,"url":null,"abstract":"<p><strong>Purpose: </strong>In radiotherapy, dose distribution conformity and compactness are critical to patient outcomes. Advanced techniques like 4π radiotherapy leverage non-coplanar beams for superior dosimetry by exploring additional degrees of freedom. However, 4π planning is computationally intensive due to large dose-loading matrices for candidate beams. This work presents an ultra-high performance parallel (UHPP) framework to accelerate high-dimensional treatment planning.

Methods. For dose calculation, we developed: 1) A two-step TERMA computation module calculating the TERMA array once per beam, enabling reuse across convolution directions; 2) A synchronized dose calculation module based on collapsed-cone convolution superposition (CCCS), arranging rays in dedicated sequences to preserve thread efficiency and minimize memory access; 3) A scattering-based coordinate transformation mapping dose from beamlet to patient Cartesian coordinates, eliminating aliasing without atomic operations. The framework includes CCCS exponential kernel calculation for varying LINAC spectra. For beam orientation optimization, we employed fast iterative shrinkage-thresholding algorithm (FISTA) with group sparsity regularization, accelerated using cuSPARSE library on GPUs. We benchmarked against Monte Carlo (MC) simulations for dose accuracy and compared computational performance to state-of-the-art (SOTA) methods. Plan quality was evaluated across four approaches: UHPP, SOTA, clinical VMAT plans, and MC calculations based on UHPP plans.

Results. Compared to MC simulations, UHPP achieved minimum 98% gamma passing rates under 1.5%/1.5mm criterion for water and slab phantoms, and average 97.35% and 92.18% under 3%/3mm criterion for pancreas and head-and-neck patients, respectively. UHPP delivered 8.86× and 6.99× speedups in dose calculation and plan optimization while maintaining comparable or superior plan quality. Both UHPP and SOTA consistently produced 4π plans outperforming clinical VMAT plans in organ-at-risk sparing and target coverage.

Conclusion. The UHPP framework delivers high dose accuracy and substantial computational speedup without sacrificing 4π planning's dosimetric advantages, supporting practical adoption of advanced 4π radiotherapy in clinical workflows.</p>","PeriodicalId":20185,"journal":{"name":"Physics in medicine and biology","volume":" ","pages":""},"PeriodicalIF":3.3000,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"An ultra-high performance parallel (UHPP) framework for complex 4π radiotherapy planning.\",\"authors\":\"Qifan Xu, Qihui Lyu, Lu Jiang, Shusen Jing, Dan Ruan, Ke Sheng\",\"doi\":\"10.1088/1361-6560/adf2f2\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<p><strong>Purpose: </strong>In radiotherapy, dose distribution conformity and compactness are critical to patient outcomes. Advanced techniques like 4π radiotherapy leverage non-coplanar beams for superior dosimetry by exploring additional degrees of freedom. However, 4π planning is computationally intensive due to large dose-loading matrices for candidate beams. This work presents an ultra-high performance parallel (UHPP) framework to accelerate high-dimensional treatment planning.

Methods. For dose calculation, we developed: 1) A two-step TERMA computation module calculating the TERMA array once per beam, enabling reuse across convolution directions; 2) A synchronized dose calculation module based on collapsed-cone convolution superposition (CCCS), arranging rays in dedicated sequences to preserve thread efficiency and minimize memory access; 3) A scattering-based coordinate transformation mapping dose from beamlet to patient Cartesian coordinates, eliminating aliasing without atomic operations. The framework includes CCCS exponential kernel calculation for varying LINAC spectra. For beam orientation optimization, we employed fast iterative shrinkage-thresholding algorithm (FISTA) with group sparsity regularization, accelerated using cuSPARSE library on GPUs. We benchmarked against Monte Carlo (MC) simulations for dose accuracy and compared computational performance to state-of-the-art (SOTA) methods. Plan quality was evaluated across four approaches: UHPP, SOTA, clinical VMAT plans, and MC calculations based on UHPP plans.

Results. Compared to MC simulations, UHPP achieved minimum 98% gamma passing rates under 1.5%/1.5mm criterion for water and slab phantoms, and average 97.35% and 92.18% under 3%/3mm criterion for pancreas and head-and-neck patients, respectively. UHPP delivered 8.86× and 6.99× speedups in dose calculation and plan optimization while maintaining comparable or superior plan quality. Both UHPP and SOTA consistently produced 4π plans outperforming clinical VMAT plans in organ-at-risk sparing and target coverage.

Conclusion. The UHPP framework delivers high dose accuracy and substantial computational speedup without sacrificing 4π planning's dosimetric advantages, supporting practical adoption of advanced 4π radiotherapy in clinical workflows.</p>\",\"PeriodicalId\":20185,\"journal\":{\"name\":\"Physics in medicine and biology\",\"volume\":\" \",\"pages\":\"\"},\"PeriodicalIF\":3.3000,\"publicationDate\":\"2025-07-22\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Physics in medicine and biology\",\"FirstCategoryId\":\"5\",\"ListUrlMain\":\"https://doi.org/10.1088/1361-6560/adf2f2\",\"RegionNum\":3,\"RegionCategory\":\"医学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"ENGINEERING, BIOMEDICAL\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Physics in medicine and biology","FirstCategoryId":"5","ListUrlMain":"https://doi.org/10.1088/1361-6560/adf2f2","RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"ENGINEERING, BIOMEDICAL","Score":null,"Total":0}
An ultra-high performance parallel (UHPP) framework for complex 4π radiotherapy planning.
Purpose: In radiotherapy, dose distribution conformity and compactness are critical to patient outcomes. Advanced techniques like 4π radiotherapy leverage non-coplanar beams for superior dosimetry by exploring additional degrees of freedom. However, 4π planning is computationally intensive due to large dose-loading matrices for candidate beams. This work presents an ultra-high performance parallel (UHPP) framework to accelerate high-dimensional treatment planning.
Methods. For dose calculation, we developed: 1) A two-step TERMA computation module calculating the TERMA array once per beam, enabling reuse across convolution directions; 2) A synchronized dose calculation module based on collapsed-cone convolution superposition (CCCS), arranging rays in dedicated sequences to preserve thread efficiency and minimize memory access; 3) A scattering-based coordinate transformation mapping dose from beamlet to patient Cartesian coordinates, eliminating aliasing without atomic operations. The framework includes CCCS exponential kernel calculation for varying LINAC spectra. For beam orientation optimization, we employed fast iterative shrinkage-thresholding algorithm (FISTA) with group sparsity regularization, accelerated using cuSPARSE library on GPUs. We benchmarked against Monte Carlo (MC) simulations for dose accuracy and compared computational performance to state-of-the-art (SOTA) methods. Plan quality was evaluated across four approaches: UHPP, SOTA, clinical VMAT plans, and MC calculations based on UHPP plans.
Results. Compared to MC simulations, UHPP achieved minimum 98% gamma passing rates under 1.5%/1.5mm criterion for water and slab phantoms, and average 97.35% and 92.18% under 3%/3mm criterion for pancreas and head-and-neck patients, respectively. UHPP delivered 8.86× and 6.99× speedups in dose calculation and plan optimization while maintaining comparable or superior plan quality. Both UHPP and SOTA consistently produced 4π plans outperforming clinical VMAT plans in organ-at-risk sparing and target coverage.
Conclusion. The UHPP framework delivers high dose accuracy and substantial computational speedup without sacrificing 4π planning's dosimetric advantages, supporting practical adoption of advanced 4π radiotherapy in clinical workflows.
期刊介绍:
The development and application of theoretical, computational and experimental physics to medicine, physiology and biology. Topics covered are: therapy physics (including ionizing and non-ionizing radiation); biomedical imaging (e.g. x-ray, magnetic resonance, ultrasound, optical and nuclear imaging); image-guided interventions; image reconstruction and analysis (including kinetic modelling); artificial intelligence in biomedical physics and analysis; nanoparticles in imaging and therapy; radiobiology; radiation protection and patient dose monitoring; radiation dosimetry