Squared polynomial extrapolation methods with cycling: An application to the positron emission tomography problem

Research output: Contribution to journalArticle

Abstract

Roland and Varadhan (Appl. Numer. Math., 55:215-226, 2005) presented a new idea called "squaring" to improve the convergence of Lemaréchal's scheme for solving nonlinear fixed-point problems. Varadhan and Roland (Squared extrapolation methods: A new class of simple and efficient numerical schemes for accelerating the convergence of the EM algorithm, Department of Biostatistics Working Paper. Johns Hopkins University, http://www.bepress.com/jhubiostat/paper63, 2004) noted that Lemaréchal's scheme can be viewed as a member of the class of polynomial extrapolation methods with cycling that uses two fixed-point iterations per cycle. Here we combine these two ideas, cycled extrapolation and squaring, and construct a new class of methods, called squared polynomial methods (SQUAREM), for accelerating the convergence of fixed-point iterations. Our main goal is to evaluate whether the squaring device is effective in improving the rate of convergence of cycled extrapolation methods that use more than two fixed-point iterations per cycle. We study the behavior of the new schemes on an image reconstruction problem for positron emission tomography (PET) using simulated data. Our numerical experiments show the effectiveness of first- and higher-order squared polynomial extrapolation methods in accelerating image reconstruction, and also their relative superiority compared to the classical, "unsquared" vector polynomial methods.

Original languageEnglish (US)
Pages (from-to)159-172
Number of pages14
JournalNumerical Algorithms
Volume44
Issue number2
DOIs
StatePublished - Feb 2007

Keywords

  • Fixed-point methods
  • Linear systems
  • Nonlinear systems
  • Polynomial extrapolation methods
  • Squaring

ASJC Scopus subject areas

  • Applied Mathematics

Fingerprint Dive into the research topics of 'Squared polynomial extrapolation methods with cycling: An application to the positron emission tomography problem'. Together they form a unique fingerprint.

  • Cite this