A vectorized Monte Carlo code for modeling photon transport in SPECT.

Published

Journal Article

A vectorized Monte Carlo computer code has been developed for modeling photon transport in single photon emission computed tomography (SPECT). The code models photon transport in a uniform attenuating region and photon detection by a gamma camera. It is adapted from a history-based Monte Carlo code in which photon history data are stored in scalar variables and photon histories are computed sequentially. The vectorized code is written in FORTRAN77 and uses an event-based algorithm in which photon history data are stored in arrays and photon history computations are performed within DO loops. The indices of the DO loops range over the number of photon histories, and these loops may take advantage of the vector processing unit of our Stellar GS1000 computer for pipelined computations. Without the use of the vector processor the event-based code is faster than the history-based code because of numerical optimization performed during conversion to the event-based algorithm. When only the detection of unscattered photons is modeled, the event-based code executes 5.1 times faster with the use of the vector processor than without; when the detection of scattered and unscattered photons is modeled the speed increase is a factor of 2.9. Vectorization is a valuable way to increase the performance of Monte Carlo code for modeling photon transport in SPECT.

Full Text

Duke Authors

Cited Authors

  • Smith, MF; Floyd, CE; Jaszczak, RJ

Published Date

  • July 1993

Published In

Volume / Issue

  • 20 / 4

Start / End Page

  • 1121 - 1127

PubMed ID

  • 8413021

Pubmed Central ID

  • 8413021

Electronic International Standard Serial Number (EISSN)

  • 2473-4209

International Standard Serial Number (ISSN)

  • 0094-2405

Digital Object Identifier (DOI)

  • 10.1118/1.597148

Language

  • eng