Fixed Point Acceleration in R

Abstract:

A fixed point problem is one where we seek a vector, X, for a function, f, such that f(X) = X. The solution of many such problems can be accelerated by using a fixed point acceleration algorithm. With the release of the FixedPoint package there is now a number of algorithms available in R that can be used for accelerating the finding of a fixed point of a function. These algorithms include Newton acceleration, Aitken acceleration and Anderson acceleration as well as epsilon extrapolation methods and minimal polynomial methods. This paper demonstrates the use of fixed point accelerators in solving numerical mathematics problems using the algorithms of the FixedPoint package as well as the squarem method of the SQUAREM package.

Cite PDF Tweet

Authors

Affiliations

Stuart Baumann

 

Margaryta Klymak

 

Published

Aug. 19, 2019

Received

May 29, 2018

DOI

10.32614/RJ-2019-037

Volume

Pages

11/1

359 - 375

Supplementary materials

Supplementary materials are available in addition to this article. It can be downloaded at RJ-2019-037.zip

Footnotes

    Reuse

    Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".

    Citation

    For attribution, please cite this work as

    Baumann & Klymak, "The R Journal: Fixed Point Acceleration in R", The R Journal, 2019

    BibTeX citation

    @article{RJ-2019-037,
      author = {Baumann, Stuart and Klymak, Margaryta},
      title = {The R Journal: Fixed Point Acceleration in R},
      journal = {The R Journal},
      year = {2019},
      note = {https://doi.org/10.32614/RJ-2019-037},
      doi = {10.32614/RJ-2019-037},
      volume = {11},
      issue = {1},
      issn = {2073-4859},
      pages = {359-375}
    }