A study of potential numerical pitfalls in GPU-based Monte Carlo dose calculation.

Abstract:

:The purpose of this study was to evaluate the impact of numerical errors caused by the floating point representation of real numbers in a GPU-based Monte Carlo code used for dose calculation in radiation oncology, and to identify situations where this type of error arises. The program used as a benchmark was bGPUMCD. Three tests were performed on the code, which was divided into three functional components: energy accumulation, particle tracking and physical interactions. First, the impact of single-precision calculations was assessed for each functional component. Second, a GPU-specific compilation option that reduces execution time as well as precision was examined. Third, a specific function used for tracking and potentially more sensitive to precision errors was tested by comparing it to a very high-precision implementation. Numerical errors were found in two components of the program. Because of the energy accumulation process, a few voxels surrounding a radiation source end up with a lower computed dose than they should. The tracking system contained a series of operations that abnormally amplify rounding errors in some situations. This resulted in some rare instances (less than 0.1%) of computed distances that are exceedingly far from what they should have been. Most errors detected had no significant effects on the result of a simulation due to its random nature, either because they cancel each other out or because they only affect a small fraction of particles. The results of this work can be extended to other types of GPU-based programs and be used as guidelines to avoid numerical errors on the GPU computing platform.

journal_name

Phys Med Biol

authors

Magnoux V,Ozell B,Bonenfant É,Després P

doi

10.1088/0031-9155/60/13/5007

subject

Has Abstract

pub_date

2015-07-07 00:00:00

pages

5007-18

issue

13

eissn

0031-9155

issn

1361-6560

journal_volume

60

pub_type

