SUGTeP115: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation
Abstract
Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrixvector and vectorvector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order selfadjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positivedefinite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDAC was benchmarked on an Nvidia GeForce TITANX GPU against an Intel i76700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memorymore »
 Authors:
 University of Alberta, Edmonton, AB (Canada)
 (Canada)
 Publication Date:
 OSTI Identifier:
 22649355
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Medical Physics; Journal Volume: 43; Journal Issue: 6; Other Information: (c) 2016 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 97 MATHEMATICAL METHODS AND COMPUTING; 60 APPLIED LIFE SCIENCES; BOLTZMANN EQUATION; FINITE ELEMENT METHOD; LEAST SQUARE FIT; LINEAR ACCELERATORS; MATHEMATICAL SOLUTIONS; MATRICES; PRODUCTIVITY
Citation Formats
Yang, R, Fallone, B, Cross Cancer Institute, Edmonton, AB, MagnetTx Oncology Solutions, Edmonton, AB, St Aubin, J, and Cross Cancer Institute, Edmonton, AB. SUGTeP115: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation. United States: N. p., 2016.
Web. doi:10.1118/1.4957005.
Yang, R, Fallone, B, Cross Cancer Institute, Edmonton, AB, MagnetTx Oncology Solutions, Edmonton, AB, St Aubin, J, & Cross Cancer Institute, Edmonton, AB. SUGTeP115: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation. United States. doi:10.1118/1.4957005.
Yang, R, Fallone, B, Cross Cancer Institute, Edmonton, AB, MagnetTx Oncology Solutions, Edmonton, AB, St Aubin, J, and Cross Cancer Institute, Edmonton, AB. 2016.
"SUGTeP115: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation". United States.
doi:10.1118/1.4957005.
@article{osti_22649355,
title = {SUGTeP115: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation},
author = {Yang, R and Fallone, B and Cross Cancer Institute, Edmonton, AB and MagnetTx Oncology Solutions, Edmonton, AB and St Aubin, J and Cross Cancer Institute, Edmonton, AB},
abstractNote = {Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrixvector and vectorvector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order selfadjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positivedefinite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDAC was benchmarked on an Nvidia GeForce TITANX GPU against an Intel i76700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memory mapping for the kronecker product between the sparse spatial and angular matrices, including Dirichlet boundary conditions, was created. Atomicity is preserved by graphcoloring overlapping nodes into separate kernel launches. The onetime mapping calculations for matrix assembly, kronecker product, and boundary condition application took 452±1ms on GPU. Matrix assembly for 16 energy groups took 556±3s on CPU, and 358±2ms on GPU using the mappings developed. The CG solver took 93±1s on CPU, and 468±2ms on GPU. Conclusion: Three computationally intensive subroutines in deterministically solving the LBTE have been formulated on GPU, resulting in two orders of magnitude speedup. Funding support from Natural Sciences and Engineering Research Council and Alberta Innovates Health Solutions. Dr. Fallone is a cofounder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta biplanar linac MR for commercialization).},
doi = {10.1118/1.4957005},
journal = {Medical Physics},
number = 6,
volume = 43,
place = {United States},
year = 2016,
month = 6
}

Purpose: Accurate radiotherapy dose calculation algorithms are essential to any successful radiotherapy program, considering the high level of dose conformity and modulation in many of today’s treatment plans. As technology continues to progress, such as is the case with novel MRIguided radiotherapy systems, the necessity for dose calculation algorithms to accurately predict delivered dose in increasingly challenging scenarios is vital. To this end, a novel deterministic solution has been developed to the first order linear Boltzmann transport equation which accurately calculates xray based radiotherapy doses in the presence of magnetic fields. Methods: The deterministic formalism discussed here with the inclusionmore »

Accelerated solution of discrete ordinates approximation to the Boltzmann transport equation via model reduction
Radiation heat transfer is an important phenomenon in many physical systems of practical interest. When participating media is important, the radiative transfer equation (RTE) must be solved for the radiative intensity as a function of location, time, direction, and wavelength. In many heattransfer applications, a quasisteady assumption is valid, thereby removing time dependence. The dependence on wavelength is often treated through a weighted sum of gray gases (WSGG) approach. The discrete ordinates method (DOM) is one of the most common methods for approximating the angular (i.e., directional) dependence. The DOM exactly solves for the radiative intensity for a finite numbermore » 
Accelerated solution of the Boltzmann equation
Methods are presented for accelerating a numerical procedure for selfconsistent solution of a kinetic equation and Poisson's equation in Plasma simulation. The kinetic equation is solved using a propagator technique, although other approaches would also benefit from the accelerated procedure. The kinetic equation is solved in a phase space of at least one spatial variable and two velocity coordinates, (Z, V[sub z], V[sub p]). V[sub p] is in the direction perpendicular to Z. In these variables it is possible to advance the reduced distribution function g(Z, V[sub z]) in time very efficiently using several [open quotes]short[close quotes] time steps withinmore » 
An asymptoticpreserving scheme for the semiconductor Boltzmann equation toward the energytransport limit
We design an asymptoticpreserving scheme for the semiconductor Boltzmann equation which leads to an energytransport system for electron mass and energy as mean free path goes to zero. As opposed to the classical driftdiffusion limit where the stiff collisions are all in one scale, new difficulties arise in the twoscale stiff collision terms because the simple BGK penalization [15] fails to drive the solution to the correct limit. We propose to set up a spatially dependent threshold on the penalization of the stiffer collision operator such that the evolution of the solution resembles a Hilbert expansion at the continuous level.more »