Multicore Performance of Block Algebraic Iterative Reconstruction Methods

Research output: Contribution to journalJournal articleResearchpeer-review

349 Downloads (Pure)


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
Issue number5
Pages (from-to)C524-C546
Publication statusPublished - 2014


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


Dive into the research topics of 'Multicore Performance of Block Algebraic Iterative Reconstruction Methods'. Together they form a unique fingerprint.

Cite this