Fixed Point Acceleration in R

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.

Stuart Baumann , Margaryta Klymak
2019-08-20

Supplementary materials

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

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, "Fixed Point Acceleration in R", The R Journal, 2019

BibTeX citation

@article{RJ-2019-037,
  author = {Baumann, Stuart and Klymak, Margaryta},
  title = {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}
}