Multicore Performance of Block Algebraic Iterative Reconstruction Methods

Research output: Contribution to journalJournal articleResearchpeer-review

327 Downloads (Pure)

Abstract

Algebraic iterative methods are routinely used for solving the ill-posed sparse linear systems arising in tomographic image reconstruction. Here we consider the algebraic reconstruction technique (ART) and the simultaneous iterative reconstruction techniques (SIRT), both of which rely on semiconvergence. Block versions of these methods, based on a partitioning of the linear system, are able to combine the fast semiconvergence of ART with the better multicore properties of SIRT. These block methods separate into two classes: those that, in each iteration, access the blocks in a sequential manner, and those that compute a result for each block in parallel and then combine these results before the next iteration. The goal of this work is to demonstrate which block methods are best suited for implementation on modern multicore computers. To compare the performance of the different block methods, we use a fixed relaxation parameter in each method, namely, the one that leads to the fastest semiconvergence. Computational results show that for multicore computers, the sequential approach is preferable.
Original languageEnglish
JournalS I A M Journal on Scientific Computing
Volume36
Issue number5
Pages (from-to)C524-C546
ISSN1064-8275
DOIs
Publication statusPublished - 2014

Keywords

  • algebraic iterative reconstruction
  • ART
  • SIRT
  • block methods
  • relaxation parameter
  • semiconvergence
  • tomographic imaging

Cite this