杂志文章
  • Dose uncertainty in radiotherapy of patients with head and neck cancer measured by in vivo ESR/alanine dosimetry using a mouthpiece.

    abstract::In order (i) to evaluate the dose uncertainty of the mouthpiece in daily use during intensity-modulated radiotherapy of patients with head and neck cancer, and (ii) to present a system for in vivo dosimetry of the oral mucosa, we equipped the mouthpiece with alanine dosimeter probes for in vivo dosimetry. The aim was ...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/56/5/010

    authors: Wagner D,Anton M,Vorwerk H

    更新日期:2011-03-07 00:00:00

  • Alanine dosimetry in strong magnetic fields: use as a transfer standard in MRI-guided radiotherapy.

    abstract::Reference dosimetry in the presence of a strong magnetic field is challenging. Ionisation chambers have shown to be strongly affected by magnetic fields. There is a need for robust and stable detectors in MRI-guided radiotherapy (MRIgRT). This study investigates the behaviour of the alanine dosimeter in magnetic field...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/ab8148

    authors: Billas I,Bouchard H,Oelfke U,Shipley D,Gouldstone C,Duane S

    更新日期:2020-06-03 00:00:00

  • Deep learning on mammary glands distribution for architectural distortion detection in digital breast tomosynthesis.

    abstract::Computer aided detection (CADe) for breast lesion can provide an important reference for radiologists in breast cancer screening. Architectural distortion (AD) is a type of breast lesion that is difficult to detect. A majority of CADe methods focus on detecting the radial pattern, which is a main characteristic of typ...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/ab98d0

    authors: Li Y,He Z,Lu Y,Ma X,Guo Y,Xie Z,Xu Z,Chen W,Chen H

    更新日期:2020-06-02 00:00:00

  • Perturbation correction for alanine dosimeters in different phantom materials in high-energy photon beams.

    abstract::In modern radiotherapy the verification of complex treatments plans is often performed in inhomogeneous or even anthropomorphic phantoms. For dose verification small detectors are necessary and therefore alanine detectors are most suitable. Though the response of alanine for a wide range of clinical photon energies in...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/61/3/N70

    authors: von Voigts-Rhetz P,Anton M,Vorwerk H,Zink K

    更新日期:2016-02-07 00:00:00

  • A theoretical study of quantitative flow measurements with constant infusion of short-lived isotopes.

    abstract::Recent developments in the technique of constant infusion and in emission computed tomography have provided a new method for quantitative regional flow measurements. This paper re-examines critically the theoretical relationship between flow and radioactivity during constant infusion of short-lived isotopes. Sensitivi...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/24/6/005

    authors: Huang SC,Phelps ME,Hoffman EJ,Kuhl DE

    更新日期:1979-11-01 00:00:00

  • Analysis on reversible/irreversible electroporation region in lung adenocarcinoma cell modelin vitrowith electric pulses delivered by needle electrodes.

    abstract::Irreversible electroporation (IRE) is a minimally invasive tumor therapy using pulsed electric field with high intensity while the important tissues such as blood vessel, bile duct, and nerve are preserved. In addition to ablation area, reversible electroporation (RE) region is also generated using needle electrodes f...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/abc12e

    authors: Lv Y,Tang X,Peng W,Cheng X,Chen S,Yao C

    更新日期:2020-10-14 00:00:00

  • Optimisation of the imaging and dosimetric characteristics of an electronic portal imaging device employing plastic scintillating fibres using Monte Carlo simulations.

    abstract::A Monte Carlo model of a novel electronic portal imaging device (EPID) has been developed using Geant4 and its performance for imaging and dosimetry applications in radiotherapy has been characterised. The EPID geometry is based on a physical prototype under ongoing investigation and comprises an array of plastic scin...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/59/22/6827

    authors: Blake SJ,McNamara AL,Vial P,Holloway L,Kuncic Z

    更新日期:2014-11-21 00:00:00

  • Stochastic versus deterministic kernel-based superposition approaches for dose calculation of intensity-modulated arcs.

    abstract::Dose calculations for radiation arc therapy are traditionally performed by approximating continuous delivery arcs with multiple static beams. For 3D conformal arc treatments, the shape and weight variation per degree is usually small enough to allow arcs to be approximated by static beams separated by 5 degrees -10 de...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/53/17/018

    authors: Tang G,Earl MA,Luan S,Wang C,Cao D,Yu CX,Naqvi SA

    更新日期:2008-09-07 00:00:00

  • Monitor unit calculation on the beam axis of open and wedged asymmetric high-energy photon beams.

    abstract::An ESTRO booklet and a report of the Netherlands Commission on Radiation Dosimetry have been published recently describing empirical methods for monitor unit (MU) calculations in symmetrical high-energy photon beams. Both documents support the same basic ideas; firstly the separation of head scatter and volume scatter...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/44/12/310

    authors: Georg D

    更新日期:1999-12-01 00:00:00

  • Novel adaptive beam-dependent margins for additional OAR sparing.

    abstract::Margins are employed in radiotherapy treatment planning to mitigate the dosimetric effects of geometric uncertainties for the clinical target volume (CTV). Unfortunately, whilst the use of margins can increase the probability that sufficient dose is delivered to the CTV, it can also result in delivering high dose of r...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/aae658

    authors: Tsang HS,Kamerling CP,Ziegenhein P,Nill S,Oelfke U

    更新日期:2018-10-29 00:00:00

  • Investigating the effect of longitudinal micro-CT imaging on tumour growth in mice.

    abstract::The aim of this study is to determine the impact of longitudinal micro-CT imaging on the growth of B16F1 tumours in C57BL/6 mice. Sixty mice received 2 × 10(5) B16F1 cells subcutaneously in the hind flank and were divided into control (no scan), 'low-dose' (80 kVp, 70 mA, 8 s, 0.07 Gy), 'medium-dose' (80 kVp, 50 mA, 3...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/56/2/002

    authors: Foster WK,Ford NL

    更新日期:2011-01-21 00:00:00

  • Dosimetry comparison and characterisation of an Al K ultrasoft x-ray beam from an MRC cold-cathode source.

    abstract::Ultrasoft x-rays of 0.3-5 keV have provided a unique tool for the investigation of intracellular mechanisms of radiation action in biological organisms, including mammalian cells. However, their use presents unique practical problems in dosimetry and experimental design. Detailed interpretation of the biological resul...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/30/10/002

    authors: Hoshi M,Goodhead DT,Brenner DJ,Bance DA,Chmielewski JJ,Paciotti MA,Bradbury JN

    更新日期:1985-10-01 00:00:00

  • A multi-scale framework with unsupervised joint training of convolutional neural networks for pulmonary deformable image registration.

    abstract::To achieve accurate and fast deformable image registration (DIR) for pulmonary CT, we proposed a Multi-scale DIR framework with unsupervised Joint training of Convolutional Neural Network (MJ-CNN). MJ-CNN contains three models at multi-scale levels for a coarse-to-fine DIR to avoid being trapped in a local minimum. It...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/ab5da0

    authors: Jiang Z,Yin FF,Ge Y,Ren L

    更新日期:2020-01-13 00:00:00

  • Experimental verification of a model for predicting transient temperature distributions by focused ultrasound.

    abstract::A thermal model for predicting time-dependent temperature distributions during ultrasound heating was tested quantitatively. The relevant thermal processes incorporated in the model are heat conduction and ultrasound power absorption, and the required input parameters include the absolute ultrasound power, the shape a...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/36/9/001

    authors: Chin RB,Zagzebski JA,Madsen EL

    更新日期:1991-09-01 00:00:00

  • Separating components of variation in measurement series using maximum likelihood estimation. Application to patient position data in radiotherapy.

    abstract::Maximum likelihood estimation (MLE) is presented as a statistical tool to evaluate the contribution of measurement error to any measurement series where the same quantity is measured using different independent methods. The technique was tested against artificial data sets; generated for values of underlying variation...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/59/20/6019

    authors: Sage JP,Mayles WP,Mayles HM,Syndikus I

    更新日期:2014-10-21 00:00:00

  • Independent Monte-Carlo dose calculation for MLC based CyberKnife radiotherapy.

    abstract::This work aims to develop, implement and validate a Monte Carlo (MC)-based independent dose calculation (IDC) framework to perform patient-specific quality assurance (QA) for multi-leaf collimator (MLC)-based CyberKnife® (Accuray Inc., Sunnyvale, CA) treatment plans. The IDC framework uses an XML-format treatment plan...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/aa97f8

    authors: Mackeprang PH,Vuong D,Volken W,Henzen D,Schmidhalter D,Malthaner M,Mueller S,Frei D,Stampanoni MFM,Dal Pra A,Aebersold DM,Fix MK,Manser P

    更新日期:2017-12-19 00:00:00

  • Experimental 3D dosimetry around a high-dose-rate clinical 192Ir source using a polyacrylamide gel (PAG) dosimeter.

    abstract::It is well known that the experimental dosimetry of brachytherapy sources presents a challenge. Depending on the particular-dosimeter used, measurements can suffer from poor spatial resolution (ion chambers), lack of 3D information (film) or errors due to the presence of the dosimeter itself distorting the radiation f...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/44/10/305

    authors: McJury M,Tapper PD,Cosgrove VP,Murphy PS,Griffin S,Leach MO,Webb S,Oldham M

    更新日期:1999-10-01 00:00:00

  • Investigation of ultrasonic properties of PAG and MAGIC polymer gel dosimeters.

    abstract::Ultrasonic speed of propagation and attenuation were investigated as a function of absorbed radiation dose in PAG and MAGIC polymer gel dosimeters. Both PAG and MAGIC gel dosimeters displayed a dependence of ultrasonic parameters on absorbed dose with attenuation displaying significant changes in the dose range invest...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/47/24/307

    authors: Mather ML,De Deene Y,Whittaker AK,Simon GP,Rutgers R,Baldock C

    更新日期:2002-12-21 00:00:00

  • Detector arrangement and sampling characteristics in rotary positron-emission computed tomography.

    abstract::A rotary positron-emission computed tomography system has been proposed in which detectors are mounted unequally spaced on a ring, so as to provide fine and uniform linear sampling during a 360 degree rotation of the ring around its centre. The usefulness of a simple iterative method for searching for a suitable arran...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/26/3/011

    authors: Yamamoto M

    更新日期:1981-05-01 00:00:00

  • Experimental characterization of lateral profiles of scanned proton and carbon ion pencil beams for improved beam models in ion therapy treatment planning.

    abstract::Scanned ion pencil beams carry a low-dose envelope which can extend up to several centimeters from the individual beam central axis. Depending on the energy and species of the beam, this halo consists mainly of secondary particles produced by nuclear interactions in the target or of particles undergoing multiple Coulo...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/56/24/009

    authors: Schwaab J,Brons S,Fieres J,Parodi K

    更新日期:2011-12-21 00:00:00

  • A new approach to develop computer-aided detection schemes of digital mammograms.

    abstract::The purpose of this study is to develop a new global mammographic image feature analysis based computer-aided detection (CAD) scheme and evaluate its performance in detecting positive screening mammography examinations. A dataset that includes images acquired from 1896 full-field digital mammography (FFDM) screening e...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/60/11/4413

    authors: Tan M,Qian W,Pu J,Liu H,Zheng B

    更新日期:2015-06-07 00:00:00

  • Evaluation of high performance data acquisition boards for simultaneous sampling of fast signals from PET detectors.

    abstract::Detectors used for positron emission tomography (PET) provide fast, randomly distributed signals that need to be digitized for further processing. One possibility is to sample the signals at the peak initiated by a trigger from a constant fraction discriminator (CFD). For PET detectors, simultaneous acquisition of man...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/50/1/003

    authors: Judenhofer MS,Pichler BJ,Cherry SR

    更新日期:2005-01-07 00:00:00

  • Assessment of noise in a digital image using the join-count statistic and the Moran test.

    abstract::We assume that the data bits of a pixel in digital images can be divided into signal and noise bits. The signal bits occupy the most significant part of the pixel and the noise bits the least significant part. The signal part of each pixel are correlated while the noise parts are uncorrelated. Two statistical methods,...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/37/2/004

    authors: Chuang KS,Huang HK

    更新日期:1992-02-01 00:00:00

  • Internal dosimetry with the Monte Carlo code GATE: validation using the ICRP/ICRU female reference computational model.

    abstract::The purpose of this work was to validate GATE-based clinical scale absorbed dose calculations in nuclear medicine dosimetry. GATE (version 6.2) and MCNPX (version 2.7.a) were used to derive dosimetric parameters (absorbed fractions, specific absorbed fractions and S-values) for the reference female computational model...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/62/5/1885

    authors: Villoing D,Marcatili S,Garcia MP,Bardiès M

    更新日期:2017-03-07 00:00:00

  • Laser shaping of corneal transplants in vitro: area ablation with small overlapping laser spots produced by a pulsed scanning laser beam using an optimizing ablation algorithm.

    abstract::Area laser lathing and trephination of donor corneas is used to produce different types of grafts for human transplantation. 193 nm (ArF excimer) laser radiation is used, since this is known to give a non-thermal laser-tissue interaction with a minimal zone of tissue damage. To guarantee the highest degree of flexibil...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/44/5/005

    authors: Homolka P,Biowski R,Kaminski S,Barisani T,Husinsky W,Bergmann H,Grabner G

    更新日期:1999-05-01 00:00:00

  • Pitch dependence of longitudinal sampling and aliasing effects in multi-slice helical computed tomography (CT).

    abstract::In this work, we investigate longitudinal sampling and aliasing effects in multi-slice helical CT. We demonstrate that longitudinal aliasing can be a significant, complicated, and potentially detrimental effect in multi-slice helical CT reconstructions. Multi-slice helical CT scans are generally undersampled longitudi...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/47/15/318

    authors: La Rivière PJ,Pan X

    更新日期:2002-08-07 00:00:00

  • Pulse shape discrimination and classification methods for continuous depth of interaction encoding PET detectors.

    abstract::In previous work we demonstrated the potential of positron emission tomography (PET) detectors with depth-of-interaction (DOI) encoding capability based on phosphor-coated crystals. A DOI resolution of 8 mm full-width at half-maximum was obtained for 20 mm long scintillator crystals using a delayed charge integration ...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/57/20/6571

    authors: Roncali E,Phipps JE,Marcu L,Cherry SR

    更新日期:2012-10-21 00:00:00

  • Derivation of mean dose tolerances for new fractionation schemes and treatment modalities.

    abstract::Avoiding toxicities in radiotherapy requires the knowledge of tolerable organ doses. For new, experimental fractionation schemes (e.g. hypofractionation) these are typically derived from traditional schedules using the biologically effective dose (BED) model. In this report we investigate the difficulties of establish...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/1361-6560/aa9836

    authors: Perkó Z,Bortfeld T,Hong T,Wolfgang J,Unkelbach J

    更新日期:2018-02-05 00:00:00

  • Cross-sections for the inactivation of ribonuclease by slow heavy ions.

    abstract::Average cross-sections for inactivation of dry ribonuclease by H, He and N projectiles with energies between 100 eV and 10 keV are presented. The trend of the damage cross-sections with particle energy is consistent with the trend of the molecular quasi-elastic scattering and ionization cross-sections. A tentative rad...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/20/6/005

    authors: Watt DE,Sutcliffe JF

    更新日期:1975-11-01 00:00:00

  • The modulation transfer function of an optical coherence tomography imaging system in turbid media.

    abstract::In this paper we describe measurements of the contrast transfer function, modulation transfer function and point-spread function of an optical coherence tomography (OCT) imaging system through scattering layers having a dimension-less scattering depth over the range 0.2-6.9. The results were found to be insensitive to...

    journal_title:Physics in medicine and biology

    pub_type: 杂志文章

    doi:10.1088/0031-9155/56/9/014

    authors: Woolliams PD,Tomlins PH

    更新日期:2011-05-07 00:00:00